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.     

Preferred longitudes of sunspot groups and high-speed solar wind streams: Evidence for a ‘solar memory’

Correlation analysis of the mean longitude distribution of sunspot groups (taken from the Greenwich Photoheliographic Results) and high-speed solar wind streams (inferred from the C9 index for geomagnetic disturbances) with the Bartels rotation period P = 27.0 days shows anti-correlation for individual cycles.In particular, the longitudes of post-maximum stable streams of cycle 18 and 19 are well anticorrelated with the preferred longitudes of sunspot groups during the maximum activity periods of these cycles. This is further analyzed using the daily Zürich sunspot number, R, between 1932 and 1980, which reveals a conspicuous similarity of cycle 18 and 19 as well as cycle 20 and 21.We conclude that there is a solar memory for preferred longitudes of activity extending at least over one, probably two cycles (i.e. one magnetic cycle of 22 years). We conjecture that this memory extends over longer intervals of time as a long-term feature of solar activity.     

A Relationship Between the Solar Rotation and Activity Analysed by Tracing Sunspot Groups

The sunspot position published in the data bases of the Greenwich Photoheliographic Results (GPR), the US Air Force Solar Optical Observing Network and National Oceanic and Atmospheric Administration (USAF/NOAA), and of the Debrecen Photoheliographic Data (DPD) in the period 1874 to 2016 were used to calculate yearly values of the solar differential-rotation parameters \(A\) and \(B\). These differential-rotation parameters were compared with the solar-activity level. We found that the Sun rotates more differentially at the minimum than at the maximum of activity during the epoch 1977?–?2016. An inverse correlation between equatorial rotation and solar activity was found using the recently revised sunspot number. The secular decrease of the equatorial rotation rate that accompanies the increase in activity stopped in the last part of the twentieth century. It was noted that when a significant peak in equatorial rotation velocity is observed during activity minimum, the next maximum is weaker than the previous one.     

Solar-Cycle-Related Variations in the Solar Differential Rotation and Meridional Flow: A Comparison

We have analysed a large set of sunspot group data (1874 – 2004) and find that the meridional flow strongly varies with the phase of the solar cycle, and the variation is quite different in the northern and the southern hemispheres. We also find the existence of considerable cycle-to-cycle variation in the meridional velocity, and about a 11-year difference between the phases of the corresponding variations in the northern and the southern hemispheres. In addition, our analysis also indicates the following: (i) the existence of a considerable difference (about 180^°) between the phases of the solar-cycle variations in the latitude-gradient terms of the northern and the southern hemispheres’ rotations; (ii) the existence of correlation (good in the northern hemisphere and weak in the southern hemisphere) between the mean solar-cycle variations of meridional flow and the latitude-gradient term of solar rotation; (iii) in the northern hemisphere, the cycle-to-cycle variation of the mean meridional velocity leads that of the equatorial rotation rate by about 11 years, and the corresponding variations have approximately the same phase in the southern hemisphere; and (iv) the directions of the mean meridional velocity is largely toward the pole in the longer sunspot cycles and largely toward the equator in the shorter cycles.     

Differential rotation and giant cell circulation of solar Ca+-network

We report new results obtained from high precision computer controlled tracings of ca. 400 bright Ca⁺-mottles made during summer 1975 in continuation of our 1974 program (Schröter and Wöhl, 1975). In particular, we looked in 1975 for the existence of a giant circulation pattern in the equatorial zone. We find for the differential rotation: = 13.93 – 2.90 sin² B (deg/day, sidereal) when combining the new measurements with those obtained in 1974. Observations from 26th April until June 19th give strong evidence that at that time four giant circulation cells, crossing the solar equator, (i.e. a nonaxisymmetric velocity field pattern with respect to the solar equator) did exist. This yields two more rapid and two slower rotating sectors with v = ±80 m s^–1. These giant cells transport angular momentum towards the equator.     

The 100—year periodic modulation of solar rotation

A periodic long-term modulation of the solar surface rotation with a time scale on the order of 100 years is found in the sunspot data from 1874 to 1992 obtained by combinig the Greenwich Photoheliographic Results from cycle 11 to cycle 20 analysed by Balthasar, Vázquez, and Wöhl and the Mitaka sunspot sketch data from cycle 18 to 22 of the National Astronomical Observatory of Japan which was the Tokyo Astronomical Observatory of the University of Tokyo until 1988. A new index of the solar rotation M defined by integrating the angular momentum density over the whole surface, which we call the angular momentum surface layer density, reached a maximum at solar cycle 14, decreased to a minimum at cycle 17, and then increased to reach another maximum at cycle 21. The increase of M means acceleration of the surface layer as a whole by transport of angular momentum from the deeper layer. This implies an decrease (increase) of the radial gradient of the differential rotation if the basic radial gradient of the differential rotation increases (decreaes) inward. The decrease of M means deceleration of the surface layer and implies an increase (decrease) of the radial gradient. The degree of the equatorial acceleration of the surface differential rotation is also found to have undergone the same 100 year periodic modulation during the same interval, reaching a minimum at cycle 14, a maximum at cycle 17, and a minimum at cycle 21 in antiphase with the modulation of M. Thus both radial and latitudinal gradients of the differential rotation increased and decreased in phase (in anti-phase) if the basic radial gradient increases (decreases) inward.     

Short-Term Variations in the Equatorial Rotation Rate of Sunspot Groups

We have detected several periodicities in the solar equatorial rotation rate of sunspot groups in the catalog Greenwich Photoheliographic Results (GPR) during the period 1931?–?1976, the Solar Optical Observing Network (SOON) during the period 1977?–?2014, and the Debrecen Photoheliographic Data (DPD) during the period 1974?–?2014. We have compared the results from the fast Fourier transform (FFT), the maximum entropy method (MEM), and the Morlet wavelet power-spectra of the equatorial rotation rates determined from SOON and DPD sunspot-group data during the period 1986?–?2007 with those of the Mount Wilson Doppler-velocity data during the same period determined by Javaraiah et al. (Solar Phys. 257, 61, 2009). We have also compared the power-spectra computed from the DPD and the combined GPR and SOON sunspot-group data during the period 1974?–?2014 to those from the GPR sunspot-group data during the period 1931?–?1973. Our results suggest a ～?250-day period in the equatorial rotation rate determined from both the Mt. Wilson Doppler-velocity data and the sunspot-group data during 1986?–?2007. However, a wavelet analysis reveals that this periodicity appears mostly around 1991 in the velocity data, while it is present in most of the solar cycles covered by the sunspot-group data, mainly near the minimum epochs of the solar cycles. We also found the signature of a period of ～?1.4 years in the velocity data during 1990?–?1995, and in the equatorial rotation rate of sunspot groups mostly around the year 1956. The equatorial rotation rate of sunspot groups reveals a strong ～?1.6-year periodicity around 1933 and 1955, a weaker one around 1976, and a strong ～?1.8-year periodicity around 1943. Our analysis also suggests periodicities of ～?5 years, ～?7 years, and ～?17 years, as well as some other short-term periodicities. However, short-term periodicities are mostly present at the time of solar minima. Hence, short-term periodicities cannot be confirmed because of the larger uncertainty in the data.     

Stable Sunspot Area Level of Debrecen Photoheliographic Data and Multivariate Correction Factor of SOON Data

We investigate the spatial and temporal variation of sunspot group areas reported by the Greenwich Photoheliographic Results (GPR), the Solar Optical Observing Network (SOON), the Kislovodsk Mountain Astronomical Station (KMAS), and the Debrecen Photoheliographic Data (DPD) databases. We identify improved correction factors for reconciling these individual records to a common scale. Our results show that the DPD sunspot group areas are stable over the studied interval (1974?–?2014). We find an improved fit between GPR and DPD sunspot group areas when using a correction factor such that \(\mathrm{GPR} = 0.975(\pm 0.006) \times \mathrm{DPD}\), independent of the position of the sunspot group on the solar disk. We also find that the scale of KMAS sunspot group areas fits that of DPD well, but has a small position-dependent trend near the limb. However, in order to set SOON sunspot group area records onto the scale of DPD, we find that there is a need for a multivariate correction factor. This multivariate correction factor has a value ranging between 1.1 and 1.9 and is dependent upon the time of the SOON observation, the distance of the group from disk center, and the observatory within the SOON network. Finally, we provide further context to the systematic bias in SOON sunspot group area observations toward lower values relative to those recorded in the GPR and DPD databases that has previously been reported in the literature. We have identified the two main contributors to the SOON area deficit; some penumbral parts are unobserved, and the spot areas are underestimated. Our analysis is vital for studies that require stable, long-term solar activity records such as solar irradiance models that estimate irradiance reduction from records of sunspot group numbers, areas, and locations.     

On the large-scale diffuse magnetic field of the sun

Although the sunspots migrate towards the equator, the large-scale weak diffuse magnetic fields of the Sun migrate poleward with the solar cycle, the polar field reversing at the time of the sunspot maxima. We apply the vector model of Dikpati and Choudhuri (1994, Paper I) to fit these observations. The dynamo layer at the base of the convection zone is taken to be the source of the diffuse field, which is then evolved in the convection zone subject to meridional circulation and turbulent diffusion. We find that the longitudinally averaged observational data can be fitted reasonably well both for positive and negative values of the-effect by adjusting the subsurface meridional flow suitably. The model will be extended in a future paper to include the decay of active regions as an extra source of the diffuse field, which may be necessary to explain the probable phase lag betweenB _r andB at lower latitudes.     

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)     

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.     

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.     

The Rotation of Sunspots During Activity Cycle 22

The solar rotation rate derived from the measurements of daily sunspot positions during activity cycle 22 (1988–1996) are presented. The averaged differential rotation is nearly symmetric with respect to the equator. The comparisons of our measurement with those for previous cycles show that the equatorial rotation rate gradually increases from cycle 18 to 22.     

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)     

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.     

Increasing Lifetime of Recurrent Sunspot Groups Within the Greenwich Photoheliographic Results

Long-lived (>20 days) sunspot groups extracted from the Greenwich Photoheliographic Results (GPR) are examined for evidence of decadal change. The problem of identifying sunspot groups that are observed on consecutive solar rotations (recurrent sunspot groups) is tackled by first constructing manually an example dataset of recurrent sunspot groups and then using machine learning to generalise this subset to the whole GPR. The resulting dataset of recurrent sunspot groups is verified against previous work by A. Maunder and other Royal Greenwich Observatory (RGO) compilers. Recurrent groups are found to exhibit a slightly larger value for the Gnevyshev?–?Waldmeier Relationship than the value found by Petrovay and van Driel-Gesztelyi (Solar Phys. 51, 25, 1977), who used recurrence data from the Debrecen Photoheliographic Results. Evidence for sunspot-group lifetime change over the previous century is observed within recurrent groups. A lifetime increase of a factor of 1.4 between 1915 and 1940 is found, which closely agrees with results from Blanter et al. (Solar Phys. 237, 329, 2006). Furthermore, this increase is found to exist over a longer period (1915 to 1950) than previously thought and provisional evidence is found for a decline between 1950 and 1965. Possible applications of machine-learning procedures to the analysis of historical sunspot observations, the determination of the magnetic topology of the solar corona and the incidence of severe space–weather events are outlined briefly.     

Sunspot nests

For the period August 1959–December 1964 the Greenwich Photoheliographic Results were searched for sunspot nests. Such a nest is a sequence of sunspot groups that appear within a small area on the solar surface and that last for several months. The search procedure is described and data for 41 probable sunspot nests are given. At least three quarters of these nests appear to be real, and not chance clusters.The nests are the same type of activity sequences as the Fleckenherde discovered by Becker (1955) and the complexes of activity pointed out by Gaizauskas et al. (1983). The complexes of activity as defined by Bumba and Howard (1965) are different patterns, however; the relation between complexes and the nests is shown.Some properties of the nests are: (i) many nests appear as double structures; (ii) single nests and components of double nests are quite compact: the effective areas are comparable to those of medium-large sunspot groups; (iii) each nest rotates at its own steady rate about the Sun; (iv) the intrinsic scatter in the rotation rates is much larger than the trend in the differential rotation; (v) displacements in latitude are less than a few meters per second; (vi) many nests live for 6 to 15 Carrington rotation periods, the minimum lifetime is not yet determined; (vii) the fraction of the sunspot groups that are members of nests is large (at least 30%).     

The covariance of latitudinal and longitudinal motions of small magnetic features

We study the covariance of longitudinal and latitudinal motions of small magnetic features after subtracting long-term averages of differential rotation and meridional flow. The covariance is generally interpreted as Reynolds stress and linked to the equatorward transport of angular momentum. Using high-resolution magnetograms taken daily with the NSO Vacuum Telescope on Kitt Peak, we determine large-scale motions by a two-dimensional crosscorrelation analysis of pairs of consecutive daily observations from which active regions are excluded, i.e., we analyze the motions of small magnetic features. In the present work, we focus on 107 day pairs obtained during the year 1988 and on 472 day pairs taken in selected intervals from 1978 to 1990. We find that all covariance values are very small (below 250 m² s⁻²), which is about one to two orders of magnitude smaller than the values from sunspot measurements derived by other authors. At active region latitudes, the masking process increases the noise, which increases the chance that the covariances at these latitudes are not significantly different from zero. We find that the results depend strongly on the temporal averaging involved. Daily unaveraged crosscorrelations lead to no apparent correlation between the residual velocities, while in the monthly averages of the 1988 data, we find a covariance of −37 ± 15 m² s⁻² at 45° with a linear correlation of −0.59, which is significantly different from zero and has the right sign for an equatorial transport of angular momentum. When we average over longer time periods, the covariance values decrease again. The annual averages of the 1978–1990 data show both no significant covariances and the smallest errors. These small covariances imply that the motions of small magnetic features do not reflect the transport of angular momentum via the mechanism of Reynolds stress. Operated by the Association of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science Foundation.     

Meridional motions of sunspot groups during eleven activity cycles

Greenwich data (1874–1976) are used for a time-dependent analysis of meridional motions of sunspot groups. We obtain the latitude-dependence of meridional motions of sunspot groups with respect to a mean latitude determined for half-year intervals. The daily meridional motions of groups are also given separately for growing and decaying sunspot groups. The development is determined from changes of sunspot areas. Our results are compared with the reductions performed by Howard (1991b) using the Mt. Wilson sunspot data from 1917 until 1985: Although we have smaller errors, we do not find any significant drift. We also do not find different trends in the meridional motions of growing as compared to decreasing sunspots.