Secular acceleration of solar rotation from 1943 to 1986

The rotation of the surface layer of the Sun is found to have been accelerated secularly from the sunspot data of 1943 to 1986. To represent the overall state of rotation of the differentially rotating Sun, we define an indexM, by integrating the angular momentum density over the whole surface of the Sun, and call it the angular momentum layer density. The indexM increased monotonically or secularly from 1943 to 1986. This period corresponds to solar cycles 18, 19, 20, and 21. The monotonic increase ofM indicates that a net angular momentum must have steadily been coming from the layer down below the surface. The differential rotation latitudinal dependence profile did not change much from cycle 18 to cycle 20, but at cycle 21 the degree of equatorial acceleration dropped. This aspect is discussed in the context of the 55-year grand cycle. Cycle 21 is the start of grand cycle VI. The latitudinal dependence is less steep at cycle 21. The time scale of secular change of the indexM reflects the time scale of change of linkage of the surface and the deep layer in form of the angular momentum transfer, and that the time scale of the profile change of the differential rotation reflects the time scale of the angular momentum transfer within the surface layer.     

The secular modulation of solar rotation from 1943 to 1992 and its time-delayed correlation with the 55-year grand cycle of the 11-year solar cycle

We devised a new method, which we call the running-segment method, to achieve high-resolution time series of indices of solar rotation for determining the latitude dependence of the differential rotation by a least-squares fitting of the daily translation of positions of sunspot groups during a fixed time segment of 11 years. The segment is moved by an amount of one year to determine the differential profile of the next point of time. Time of the determined rotation data is defined by an arithmetic mean of the beginning and ending years of the segment. The rotation underwent an acceleration from 1948 to 1974 and a deceleration from 1974 to 1987. We found that the time profile of the indexM, the angular momentum surface layer density defined by integration of the angular momentum volume density over the whole surface, follows almost exactly the time profile of the 11-year running mean of the yearly mean of the sunspot relative number with a delay time of about 20 years. The acceleration (deceleration) phase corresponds to the ascending (descending) phase of amplitude of the 11-year solar cycle of cycle 16 (19) to cycle 19 (20) with a delay time of about 20 years. The cycles 15–20 correspond to the 55-year grand cycle V of the 11-year cycle. The delay time of about 20 years agrees well with the delay time predicted by a nonlinear dynamo theory of the solar cycle for driving the 55-year modulation of the 11-year solar cycle. The agreement suggests that the Lorentz force of the magnetic field of the solar cycle during grand cycle V drives the solar rotation modulation from 1948 to 1987 and that the force needed about 20 years to modify the rotation during 1943–1992.     

Solar differential rotation derived from sunspot observations

Sunspot drawings obtained at National Astronomical Observatory of Japan during the years 1954–1986 were used to determine the differential rotation of the Sun. From the limited data set of three solar cycles it was found that three factors (the level of cycle activity, the cycle phase, and sunspot type) affect the solar rotation rate. The differential rotation varies from cycle to cycle in such a way that the rotation velocity in the low activity cycle (cycle 20) is higher than in the high-activity cycle (cycle 19). The equatorial rotation rate shows a systematic variation within each cycle. The rate is higher at the beginning of the cycle and decreases subsequently. Although quite small, the variation of solar differential rotation with respect to Zürich sunspot type was found. The H and J types show the slowest rotation among all the sunspot types.     

Differential rotation of the sun determined tracing sunspots and oscillations of sunspot tilt angle

The time and spatial characteristics of 324 large sunspots (S50 millionths of the solar hemisphere) selected from the Abastumani Astrophysical Observatory photoheliogram collection (1950–1990) have been studied. The variations of sunspot angular rotation velocity residuals and oscillations of sunspot tilt angle were analyzed. It has been shown that the differential rotation rate of selected sunspots correlates on average with the solar cycle. The deceleration of differential rotation of large sunspots begins on the ascending arm of the activity curve and ends on the descending arm reaching minimum near the epochs of solar activity maxima. This behavior disappears during the 21st cycle. The amplitudes and periods of sunspot tilt-angle oscillations correlate well with the solar activity cycle. Near the epochs of activity maximum there appear sunspots with large amplitudes and periods showing a significant scatter while the scatter near the minimum is rather low. We also found evidence of phase difference between the sunspot angular rotation velocity and the amplitudes and periods of tilt-angle oscillations.     

On the solar rotation and activity

The interaction between differential rotation and magnetic fields in the solar convection zone was recently modelled by Brun (2004). One consequence of that model is that the Maxwell stresses can oppose the Reynolds stresses, and thus contribute to the transport of the angular momentum towards the solar poles, leading to a reduced differential rotation. So, when magnetic fields are weaker, a more pronounced differential rotation can be expected, yielding a higher rotation velocity at low latitudes taken on the average. This hypothesis is consistent with the behaviour of the solar rotation during the Maunder minimum. In this work we search for similar signatures of the relationship between the solar activity and rotation determined tracing sunspot groups and coronal bright points. We use the extended Greenwich data set (1878–1981) and a series of full-disc solar images taken at 28.4 nm with the EIT instrument on the SOHO spacecraft (1998–2000). We investigate the dependence of the solar rotation on the solar activity (described by the relative sunspot number) and the interplanetary magnetic field (calculated from the interdiurnal variability index). Possible rotational signatures of two weak solar activity cycles at the beginning of the 20th century (Gleissberg minimum) are discussed. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Meridional Motions and Reynolds Stress Determined by Using Kanzelhöhe Drawings and White Light Solar Images from 1964 to 2016

Sunspot position data obtained from Kanzelhöhe Observatory for Solar and Environmental Research (KSO) sunspot drawings and white light images in the period 1964 to 2016 were used to calculate the rotational and meridional velocities of the solar plasma. Velocities were calculated from daily shifts of sunspot groups and an iterative process of calculation of the differential rotation profiles was used to discard outliers. We found a differential rotation profile and meridional motions in agreement with previous studies using sunspots as tracers and conclude that the quality of the KSO data is appropriate for analysis of solar velocity patterns. By analyzing the correlation and covariance of meridional velocities and rotation rate residuals we found that the angular momentum is transported towards the solar equator. The magnitude and latitudinal dependence of the horizontal component of the Reynolds stress tensor calculated is sufficient to maintain the observed solar differential rotation profile. Therefore, our results confirm that the Reynolds stress is the dominant mechanism responsible for transport of angular momentum towards the solar equator.     

Internal rotation of the Sun as inferred from GONG observations

Helioseismology is a direct and most informative method of studying the structure and dynamics of the Sun. Determining the internal differential rotation of the Sun requires that the frequencies of its eigentones be estimated with a high accuracy, which is possible only on the basis of continuous long-term observations. The longest quasi-continuous series of data have been obtained by the Global Oscillation Network Group (GONG). The parameters of each individual mode of solar acoustic oscillations with low spherical degrees l=0, 1, 2, 3, 4, 5, 6 are determined by using 1260-day-long series of GONG observations. The mean frequency splitting by rotation for the modes of each radial order n is calculated by using all possible combinations between the eigenfrequencies in multiplets. As a result, it has become possible to statistically estimate the splitting and its measurement errors for the modes of each radial order. The mean splitting for each given degree l=1–6 is presented under the assumption of its independence of oscillation frequency, which holds for the achieved accuracy. The frequencies and splittings for the modes with low spherical degrees l, together with the MDI group results for higher degrees l, are used to invert the radial profile of solar angular velocity. Using the SOLA method to solve the inverse problem of restoring the rotation profile has yielded solutions sensitive to the deepest stellar interiors. Our results indicate that the solar core rotates faster than the surface, and there may be a local minimum in angular velocity at its boundary.     

A large-eddy simulation of turbulent compressible convection: differential rotation in the solar convection zone

We present the results of two simulations of the convection zone, obtained by solving the full hydrodynamic equations in a section of a spherical shell. The first simulation has cylindrical rotation contours (parallel to the rotation axis) and a strong meridional circulation, which traverses the entire depth. The second simulation has isorotation contours about mid-way between cylinders and cones, and a weak meridional circulation, concentrated in the uppermost part of the shell.
We show that the solar differential rotation is directly related to a latitudinal entropy gradient, which pervades into the deep layers of the convection zone. We also offer an explanation of the angular velocity shear found at low latitudes near the top. A non-zero correlation between radial and zonal velocity fluctuations produces a significant Reynolds stress in that region. This constitutes a net transport of angular momentum inwards, which causes a slight modification of the overall structure of the differential rotation near the top. In essence, the thermodynamics controls the dynamics through the Taylor–Proudman momentum balance . The Reynolds stresses only become significant in the surface layers, where they generate a weak meridional circulation and an angular velocity 'bump'.     

22-Year periodicity in the solar differential rotation

Using the data on sunspot groups compiled during 1879–1975, we determined variations in the differential rotation coefficientsA andB during the solar cycle. The variation in the equatorial rotation rateA is found to be significant only in the odd numbered cycles, with an amplitude ∼ 0.01 μ rads^-1. There exists a good anticorrelation between the variations of the differential rotation rateB derived from the odd and even numbered cycles, suggesting existence of a ‘22-year’ periodicity inB. The amplitude of the variation ofB is ∼ 0.05 μ rad s^-1.     

Meridional Motion and Reynolds Stress from Debrecen Photoheliographic Data

The Debrecen Photoheliographic Data catalogue is a continuation of the Greenwich Photoheliographic Results providing daily positions of sunspots and sunspot groups. We analyse the data for sunspot groups focussing on meridional motions and transfer of angular momentum towards the solar equator. Velocities are calculated with a daily shift method including an automatic iterative process of removing the outliers. Apart from the standard differential rotation profile, we find meridional motion directed towards the zone of solar activity. The difference in measured meridional flow in comparison to Doppler measurements and some other tracer measurements is interpreted as a consequence of different flow patterns inside and outside of active regions. We also find a statistically significant dependence of meridional motion on rotation velocity residuals confirming the transfer of angular momentum towards the equator. Analysis of horizontal Reynolds stress reveals that the transfer of angular momentum is stronger with increasing latitude up to about \(40^{\circ}\), where there is a possible maximum in absolute value.     

Chromospheric activity of late-type giants and supergiants: Reappearance of dynamo activity in the interior due to the spin-up of the core in evolution

The renewal of chromospheric activity in red giants and supergiants is interpreted in terms of the reappearance of dynamo activity in the interior due to the spin-up of the core caused by its contraction in the course of evolution from the main sequence to the giant stage. A region of very high rotational shear (differential rotation) develops between the core, which spins up by a large factor through the drastic contraction, and the envelope, which spins down in contrast by virtue of expansion. Mechanisms of angular momentum transfer may operate to smear this large shear, and bring the inner part of the envelope into sheared rotation. A convective layer, on the other hand, develops in the envelope from the surface inwards, when the envelope expands and the temperature is lowered. A dynamo layer, or a layer in which the sheared rotation co-exists with the convection (the presence of a remnant magnetic field being postulated), will thus reappear in the inner part of the envelope when the envelope-convection reaches down and invades the layer of sheared rotation. Surface chromospheric activity due to the magnetic field is thus renewed when the regenerated magnetic field is brought up to the surface by the envelope-convection. These phenomena occur as the star evolves into the giant stage and hence explain the observed characteristic of gradual revival of chromospheric activity from the subgiant to the giant stage. Visiting Scientist at the Tokyo Astronomical Observatory under the support of the Japan Society for the Promotion of Sciences. On leave from the Indian Institute of Astrophysics. Deceased on 1982 August 19     

The subsurface radial gradient of solar angular velocity from MDI f-mode observations

We report quantitative analysis of the radial gradient of solar angular velocity at depths down to about 15 Mm below the solar surface for latitudes up to 75° using the Michelson Doppler Imager (MDI) observations of surface gravity waves (fmodes) from the Solar and Heliospheric Observatory (SOHO). A negative outward gradient of around –400 nHz/R , equivalent to a logarithmic gradient of the rotation frequency with respect to radius which is very close to –1, is found to be remarkably constant between the equator and 30° latitude. Above 30° it decreases in absolute magnitude to a very small value at around 50°. At higher latitudes the gradient may reverse its sign: if so, this reversal takes place in a thin layer extending only 5 Mm beneath the visible surface, as evidenced by the most superficial modes (with degrees l>250). The signature of the torsional oscillations is seen in this layer, but no other significant temporal variations of the gradient and value of the rotation rate there are found.     

Long-Term Variations in Solar Differential Rotation and Sunspot Activity

The solar equatorial rotation rate, determined from sunspot group data during the period 1879–2004, decreased over the last century, whereas the level of activity has increased considerably. The latitude gradient term of the solar rotation shows a significant modulation of about 79 year, which is consistent with what is expected for the existence of the Gleissberg cycle. Our analysis indicates that the level of activity will remain almost the same as the present cycle during the next few solar cycles (i.e., during the current double Hale cycle), while the length of the next double Hale cycle in sunspot activity is predicted to be longer than the current one. We find evidence for the existence of a weak linear relationship between the equatorial rotation rate and the length of sunspot cycle. Finally, we find that the length of the current cycle will be as short as that of cycle 22, indicating that the present Hale cycle may be a combination of two shorter cycles. Presently working for the Mt. Wilson Solar Archive Digitization Project at UCLA.     

The Mechanism involved in the Reversals of the Sun's Polar Magnetic Fields

Models of the polarity reversals of the Sun's polar magnetic fields based on the surface transport of flux are discussed and are tested using observations of the polar fields during Cycle 23 obtained by the National Solar Observatory at Kitt Peak. We have extended earlier measurements of the net radial flux polewards of ±60° and confirm that, despite fluctuations of 20%, there is a steady decline in the old polarity polar flux which begins shortly after sunspot minimum (although not at the same time in each hemisphere), crosses the zero level near sunspot maximum, and increases, with reversed polarity during the remainder of the cycle. We have also measured the net transport of the radial field by both meridional flow and diffusion across several latitude zones at various phases of the Cycle. We can confirm that there was a net transport of leader flux across the solar equator during Cycle 23 and have used statistical tests to show that it began during the rising phase of this cycle rather than after sunspot maximum. This may explain the early decrease of the mean polar flux after sunspot minimum. We also found an outward flow of net flux across latitudes ±60° which is consistent with the onset of the decline of the old polarity flux. Thus the polar polarity reversals during Cycle 23 are not inconsistent with the surface flux-transport models but the large empirical values required for the magnetic diffusivity require further investigation.     

Chromospheric rotation during 1972–73, years of declining activity

The rotational behaviour of the chromosphere, observed in the Ca ii K₃ line, is studied during 1972–1973, years of decreasing solar activity. Daily chromospheric filtergrams, detected at the Anacapri Observatory, are digitized by means of a flying-spot photometer, controlled by computer. The time series of the daily chromospheric data detected at central meridian, relative to 30 consecutive latitude zones, are analyzed to determine the recurrence tendency due to the rotation of long-lived chromospheric features. The computed rotation rate is independent of latitude, in agreement with the results obtained for the green corona during the years before sunspot minimum. Namely both chromospheric and coronal features, with lifetime exceeding one solar rotation, are almost not affected by differential rotation before sunspot minimum.     

An interpretation of the broad-band circular polarization of sunspots

The broad-band circular polarization of sunspots is discussed on the basis of the observations made in the Okayama Astrophysical Observatory. The observation with the spectrograph proves that it is the integrated polarization of spectral lines in the observed spectral range. A velocity gradient in the line-of-sight can produce this integrated polarization due to the differential saturation between Zeeman components of magnetically sensitive lines. The observed degree of polarization and its spatial distribution in sunspots is explained when we introduce a differentially twisted magnetic field in addition to the velocity gradient. The differential twist has the azimuth rotation of the magnetic field along the line-of-sight and generates the circular polarization from the linear polarization due to the magneto-optical effect. The required azimuth rotation is reasonable and amounts at most to 30°. The required velocity gradient is compatible with the line asymmetry and its spatial distribution observed in sunspots. The observed polarity rule leads to the conclusion that the sunspot magnetic field has the differential twist with the right-handed azimuth rotation relative to the direction of the main magnetic field, without regard to the magnetic polarity and to the solar cycle. The twist itself is left-handed under the photosphere, when the sunspot is assumed to be a unwinding emerging magnetic field.     

The peculiarities of rotation of the solar polar regions

The sidereal rotation rate of the high-latitude solar regions is examined using long-lived photospheric polar faculae. The observations were carried out with the photoheliograph of Kislovodsk Mountain Station of the Pulkovo Observatory from 1982 to 1986. The following facts have been established: (a) There is a differential rotation of the polar faculae close to the maximum of solar activity, while the amount of latitude gradient of solar rotation decreases towards the sunspot minimum; (b) small differences of rotation in the northern and southern hemispheres of the Sun are observed; (c) some deviations of differential rotation curves constructed for each Carrington rotation from the mean curve of differential rotation are revealed. The total amplitude of the maximum positive and negative excesses is about 40–50 m s^–1. The positive surplus velocities of solar rotation (the amplitude of which is about 20–25 m s^–1) move in the form of a wave from heliographic latitudes 40° with a velocity of 1.6 m s^–1. The latitude width of this flow is B 15°. This wave of abnormally high velocity starts in the year of minimum solar activity and reaches the pole 11 years later. The picture is symmetrical relative to the equator.     

Sunspot Catalogue of the Observatory of the University of Coimbra (1929 – 1941)

A sunspot catalogue was published by the Coimbra Astronomical Observatory (Portugal), which is now called the Geophysical and Astronomical Observatory of the University of Coimbra, for the period 1929?–?1941. We digitalised data included in that catalogue and provide a machine-readable version. We show the reconstructions for the (total and hemispheric) sunspot number index and sunspot area according to this catalogue and compare it with the sunspot number index (version 2) and the Balmaceda sunspot area series (Balmaceda et al. in J. Geophys. Res.114, A07104, 2009). Moreover, we also compared the Coimbra catalogue with records made at the Royal Greenwich Observatory. The results demonstrate that the historical catalogue compiled by the Coimbra Astronomical Observatory contains reliable sunspot data and can therefore be considered for studies about solar activity.     

The Relative Phase Asynchronization between Sunspot Numbers and Polar Faculae

The monthly sunspot numbers compiled by Temmer et al. and the monthly polar faculae from observations of the National Astronomical Observatory of Japan, for the interval of March 1954 to March 1996, are used to investigate the phase relationship between polar faculae and sunspot activity for total solar disk and for both hemispheres in solar cycles 19, 20, 21 and 22. We found that (1) the polar faculae begin earlier than sunspot activity, and the phase difference exhibits a consistent behaviour for different hemispheres in each of the solar cycles, implying that this phenomenon should not be regarded as a stochastic fluctuation; (2) the inverse correlation between polar faculae and sunspot numbers is not only a long-term behaviour, but also exists in short time range; (3) the polar faculae show leads of about 50–71 months relative to sunspot numbers, and the phase difference between them varies with solar cycle; (4) the phase difference value in the northern hemisphere differs from that in the southern hemisphere in a solar cycle, which means that phase difference also existed between the two hemispheres. Moreover, the phase difference between the two hemispheres exhibits a periodical behaviour. Our results seem to support the finding of Hiremath (2010).     

Differential rotation and meridional flow of Arcturus

The spectroscopic variability of Arcturus hints at cyclic activity cycle and differential rotation. This could provide a test of current theoretical models of solar and stellar dynamos. To examine the applicability of current models of the flux transport dynamo to Arcturus, we compute a mean‐field model for its internal rotation, meridional flow, and convective heat transport in the convective envelope. We then compare the conditions for dynamo action with those on the Sun. We find solar‐type surface rotation with about 1/10th of the shear found on the solar surface. The rotation rate increases monotonically with depth at all latitudes throughout the whole convection zone. In the lower part of the convection zone the horizontal shear vanishes and there is a strong radial gradient. The surface meridional flow has maximum speed of 170 m/s and is directed towards the equator at high and towards the poles at low latitudes. Turbulent magnetic diffusivity is of the order 10¹⁵–10¹⁶ cm²/s. The conditions on Arcturus are not favorable for a circulation‐dominated dynamo (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)