Free Access
Issue
A&A
Volume 603, July 2017
Article Number A134
Number of page(s) 8
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201630377
Published online 20 July 2017

© ESO, 2017

1. Introduction

The investigation of rotation profiles of different layers of the solar interior and its atmosphere is important for our understanding of the physical processes driving the solar activity and the connections between the Sun and heliosphere. It is one of the goals of the forthcoming Solar Orbiter mission and other observational missions as well. In fact, studies of the solar rotation rates already began in the 17th century after the first discovery of sunspots. Despite a large amount of accumulated data, we still need more observations and improved mathematical modelling to get more insight into the interconnection of the various physical processes operating at different radial distances inside and outside the Sun.

An important dynamical indicator that is connected to both the solar interior and its atmosphere is the solar wind. The solar wind, a global plasma flow from the Sun into the interplanetary space, is characterised by a variety of physical properties. Regarding characteristic flow patterns, the solar wind is bimodal and consists of the so-called fast and slow winds, and includes co-rotating interaction regions (CIRs) that represent a transitional flow pattern with a significant cross gradient of the speed. The characteristic velocities, location, morphology, and other plasma parameters of this flow pattern vary with the solar cycle (Miralles et al. 2001, 2002, 2004). For instance, the speed of the fast solar wind is about 600−800 km s-1 along “open” magnetic field lines that are identified as coronal holes (CHs; Krieger et al. 1973; Pneuman 1973; Noci 1973). Coronal holes are large, magnetically unipolar (Wilcox 1968; Altschuler et al. 1972; Timothy et al. 1975) low-density plasma (Munro & Withbroe 1972) regions in the solar corona. In epochs of the solar minimum, large CHs are mainly concentrated at high latitudes around the polar regions, while the main streamer belt and the heliospheric current sheet are located near the equator, enabling the formation of the slow solar wind flow pattern with approximately 400 km s-1 velocity. During solar maxima, such a strict dichotomy of the fast and slow wind patterns disappears as the solar atmosphere becomes populated with many small fragmentised CHs and streamers. This situation results in a more or less homogenous mixture of the fast, slow, and CIR flow patterns both in latitude and longitude. Therefore, the study of the rotational and statistical properties of CHs (as containers of fast solar wind patterns) is very important. As a matter of fact, their distributions are of a large-scale nature and seemingly are linked to the global solar magnetic field. Consequently, they are often considered as indicators of the large-scale magnetic field in the solar corona (Obridko & Shelting 1989; Bumba et al. 1995). Recently, Bilenko & Tavastsherna (2016) investigated some temporal and spatial regularities in CH distributions during several solar cycles. Their results indicate that the longitudinal distribution of non-polar CHs follows all structural changes of the global magnetic field. In general, the formation and acceleration mechanisms of various flow patterns are linked to the understanding of the near Sun physical processes and their understanding requires investigation of the heliosphere as a global system (for instance see recently reported analytical modelling; Shergelashvili & Fichtner 2012). Currently, the attention is stimulated by the forthcoming Solar Orbiter space mission.

The first attempts to study rotation rate distributions of CHs over different latitudes indicated rigid rotational profiles (Wagner 1975, 1976; Timothy et al. 1975; Adams 1976; Bohlin 1977; Hiremath & Hegde 2013; Japaridze et al. 2015). However, more recent studies showed the presence of pronounced differential rotations (Shelke & Pande 1985; Obridko & Shel’Ting 1988; Insley et al. 1995; Navarro-Peralta & Sanchez-Ibarra 1994). Despite several improvements of the observational methods and the automation of search algorithms, final conclusions regarding this issue are currently not available. Some researchers report that the CH rotation rates vary with the solar cycle (Harvey & Sheeley 1987; Obridko & Shel’Ting 1988; Navarro-Peralta & Sanchez-Ibarra 1994). During solar minimum CHs rotate more rigidly unlike to the maximum. In this regard Nash et al. (1988) concluded that this difference may be a consequence of connection between CHs and randomly emerging photospheric magnetic fields. In addition, Mancuso & Giordano (2011) reported that coronal magnetic structures rotate less differentially and much faster than the lower layers of the solar atmosphere, especially above of 40° latitude.

In the current paper, we investigate the CH rotational dynamics covering the period from 1 January 2013 to 20 April 2015, which included the prolonged peak of solar cycle 24. The results are analysed regarding the distributions over latitudinal locations and effective areas of CHs separately. Finally, we compare the CH rotational rates to the photospheric rotation and rotational profiles of the solar interior obtained from helioseismic measurements (Schou et al. 1998).

2. Data and method

2.1. Observational data

We have retrieved data of the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) from the Heliophysics Events Knowledgebase (HEK) via the SolarSoft (SSW) IDL tools for the period from 1 January 2013 to 20 April 2015, which corresponds to the maximum of solar cycle 24. Using the database downloaded from HEK, we created catalogues of coronal holes from 2010 up to date. This catalogue is created in the framework of the SOLSPANET project (“Solar and Space Weather Network of Excellence”) and is publicly available1.

The AIA is equipped with four 4096 × 4096 detectors with a pixel size of 0.6 arcsec. The AIA records a full set of near-simultaneous images with a 12 s cadence. The instrument contains 10 different wavelength channels; 7 extreme ultraviolet (EUV) channels (131, 171, 193, 211, 335, 94, and 304 Å), which produce the full-disc images of the corona and transition region; and 3 channels in white light and UV bands (Lemen et al. 2012; Boerner et al. 2011).

Even though CHs are clearly visible in X-ray and EUV images as dark cavities, it is difficult to detect these objects using a single EUV wavelength. Filaments also resemble dark structures in EUV images, and in some cases, it is hard to distinguish and identify CHs. During the first studies of CHs, they were identified manually by observers. In past years, various automatic detection methods have been developed, most of which are based on the intensity threshold method (Delouille et al. 2007; Kirk et al. 2009; Krista & Gallagher 2009; Rotter et al. 2012). Recently, geometrical classification (Reiss et al. 2014) and supervised machine learning methods (Reiss et al. 2015) of these coronal features have been developed, which represent more advanced CH detection approaches. In this work, we used data from the spatial possibilistic clustering algorithm (SPoCA), which is based on fuzzy clustering of intensity values (Verbeeck et al. 2013, 2014). This algorithm provides the exact localisation and characterisation of solar active regions and coronal holes on SDO/AIA193 Å wavelength images. This channel provides images of the corona at approximately 3700 ± 700 km above the photosphere (Aschwanden et al. 2013). The SPoCA yields CHs data every four hours and gives us the possibility to track the motion of individual CHs as described in Verbeeck et al. (2013).

2.2. Method of analysis

The rotational velocities of particular CHs were calculated using the position of the geometric centre of a polygon for each object. This method is justified because we measure the velocities of CHs every four hours. This time interval is too short for significant modifications of the shape and area of the CHs. Therefore, the velocity of the geometric centre of the polygon corresponds to the velocity of the CH itself.

Obviously, there are exceptions, in particular, significant changes might occur in the surface area or shape of the CHs, and then the velocity of the geometric centre may no longer reflect the movement of the CH. The CHs might also be broken in a few parts or merge with other CHs. Such cases are excluded from our statistical analysis.

To investigate the probable differential motion of coronal holes, we divided the solar surface into 32 latitudinal zones, each 3.75 degrees wide with the following intervals: from 0 to ± 3.75°, ± 3.75−7.5°, ± 7.5−11.25°, ..., ± 56.25−60°. We then classified the CHs into three groups, viz. relatively small-, medium-, and large-sized objects, namely with an area less than 10 000 Mm2, between 10 000−40 000 Mm2 and larger than 40 000 Mm2, respectively. The method of CH area estimation is standard (Verbeeck et al. 2014) and the values are given in related catalogues (i.e. SPoCA and SOLSPANET).

thumbnail Fig. 1

Distribution of 3056 CHs with latitude and areas during the years 2013–2015. Black crosses represent the individual CHs.

According to these criteria, we developed a special programme using the Interactive Data Language (IDL). This code uses the heliographic coordinates of the CHs geometric centres to calculate the displacement of each object after every four hours. Finally, it provides a measurement of the angular velocities in the different latitudinal zones and for the three area groups. Using this method, we calculated in total angular velocities of 540 objects during the studied time span.

3. Results

3.1. Distribution of coronal holes over latitude and areas

We studied the CH distributions with respect to latitudes and CH areas. The results are presented in Figs. 1 and 2. The horizontal axis corresponds to the latitude; the vertical axis denotes the average area of the CHs. Black crosses correspond to individual CHs, for each graph. Figure 1 shows the 2D distribution for all detected CHs. Figure 2 shows the distribution separately for small (panel a), medium (panel b), and large (panel c) CHs.

As we can see, in this epoch (solar maximum) CHs are observed at all latitudes. This is in contrast to solar minimum (de Toma 2011), when they are only found near high latitudes and at the poles. The prime characteristic of this distribution is that small CHs, with an area less than 2000 Mm2, are not observed in the equatorial zone, as is clear from Fig. 2, panel a. According to panel b of Fig. 2, extremely large CHs are not found is this region either.

We noticed a north-south asymmetry in the distribution of coronal holes during the studied time window. In the northern hemisphere, coronal holes were observed about 1.48 times more frequently (1822, 59, 62%) than in the southern hemisphere (1234, 40, 38%). Accordingly, the total number of detected CHs was 3056.

In general, the sizes of CHs vary within a quite wide range: the maximum size can be one-third or more of the solar disc area. In our case, the smallest detected coronal hole was 71 Mm2 (0.0012% of the solar disc surface area). The largest coronal hole occupied about 203.257 Mm2 (9.24%).

Finally, as mentioned above, for the study of the distribution of CHs over the solar disc we had registered in total 3056 CHs. However, we selected only 540 CHs for the calculation of the rotation rates, using the following criteria: (i) they must be located on the disc within the 120° latitude and 120° longitude box; (ii) the lifetime must be at least 12 h (i.e. they must appear in at least three snapshots); and (iii) the form and the area of the CH polygon must not change significantly during the observational time. The latitudinal distribution of these selected objects is not homogenous (Fig. 3) and is similar to the distribution of the whole CH set shown in Fig. 1. Although we consider 3.75° wide latitudinal zones, at least four CHs were studied in each zone, while less than 10 objects were present in only six latitudinal zones.

The latitudinal distribution of CHs divided by the area groups is presented in Fig. 4. The distribution of CHs having areas less than 10 000 Mm2 is similar to the distribution of all CHs because they constitute the majority of the observed objects. As a matter of fact, approximately 53% of the total number of CHs is in this category (i.e. 285 CHs, Fig. 4 panel (a)). Also, 195 average-sized objects were studied in all latitudinal zones (Fig. 4 panel b). Coronol holes with areas larger than 40 000 Mm2 are the rarest; only 60 objects fall in this category (Fig. 4 panel c).

3.2. Rotation rates of coronal holes

After the calculation of the CH rotation rates, we concluded that their rotation profiles differ significantly from the photospheric differential rotation (Fig. 5). These results indicate that the latitudinal gradients of the rotation rate are noticeably smaller than those observed in the photosphere (see Fig. 5). The average rotation rate for all CHs is 13.86(±0.05)°/ d. In the vicinity of the equator the velocity is approximately 14.15(±0.07)°/ d (i.e. a rotation period of 25.5 days). This value gradually decreases towards the poles, and in the vicinity of ± 55−60° latitude the rotation rates are 12.82(±0.07)°/ d in the southern hemisphere and 12.9( ± 0.07)°/ d in the northern hemisphere (Fig. 5, panel d).

The four panels in Fig. 5 present the rotation profiles for all, small-, medium-, and large-sized CHs, respectively. We have overlayed the photospheric rotation profiles estimated by two different methods: first, from the helioseismology data (Schou et al. 1998; black dashed line) and, second, the Doppler shift measurement result (Snodgrass & Ulrich 1990; black dot-dashed line). The black asterisks denote the observed values of the rotation rates of CHs. The red and blue curves are obtained by fitting the parabolic law A + Bsin2(θ) + Csin4(θ), representing the latitudinal profile of the rotation rate of CHs in the northern and southern hemispheres, respectively (Fig. 5).

Each of the plots in Fig. 5 indicates that the rotation profiles of CHs differ from both estimates of the photospheric rotation. In particular, the sidereal angular velocities of CHs, measured with the low corona EUV data, are higher than those rates in the photosphere above 35–40 degrees latitude; we used a constant for transformation from synodic to sidereal rotation rate given in Wittmann (1996). These results are in good agreement with other measurements of coronal rotation rates, for instance, those reported by Mancuso & Giordano (2011). The values of the coefficients that we obtained after the fitting are presented in Table 1.

thumbnail Fig. 2

Distribution of CHs in years 2013–2015 by area groups. Black crosses represent individual CHs. Panel a) represents CHs with an average area less than 10 000 Mm2; panel b) shows CHs between 10 000−40 000 Mm2; and panel c) represents CHs larger than 40 000 Mm2.

thumbnail Fig. 3

Distribution of the 540 selected CHs over latitude and area. Black crosses represent the individual CHs.

thumbnail Fig. 4

Distribution of investigated CHs by latitudinal zones. Panel a) represents coronal holes with average area less than 10 000 Mm2; panel b) shows CHs with area in the range 10 000−40 000 Mm2; and panel c) shows CHs larger than 40 000 Mm2.

thumbnail Fig. 5

CH sidereal rotation rates for 32 latitudinal zones 3.75° wide, and their corresponding curves fitted for each hemisphere separately. The black dot-dashed line indicates the photosphere rotation rate obtained from Doppler shift measurements (Snodgrass & Ulrich 1990). The black dashed line indicates the photosphere rotation rate obtained from the helioseismology (Schou et al. 1998). The black crosses are the rotation rates of CHs for the corresponding latitudes. The red curve corresponds to the rotation profile of CHs in the southern hemisphere, while the blue curve shows the rotation profile of CHs in the northern hemisphere. Panel a) provides information about the rotation rates of CHs with areas less than 10 000 Mm2; the corresponding error bar is 0.07°/ d (2.25 nHz); panel b) shows results for CHs with areas from 10 000 to 40 000 Mm2, error = 0.1°/ d (3.22 nHz); panel c) shows results for CHs with areas greater than 40 000 Mm2, error = 0.09°/ d (2.89 nHz); and panel d) shows the results for all CHs, error = 0.07°/ d (2.25 nHz). Standard errors are calculated as σ/N\hbox{$\sqrt N$} for each area group.

4. Discussion

4.1. Probable connection to the tachocline area

After several long-term investigations of solar photosphere rotation features, various authors have obtained different rotation profiles. We presented two illustrative examples for comparison in Fig. 5. In fact, different observations can give slightly different profiles. However, it is most likely that all of these rotation profiles are located between the two curves obtained from observations by Snodgrass & Ulrich (1990), during the years 1967–1987, and from the helioseismological observation.

We compare the average angular velocities of CHs calculated in the 3.75-degree width latitudinal zones (in total 16 latitudinal zones) with the solar interior rotation rate values in the same latitudinal zones from 0.6 to 1R. For illustration, we presented in Fig. 6 only four examples for the latitudinal zones at 0°, 20°, 40°, and 60°. Now, in Fig. 7 we show which depth of the solar interior matches the angular velocity of CHs. The asterisks in Fig. 7 represent the region in which the solar interior rotates with the same rotation rate as the CHs in the corresponding latitudinal zone. These values vary from 0.62 to 0.85 R; while on average, it is 0.71( ± 0.014) R. Figure 8 represents two different photospheric rotation profiles (dot-dashed line and dashed line), the solar interior rotation profile at 0.71 R (dotted line) and the angular velocity values of CHs combined for both hemispheres (solid line) and separately for the northern and southern hemispheres (plus and cross signs). The plots clearly show that the rotation rates of CHs directly follow the curve that corresponds to the solar internal rotation at the radius around 0.71 R. Moreover, in the vicinity of the equator the rotation rates of CHs are quantitatively similar to the photospheric rotation values obtained from the helioseismological observations; above 30° latitude the angular velocities of the photosphere decrease significantly, while the solar interior and CHs rotation profiles almost identically coincide with each other. Figure 8 shows that the interior and CHs rotate faster than the photosphere above this latitude, and this is in accordance with other measurements (Mancuso & Giordano 2011).

Table 1

Values of A, B, and C coefficients for the sidereal rotation of CHs after fitting the data on the pattern A + Bsin2(θ) + Csin4(θ).

thumbnail Fig. 6

Comparison between the CHs and solar interior rotation rates. The horizontal line represents the coronal hole sidereal rotation rate for the corresponding latitudinal zone. The curve represents the solar interior rotation profile for the same latitudinal zones. The crossing point indicates the depth where the CHs and the solar interior rotate with the same angular velocity. The vertical line corresponds to the bottom of the convective zone.

Our findings indicate that CHs could be linked to the deep layers of the Sun below the convective zone, in the tachocline, where the solar global magnetic field presumably originates. It turns out that although the rotation rates and photosphere angular velocities of CHs (helioseismological results) are quantitatively close to each other in the equatorial region, in general, the latitudinal characteristics of CHs rotation do not match any of the photospheric rotation profiles. These results may indicate that the magnetic field lines of the CHs are connected to the tachocline and the lower layers of the convection zone, approximately located at 0.71 R. We give a qualitative description of the possible physical mechanism that might provide link between CHs with deep subphotospheric layers in the conclusions.

Recently, similar results were obtained by Hiremath & Hegde (2013), who investigated CHs rotation rates based on SOHO/EIT data during the years 2001–2008. Their average angular velocity (438 nHz) corresponds to the solar interior rotation rate at 0.62( ± 0.10) solar radii. This value is slightly smaller compared to the depth obtained by us, although the main conclusion that CHs magnetic field lines are connected to the inner layers of the Sun still remains.

5. Conclusions

Our goal was to evaluate the solar CH sidereal rotation rates and to study the differences in the rotation characteristics among CHs of different sizes. To achieve these objectives, we used the CH time series with data from the SPoCA, which uses the SDO/AIA data in the 193 Å wavelength. In total, we examined 540 objects in a 120 × 120-degree area on the solar disc from 1 January 2013 to 20 April 2015.

The analysis indicates that there was a north-south asymmetry in the distribution of CHs during the years 2013–2015. Approximately 60% of the detected objects were located in the northern hemisphere. In this epoch, CHs were presented at all latitudes. However, the largest and smallest objects were found only at high latitudes (cf. Fig. 2).

The average angular velocity of the 540 investigated objects is 13.86( ± 0.05)°/ d. In the equatorial area, the rotation rate is 14.15( ± 0.07)°/ d, while in the ± 55−60° latitudinal zone it decreases to approximately 12.9( ± 0.07)°/ d. We found that the CH rotation rate is drastically different from the photospheric rotation rate.

The CH rotation profile almost perfectly coincides with the rotation profile of the tachocline and the lower layers of the convection zone at around 0.71 R. This fact gives an indication that CHs may be linked to the solar global magnetic field, which originates in the tachocline at this depth.

We examined objects only in the ± 60° latitudinal zone. Further study is needed to prove similar conclusions for latitudes above ± 60°. Also, mathematical modelling, both analytical and numerical, is desirable to get further insight into the connection of the CH roots with the deep layers of the solar interior and global solar magnetic field.

thumbnail Fig. 7

Connection between the CHs and the solar interior.

thumbnail Fig. 8

Sidereal rotation profile comparison between different layers of the Sun. The dot-dashed line indicates the photosphere rotation rate (Snodgrass & Ulrich 1990); the dashed line indicates the photosphere rotation rate according to Schou et al. (1998); the dotted line denotes the solar interior rotation rate at 0.71 R; the solid line is the rotation rate of the CHs; the plus signs are the rotation rates of the CHs in the northern hemisphere; and the cross signs are the rotation rates of the CHs in the southern hemisphere.

In this work, we report the coincidence of the CH rotation rates with those in the lower convection zone, which allows us to propose the hypothesis that CHs may be linked to the deep solar layers through the configuration of the global magnetic field. Of course, a statement of such causality requires rigorous modelling of the mentioned field configurations and related connections. Pursuing of such modelling is not the primary aim of the current work and is somewhat away from its scope. However, we can make here some indicative remarks about a possible explanation of found coincidences. In a recent study, Ruždjak et al. (2004) reported the similar analysis of rotation of sunspot groups with different age. The main conclusion of this work was that in the beginning their lifetime the young sunspot groups are plausibly anchored at higher depths (about 0.93 R), while the anchoring depth becomes more and deeper and may reach lower convection zone for those sunspot groups that are recursive. This may be related to the gradual emergence of magnetic flux from tachocline and its connection to the magnetic features formed on the solar surface. On the other hand, there is the observation of the fact that CHs may originate from coronal funnels (Tu et al. 2005) meaning that coronal holes and the fast wind patterns are tightly linked with the lower solar atmosphere and magnetic field topology there. Based on these reports, one can suppose that CHs may be formed and connected with the global residual magnetic field emerged in long-living recursive sunspot groups and thus provided a link between CHs and solar deeper layers. A recent work by McIntosh et al. (2014, 2015) has demonstrated that the sunspot cycle is driven by the rotation of the deep interior of the Sun. These authors have explained it in terms of the intra- and extra-hemispheric interaction between the overlapping activity bands of the 22-year magnetic polarity cycle, which are driven by the rotation of the deep interior of the Sun. At the same time there might be rather small CHs, which presumably belong to the group with an area less 10 000 Mm2, with short lifetimes (of the order of few days). These CHs can be linked with a residual magnetic field of weaker non-recursive sunspot groups, which are linked to the CHs anchored at higher layers. This intuitively drawn scenario is also in agreement with conclusions given by Nash et al. (1988). Possibly this can be an indirect explanation of the presence of two intersection points between the helioseismic radial rotation profile and CH mean rotation rates. The proof or denying of this stated hypothesis needs a proper analytical and/or numerical modelling, which can become the subject of future studies.


Acknowledgments

The work was supported by European FP7-PEOPLE-2010-IRSES-269299 project – SOLSPANET. The work was also supported by Shota Rustaveli National Science Foundation grant DI-2016-52. Work of S.R.B. was supported under Shota Rustaveli National Science Foundation grants for doctoral students – PhDF2016_204 and grant for young scientists for scientific research internships abroad IG/50/1/16. B.M.S. and M.L.K. acknowledge the support by the Austrian Fonds zur Foerderung der Wissenschaftlichen Forschung within the project P25640-N27 and Leverhulme Trust grant IN-2014-016. The work of T.Z. was supported by the Austrian Fonds zur Forderung der wissenschaftlichen Forschung under projects P25640-N27 and P26181-N27 and from the Georgian Shota Rustaveli National Science Foundation projects DI-2016-17 and 217146. A.K. was supported by NASA grant NNX14AB70G. M.L.K. additionally acknowledges the support of the FWF projects I2939-N27 and P25587-N27, and the grants Nos. 16-52-14006, 14-29-06036 of the Russian Fund for Basic Research. We are thankful to the anonymous referee for constructive comments that led to significant improvement of the manuscript content.

References

  1. Adams, W. M. 1976, Sol. Phys., 47, 601 [NASA ADS] [CrossRef] [Google Scholar]
  2. Altschuler, M. D., Trotter, D. E., & Orrall, F. Q. 1972, Sol. Phys., 26, 354 [Google Scholar]
  3. Aschwanden, M. J., Boerner, P., Schrijver, C. J., & Malanushenko, A. 2013, Sol. Phys., 283, 5 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bilenko, I. A., & Tavastsherna, K. S. 2016, Sol. Phys., 291, 2329 [NASA ADS] [CrossRef] [Google Scholar]
  5. Boerner, P., Warren, H. P., Testa, P., Weber, M., & Schrijver, C. J. 2011, AGU Fall Meeting Abstracts [Google Scholar]
  6. Bohlin, J. D. 1977, Sol. Phys., 51, 377 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bumba, V., Klvana, M., & Sykora, J. 1995, A&A, 298, 923 [NASA ADS] [Google Scholar]
  8. Delouille, V., Barra, V., & Hochedez, J. 2007, AGU Fall Meeting Abstracts [Google Scholar]
  9. de Toma, G. 2011, Sol. Phys., 274, 195 [NASA ADS] [CrossRef] [Google Scholar]
  10. Harvey, J., & Sheeley, Jr., N. R. 1987, BAAS, 19, 935 [NASA ADS] [Google Scholar]
  11. Hiremath, K. M., & Hegde, M. 2013, ApJ, 763, 137 [NASA ADS] [CrossRef] [Google Scholar]
  12. Insley, J. E., Moore, V., & Harrison, R. A. 1995, Sol. Phys., 160, 1 [NASA ADS] [CrossRef] [Google Scholar]
  13. Japaridze, D. R., Bagashvili, S. R., Shergelasvili, B. M., & Chargeishvili, B. B. 2015, Astrophysics, 58, 575 [NASA ADS] [CrossRef] [Google Scholar]
  14. Kirk, M. S., Pesnell, W. D., Young, C. A., & Hess Webber, S. A. 2009, Sol. Phys., 257, 99 [Google Scholar]
  15. Krieger, A. S., Timothy, A. F., & Roelof, E. C. 1973, Sol. Phys., 29, 505 [NASA ADS] [CrossRef] [Google Scholar]
  16. Krista, L. D., & Gallagher, P. T. 2009, Sol. Phys., 256, 87 [Google Scholar]
  17. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
  18. Mancuso, S., & Giordano, S. 2011, ApJ, 729, 79 [NASA ADS] [CrossRef] [Google Scholar]
  19. McIntosh, S. W., Wang, X., Leamon, R. J., et al. 2014, ApJ, 792, 12 [NASA ADS] [CrossRef] [Google Scholar]
  20. McIntosh, S. W., Leamon, R. J., Krista, L. D., et al. 2015, Nature Commun., 6, 6491 [Google Scholar]
  21. Miralles, M. P., Cranmer, S. R., Panasyuk, A. V., Romoli, M., & Kohl, J. L. 2001, ApJ, 549, L257 [NASA ADS] [CrossRef] [Google Scholar]
  22. Miralles, M. P., Cranmer, S. R., & Kohl, J. L. 2002, in From Solar Min to Max: Half a Solar Cycle with SOHO, ed. A. Wilson, ESA SP, 508, 351 [NASA ADS] [Google Scholar]
  23. Miralles, M. P., Cranmer, S. R., & Kohl, J. L. 2004, Adv. Space Res., 33, 696 [NASA ADS] [CrossRef] [Google Scholar]
  24. Munro, R. H., & Withbroe, G. L. 1972, ApJ, 176, 511 [NASA ADS] [CrossRef] [Google Scholar]
  25. Nash, A. G., Sheeley, Jr., N. R., & Wang, Y.-M. 1988, Sol. Phys., 117, 359 [NASA ADS] [CrossRef] [Google Scholar]
  26. Navarro-Peralta, P., & Sanchez-Ibarra, A. 1994, Sol. Phys., 153, 169 [NASA ADS] [CrossRef] [Google Scholar]
  27. Noci, G. 1973, Sol. Phys., 28, 403 [NASA ADS] [CrossRef] [Google Scholar]
  28. Obridko, V. N., & Shel’Ting, B. D. 1988, Byulletin Solnechnye Dannye Akademie Nauk SSSR, 1988, 89 [NASA ADS] [Google Scholar]
  29. Obridko, V. N., & Shelting, B. D. 1989, Sol. Phys., 124, 73 [NASA ADS] [CrossRef] [Google Scholar]
  30. Pneuman, G. W. 1973, Sol. Phys., 28, 247 [NASA ADS] [CrossRef] [Google Scholar]
  31. Reiss, M., Temmer, M., Rotter, T., Hofmeister, S. J., & Veronig, A. M. 2014, Central European Astrophysical Bulletin, 38, 95 [NASA ADS] [Google Scholar]
  32. Reiss, M. A., Hofmeister, S. J., De Visscher, R., et al. 2015, J. Space Weather and Space Climate, 5, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Rotter, T., Veronig, A. M., Temmer, M., & Vršnak, B. 2012, Sol. Phys., 281, 793 [NASA ADS] [CrossRef] [Google Scholar]
  34. Ruždjak, D., Ruždjak, V., Brajša, R., & Wöhl, H. 2004, Sol. Phys., 221, 225 [NASA ADS] [CrossRef] [Google Scholar]
  35. Schou, J., Antia, H. M., Basu, S., et al. 1998, ApJ, 505, 390 [NASA ADS] [CrossRef] [Google Scholar]
  36. Shelke, R. N., & Pande, M. C. 1985, Sol. Phys., 95, 193 [NASA ADS] [CrossRef] [Google Scholar]
  37. Shergelashvili, B. M., & Fichtner, H. 2012, ApJ, 752, 142 [NASA ADS] [CrossRef] [Google Scholar]
  38. Snodgrass, H. B., & Ulrich, R. K. 1990, ApJ, 351, 309 [Google Scholar]
  39. Timothy, A. F., Krieger, A. S., & Vaiana, G. S. 1975, Sol. Phys., 42, 135 [NASA ADS] [CrossRef] [Google Scholar]
  40. Tu, C.-Y., Zhou, C., Marsch, E., et al. 2005, Science, 308, 519 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  41. Verbeeck, C., Higgins, P. A., Colak, T., et al. 2013, Sol. Phys., 283, 67 [NASA ADS] [CrossRef] [Google Scholar]
  42. Verbeeck, C., Delouille, V., Mampaey, B., & De Visscher, R. 2014, A&A, 561, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Wagner, W. J. 1975, ApJ, 198, L141 [NASA ADS] [CrossRef] [Google Scholar]
  44. Wagner, W. J. 1976, in Basic Mechanisms of Solar Activity, eds. V. Bumba, & J. Kleczek, IAU Symp., 71, 41 [Google Scholar]
  45. Wilcox, J. M. 1968, Space Sci. Rev., 8, 258 [Google Scholar]
  46. Wittmann, A. D. 1996, Sol. Phys., 168, 211 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Values of A, B, and C coefficients for the sidereal rotation of CHs after fitting the data on the pattern A + Bsin2(θ) + Csin4(θ).

All Figures

thumbnail Fig. 1

Distribution of 3056 CHs with latitude and areas during the years 2013–2015. Black crosses represent the individual CHs.

In the text
thumbnail Fig. 2

Distribution of CHs in years 2013–2015 by area groups. Black crosses represent individual CHs. Panel a) represents CHs with an average area less than 10 000 Mm2; panel b) shows CHs between 10 000−40 000 Mm2; and panel c) represents CHs larger than 40 000 Mm2.

In the text
thumbnail Fig. 3

Distribution of the 540 selected CHs over latitude and area. Black crosses represent the individual CHs.

In the text
thumbnail Fig. 4

Distribution of investigated CHs by latitudinal zones. Panel a) represents coronal holes with average area less than 10 000 Mm2; panel b) shows CHs with area in the range 10 000−40 000 Mm2; and panel c) shows CHs larger than 40 000 Mm2.

In the text
thumbnail Fig. 5

CH sidereal rotation rates for 32 latitudinal zones 3.75° wide, and their corresponding curves fitted for each hemisphere separately. The black dot-dashed line indicates the photosphere rotation rate obtained from Doppler shift measurements (Snodgrass & Ulrich 1990). The black dashed line indicates the photosphere rotation rate obtained from the helioseismology (Schou et al. 1998). The black crosses are the rotation rates of CHs for the corresponding latitudes. The red curve corresponds to the rotation profile of CHs in the southern hemisphere, while the blue curve shows the rotation profile of CHs in the northern hemisphere. Panel a) provides information about the rotation rates of CHs with areas less than 10 000 Mm2; the corresponding error bar is 0.07°/ d (2.25 nHz); panel b) shows results for CHs with areas from 10 000 to 40 000 Mm2, error = 0.1°/ d (3.22 nHz); panel c) shows results for CHs with areas greater than 40 000 Mm2, error = 0.09°/ d (2.89 nHz); and panel d) shows the results for all CHs, error = 0.07°/ d (2.25 nHz). Standard errors are calculated as σ/N\hbox{$\sqrt N$} for each area group.

In the text
thumbnail Fig. 6

Comparison between the CHs and solar interior rotation rates. The horizontal line represents the coronal hole sidereal rotation rate for the corresponding latitudinal zone. The curve represents the solar interior rotation profile for the same latitudinal zones. The crossing point indicates the depth where the CHs and the solar interior rotate with the same angular velocity. The vertical line corresponds to the bottom of the convective zone.

In the text
thumbnail Fig. 7

Connection between the CHs and the solar interior.

In the text
thumbnail Fig. 8

Sidereal rotation profile comparison between different layers of the Sun. The dot-dashed line indicates the photosphere rotation rate (Snodgrass & Ulrich 1990); the dashed line indicates the photosphere rotation rate according to Schou et al. (1998); the dotted line denotes the solar interior rotation rate at 0.71 R; the solid line is the rotation rate of the CHs; the plus signs are the rotation rates of the CHs in the northern hemisphere; and the cross signs are the rotation rates of the CHs in the southern hemisphere.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.