On the possibility of constructing a radiative sunspot model in magnetohydrostatic equilibrium

It is currently believed that it is impossible to construct a radiative sunspot model in magnetohydrostatic equilibrium unless magnetic fields below the surface are excessively large (> 100 kG). This belief is based on results obtained using the mixing length theory of convection. We wish to point out that by using a different theory of convection, due to Öpik (1950), it is possible to compute a radiative sunspot model in which the field becomes no greater than 9000 G. By applying two boundary conditions, (i) depth of spot equals depth of convection zone, (ii) magnetic field has zero gradient at the base of the spot, we show that a radiative spot has a unique effective temperature for a given Wilson depression, . For = 650 km, we find T _e = 3800K ; for = 150 km, T _e = 3950K. According to our model, spots having T _e cooler than these values should not exist.     

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.     

A new empirical model of a sunspot umbra

H and K-line spectra of quiescent prominences, taken with the slit placed normal to the limb, commonly reveal a gas streaming (5–50 km/s) that is peculiar to the upper edge of these objects. On the average this streaming is uni-directional and consistent with a hypothetical east-west wind.Operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation.     

A refined measurement of the sunspot radiative flux deficit

New measurements of the radiative flux deficits of two large sunspots are presented, based on detailed isophotometric maps. Results are given separately for umbrae and penumbrae. The umbral and penumbral deficits are 4–5 × 10¹⁰ and 1–1.5 × 10¹⁰ erg cm^–2 s^–1 respectively, the larger figures ref to the larger spot. Over limited areas centered on the umbral cores the deficits for the two spots amount to 76 and 86% of the photospheric flux.     

Preliminary observation of missing energy flux of sunspot

We have searched for the missing energy flux of the sunspot by measuring white-light photographs over the region around the sunspots exceeding ten times of the diameter of spots. It was found that the excess brightness amounts to 0.3–0.8% of the local continuum intensity immediately beyond the penumbra and smoothly decrease to zero at the edge of our inspected radii of 130 000–160000km. This result suggests that the flux deficit in sunspots of the order of 30% over the umbra-penumbra is roughly the same as the excess flux found in the very large surrounding region which includes photospheric facular contributions. Among nineteen spots fourteen showed this extended bright ring. Narrow band photoelectric scans revealed that the excess radiation originates in the continuum, not from weakened absorption lines.     

Studies of grand minima in sunspot cycles by using a flux transport solar dynamo model

We propose that grand minima in solar activity are caused by simultane- ous fluctuations in the meridional circulation and the Babcock-Leighton mechanism for the poloidal field generation in the flux transport dynamo model. We present the following results： （a） fluctuations in the meridional circulation are more effective in producing grand minima; （b） both sudden and gradual initiations of grand minima are possible; （c） distributions of durations and waiting times between grand minima seem to be exponential; （d） the coherence time of the meridional circulation has an effect on the number and the average duration of grand minima, with a coherence time of about 30 yr being consistent with observational data. We also study the occurrence of grand maxima and find that the distributions of durations and waiting times between grand maxima are also exponential, like the grand minima. Finally we address the question of whether the Babcock-Leighton mechanism can be operative during grand minima when there are no sunspots. We show that an a-effect restricted to the upper portions of the convection zone can pull the dynamo out of the grand minima and can match various observational requirements if the amplitude of this a-effect is suitably fine-tuned.     

A model of the sunspot umbra

Equations governing the structure of the umbra of a single spot have been integrated on the spot-axis. It is shown that a consistent umbral model can be obtained only for a narrow range of the electron pressures at the spot surface. The spot center is found to be at a depth of about 400 km below the normal solar surface. The reduced energy flux observed at the surface is assumed to be due to the effect of the spot-magnetic field on subphotospheric convection and an empirical factor is introduced to take into account the reduction in the convective energy flux. With an inferred expression for as a function of the internal and magnetic-energy density it is shown that in a consistent model the physical variables in the spot approach their ambient values at about 2330 km below the undisturbed solar surface and at the same depth the energy flux approaches the normal solar value and the magnetic interference with convection vanishes. The present investigation is a refinement of an earlier paper by Chitre (1963).Supported in part by the Office of Naval Research [Nonr-220(47)] and the National Science Foundation [GP-5391].     

Acoustic wave ducting along magnetic flux tubes in the sunspot regions

The acoustic waves generated in the solar atmosphere propagate globally as well as upwards. These waves interact with the solar magnetic field structures and are ducted upwards. The velocity of these modified acoustic waves is shown to vary in a modelled solar atmosphere. The solar plasma propagating upwards with these waves are likely to alter the observed features of spicules, granules, and supergranules during changing phases of sunspot regions.     

Oscillations of flux tube bundles and the structure of sunspot magnetic fields

Observational aspects of the previously found quasi-hourly oscillations of magnetic fragments in sunspot polar coordinates are investigated. The orientation of the oscillations is shown to be azimuthally anisotropic, with their amplitude reaching a maximum in penumbra at a distance of ～0.8 sunspot radius (the maximum amplitude is estimated to be 3700 km). Based on the detected deviations of the oscillations from the radial direction, we numerically simulate the horizontal configuration of field lines in the region of the major spots in bipolar groups.     

Rotational modulation of Ca K flux ratio and sunspot number

The results of daily monitoring of the Ca K flux variations of integrated sunlight during the period 1 July 1979–18 February 1980 are presented, and the data analyzed by both power spectrum and autocorrelation methods. The latter method is found to be superior; in it the autocorrelation function is fit to a model having three parameters, viz. the maximum long-term autocorrelation c ₁, the rotational period ₀ and an effective lifetime ₁. Values of c ₁ = 0.26, ₀ = 25.9 days, and ₁ = 500 days are obtained for the Ca K flux ratio; somewhat different values of c ₁ = 0.77, ₀ = 28.0 days, and ₁ = 89 days are obtained for the sunspot number for the same period. Rotational modulation is clearly evident for the relatively weak chromosphere of the Sun and should be even more easily observed for more active stars.     

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.     

The influence of the sunspot model on the Li-abundance

The dependence of the Li-abundance on the equivalent width of the Li-resonance doublet at 6708 Å is calculated for different umbral models. The choice of the model strongly influences the deduced Li-abundance (up to a factor 8 or log _Li = 0.9) even when using recent umbral models. Detailed discussion of the observations and reduction with the most suitable umbral model (Stellmacher and Wiehr, 1970) leads to an abundance of log _Li = 1.1±0.05 (in the log _H = 12.0 scale).     

prior to surface appearance of sunspot flux tubes and magneto-thermal pulsation of the Sun

We found an evidence that the luminosity of the Sun systematically decreased about 20 days before sunspot surface appearance by analysing time-lag correlation of time derivatives of running mean time profiles of the data of the Active Cavity Radiometer Irradiance Monitor (ACRIM) I experiment on board of Solar Maximum Mission (SMM) and of the data of the daily sunspot number. This indicates that sunspot flux tube cooling and heat transport blocking by the flux tubes start to take place in the interior of the solar convection zone well before the sunspot surface appearance. From this finding and our previous finding that the luminosity of the Sun systematically increased and the blocked heat appeared on the surface about 50 days after the sunspot surface appearance, a new view of sunspot formation and dynamics and a new view of the luminosity modulation emerged. (i) Sunspots of a solar cycle are formed from clusters of flux tubes which can be seen in the running mean time profile of the sunspot number as a peak with duration on the order of 100 to 200 days. (ii) Heat flow is blocked by the cluster of sunspot flux tubes inside the convection zone to decrease the luminosity about 20 days before the surface emergence of the sunspot cluster. (iii) The blocked heat appears on the surface about 50 days after the surface emergence of the cluster of sunspot flux tubes to heat up the surface. This appears as a thermal pulse in the running mean time profile of the ACRIM dat in between the peaks of the sunspot running mean time profile. This process of heating the surface makes the temperature gradient less steep and weakens the buoyancy of sunspot flux tubes below the surface. (vi) The radiative cooling of the surface layer by the excess heat release steepens the temperature gradient so that the buoyancy of the sub-surface magnetic flux tubes becomes stronger to cause the next surge of emergence of a cluster of sunspots and other magnetic activities, which creates a peak in the time profile of the sunspot number. We call this peak a magnetic pulse of the Sun and the coupled process of alternating pulsed appearance of heat and sunspots the magneto-thermal pulsation of the Sun.     

Polar coronal holes and the sunspot cycle. A new method to predict sunspot numbers

Twenty years ago, Ohl (1966, 1968) found a correlation between geomagnetic activity around the minimum of the solar cycle and the Wolf sunspot number in the maximum of the following solar cycle. In this paper we shall show that such a relation means indeed a relation between the polar coronal holes area around the minimum of the solar cycle and the sunspot number in the maximum of the next. In fact, a very high positive correlation exists between the temporal evolution of the size of polar coronal holes and the Wolf sunspot number 6.3. years later.     

A new method to calculate the resonant response of a sunspot model atmosphere to magneto-atmospheric waves

In order to understand the observed oscillations in sunspots we present a new method for calculating the resonant response of a realistic semi-empirical model of the sunspot umbral atmosphere and subphotosphere to magneto-atmospheric waves in a vertical magnetic field. The depth dependence of both the adiabatic coefficient and the turbulent pressure is taken into account. This requires an extension of the wave equations by Ferraro & Plumpton (1958). We compare the coefficients of wave transmission, re flection, and conversion between fast mode and slow mode waves for different assumptions, compare the results with those from earlier modelling efforts, and point out possible sources of mistakes. The depth dependence of the adiabatic coefficient strongly influences the resulting spectrum of resonance frequencies. The condition of a conservation of wave flux is violated if the depth dependence of the turbulent pressure is not properly considered.     

To the description of long-term variations in the solar magnetic flux: The sunspot area index

We show that the Wolf sunspot numbers W and the group sunspot numbers GSN are physically different indices of solar activity and that it is improper to compare them. Based on the approach of the so-called “primary” indices from the observational series of W(t) and GSN(t), we suggest series of yearly mean sunspot areas beginning in 1610 and monthly mean sunspot areas beginning in 1749.     

The spectrum of umbral oscillations of the Scheuer-Thomas sunspot model

Proper periods are calculated for oscillations of an isothermal gaseous layer with a homogeneous vertical magnetic field, limited by two rigid horizontal boundaries. The obtained results are in good agreement with the observed umbral oscillations in a sunspot.