Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Variations of the low‐degree solar p‐mode frequencies for 10 years of GOLF observations

We experiment with a method of measuring the frequency of solar p modes, intended to extend the passband for the variations of the frequency spectrum as high as possible. So far this passband is limited to a fraction of μ Hz for the classical analysis based on numerical fits of a theoretical line profile to a power spectrum averaged over periods lasting at least several weeks. This limit for the present analysis can be shifted to the mHz range, corresponding to some of the “5 min” oscillations, but in this range we use a lower resolution which allows us to separate odd and even p modes. We show an example of the results for long term variations and apply this analysis to search for a modulation of the p‐mode frequency spectrum by asymptotic series of solar g modes. A faint signal is found in the analysis of 10 years of GOLF data. This very preliminary result possibly indicates the detection of a small number of g modes of degree l = 1. A tentative determination of an observational value of the parameter P₀ follows. P₀ is the scaling factor of the asymptotic series of g modes and is a key data for solar core physics. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Diamagnetic pumping near the base of a stellar convection zone

The property of inhomogeneous turbulence in conducting fluids to expel large‐scale magnetic fields in the direction of decreasing turbulence intensity is shown as important for the magnetic field dynamics near the base of a stellar convection zone. The downward diamagnetic pumping confines a fossil internal magnetic field in the radiative core so that the field geometry is appropriate for formation of the solar tachocline. For the stars of solar age, the diamagnetic confinement is efficient only if the ratio of turbulent magnetic diffusivity η_T of the convection zone to the (microscopic or turbulent) diffusivity η_in of the radiative interior is η_T/η_in ≥ 10⁵. Confinement in younger stars requires larger η_T/η_in. The observation of persistent magnetic structures on young solar‐type stars can thus provide evidence for the nonexistence of tachoclines in stellar interiors and on the level of turbulence in radiative cores. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spectroscopy of solar prominences from space and ground

Two quiescent solar prominences were observed in July 2000 from SUMER aboard SOHO and from the two German solar telescopes at Tenerife. Two‐dimensional images taken at the VTT simultaneously in the spectral lines Hβ at 4862 Å and Ca II at 8542 Å show no significant spatial variation of their pressure‐sensitive emission ratio. Slit spectra of the Ca II 8542 Å and He I 10830 Å lines obtained at the Gregory‐Coudé telescope yield 8000 K < T_kin < 9000 K and 3 km/s < V_n–th < 8 km/s. Among the various spectral ranges observed with SUMER, we first investigate the Lyman emission lines, which were fitted by Gaussians yielding reliable spectral radiances and line widths for the series members 5 < k < 18. A determination of the level population gives for the lower series members a Boltzmann temperature of 60 000 K, the higher members being over‐populated. This temperature indicates an origin of the Lyman lines from hot surroundings of the cool prominence body seen in the ground‐based data; this also holds for the ‘hotter’ SUMER lines.     

Structure of a simple sunspot from the inversion of IR spectral data

Analysis of spectral data of two neighboring infrared lines, Fe I 15648.5 Å (g = 3) and FeI 15652.9 Å (g_eff = 1.53) are carried out for a simple sunspot when it was near the solar disk center (μ = 0.92), to understand the basic structure of sunspot magnetic field. Inversions of Stokes profiles are carried out to derive different atmospheric parameters both as a function of location within the sunspot and height in the atmosphere. As a result of the inversion we have obtained maps of magnetic field strength, temperature, line‐of‐sight velocity, field inclination and azimuth for different optical depth layers between log(τ₅) = 0 and log(τ₅) = –2.0. In this paper we present few results from our inversion for a layer averaged between log(τ5) from 0.0 to –0.5.     

Statistical Assessment of Photospheric Magnetic Features in Imminent Solar Flare Predictions

In this study we use the ordinal logistic regression method to establish a prediction model, which estimates the probability for each solar active region to produce X-, M-, or C-class flares during the next 1-day time period. The three predictive parameters are (1) the total unsigned magnetic flux T _flux, which is a measure of an active region’s size, (2) the length of the strong-gradient neutral line L _gnl, which describes the global nonpotentiality of an active region, and (3) the total magnetic dissipation E _diss, which is another proxy of an active region’s nonpotentiality. These parameters are all derived from SOHO MDI magnetograms. The ordinal response variable is the different level of solar flare magnitude. By analyzing 174 active regions, L _gnl is proven to be the most powerful predictor, if only one predictor is chosen. Compared with the current prediction methods used by the Solar Monitor at the Solar Data Analysis Center (SDAC) and NOAA’s Space Weather Prediction Center (SWPC), the ordinal logistic model using L _gnl, T _flux, and E _diss as predictors demonstrated its automatic functionality, simplicity, and fairly high prediction accuracy. To our knowledge, this is the first time the ordinal logistic regression model has been used in solar physics to predict solar flares.     

p‐mode power variation with solar atmosphere as observed in the Na D1 and K spectral lines

In this work we investigate p‐mode power variation with solar atmosphere. To this aim, we use THÉMIS observations of the Na D1 (λ 5896 Å) and K (λ 7699 Å) spectral lines. While the formation heights of the K spectral line are essentially located in the photospheric layer, the formation heights of the Na D1 line span a much wider region: from photosphere up to chromosphere. Hence, we had the opportunity to infer p‐mode power variation up to the chromospheric layer. By analyzing power spectra obtained by temporal series at different points of the Na D1 and K spectral lines, we confirm and quantify the increase in p‐mode power towards higher atmospheric layers. Furthermore, the large span in formation heights of the Na D1 line induces a larger enhancement of p‐mode power with solar atmosphere compared to the K spectral line. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The hemispheric variation of the flare index during solar cycles 20–23

In this paper, the monthly counts of flare index in the northern and southern hemispheres are used to investigate the hemispheric variation of the flare index in each of solar cycles 20–23. It is found that, (1) the flare index is asymmetrically distributed in each solar cycle and its asymmetry is a real phenomenon; (2) the flare index in the northern hemisphere begins earlier than that in the southern hemisphere in each of solar cycles 20–23, and the phase shifts between the two hemispheres show an odd‐even pattern; (3) although the flare index dominating in a hemisphere does not mean that it leads in phase in this hemisphere in individual solar cycle, these two features have an intrinsic relationship. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Detection and temporal coherence of p‐modes below 1.4 mHz

Data collected recently by the helioseismic experiments aboard the SOHO spacecraft have allowed the detection of low degree p‐modes with increasingly lower order n. In particular, the GOLF experiment is currently able to unambiguously identify low degree modes with frequencies as low as 1.3 mHz. The detection of p‐modes with very low frequency (i.e., low n), is difficult due to the low signal‐to‐noise ratio in this spectral region and its contamination by solar signals that are not of acoustic origin. To address this problem without using any theoretical a priory, we propose a methodology that relies only on the inversion of observed values to define a spectral window for the expected locations of these low frequency modes. The application of this method to 2920‐day‐long GOLF observations is presented and its results discussed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Cyclical behavior of solar filaments

The Carte Synoptique catalogue of solar filaments from 1919 March to 1957 July, corresponding to complete cycles 16‐18, is utilized to show the latitudinal migrations of solar filaments at low (≤50°) and high (>50°) latitudes and the latitudinal distributions of solar filaments for all solar filaments, solar filaments whose maximum lengths during solar disk passage are less than or equal to 70° and solar filaments whose maximum lengths during solar disk passage are larger than 70°. The results show the following. (1) The latitudinal migrations of all low‐latitude solar filaments and low‐latitude solar filaments whose maximum lengths during solar disk passage are less than or equal to 70° follow the Spörer sunspot law. However, the latitudinal migration of low‐latitude solar filaments whose maximum lengths during solar disk passage are larger than 70° do not follow the Spörer sunspot law: there is no equatorward and no poleward drift. The latitudinal migration of high‐latitude solar filaments whose maximum lengths during solar disk passage are larger than 70° is more significant than those of all high‐latitude solar filaments and high‐latitude solar filaments whose maximum lengths during solar disk passage are less than or equal to 70°: there is a poleward migration from the latitude of about 50° to 70° and an equatorward migration from the latitude of about 70° to 50° of all high‐latitude solar filaments and high‐latitude solar filaments whose maximum lengths during solar disk passage are less than or equal to 70° and there is a poleward migration from the latitude of about 50° to 80° and an equatorward migration from the latitude of about 80° to 50° of high‐latitude solar filaments whose maximum lengths during solar disk passage are larger than 70°. (2) The statistical characteristics of latitudinal distribution of solar filaments whose maximum lengths during solar disk passage are larger than 70° is different from those of all solar filaments and solar filaments whose maximum lengths during solar disk passage are less than or equal to 70° (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Effect of Abnormal Granulation on Acoustic Wave Travel Times and Mode Frequencies

Observations indicate that in plage areas (i.e. in active regions outside sunspots) acoustic waves travel faster than in the quiet Sun, leading to shortened travel times and higher p-mode frequencies. Coupled with the 11-year variation of solar activity, this may also explain the solar cycle variation of oscillation frequencies. While it is clear that the ultimate cause of any difference between the quiet Sun and plage is the presence of magnetic fields of order 100 G in the latter, the mechanism by which the magnetic field exerts its influence has not yet been conclusively identified. One possible such mechanism is suggested by the observation that granular motions in plage areas tend to be slightly “abnormal”, dampened compared to the quiet Sun. In this paper we consider the effect that abnormal granulation observed in active regions should have on the propagation of acoustic waves. Any such effect is found to be limited to a shallow surface layer where sound waves propagate nearly vertically. The magnetically suppressed turbulence implies higher sound speeds, leading to shorter travel times. This time shift Δ τ is independent of the travel distance, while it shows a characteristic dependence on the assumed plage field strength. As a consequence of the variation of the acoustic cutoff with height, Δ τ is expected to be significantly higher for higher frequency waves within the observed regime of 3 – 5 mHz. The lower group velocity near the upper reflection point further leads to an increased envelope time shift, as compared to the phase shift. p-mode frequencies in plage areas are increased by a corresponding amount, Δ ν/ν=ν Δ τ. These characteristics of the time and frequency shifts are in accordance with observations. The calculated overall amplitudes of the time and frequency shifts are comparable to, but still significantly less than (by a factor of 2 to 5), those suggested by measurements.     

IR spectroscopy of COmosphere dynamics with the CO ?rst overtone band

We discuss observations of the weak ?rst overtone (Δν = 2) CO absorption band near 2300 nm with the U.S. National Solar Observatory Array Camera (NAC), a modern mid‐infrared detector. This molecular band provides a thermal diagnostic that forms lower in the atmosphere than the stronger fundamental band near 4600 nm. The observed center‐to‐limb increase in CO line width qualitatively agrees with the proposed higher temperature shocks or faster plasma motions higher in the COmosphere. The spatial extent of chromospheric shock waves is currently at or below the diffraction limit of the available CO lines at existing telescopes. Five minute period oscillations in line strength and measured Doppler shifts are consistent with the p‐mode excitation of the photospheric gas. We also show recent efforts at direct imaging at 4600 nm. We stress that future large‐aperture solar telescopes must be teamed with improved, dynamic mid‐infrared instruments, like the NAC, to capitalize on the features that motivate such facilities (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A solar super‐flare as cause for the 14C variation in AD 774/5?

We present further considerations regarding the strong ¹⁴C variation in AD 774/5. For its cause, either a solar super‐flare or a short gamma‐ray burst were suggested. We show that all kinds of stellar or neutron star flares would be too weak for the observed energy input at Earth in AD 774/5. Even though Maehara et al. (2012) present two super‐flares with ∼10³⁵ erg of presumably solar‐type stars, we would like to caution: These two stars are poorly studied and may well be close binaries, and/or having a M‐type dwarf companion, and/or may be much younger and/or much more magnetic than the Sun – in any such case, they might not be true solar analog stars. From the frequency of large stellar flares averaged over all stellar activity phases (maybe obtained only during grand activity maxima), one can derive (a limit of) the probability for a large solar flare at a random time of normal activity: We find the probability for one flare within 3000 years to be possibly as low as 0.3 to 0.008 considering the full 1σ error range. Given the energy estimate in Miyake et al. (2012) for the AD 774/5 event, it would need to be ∼2000 stronger than the Carrington event as solar super‐flare. If the AD 774/5 event as solar flare would be beamed (to an angle of only ∼24°), 100 times lower energy would be needed. A new AD 774/5 energy estimate by Usoskin et al. (2013) with a different carbon cycle model, yielding 4 ot 6 time lower ¹⁴C production, predicts 4–6 times less energy. If both reductions are applied, the AD 774/5 event would need to be only ∼4 times stronger than the Carrington event in 1859 (if both had similar spectra). However, neither ¹⁴C nor ¹⁰Be peaks were found around AD 1859. Hence, the AD 774/5 event (as solar flare) either was not beamed that strongly, and/or it would have been much more than 4‐6 times stronger than Carrington, and/or the lower energy estimate (Usoskin et al. 2013) is not correct, and/or such solar flares cannot form (enough) ¹⁴C and ¹⁰Be. The 1956 solar energetic particle event was followed by a small decrease in directly observed cosmic rays. We conclude that large solar super‐flares remain very unlikely as the cause for the ¹⁴C increase in AD 774/5. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multi‐epoch time‐resolved photometry of the eclipsing polar CSS081231:071126+440405

The eclipsing polar CSS081231:071126+440405 turned bright (V_max ∼ 14.5) in late 2008 and was subsequently observed intensively with small and medium‐sized telescopes. A homogeneous analysis of this comprehensive dataset comprising 109 eclipse epochs is presented and a linear ephemeris covering the five years of observations, about 24000 orbital cycles, is derived. Formally this sets rather tight constraints on the mass of a hypothetical circumbinary planet, M_pl ≤ 2 M_Jup. This preliminary result needs consolidation by long‐term monitoring of the source. The eclipse lasts 433.08 ± 0.65 s, and the orbital inclination is found to be i = 79.3°–83.7°. The centre of the bright phase displays accretion‐rate dependent azimuthal shifts. No accretion geometry is found that explains all observational constraints, suggesting a complex accretion geometry with possible pole switches and a likely non‐dipolar field geometry. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analytical theory of tachocline transport at mid latitudes

We provide a theory of magnetic diffusion, momentum transport, and mixing in the solar tachocline by considering magnetohydrodynamics (MHD) turbulence on a β plane subject to a large scale shear (provided by the latitudinal differential rotation). In the strong magnetic field regime, we find that the turbulent viscosity and diffusivity are reduced by magnetic fields only, similarly to the two-dimensional MHD case (without Rossby waves). In the weak magnetic field regime, we find a crossover scale (L_R) from a Alfvén dominated regime (on small scales) to a Rossby dominated regime (on large scales). For parameter values typical of the tachocline, L_R is larger than the solar radius so that Rossby waves are unlikely to play an important role in the transport of magnetic field and angular momentum. This is mainly due to the enhancement of magnetic back-reaction by shearing which efficiently generates small scales, thus strong currents. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent dynamo near tachocline and reconstruction of azimuthal magnetic field in the solar convection zone

In order to extend the abilities of the αΩ dynamo model to explain the observed regularities and anomalies of the solar magnetic activity, the negative buoyancy phenomenon and the magnetic quenching of the α effect were included in the model, as well as newest helioseismically determined inner rotation of the Sun were used. Magnetic buoyancy constrains the magnitude of toroidal field produced by the Ω effect near the bottom of the solar convection zone (SCZ). Therefore, we examined two “antibuoyancy” effects: i) macroscopic turbulent diamagnetism and ii) magnetic advection caused by vertical inhomogeneity of fluid density in the SCZ, which we call the ∇ρ effect. The Sun's rotation substantially modifies the ∇ρ effect. The reconstruction of the toroidal field was examined assuming the balance between mean‐field magnetic buoyancy, turbulent diamagnetism and the rotationally modified ∇ρ effect. It is shown that at high latitudes antibuoyancy effects block the magnetic fields in the deep layers of the SCZ, and so the most likely these deep‐rooted fields could not become apparent at the surface as sunspots. In the near‐equatorial region, however, the upward ∇ρ effect can facilitate magnetic fields of about 3000 – 4000 G to emerge through the surface at the sunspot belt. Allowance for the radial inhomogeneity of turbulent velocity in derivations of the helicity parameter resulted in a change of sign of the α effect from positive to negative in the northern hemisphere near the bottom of the SCZ. The change of sign is very important for direction of the Parker's dynamo‐waves propagation and for parity of excited magnetic fields. The period of the dynamo‐wave calculated with allowance for the magnetic quenching is about seven years, that agrees by order of magnitude with the observed mean duration of the sunspot cycles. Using the modern helioseismology data to define dynamo‐parameters, we conclude that north‐south asymmetry should exist in the meridional field. At low latitudes in deep layers of the SCZ, the αΩ dynamo excites most efficiency the dipolar mode of the meridional field. Meanwhile, in high‐latitude regions a quadrupolar mode dominates in the meridional field. The obtained configuration of the net meridional field is likely to explain the magnetic anomaly of polar fields (the apparent magnetic “monopole”) observed near the maxima of solar cycles. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Relationship between group sunspot numbers and Wolf sunspot numbers

Continuous wavelet transform and cross‐wavelet transform have been used to investigate the phase periodicity and synchrony of the monthly mean Wolf (R_z) and group (R_g) sunspot numbers during the period of June 1795 to December 1995. The Schwabe cycle is the only one common period in R_g and R_z, but it is not well‐defined in case of cycles 5–7 of R_g and in case of cycles 5 and 6 of R_z. In fact, the Schwabe period is slightly different in R_g and R_z before cycle 12, but from cycle 12 onwards it is almost the same for the two time series. Asynchrony of the two time series is more obviously seen in cycles 5 and 6 than in the following cycles, and usually more obviously seen around the maximum time of a cycle than during the rest of the cycle. R_g is found to fit R_z better in both amplitudes and peak epoch during the minimum time time of a solar cycle than during the maximum time of the cycle, which should be caused by their different definition, and around the maximum time of a cycle, R_g is usually less than R_z. Asynchrony of R_g and R_z should somewhat agree with different sunspot cycle characteristics exhibited by themselves (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

27‐day variations of solar indices and cosmic ray neutron monitor intensities

Solar activity indices (coronal, chromospheric as well as photospheric) and cosmic ray neutron monitor rates (different cut‐off rigidity) have been used to study 27‐day variations in the years from 1957 to 2004. Daily data were employed for this purpose, analysed by the FFT and wavelet techniques. To work with a continuous data set for the cosmic rays (CR), the ‘Composite Cosmic Ray’ (CCR) set was first created from the observations carried out at different neutron monitor stations. The CCR frequency analysis shows significant 27‐day variations in the intensity of CR, with its amplitude's values very sensitive to the sign of the quantity qA. The most significant 27‐day variations of CR were found not to correlate with those of other solar indices. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Visual and photoelectric measurements of the solar diameter (1972‐2002): Methods and results

Measurements of the solar diameter using both visual and photoelectric drift scan techniques have been made since 1972 using two almost identical 45‐cm Gregory‐Coudé telescopes at Locarno/Switzerland and Izaña/Tenerife. The method, in which a time measurement substitutes an angular measurement, is especially suited to obtain about 30 measurements of the absolute solar semidiameter per day. During the years 1972–2002 a total of 10996 visual timing measurements have been made on 320 observing days, an additional 1373 photoelectric recordings have been obtained on 117 observing days. The data were used to study the long‐term behaviour of the solar semidiameter R at unit distance and its possible variations. No fluctuations dR in excess of about ±0.05″ have been found, neither long‐term nor short‐term. The photoelectric semidiameter, which refers to the continuum at λ ≈ 585 nm, is Rphot = (959.89 ± 0.12)″. The visual semidiameter, which refers to the footpoint of the limb intensity profile at λ ≈ 550 nm, is Rvis = (960.62 ± 0.02)″.