Sunspot and stray light observations during the 1971, February 25 partial solar eclipse

Observations of a sunspot during and after a partial solar eclipse are described. The amount of scattered light confirms the existence of a spread function component with a half width of 10″. The observations also indicate the possibility of severely underestimating this component by aureole measurements. Umbral continuum intensities of 0.10 I _⊙ in the red spectral region were directly measured, the correction for scattered light amounts to 0.02 I _⊙. Intensities calculated with four umbral models are larger than the observed values, indicating this sunspot to be cooler by some 100 K. The wings of two strong Ca i lines are equally explained by the models of Henoux, Kneer, and Stellmacher/Wiehr. Yun's model can be ruled out because of too high a temperature.     

Basic properties of sunspots: equilibrium,stability and long-term eigen oscillations

A number of fundamental questions as regards the physical nature of sunspots are formulated. In order to answer these questions, we apply the model of a round-shaped unipolar sunspot with a lower boundary consisting of cool plasma and with strong magnetic field at the depth of about 4 Mm beneath the photosphere, in accordance with the data of local helioseismology and with certain physically sound arguments (the shallow sunspot model). The magnetic configuration of a sunspot is assumed to be close to the observed one and similar to the magnetic field of a round solenoid of the appropriate size. The transverse (horizontal) and longitudinal (vertical) equilibria of a sunspot were calculated based on the thermodynamic approach and taking into account the magnetic, gravitational, and thermal energy of the spot and the pressure of the environment. The dependence of the magnetic field strength in the sunspot center, B ₀, on the radius of the sunspot umbra a is derived theoretically for the first time in the history of sunspot studies. It shows that the magnetic field strength in small spots is about 700 Gauss (G) and then increases monotonically with a, tending asymptotically to a limit value of about 4000 G. This dependence, B ₀(a) includes, as parameters, the gravity acceleration on the solar surface, the density of gas in the photosphere, and the ratio of the radius of the spot (including penumbra), a _p, to the radius of its umbra a. It is shown that large-scale subsurface flows of gas in the sunspot vicinity, being the consequence but not the cause of sunspot formation, are too weak to contribute significantly to the pressure balance of the sunspot. Stability of the sunspot is provided by cooling of the sunspot plasma and decreasing of its gravitational energy due to the vertical redistribution of the gas density when the geometric Wilson depression of the sunspot is formed. The depth of a depression grows linearly with B ₀, in contrast to the quadratic law for the magnetic energy. Therefore, the range of stable equilibria turns out to be limited: large spots, with radius a larger than some limit value (about 12–18 Mm, depending on the magnetic field configuration), are unstable. It explains the absence of very large spots on the Sun and the appearance of light bridges in big spots that divide the spot into a few parts. The sunspots with B ₀≈2.6÷2.7 kilogauss (kG) and a≈5 Mm are most stable. For these spots, taken as a single magnetic structure, the period of their vertical eigen oscillations is minimal and amounts, according to the model, to 10–12 hours. It corresponds well to the period derived from the study of long-term oscillations of sunspots using SOHO/MDI data.     

Sunspot models with twisted magnetic field

The return-flux sunspot model is generalized by including azimuthal magnetic field B _Φ . The basic equation is obtained and numerical solutions are compared with the analogous solutions for the Schlüter-Temesvary sunspot theory for two cases: B _Φ ～ B _r and B _Φ ～ rB _r . The solutions demonstrate that the twisting of the sunspot magnetic field decreases with height. Our models confirm Yun's early statement: the azimuthal field only slightly influences sunspot structure.     

The Relative Umbral Area in Spot Groups as an Index of Cyclic Variation of Solar Activity

The Greenwich series of data was used to study the ratio [q] of the total umbra area to the total area of the sunspot group (for brevity “relative umbral area”) for the period 1874?–?1976. It was revealed that the annual mean value of q varied in time from 0.15 to 0.28 and reached its maximum in the early 1930s. The dependence of q on the sunspot group area [S] was considered to show that the smallest groups, of area less than 100 m.v.h. (millionths of the visible hemisphere), contributed most significantly to the temporal variation of q. In contrast to the earlier results, the dependence obtained proved to be rather complicated. The coefficients of the linear expansion q(S) are themselves dependent on the sunspot-group area and time [t]; i.e. the relation of q to both S and t is nonlinear. Only in sunspot groups with a large area does dependence disappear, and q becomes constant, equal to 0.18. This is the value given in textbooks. The relations obtained show that the relive umbral area and the relative number of small groups are important parameters of the secular variation of solar activity. In particular, they may account for variations in the mean magnetic field in active regions, the complexity of a group according to the magnetic classification, the flare activity of a sunspot group, and its geophysical impact. It is conjectured that the parameter q describes the time-varying relative contribution from the interior and subsurface dynamo mechanisms.     

Comparison of three solar activity indices based on sunspot observations

The following sunspot formation indices are analyzed: the relative sunspot number R _z, the normalized sunspot group number R _g, and the total sunspot area A. Six empirical formulas are derived to describe the relations among these indices after 1908. The earlier data exhibit systematic deviations from these formulas, which can be attributed to systematic errors of the indices. The Greenwich data on the sunspot total area A and the sunspot group number in 1874–1880 are found to be doubtful. Erroneous data at the beginning of the Greenwich series must spoil the values of the index R _g in the XVII–XIX centuries. The Hoyt-Schatten series of R _g may be less reliable than the well-known Wolf number series R _z.     

Study of Differences Between Sunspot and White Light Facular Area Data Determined from SDO/HMI and SOHO/MDI Observations

Sunspot and white light facular areas are important data for solar activity and are used, for example, in the study of the evolution of sunspots and their effect on solar irradiance. Solar Dynamic Observatory??s Helioseismic and Magnetic Imager (SDO/HMI) solar images have much higher resolution (??0.5????pixel^?1) than Solar and Heliospheric Observatory??s Michelson Doppler Imager (SOHO/MDI) solar images (??2????pixel^?1). This difference in image resolution has a significant impact on the sunspot and white light facular areas measured in the two datasets. We compare the area of sunspots and white light faculae derived from SDO/HMI and SOHO/MDI observations. This comparison helps the calibration of the SOHO sunspot and facular area to those in SDO observations. We also find a 0.22 degree difference between the North direction in SDO/HMI and SOHO/MDI images.     

SID flares and sunspot morphology

The degree of association between geoeffective (SID producing) flares (hereafter called SID flares) and sunspot morphology is examined. It is found that: (1) the frequency of SID flares associated with sunspot groups is linear function of sunspot area and rate of change in area; (2) the SID flare intensity is dependent on the sunspot area and on the magnetic morphology (field geometry); (3) the probability of a sunspot group being magnetically complex (henceforth called complex ratio) is a linear function of spot area, the larger this area the more likely a group is in the βγ or δ magnetic class; (4) the complex ratio exhibits the greatest degree of association to SID flare frequency. We conclude from these results that a higher frequency of D-region ionizing flares (emitting a soft X-ray flux >2 × 10^?3 erg cm^?2 s^?1) is likely to accompany the disk transit of large area, complex spot groups. This combination of morphological factors reflects a shearing of the associated force-free magnetic field, with accumulation of free magnetic energy to power SID flares. Mutual polarity intrusion would be one observational signature of the pre-flare energy storing process.     

Sunspot as an isolated magnetic structure: Stability and oscillations

A simple energy model of a sunspot as a compact magnetic feature is described where the main energy contribution is provided by the coolest and most compressed part of the magnetic force tube of the spot at depths ranging from Wilson’s depression level (300–500 km) down to 2–3 thousand km. The equilibrium and stability conditions for such a system are analyzed using the variation principle, and oscillations of the system as a whole about the inferred equilibrium position are studied. The sunspot is shown to be stable in the magnetic field strength interval from 0.8–1 to 4–5 kG. The dependence of the eigenfrequency on magnetic field strength ω(B) is computed for the main oscillatory mode, where only the umbra of the sunspot takes part in oscillations, ω = ω ₁ (B). Lower subharmonics may appear in the case where penumbra too becomes involved in the oscillatory process: ω ₂ = ω ₁/2, ω ₃ = ω ₁/3. Theoretical curves agree well with the observational data obtained in Pulkovo using various independent methods: from temporal variations of sunspot magnetic field and from line-of-sight-velocity measurements. The periods of oscillations found range from 40 to 200 minutes.     

Total sunspot area before 1875 (reconstruction)

Empirical functions approximating the dependences of total sunspot area A on relative sunspot number W and group sunspot number GN have been found. In the function A(W), allowance for its dependence on the secular activity cycle has been made; it is shown that this allowance is not needed for the function A(GN). The yearly mean A for 1700–1874 have been reconstructed using these functions and the available W and GN time series. Having supplemented the original data with archival observations, we have been able to reconstruct the monthly mean A _W since 1821. We discuss the causes of the systematic difference between the reconstructions using W and GN.     

Recurrent magnetic activity,sunspot number and its rate of decline

Power spectral densities computed from low-latitude horizontal intensity of the Earth's magnetic field over two-year periods of declining phases of solar cycles 16 to 19 show a close relationship with the maximum relative sunspot number of the following solar cycles. The maximum sunspot number shows an exponential rise with the power density near 1/27 cd^?1; maximum R _z,however, increases linearly with power density near 1/14 cd^?1. It is also shown that the rate of decline of sunspot number in a solar cycle is almost exactly related, linearly, to power spectral density for the preceding solar cycle. Power densities near 1/27 and 1/14 cd^?1 in declining phase of solar cycle appear to be satisfactory indices for the maximum relative sunspot number of the following cycle and its rate of decline thereafter.     

The Behavior of Sunspot Contrast during Cycle 23

Results are presented from a study of various sunspot contrast parameters in broadband red (672.3 nm) Cartesian full-disk digital images taken at the San Fernando Observatory (SFO) over eight years, 1997 – 2004, of the twenty-third sunspot cycle. A subset of over 2700 red sunspots was analyzed and values of average and maximum sunspot contrast as well as maximum umbral contrast were compared to various sunspot parameters. Average and maximum sunspot contrasts were found to be significantly correlated with sunspot area (r _s=− 0.623 and r _s=− 0.714, respectively). Maximum umbral contrast was found to be significantly correlated with umbral area (r _s=− 0.535). These results are in agreement with the works of numerous other authors. No significant dependence was detected between average contrast, maximum contrast, or maximum umbral contrast during the rising phase of the solar cycle (r _s=0.024, r _s=0.033, and r _s=0.064, respectively). During the decay phase, no significant correlation was found between average contrast or maximum contrast and time (r _s=− 0.057 and r _s=0.009, respectively), with a weak dependence seen between maximum umbral contrast and cycle (r _s=0.102).     

Interpretation of vector magnetograph data including magneto-optic effects

In this paper, the presence of Faraday rotation in measurements of the orientation of a sunspot's transverse magnetic field is investigated. Using observations obtained with the Marshall Space Flight Center's (MSFC) vector magnetograph, the derived vector magnetic field of a simple, symmetric sunspot is used to calculate the degree of Faraday rotation in the azimuth of the transverse field as a function of wavelength from analytical expressions for the Stokes parameters. These results are then compared with the observed rotation of the field's azimuth which is derived from observations at different wavelengths within the Fei 5250 Å spectral line. From these comparisons, we find: the observed rotation of the azimuth is simulated to a reasonable degree by the theoretical formulations if the line-formation parameter η _o is varied over the sunspot; these variations in η _o are substantiated by the line-intensity data; for the MSFC system, Faraday rotation can be neglected for field strengths less than 1800 G and field inclinations greater than 45°; to minimize the effects of Faraday rotation in sunspot umbrae, MSFC magnetograph measurements must be made in the far wings of the Zeeman-sensitive spectral line.     

The cooling of a sunspot

This paper considers the recent criticism by Mullan (1973) of sunspot models and the cooling mechanism which I have proposed in Papers I, II and III of this series. The discussion of the cooling produced by an idealized flow cycle has been extended to include vertical temperature gradients which are consistent with a convectively unstable atmosphere. This leads to an expression for Mullan's parameter f (the ratio in which estimates of the energy flux based on an idealized Carnot cycle should be reduced) which is appropriate to this situation. It is shown that, for a cycle similar to that of Paper III, f = 0.82, while for one which has a vertical extent of order 5 Mm, f= 0.4. Hence the energy flux which, in principle, can be transported away from a sunspot by such a cycle is conservatively estimated to be 1.1 × 10²⁹ erg s^?1 compared with a typical sunspot energy deficit of 2.2 × 10²⁹ erg s^?1. Other criticisms relating to the magnetic field amplification and the ‘cool one’ model are discussed. It is concluded that the essential features of these models remain valid and that the modifications suggested by Mullan's criticism greatly increase their applicability to the sunspot problem.     

The use of abnormal quiet days in Sq(H) for predicting the magnitude of sunspot maximum at the time of preceding sunspot minimum

The previously established connection between the occurence of AQDs (“abnormal quiet days” when the phase of the solar diurnal variation of horizontal magnetic field, Sq(H), at a mid-latitude northern hemisphere station is anomalous) at sunspot minimum and the magnitude of the following sunspot maximum is examined in the light of our recent improved understanding of the nature and cause of AQDs. A small contribution to the relationship is found to arise from variations from cycle to cycle in the additional northward field which is characteristic of AQDs and leads to a reduced Sq(H) amplitude at stations poleward of the Sq focus. However, the main factor which determines the connection is a variation from one sunspot minimum to another of the amplitude of the small southward bay-like field perturbations which constitute the AQD events, and evidence is presented which suggests that this parameter may be quantitatively related to the extent of southward swing of the B_z component of the interplanetary magnetic field which determines the energy transfer from the solar wind into the magnetospheric tail. It thus appears that the magnitude of southward swing in B_z might be another solar parameter which anticipates the size of a forthcoming sunspot cycle during its build-up over the declining phase of the previous cycle and at the minimum.     

An Estimate for the Size of Sunspot Cycle 24

For the sunspot cycles in the modern era (cycle?10 to the present), the ratio of R _Z(max)/R _Z(36th month) equals 1.26±0.22, where R _Z(max) is the maximum amplitude of the sunspot cycle?using smoothed monthly mean sunspot number and R _Z(36th month) is the smoothed monthly mean sunspot number 36 months after cycle?minimum. For the current sunspot cycle?24, the 36th month following the cycle?minimum occurred in November 2011, measuring?61.1. Hence, cycle?24 likely will have a maximum amplitude of about 77.0±13.4 (the one-sigma prediction interval), a value well below the average R _Z(max) for the modern era sunspot cycles (about 119.7±39.5).     

Long-term oscillations of sunspots and a special class of artifacts in SOHO/MDI and SDO/HMI data

A specific type of artifacts (named as “p2p”), that originate due to displacement of the image of a moving object along the digital (pixel) matrix of receiver are analyzed in detail. The criteria of appearance and the influence of these artifacts on the study of long-term oscillations of sunspots are deduced. The obtained criteria suggest us methods for reduction or even elimination of these artifacts. It is shown that the use of integral parameters can be very effective against the “p2p” artifact distortions. The simultaneous observations of sunspot magnetic field and ultraviolet intensity of the umbra have given the same periods for the long-term oscillations. In this way the real physical nature of the oscillatory process, which is independent of the artifacts have been confirmed again. A number of examples considered here confirm the dependence between the periods of main mode of the sunspot magnetic field long-term oscillations and its strength. The dependence was derived earlier from both the observations and the theoretical model of the shallow sunspot. The anti-phase behavior of time variations of sunspot umbra area and magnetic field of the sunspot demonstrates that the umbra of sunspot moves in long-term oscillations as a whole: all its points oscillate with the same phase.     

A note on the values of the vertical gradient of the magnetic field in the return flux sunspot model

The Return Flux (RF) sunspot model (Osherovich, 1982) imposes a restriction on the value of the vertical gradient of the magnetic field, dB/dz, analogous to a restriction implied by the self-similar sunspot model of Schlüter and Temesvary (ST). The maximum value of the gradient, (dB/dz)_max, is shown to be 10% smaller in the RF model than in the ST model. The dependence of (dB/dz)_max on the sunspot radius is predicted.     

Lithium abundance in sunspots

The observations of lithium were carried out with the TST-2 telescope at CrAO on August 15–20, 2006. A sunspot model was calculated for the dates of observations. The lithium abundance in a sunspot and in the undisturbed photosphere was determined. It is log(N _Li) = 1.35 for the sunspot and log (N _Li) = 1.05 for the undisturbed photosphere.     

Solar cycle variation of heliosphepic parameters

This paper makes a statistical analysis of the solar cycle variation of heliospheric quantities observed at 1 AU. Two kinds of solar cycle variation with different characters have been identified, i.e. the sunspot and coronal-hole cycles, the latter is characterized by the coronal structure lifetime L_C. The kinetic and internal energy parameters of solar wind particles follow the coronal-hole cycle, reaching a maximum at 1973 or 1974. Certain parameter combinations involving IMF quantities are found to be the looked-for heliospheric quantities that follow the sunspot cycle. Among them the ratio of magnetic to kinetic energy density μ_k and the ratio of magnetic to thermal pressure μ_p show a positive correlation with the sunspot number R while the plasma parameter and the Alfvenic Mach number M_a correlate negatively with the sunspot number. The solar wind temperature T, velocity V as well as the adiabatic sound speed C_s vary basically with the coronal-hole cycle as compared to the Alfvenic speed C_a following the sunspot cycle, while the magnetic sonic speed C_as possesses a dual nature. The implication of the above results is that the solar winds observed during the sunspot maximum and the coronal-hole maximum years differ basically in their characters. The former, on the average, is a stream with a low speed, temperature, and density, but under highly magnetic control; while the latter is one with high speed, temperature and density, but under weakly magnetic control. The output of the mass, kinetic and thermal energy fluxes in the former is much less than in the latter.