Next Article in Journal
Assessment of Climate Variability among Seasonal Trends Using In Situ Measurements: A Case Study of Punjab, Pakistan
Next Article in Special Issue
Tracking the Endogenous Dynamics of the Solfatara Volcano (Campi Flegrei, Italy) through the Analysis of Ground Thermal Image Temperatures
Previous Article in Journal
Energy and Environmental Implications of Using Energy-Harvesting Speed Humps in Nablus City, Palestine
Previous Article in Special Issue
Scaling Properties and Persistence of Long-Term Solar Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stellar Turbulent Convection: The Multiscale Nature of the Solar Magnetic Signature

Physics Department, University of Rome Tor Vergata, I-00133 Rome, Italy
*
Author to whom correspondence should be addressed.
Submission received: 30 June 2021 / Revised: 15 July 2021 / Accepted: 19 July 2021 / Published: 22 July 2021

Abstract

:
The multiscale dynamics associated with turbulent convection present in physical systems governed by very high Rayleigh numbers still remains a vividly disputed topic in the community of astrophysicists, and in general, among physicists dealing with heat transport by convection. The Sun is a very close star for which detailed observations and estimations of physical properties on the surface, connected to the processes of the underlying convection zone, are possible. This makes the Sun a unique natural laboratory in which to investigate turbulent convection in the hard turbulence regime, a regime typical of systems characterized by high values of the Rayleigh number. In particular, it is possible to study the geometry of convection using the photospheric magnetic voids (or simply voids), the quasi-polygonal quiet regions nearly devoid of magnetic elements, which cover the whole solar surface and which form the solar magnetic network. This work presents the most extensive statistics, both in the spatial scales studied (1–80 Mm) and in the temporal duration (SC 23 and SC 24), to investigate the multiscale nature of solar magnetic patterns associated with the turbulent convection of our star. We show that the size distribution of the voids, in the 1–80 Mm range, for the 317,870 voids found in the 692 analyzed magnetograms, is basically described by an exponential function.

1. Introduction

The energy flux through the stellar interior is one of the fundamental quantities that govern the stellar structure. This flux is largely convective in the outermost layers of Sun-like stars (e.g., [1]) and becomes convective throughout the whole stellar structure in low-mass stars and brown dwarfs (e.g., [2]). Therefore, the heat transport and dynamics in these convection zones (CZs) are defined by organized flows of plasma. These flows are believed to contribute substantially to developing and maintaining the global stellar magnetic dynamo process, the process that is responsible for the cyclic regeneration of the stellar/solar large-scale magnetic field, which is at the basis of all those processes known as stellar activity (e.g., [3,4]).
While researchers have reached consensus on the basic ingredients needed to trigger the global stellar dynamo process (see again, [4]), how exactly this dynamo works remains unclear (e.g., [5,6,7,8]). Furthermore, it should be remembered that in stellar CZs the Rayleigh number ( R a ) typically exceeds the value of 10 15 (e.g., [9,10,11]). We recall that R a quantifies the relative magnitude of thermal forcing to dissipative forces in the fluid convective flows. Indeed, it represents the buoyancy forcing divided by the product of the rate of heat diffusion and viscous drag, at least in the rather crude Boussinesq convection (i.e., nearly incompressible plasma) approximation in the outermost stellar CZs. The expression for R a is given by:
R a = g α Δ T L 3 ν κ
with g the local gravity, α the thermal expansion coefficient, Δ T the temperature difference between top and bottom CZs, L the CZ thickness, ν the kinematic viscosity and κ the plasma thermal diffusivity. R a is estimated to be of the order of 10 23 in the outer part of the Sun (e.g., [12]). The high R a value expected in stellar CZs has two main effects. The first is that the convection becomes strongly turbulent [13]. This regime has important implications on the structuring of convective cells and on the dynamic properties of temperature distribution within the fluid (for details refer to [14,15]). The second, of paramount importance in stars, but more broadly in all physical systems with high R a , involves the efficiency in the transport of energy by convection. This efficiency can be quantified by the Nusselt number, N u , which is the non-dimensional ratio of the total energy flux to the conductive (radiative) flux allowed for the same Δ T between top and bottom CZs.
When looking for a functional relationship between N u and R a , simple scaling arguments supported by laboratory experiments have suggested that N u R a β . However, although for β there are both theoretical evaluations, based on stability arguments and dimensional considerations, and estimates from laboratory experiments under various experimental conditions, the basic properties of heat transport in the hard turbulence regime remain unclear. In particular, it is not clear whether there exists an asymptotic regime that is supposed to occur at very high R a . Indeed, as R a increases, displacements are observed from a basic assumption that N u scales as R a 1 3 and this indicates the need for further analyses and experiments to clarify the heat transfer efficiency of the turbulent convective flows in high R a regimes (e.g., [16]).
In addition to these basic uncertainties in the physics governing turbulent convection, we must be aware that in the stellar case there are further possible complications. In fact, in this case, we must take into account that the dynamics of the plasma, stratified and compressible in the outermost regions of the star, is typically modified by the fact that the plasma flows interact with or even generate magnetic fields and, possibly, by the influence of the Coriolis force due to the stellar rotation.
This theoretical and experimental difficulty of quantifying the dependence of N u on R a in physical systems characterized by stratified plasma in the magneto-convective regime and very high R a makes the Sun the ideal laboratory for these studies.
Indeed, the Sun is the only star for which detailed observations and estimations of physical properties on the surface and in the CZ are currently possible and represents a unique natural laboratory in which to study turbulent convection in the hard turbulence regime. This is especially the case since spatial, temporal and spectral high-resolution observations of its photosphere are possible from ground-based (e.g., [17,18,19,20,21,22]) or space (e.g., [23,24]) instruments.
Detailed observation of the Sun’s surface and its scaling properties related to temperature, velocity or magnetic structures raises several interesting questions about pattern formation in turbulent convection conditions at high R a (e.g., [25,26,27,28,29,30,31,32]). For many of these questions there are still no certain answers. For example, what determines the distribution of size and average life (e.g., [21,31,33]) of convective structures? Are the classic scales such as granulation, mesogranulation and supergranulation real convective scales or do they instead emerge from self-organization of surface flows (e.g., [27,34]) or from methods of observation of the sun (e.g., [28,35]), i.e., emergent structures/patterns resulting from physical processes sampled, and consequently filtered, at particular spatial and temporal scales? How is the persistence of structures organized for many tens of hours (e.g., [26,31,36]) possible in a regime of hard turbulence convection?
Photospheric magnetic voids (PMVs) offer a unique opportunity for studying convective signatures in a physical system governed by very high R a [37,38].
PMVs, or simply voids, are multiscale (Mm to tens of Mm), quasi-polygonal quiet-sun regions nearly devoid of magnetic elements, which cover the whole solar surface [39,40] and which form the well-known magnetic network typically associated with the chromospheric network (e.g., [41,42,43,44,45,46]). The relative absence of a photospheric magnetic field in these regions is due to the advective transport process generated by convective flows.
Discussing the origin of the photospheric magnetic field is beyond the scope of this paper. Then, we refer the reader to the literature that deals with the dynamo processes, global and local, which are believed to be at the origin of the emergence, dynamics and cancellation of this magnetic field (e.g., [47,48,49,50,51]). Of interest for this work are the convective transport mechanisms of magnetic elements and the consequent magnetic pattern formation that is intimately linked to the boundaries of the underlying convective structures. In particular, we refer to the pattern formed by magnetic elements that are swept at the edges of convective structures, which leads to the formation of an irregular quasi-polygonal magnetic network that outlines the boundaries of convective cells [37,38,39,40].
Once we have the numerical tools to detect the PMVs, we study statistical properties that will be used to investigate the multiscale properties of the convective structures highlighted by the associated magnetic pattern. In this work we make use of the largest dataset of quiet-sun magnetograms used so far with the PMV algorithm to extend the statistics of previous works both in terms of spatial and temporal scales. PMVs were detected in the range 2–10 Mm in [39] and in the range 10–60 Mm in [40]. In this work we access the spatial scales from 1 to 80 Mm. Furthermore, we compare PMVs from SC23 and SC24 and we apply the technique to magnetograms from the Helioseismic and Magnetic Imager instrument for the first time.

2. Materials and Methods

2.1. Datasets

With the aim of identifying PMVs representative exclusively of photospheric convective patterns we selected quiet-sun line-of-sight magnetograms acquired during periods characterized by a low activity, typically associated with periods of solar minimum. This choice means that the possible effect of extended active regions on the magnetic pattern of convective origin is negligible. Data sets of quiet-sun magnetic images used in this work, and gathered from three different instruments on board three different space missions, i.e., SOHO, SDO and HINODE, are summarized in Table 1.
In detail, the line-of-sight magnetograms were acquired from the following sources.
  • The Solar Oscillations Investigation/Michelson Doppler Imager (SOI/MDI) instrument [52], on board the SOHO satellite, collected solar images from 1996 to 2011 using the spectral line Ni I 676.78 nm (e.g., [53,54]). The images were recorded by a 1024 × 1024 CCD camera in two spatial resolution modes: full disk and high-resolution of the central part of the disk. In this work, 191 original high-resolution images acquired during the solar activity minimum between SC 22 and SC 23 are used, covering a period of 7 months from 1 August 1996 to 28 February 1997. The images have a field of view of 11 × 11 arc min 2 with a plate scale of 0.625 arcsecond per pixel and a (diffraction-limited) resolution of 1.25 arc min (see http://soi.stanford.edu/sssc/progs/mdi/calib.html (accessed on 1 July 2021) for more information about data).
  • The Helioseismic and Magnetic Imager (HMI) is part of the NASA Solar Dynamics Observatory (SDO) mission [55]. HMI provides stabilized one-arcsecond-resolution full-disk Doppler velocity, line-of-sight magnetic flux, and continuum proxy images every 45 s, and vector magnetic field maps every 90 or 135 s depending on the image frame sequence selected. The solar image nearly fills the 4096 × 4096 pixel CCD camera allowing a pixel scale equal to 0.5 arcsecond per pixel [56,57]. HMI quiet magnetograms used in this work are the central sub-array of the original 4096 × 4096 full-disk images and they cover the whole mission from May 2010 to May 2019 (see http://jsoc.stanford.edu (accessed on 1 July 2021) for more information about data).
  • The Narrowband Filter Imager (NFI) of the Solar Optical Telescope (SOT) is on board the Hinode satellite [58]. The observation was recorded on 10 March 2007 between 11:37 and 14:34 UT, and views a portion of solar photosphere of 302 × 162 arcsecond 2 at the disc center. The original FITS header declares a different plate scale of 0.1476 arcsec/pixel and 0.1585 arcsec/pixel along x and y, respectively. Assuming for the solar radius the estimate of 696,342 ± 65 km [59] and the corresponding apparent radius of 959.6 ± 0.1 arcseconds [60], this difference corresponds to ≃7.9 km, which has a negligible effect on the PMV identification algorithm and is therefore not taken into account in the analysis (see http://sdc.uio.no/sdc/ (accessed on 1 July 2021) for more information about data).

2.2. PMV-Finding Algorithm

The PMV-finding procedure we adopt in this work is described in full detail in [40] (henceforth BSD14). The procedure is based on that proposed by [61] (henceforth AM98). Compared to the algorithm proposed by AM98, our climbing algorithm has been significantly improved, making the BSD14 procedure faster. Indeed, this algorithm produces results faster then the conventional algorithm proposed by AM98, and actually we have reduced the time required for the PMV identification (about 25 times shorter than in AM98). The reduction of calculation times was useful in the testing phase allowing a greater number of analysis runs.
The steps of the PMV-finding procedure are as follows: 1. binarization of a magnetogram by means of the classification of pixels as magnetic or non-magnetic; 2. detecting empty region in the magnetic pattern; 3. recursive growth of the empty circles; 4. labeling of the circles for the definition of the single PMV.
To understand the PMV-finding procedure we start from the definition of distance field ( D F ) reported in AM98. This in turn is based on the empty sphere method [62]. Our scalar distance field D : L 2 R is analogously defined, in a 2D space, as the distance of a given point x to the nearest magnetic element M:
D ( x ) = m i n j { | x M j | }
where M j , j = 1 , 2 , , N are the sites of the magnetic elements highlighted by the binarized magnetogram. The climbing algorithm of AM98, based on the path along the monotonically increasing D F ( x , y ) is replaced by a faster algorithm in which the climbing is carried out only on the centers of circles defined by a recursive algorithm of circle packing (Figure 1, Panel c). In each step a circle is drawn centered on the maximum of the current D F and having a radius equal to the local D F value (i.e., extending the radius up to the touch with the nearest magnetic structure). In the following step we insert the last generated circle as a fictitious magnetic area in the previously used binarized image and we calculate a new D F . Likewise in the Apollonian circle packing [63,64], a sequence of circles C 1 , C 2 , C 3 , having non-crescent diameters (therefore circles whose diameters are smaller than or at most equal to the diameters of the circles of the previous iteration) and filling the regions “empty” of magnetic elements is recursively generated. The iteration is stopped when the current C i radius R is smaller than an assigned R m i n (please see BSD14 for details). The climbing algorithm is then executed only for the centers of the filling circles, saving a noteworthy amount of CPU time.
In Figure 1 an example of the complete process of the circle packing algorithm coverage relative to a typical SOHO/MDI magnetogram is shown.
The rectangular region (a) shows a sub-sample of an original magnetogram. To define the rectangular magnetic structures (region b) we used a fixed threshold T equal to 3 times the standard deviation ( σ ) of the absolute magnetic field values in the selected sub-region. In detail, the thresholds are 16 G, 15 G and 98 G for SOHO/MDI, SDO/HMI and HINODE/SOT, respectively. Please refer to BSD14 for an in-depth discussion on the projection effects of the magnetic field ( B L o S = B n o r m a l ) and on the choice of the threshold values. The rectangular region (c) is obtained from the application of the adaptive Apollonian circle packing procedure to cover the void area and the labeling of each circle center to the appropriate PMV by the climbing algorithm. Different colors identify different PMVs. Finally, the rectangular region (d) shows an estimation of PMV borders.
For the research of the PMVs, above reported, we used the same threshold for any analyzed magnetogram, i.e., T = 3 σ . We have seen that this value is valid for two of the datasets, i.e., SDO/HMI and SOHO/MDI, and the derived distributions are very similar (refer to the “Results and Discussion” section).
In the case of Hinode SOT/SP we used a single high-resolution magnetogram that was sufficiently large enough to allow a statistically accurate analysis of the PMVs present in it. We refer to [39] (henceforth BGS13) for the details regarding the magnetogram used. Here we recall only the main characteristics of the magnetogram. We underline that, unlike the void-finding algorithm used in BGS13, in this work we used for consistency the same fast procedure just described in this section that was applied to the other two datasets.
To study the effect of the threshold in the detection of PMVs in the high-resolution magnetogram, we used different thresholds. In fact, we expect that as the threshold increases, the regions emptied by magnetic elements tend to increasingly resemble the PMVs obtained from magnetograms with lower spatial resolution and magnetic sensitivity (i.e., SOHO/MDI and SDO/HMI). This is because as the threshold increases, the smaller PMVs highlighted by the lowest magnetic field values will disappear. To better explain what has just been discussed, we report in Figure 2 the result of the application to the Hinode SOT/SP magnetogram of our procedure using two different thresholds, 3 σ and 6 σ , respectively. The effect discussed above is clearly seen. As the threshold value increases, the emptied regions outlined by the weaker fields tend to disappear and aggregate into larger regions. From a point of view of the distribution of PMVs as a function of size, this effect involves a gradual decrease in the measured slope and an approach to the slope observed for the larger scales, i.e., SDO/HMI and SOHO/MDI, due to the fact that the magnetograms of these two missions are less spatially resolved and less sensitive than the magnetogram derived from Hinode observations.

3. Results and Discussion

The application of our PMV-finding procedure to 692 magnetograms found and labeled 317,870 PMVs. The equivalent diameter, i.e., the diameter of the circle of equal area of the detected PMV, was calculated and the PMV distribution is reported in Figure 3. The distribution of PMVs corresponds to what in BSD14 is called Void Size Distribution (or VSD). The present analysis confirms what is reported and discussed in BSD14. In particular, the PMV distribution is described by a near-exponential function in the observed range. Indeed, if we pass a best-fit line in the PMV counts derived from the Hinode data (for T = 6 σ ) and in the averaged values of the counts derived from the SOHO/MDI and SDO/HMI magnetograms, we obtain for the line (black dash-dot line in Figure 3) the slope S = 0.058 ± 0.002 , which gives the exponential constant in the equation describing our PMV distribution:
P M V c d 10 0.058 × E D
where P M V c d is the PMV density in C o u n t s / a r c s e c 2 and E D is the PMV equivalent diameter in Mm. For the sake of completeness, we report that the residual distribution between the observed distribution and the exponential fit shows systematics that will be the subject of a future study. Indeed, as found in BSD14, the PMV distributions from SOHO/MDI and SDO/HMI data show a break at about 35 Mm. This scale is linked to one of the sources of systematics once a single exponential fit is performed on the entire range of spatial scales. The regime at large spatial scales stands up to 80 Mm without finding any preferential scales. Furthermore, in this work we find that the smaller spatial scales (Hinode data for T = 6 σ ) seem to fit better with the trend for the larger spatial scales (>35 Mm). Nevertheless, as discussed in BSD14, the break at 35 Mm could be related to a convective instability scale able to inject downflow plumes. A possible instability mechanism that could support such a scale has been described in [65]. These plumes could become a new site for magnetic field concentrations, reducing the number of voids with equivalent diameters larger than 35 Mm. Further investigation is needed to support this hypothesis.
The statistic presented in this work is much greater than in BSD14. Indeed, present data cover two solar cycles, and the SOHO/MDI and SDO/HMI magnetograms refer to regions of quiet-sun relative to SC 23 and SC 24, respectively. Thus, this underlines how the dimensions of the large (i.e., equivalent diameter > 3 Mm) PMVs are very similar in the two cycles.
The analysis of the Hinode SOT/SP high-resolution magnetogram performed with the same PMV-finding procedure used to study large scale ranges allows us to extend the distribution up to the equivalent size of about 1 Mm. Apart from expected problems of connection between the two distributions, no particular organization scales are highlighted.
Furthermore, the near-exponential function that describes the VSD can also be used to describe other physical quantities in systems characterized by turbulent convection. For instance, the probability density of normalized temperature fluctuations is close to an exponential in the hard-turbulence regime ( R a > 10 8 ) [13], which is the case of photospheric turbulent convection ( R a 10 11 ) [66], while it transits to a Gaussian in the soft turbulent convective regime ( R a < 10 7 ).
The absence of particular features in the distribution in the 1–80 Mm range rules out the presence of special scales (i.e., supergranulation or mesogranulation) when the magnetic pattern formed by surface convective flows is used to identify convective cell borders. This result supports the multiscale hypothesis of turbulent convective flows on the solar photosphere reported by several authors [39,40,67,68].
As clearly described in [67], theoretical models had anticipated the idea that convective structures, resulting from solar plasma flows, show a smooth spectrum, although the canonical literature typically refers to distinct scales of solar convection, such as granulation, mesogranulation and supergranulation. The multiscale view has been supported over time, not only in the papers just mentioned, but also in papers that studied the properties of convective velocity or intensity patterns at different scales [69,70], or introduced the concept of ‘trees of fragmenting granules’ to study families of replicated structures at different scales starting from single granules [71].

4. Conclusions

This work presents the most extensive statistics, both in the spatial scales studied and in the temporal duration, of photospheric regions substantially “emptied” of an intense magnetic field. By intense magnetic field we mean magnetic field values higher than a suitably defined threshold (see Section 2). Our method relies on the well-established hypothesis (e.g., [72]) that these “emptied” regions, which we have called PMVs or voids, are produced by magnetic elements that are swept to the boundaries of convective cells by superficial advection flows and underline their pattern. As described in Section 2.2, we identified the PMVs in the selected quiet LoS magnetograms using a PMV finder based on a fast climbing algorithm.
For PMV statistics (number density versus effective diameter) we found remarkable agreement with previous void analysis and with theoretical predictions that suggested the absence of preferential scales in the organization of convective structures at the surface of the Sun.
Another important aspect of our methodology is that this behavior is extended in the spatial range 1–80 Mm. This range covers all the scales classically referred to as “typical” of the photospheric convective organization, i.e., granulation, mesogranulation, supergranulation and giant cells. This multiscale behavior therefore seems to be typical of the turbulent convection of our star, as anticipated on theoretical grounds from the very beginning [67] and consequently supports the multiscale hypothesis of turbulent convection in stars.
A further important aspect to underline is the difference between turbulent convection and classical turbulence (for example the turbulence developed in terms of a hierarchy of scales in which the plasma velocity fluctuations follow the Kolmogorov scaling characterized by a power spectrum with a slope of −5/3). In this regard, we underline that the present work is focused on the analysis of the dimensions of convective cells in a range of scales far from the typical scales in which a local turbulence model can be applied, i.e., a scale of the order of the pressure scale height, 150 km. We refer to the paper [73] for an in-depth discussion on the nature of solar granulation, and in general of stellar convection, in terms of turbulence models.
Conclusions on whether the distribution we found is compatible or not with the presence of particular asymptotic regimes that are supposed to occur at very high R a and which are connected to heat transfer efficiency of the turbulent convection are more difficult to draw. The arduous future task will be that of combining the uncertainties in linking N u on R a , in physical systems characterized by stratified plasma in a magneto-convective regime and very high R a , to estimate the global efficiency in convective energy transport by refining the sophisticated models currently available but which, unfortunately, have computational domains that probe a range of scales that are still too small compared to the solar case.

Author Contributions

Conceptualization, F.B. and S.S.; methodology, S.S., F.B. and D.D.M.; software, S.S.; validation, F.B., D.D.M. and L.G.; formal analysis, all; investigation, all; resources, S.S. and F.B.; data curation, S.S.; writing—original draft preparation, F.B. and S.S.; writing—review and editing, all; supervision, F.B. and S.S.; funding acquisition, F.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Italian MIUR-PRIN grant 2017APKP7T on Circumterrestrial Environment: Impact of Sun–Earth Interaction and by the European Union’s Horizon 2020 research and innovation program under grant agreement No. 824135 (SOLARNET).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

SOHO is a project of international cooperation between ESA and NASA. The authors thank the MDI team for MDI data. SDO/HMI data are courtesy of NASA/SDO and the HMI science teams. Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as a domestic partner and NASA and STFC (UK) as international partners.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ledoux, P.; Schwarzschild, M.; Spiegel, E.A. On the Spectrum of Turbulent Convection. Astrophys. J. 1961, 133, 184. [Google Scholar] [CrossRef]
  2. Ventura, P.; Zeppieri, A. Non instantaneous mixing: Deuterium burning in very low mass stars and brown dwarfs. Astron. Astrophys. 1998, 340, 77–80. [Google Scholar]
  3. Brun, A.S.; Browning, M.K. Magnetism, dynamo action and the solar-stellar connection. Living Rev. Sol. Phys. 2017, 14, 4. [Google Scholar] [CrossRef] [PubMed]
  4. Charbonneau, P. Dynamo models of the solar cycle. Living Rev. Sol. Phys. 2020, 17, 4. [Google Scholar] [CrossRef]
  5. Brandenburg, A.; Bigazzi, A.; Subramanian, K. The helicity constraint in turbulent dynamos with shear. Mon. Not. R. Astron. Soc. 2001, 325, 685–692. [Google Scholar] [CrossRef] [Green Version]
  6. Phillips, A.; Brooke, J.; Moss, D. The importance of physical structure in solar dynamo models. Astron. Astrophys. 2002, 392, 713–727. [Google Scholar] [CrossRef]
  7. Bigazzi, A.; Ruzmaikin, A. The Sun’s Preferred Longitudes and the Coupling of Magnetic Dynamo Modes. Astrophys. J. 2004, 604, 944–959. [Google Scholar] [CrossRef] [Green Version]
  8. Lehtinen, J.J.; Spada, F.; Käpylä, M.J.; Olspert, N.; Käpylä, P.J. Common dynamo scaling in slowly rotating young and evolved stars. Nat. Astron. 2020, 4, 658–662. [Google Scholar] [CrossRef] [Green Version]
  9. Berrilli, F. Spatio-temporal patterns in solar surface convection. Mem. S.A.It. 2003, 74, 576. [Google Scholar]
  10. Brun, A.; Miesch, M. Stellar convection simulations. Scholarpedia 2008, 3, 4278. [Google Scholar] [CrossRef]
  11. Nigro, G. Attracted by the fascinating magnetism of the Sun. Nuovo Cimento C Geophys. Space Phys. C 2019, 42, 2. [Google Scholar] [CrossRef]
  12. Niemela, J.J.; Skrbek, L.; Sreenivasan, K.R.; Donnelly, R.J. Turbulent convection at very high Rayleigh numbers. Nature 2000, 404, 837–840. [Google Scholar] [CrossRef]
  13. Yakhot, V. Probability distributions in high-Rayleigh number Bénard convection. Phys. Rev. Lett. 1989, 63, 1965–1967. [Google Scholar] [CrossRef] [PubMed]
  14. Libchaber, A. From Chaos to Turbulence in Benard Convection. Proc. R. Soc. Lond. Ser. A 1987, 413, 63–69. [Google Scholar] [CrossRef]
  15. Castaing, B.; Gunaratne, G.; Kadanoff, L.; Libchaber, A.; Heslot, F. Scaling of hard thermal turbulence in Rayleigh-Benard convection. J. Fluid Mech. 1989, 204, 1–30. [Google Scholar] [CrossRef] [Green Version]
  16. Urban, P.; Hanzelka, P.; Musilová, V.; Králík, T.; La Mantia, M.; Srnka, A.; Skrbek, L. Heat transfer in cryogenic helium gas by turbulent Rayleigh-Bénard convection in a cylindrical cell of aspect ratio 1. New J. Phys. 2014, 16, 053042. [Google Scholar] [CrossRef] [Green Version]
  17. Komm, R.; Mattig, W.; Nesis, A. The small-scale velocity field in the solar photosphere. Astron. Astrophys. 1991, 243, 251–262. [Google Scholar]
  18. Montagne, M.; Mueller, R.; Vigneau, J. The photosphere of the Sun: Statistical correlations between magnetic field, intensity and velocity. Astron. Astrophys. 1996, 311, 304–310. [Google Scholar]
  19. Berrilli, F.; Florio, A.; Consolini, G.; Bavassano, B.; Briand, C.; Bruno, R.; Caccin, B.; Carbone, V.; Ceppatelli, G.; Egidi, A.; et al. Dependence of the photospheric vertical flow characteristics on the granule dimension. Astron. Astrophys. 1999, 344, L29–L32. [Google Scholar]
  20. Hirzberger, J.; Koschinsky, M.; Kneer, F.; Ritter, C. High resolution 2D-spectroscopy of granular dynamics. Astron. Astrophys. 2001, 367, 1011–1021. [Google Scholar] [CrossRef] [Green Version]
  21. Berrilli, F.; Consolini, G.; Pietropaolo, E.; Caccin, B.; Penza, V.; Lepreti, F. 2-D multiline spectroscopy of the solar photosphere. Astron. Astrophys. 2002, 381, 253–264. [Google Scholar] [CrossRef] [Green Version]
  22. Hoekzema, N.M.; Rimmele, T.R.; Rutten, R.J. Small-scale topology of solar atmospheric dynamics. V. Acoustic events and internetwork grains. Astron. Astrophys. 2002, 390, 681–690. [Google Scholar] [CrossRef]
  23. Welsch, B.T. The photospheric Poynting flux and coronal heating. Pub. Astron. Soc. Jpn. 2015, 67, 18. [Google Scholar] [CrossRef] [Green Version]
  24. Utz, D.; del Toro Iniesta, J.C.; Bellot-Rubio, L.; Bodnárová, M.; Muller, R.; Bárta, M.; Thonhofer, S.; Hanslmeier, A. Long time variations of Magnetic Bright Points observed by Hinode/SOT. Cent. Eur. Astrophys. Bull. 2015, 39, 91–100. [Google Scholar]
  25. Straus, T.; Deubner, F.L.; Fleck, B. Is mesogranulation a distinct regime of convection? Astron. Astrophys. 1992, 256, 652–659. [Google Scholar]
  26. Getling, A.V.; Brandt, P.N. Regular structures of the solar photosphere. (Persistence of the granular field and trenching in the brightness relief). Astron. Astrophys. 2002, 382, L5–L8. [Google Scholar] [CrossRef]
  27. Rast, M.P. The Scales of Granulation, Mesogranulation, and Supergranulation. Astrophys. J. 2003, 597, 1200–1210. [Google Scholar] [CrossRef]
  28. Berrilli, F.; Del Moro, D.; Russo, S.; Consolini, G.; Straus, T. Spatial Clustering of Photospheric Structures. Astrophys. J. 2005, 632, 677–683. [Google Scholar] [CrossRef]
  29. Caroli, A.; Giannattasio, F.; Fanfoni, M.; Del Moro, D.; Consolini, G.; Berrilli, F. Turbulent convective flows in the solar photospheric plasma. J. Plasma Phys. 2015, 81, 495810514. [Google Scholar] [CrossRef]
  30. Cossette, J.F.; Rast, M.P. Supergranulation as the Largest Buoyantly Driven Convective Scale of the Sun. Astrophys. J. Lett. 2016, 829, L17. [Google Scholar] [CrossRef] [Green Version]
  31. Giannattasio, F.; Berrilli, F.; Consolini, G.; Del Moro, D.; Gošić, M.; Bellot Rubio, L. Occurrence and persistence of magnetic elements in the quiet Sun. Astron. Astrophys. 2018, 611, A56. [Google Scholar] [CrossRef] [Green Version]
  32. Viavattene, G.; Consolini, G.; Giovannelli, L.; Berrilli, F.; Del Moro, D.; Giannattasio, F.; Penza, V.; Calchetti, D. Testing the Steady-State Fluctuation Relation in the Solar Photospheric Convection. Entropy 2020, 22, 716. [Google Scholar] [CrossRef]
  33. Del Moro, D.; Berrilli, F.; Duvall, T.L.J.; Kosovichev, A.G. Dynamics and Structure of Supergranulation. Sol. Phys. 2004, 221, 23–32. [Google Scholar] [CrossRef]
  34. Ploner, S.R.O.; Solanki, S.K.; Gadun, A.S. Is solar mesogranulation a surface phenomenon? Astron. Astrophys. 2000, 356, 1050–1054. [Google Scholar]
  35. Roudier, T.; Rieutord, M.; Malherbe, J.M.; Vigneau, J. Determination of horizontal velocity fields at the sun’s surface with high spatial and temporal resolution. Astron. Astrophys. 1999, 349, 301–311. [Google Scholar]
  36. Pandey, A.; Scheel, J.D.; Schumacher, J. Turbulent superstructures in Rayleigh-Bénard convection. Nat. Commun. 2018, 9, 2118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Roche, P.E. The ultimate state of convection: A unifying picture of very high Rayleigh numbers experiments. New J. Phys. 2020, 22, 073056. [Google Scholar] [CrossRef]
  38. Blass, A.; Zhu, X.; Verzicco, R.; Lohse, D.; Stevens, R.J.A.M. Flow organization and heat transfer in turbulent wall sheared thermal convection. J. Fluid Mech. 2020, 897, A22. [Google Scholar] [CrossRef]
  39. Berrilli, F.; Scardigli, S.; Giordano, S. Multiscale Magnetic Underdense Regions on the Solar Surface: Granular and Mesogranular Scales. Sol. Phys. 2013, 282, 379–387. [Google Scholar] [CrossRef] [Green Version]
  40. Berrilli, F.; Scardigli, S.; Del Moro, D. Magnetic pattern at supergranulation scale: The void size distribution. Astron. Astrophys. 2014, 568, A102. [Google Scholar] [CrossRef] [Green Version]
  41. Simon, G.W.; Leighton, R.B. Velocity Fields in the Solar Atmosphere. III. Large-Scale Motions, the Chromospheric Network, and Magnetic Fields. Astrophys. J. 1964, 140, 1120. [Google Scholar] [CrossRef]
  42. Labonte, B.J.; Howard, R. The Magnetic Flux in the Quiet Sun Network. Sol. Phys. 1982, 80, 15–19. [Google Scholar] [CrossRef]
  43. Liu, Y.; Zhang, H.; Ai, G.; Wang, H.; Zirin, H. The evolution of magnetic network elements in the quiet Sun. Astron. Astrophys. 1994, 283, 215–220. [Google Scholar]
  44. Hagenaar, H.J.; Schrijver, C.J.; Title, A.M. The Distribution of Cell Sizes of the Solar Chromospheric Network. Astrophys. J. 1997, 481, 988–995. [Google Scholar] [CrossRef]
  45. Berrilli, F.; Florio, A.; Ermolli, I. On the Geometrical Properties of the Chromospheric Network. Sol. Phys. 1998, 180, 29–45. [Google Scholar] [CrossRef]
  46. Martínez González, M.J.; Bellot Rubio, L.R.; Solanki, S.K.; Martínez Pillet, V.; Del Toro Iniesta, J.C.; Barthol, P.; Schmidt, W. Resolving the Internal Magnetic Structure of the Solar Network. Astrophys. J. Lett. 2012, 758, L40. [Google Scholar] [CrossRef]
  47. Cattaneo, F. On the Origin of Magnetic Fields in the Quiet Photosphere. Astrophys. J. Lett. 1999, 515, L39–L42. [Google Scholar] [CrossRef]
  48. Tobias, S.M.; Cattaneo, F.; Brummell, N.H. On the Generation of Organized Magnetic Fields. Astrophys. J. 2011, 728, 153. [Google Scholar] [CrossRef]
  49. Okunev, O.V.; Domínguez Cerdeña, I.; Puschmann, K.G.; Kneer, F.; Sánchez Almeida, J. Quiet sun magnetic fields vs. polar faculae-local vs. global dynamo? Astron. Nachrichten 2005, 326, 205–207. [Google Scholar] [CrossRef]
  50. Wiegelmann, T.; Thalmann, J.K.; Solanki, S.K. The magnetic field in the solar atmosphere. Astron. Astrophys. Rev. 2014, 22, 78. [Google Scholar] [CrossRef]
  51. Cameron, R.; Schüssler, M. The crucial role of surface magnetic fields for the solar dynamo. Science 2015, 347, 1333–1335. [Google Scholar] [CrossRef] [Green Version]
  52. Scherrer, P.H.; Bogart, R.S.; Bush, R.I.; Hoeksema, J.T.; Kosovichev, A.G.; Schou, J.; Rosenberg, W.; Springer, L.; Tarbell, T.D.; Title, A.; et al. The Solar Oscillations Investigation—Michelson Doppler Imager. Sol. Phys. 1995, 162, 129–188. [Google Scholar] [CrossRef]
  53. Criscuoli, S.; Ermolli, I.; Del Moro, D.; Giorgi, F.; Tritschler, A.; Uitenbroek, H.; Vitas, N. Line Shape Effects on Intensity Measurements of Solar Features: Brightness Correction to SOHO MDI Continuum Images. Astrophys. J. 2011, 728, 92. [Google Scholar] [CrossRef]
  54. Chitta, L.P.; Jain, R.; Kariyappa, R.; Jefferies, S.M. Observations of the Interaction of Acoustic Waves and Small-scale Magnetic Fields in a Quiet Sun. Astrophys. J. 2012, 744, 98. [Google Scholar] [CrossRef]
  55. Pesnell, W.D.; Thompson, B.J.; Chamberlin, P.C. The Solar Dynamics Observatory (SDO). Sol. Phys. 2012, 275, 3–15. [Google Scholar] [CrossRef] [Green Version]
  56. Schou, J.; Scherrer, P.H.; Bush, R.I.; Wachter, R.; Couvidat, S.; Rabello-Soares, M.C.; Bogart, R.S.; Hoeksema, J.T.; Liu, Y.; Duvall, T.L.; et al. Design and Ground Calibration of the Helioseismic and Magnetic Imager (HMI) Instrument on the Solar Dynamics Observatory (SDO). Sol. Phys. 2012, 275, 229–259. [Google Scholar] [CrossRef] [Green Version]
  57. Scherrer, P.H.; Schou, J.; Bush, R.I.; Kosovichev, A.G.; Bogart, R.S.; Hoeksema, J.T.; Liu, Y.; Duvall, T.L.; Zhao, J.; Title, A.M.; et al. The Helioseismic and Magnetic Imager (HMI) Investigation for the Solar Dynamics Observatory (SDO). Sol. Phys. 2012, 275, 207–227. [Google Scholar] [CrossRef] [Green Version]
  58. Tsuneta, S.; Ichimoto, K.; Katsukawa, Y.; Nagata, S.; Otsubo, M.; Shimizu, T.; Suematsu, Y.; Nakagiri, M.; Noguchi, M.; Tarbell, T.; et al. The Solar Optical Telescope for the Hinode Mission: An Overview. Sol. Phys. 2008, 249, 167–196. [Google Scholar] [CrossRef]
  59. Emilio, M.; Kuhn, J.R.; Bush, R.I.; Scholl, I.F. Measuring the Solar Radius from Space during the 2003 and 2006 Mercury Transits. Astrophys. J. 2012, 750, 135. [Google Scholar] [CrossRef] [Green Version]
  60. Egidi, A.; Caccin, B.; Sofia, S.; Heaps, W.; Hoegy, W.; Twigg, L. High-Precision Measurements of the Solar Diameter and Oblateness by the Solar Disk Sextant (SDS) Experiment. Sol. Phys. 2006, 235, 407–418. [Google Scholar] [CrossRef]
  61. Aikio, J.; Mähönen, P. A Simple Void-searching Algorithm. Astrophys. J. 1998, 497, 534–540. [Google Scholar] [CrossRef] [Green Version]
  62. Einasto, J.; Einasto, M.; Gramann, M. Structure and formation of superclusters -IX. Self-similarity of voids. Mon. Not. R. Astron. Soc. 1989, 238, 155–177. [Google Scholar] [CrossRef] [Green Version]
  63. Bessis, D.; Demko, S. Generalized Apollonian packings. Commun. Math. Phys. 1990, 134, 293–319. [Google Scholar] [CrossRef]
  64. Oliveira, R.S.; Andrade, J.S., Jr.; Andrade, R.F.S. Fluid flow through Apollonian packings. Phys. Rev. E 2010, 81, 047302. [Google Scholar] [CrossRef] [PubMed]
  65. Rast, M.P. On the Nature of “Exploding” Granules and Granule Fragmentation. Astrophys. J. 1995, 443, 863. [Google Scholar] [CrossRef]
  66. Bray, R.J.; Loughhead, R.E.; Durrant, C.J. Book-Review—the Solar Granulation—ED.2. J. Br. Astron. Assoc. 1984, 94, 196. [Google Scholar]
  67. Nordlund, Å.; Stein, R.F.; Asplund, M. Solar Surface Convection. Living Rev. Sol. Phys. 2009, 6, 2. [Google Scholar] [CrossRef] [Green Version]
  68. Yelles Chaouche, L.; Moreno-Insertis, F.; Martínez Pillet, V.; Wiegelmann, T.; Bonet, J.A.; Knölker, M.; Bellot Rubio, L.R.; del Toro Iniesta, J.C.; Barthol, P.; Gandorfer, A.; et al. Mesogranulation and the Solar Surface Magnetic Field Distribution. Astrophys. J. Lett. 2011, 727, L30. [Google Scholar] [CrossRef] [Green Version]
  69. Schrijver, C.J.; Hagenaar, H.J.; Title, A.M. On the Patterns of the Solar Granulation and Supergranulation. Astrophys. J. 1997, 475, 328–337. [Google Scholar] [CrossRef] [Green Version]
  70. Berrilli, F.; Del Moro, D.; Consolini, G.; Pietropaolo, E.; Duvall, T.L., Jr.; Kosovichev, A.G. Structure Properties of Supergranulation and Granulation. Sol. Phys. 2004, 221, 33–45. [Google Scholar] [CrossRef]
  71. Roudier, T.; Muller, R. Relation between families of granules, mesogranules and photospheric network. Astron. Astrophys. 2004, 419, 757–762. [Google Scholar] [CrossRef]
  72. Zwaan, C. Elements and patterns in the solar magnetic field. Ann. Rev. Astron. Astrophys. 1987, 25, 83–111. [Google Scholar] [CrossRef]
  73. Nordlund, A.; Spruit, H.C.; Ludwig, H.G.; Trampedach, R. Is stellar granulation turbulence? Astron. Astrophys. 1997, 328, 229–234. [Google Scholar]
Figure 1. The circle packing algorithm on a SOHO/MDI magnetogram: (a) the original magnetogram, (b) the binarized version of the magnetogram, (c) the circle filling result (same color means association with the same PMV), (d) PMV borders (shown with solid black line).
Figure 1. The circle packing algorithm on a SOHO/MDI magnetogram: (a) the original magnetogram, (b) the binarized version of the magnetogram, (c) the circle filling result (same color means association with the same PMV), (d) PMV borders (shown with solid black line).
Atmosphere 12 00938 g001
Figure 2. We show the result of the application to the Hinode SOT/SP magnetogram of our procedure using two different thresholds, 3 σ (upper panels) and 6 σ (lower panels). As the threshold value increases, the emptied regions outlined by the weaker fields tend to disappear and aggregate into larger regions.
Figure 2. We show the result of the application to the Hinode SOT/SP magnetogram of our procedure using two different thresholds, 3 σ (upper panels) and 6 σ (lower panels). As the threshold value increases, the emptied regions outlined by the weaker fields tend to disappear and aggregate into larger regions.
Atmosphere 12 00938 g002
Figure 3. Distribution of PMV number versus equivalent diameter of voids. The plot shows the number density of PMVs obtained with our void detection algorithm. The different curves were obtained using different LoS magnetograms and magnetic thresholds. The circles denote the number density of “small”-scale PMVs detected in the Hinode magnetogram. The different colors refer to different magnetic thresholds. In detail, black circles refer to a magnetic threshold T = 3 σ , red circles refer to a magnetic threshold T = 3.5 σ , blue circles refer to a magnetic threshold T = 4 σ , and green circles refer to a magnetic threshold T = 6 σ . It is worth noting that as the threshold increases the number of small-sized HINODE PMVs decreases and the slope of the distributions approaches the slope derived for full-disk magnetograms. The pink diamonds and the brown squares denote the number density of “medium/large”-scale PMVs detected in SOHO/MDI (SC 22) and SDO/HMI (SC 23) magnetograms, respectively. The black dash-dot line, with slope S = 0.058 ± 0.0024 , represents the exponential behavior of PMV size distribution.
Figure 3. Distribution of PMV number versus equivalent diameter of voids. The plot shows the number density of PMVs obtained with our void detection algorithm. The different curves were obtained using different LoS magnetograms and magnetic thresholds. The circles denote the number density of “small”-scale PMVs detected in the Hinode magnetogram. The different colors refer to different magnetic thresholds. In detail, black circles refer to a magnetic threshold T = 3 σ , red circles refer to a magnetic threshold T = 3.5 σ , blue circles refer to a magnetic threshold T = 4 σ , and green circles refer to a magnetic threshold T = 6 σ . It is worth noting that as the threshold increases the number of small-sized HINODE PMVs decreases and the slope of the distributions approaches the slope derived for full-disk magnetograms. The pink diamonds and the brown squares denote the number density of “medium/large”-scale PMVs detected in SOHO/MDI (SC 22) and SDO/HMI (SC 23) magnetograms, respectively. The black dash-dot line, with slope S = 0.058 ± 0.0024 , represents the exponential behavior of PMV size distribution.
Atmosphere 12 00938 g003
Table 1. For each instrument the table reports the number of analyzed magnetograms, the beginning and end of the period used to select the quiet-sun magnetograms, the pixel scale of the magnetograms, the size of the magnetograms in pixels and the total number of identified PMVs.
Table 1. For each instrument the table reports the number of analyzed magnetograms, the beginning and end of the period used to select the quiet-sun magnetograms, the pixel scale of the magnetograms, the size of the magnetograms in pixels and the total number of identified PMVs.
SOHO/MDISDO/HMIHINODE/SOT
Number of magnetograms1915001
Start1996-08-012010-05-122007-03-10
End1997-02-282019-05-272007-03-10
Pixel scale (asec/pixel)0.6200.5040.150
Size1024 × 10241024 × 10242047 × 994
Detected PMVs133,569183,359942
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Scardigli, S.; Berrilli, F.; Del Moro, D.; Giovannelli, L. Stellar Turbulent Convection: The Multiscale Nature of the Solar Magnetic Signature. Atmosphere 2021, 12, 938. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12080938

AMA Style

Scardigli S, Berrilli F, Del Moro D, Giovannelli L. Stellar Turbulent Convection: The Multiscale Nature of the Solar Magnetic Signature. Atmosphere. 2021; 12(8):938. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12080938

Chicago/Turabian Style

Scardigli, Stefano, Francesco Berrilli, Dario Del Moro, and Luca Giovannelli. 2021. "Stellar Turbulent Convection: The Multiscale Nature of the Solar Magnetic Signature" Atmosphere 12, no. 8: 938. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12080938

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop