-
Solar activity in the past millennia can only be reconstructed from
cosmogenic radionuclide records in terrestrial archives. However, because of
the diversity of the proxy archives, it is difficult to build a homogeneous
reconstruction. Here we provide a new consistent multiproxy reconstruction of
the solar activity over the last 9000 years, using available long-span datasets
of 10Be and 14C in terrestrial archives. A new method, based on a Bayesian
approach, was applied for the first time to solar activity reconstruction. A
Monte Carlo search for the most probable value of the modulation potential was
performed to match data from different datasets for a given time. We used six
10Be series from Greenland and Antarctica, and the global 14C production
series. The 10Be series were resampled to match wiggles related to the grand
minima in the 14C reference dataset. The GRIP and the EDML 10Be series diverge
from each other during the second half of the Holocene, while the 14C series
lies between them. A likely reason for this is the insufficiently precise
beryllium transport and deposition model for Greenland. A slow 6-millennia
variability with lows at ca. 5500 BC and 1500 AD of solar activity is found.
Two components of solar activity can be statistically distinguished: the main
'normal' component and a component corresponding to grand minima. A possible
existence of a component representing grand maxima is indicated, but it cannot
be separated from the main component in a statistically significant manner. A
new consistent reconstruction of solar activity over the last nine millennia is
presented with the most probable values of decadal sunspot numbers and their
realistic uncertainties. Independent components of solar activity corresponding
to the main moderate activity and the grand-minimum state are identified; they
may be related to different operation modes of the dynamo.
-
A ground-based neutron monitor is a standard tool to measure cosmic ray
variability near Earth, and it is crucially important to know its yield
function for primary cosmic rays. Although there are several earlier
theoretically calculated yield functions, none of them agrees with experimental
data of latitude surveys of sea-level neutron monitors, thus suggesting for an
inconsistency. A newly computed yield function of the standard sea-level 6NM64
neutron monitor is presented here separately for primary cosmic ray protons and
$\alpha-$particles, the latter representing also heavier species of cosmic
rays. The computations have been done using the GEANT-4 Planetocosmics
Monte-Carlo tool and a realistic curved atmospheric model. For the first time,
an effect of the geometrical correction of the neutron monitor effective area,
related to the finite lateral expansion of the cosmic ray induced atmospheric
cascade, is considered, that was neglected in the previous studies. This
correction slightly enhances the relative impact of higher-energy cosmic rays
(energy above 5--10 GeV/nucleon) in neutron monitor count rate. The new
computation finally resolves the long-standing problem of disagreement between
the theoretically calculated spatial variability of cosmic rays over the globe
and experimental latitude surveys. The newly calculated yield function,
corrected for this geometrical factor, appears fully consistent with the
experimental latitude surveys of neutron monitors performed during three
consecutive solar minima in 1976--77, 1986--87 and 1996--97. Thus, we provide a
new yield function of the standard sea-level neutron monitor 6NM64 that is
validated against experimental data.
-
Long-term records of sunspot number and concentrations of cosmogenic
radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's
magnetic activity over hundreds and thousands of years. We identify several
clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100
years. We found that the periods of the first five spherical harmonics of the
slow magnetic Rossby mode in the presence of a steady toroidal magnetic field
of 1200-1300 G in the lower tachocline are in perfect agreement with the time
scales of observed variations. The steady toroidal magnetic field can be
generated in the lower tachocline either due to the steady dynamo magnetic
field for low magnetic diffusivity or due to the action of the latitudinal
differential rotation on the weak poloidal primordial magnetic field, which
penetrates from the radiative interior. The slow magnetic Rossby waves lead to
variations of the steady toroidal magnetic field in the lower tachocline, which
modulate the dynamo magnetic field and consequently the solar cycle strength.
This result constitutes a key point for long-term prediction of the cycle
strength. According to our model, the next deep minimum in solar activity is
expected during the first half of this century.
-
We report on the 2012 May 17 Ground Level Enhancement (GLE) event, which is
the first of its kind in Solar Cycle 24. This is the first GLE event to be
fully observed close to the surface by the Solar Terrestrial Relations
Observatory (STEREO) mission. We determine the coronal mass ejection (CME)
height at the start of the associated metric type II radio burst (i.e., shock
formation height) as 1.38 Rs (from the Sun center). The CME height at the time
of GLE particle release was directly measured from a STEREO image as 2.32 Rs,
which agrees well with the estimation from CME kinematics. These heights are
consistent with those obtained for cycle-23 GLEs using back-extrapolation. By
contrasting the 2012 May 17 GLE with six other non-GLE eruptions from
well-connected regions with similar or larger flare size and CME speed, we find
that the latitudinal distance from the ecliptic is rather large for the non-GLE
events due to a combination of non-radial CME motion and unfavorable solar B0
angle, making the connectivity to Earth poorer. We also find that the coronal
environment may play a role in deciding the shock strength.
-
We present the results of full new calculation of radiocarbon 14C production
in the Earth atmosphere, using a numerical Monte-Carlo model. We provide, for
the first time, a tabulated 14C yield function for the energy of primary cosmic
ray particles ranging from 0.1 to 1000 GeV/nucleon. We have calculated the
global production rate of 14C, which is 1.64 and 1.88 atoms/cm2/s for the
modern time and for the pre-industrial epoch, respectively. This is close to
the values obtained from the carbon cycle reservoir inventory. We argue that
earlier models overestimated the global 14C production rate because of outdated
spectra of cosmic ray heavier nuclei. The mean contribution of solar energetic
particles to the global 14C is calculated as about 0.25% for the modern epoch.
Our model provides a new tool to calculate the 14C production in the Earth's
atmosphere, which can be applied, e.g., to reconstructions of solar activity in
the past.
-
The most powerful explosions on the Sun [...] drive the most severe
space-weather storms. Proxy records of flare energies based on SEPs in
principle may offer the longest time base to study infrequent large events. We
conclude that one suggested proxy, nitrate concentrations in polar ice cores,
does not map reliably to SEP events. Concentrations of select radionuclides
measured in natural archives may prove useful in extending the time interval of
direct observations up to ten millennia, but as their calibration to solar
flare fluences depends on multiple poorly known properties and processes, these
proxies cannot presently be used to help determine the flare energy frequency
distribution. Being thus limited to the use of direct flare observations, we
evaluate the probabilities of large-energy solar explosions by combining solar
flare observations with an ensemble of stellar flare observations. We conclude
that solar flare energies form a relatively smooth distribution from small
events to large flares, while flares on magnetically-active, young Sun-like
stars have energies and frequencies markedly in excess of strong solar flares,
even after an empirical scaling with the mean activity level of these stars. In
order to empirically quantify the frequency of uncommonly large solar flares
extensive surveys of stars of near-solar age need to be obtained, such as is
feasible with the Kepler satellite. Because the likelihood of flares larger
than approximately X30 remains empirically unconstrained, we present indirect
arguments, based on records of sunspots and on statistical arguments, that
solar flares in the past four centuries have likely not substantially exceeded
the level of the largest flares observed in the space era, and that there is at
most about a 10% chance of a flare larger than about X30 in the next 30 years.
-
We present an overview of the observed properties of the GLEs and those of
the associated flares and CMEs. The solar eruptions are very intense involving
X-class flares and extreme CME speeds (average ~2000 km/s). The active regions
in which the GLE events originate are generally large: 1290 msh (median 1010
msh) compared to 934 msh (median: 790 msh) for SEP-producing active regions.
The initial acceleration of GLE-associated CMEs is much larger (by a factor of
2) than that of ordinary CMEs (2.3 km/s2 vs.1 km/s2). The GLE particle release
is delayed with respect to the onset of all electromagnetic signatures of the
eruptions: type II bursts, low frequency type III bursts, soft X-ray flares and
CMEs. The presence of metric type II radio bursts some 17 min (median: 16 min;
range: 3 to 48 min) before the GLE onset indicates shock formation well before
the particle release. The release of GLE particles occurs when the CMEs reach
an average height of ~3.09 Rs for well-connected events. For poorly connected
events, the average CME height at GLE particle release is ~66% larger (mean:
5.18 Rs). The longitudinal dependence is consistent with shock accelerations
because the shocks from poorly connected events need to expand more to cross
the field lines connecting to an Earth observer. The CME height at metric type
II burst onset is in the narrow range 1.29 to 1.8 Rs, with A mean of 1.53 Rs.
The CME heights at metric type II burst onset and GLE particle release
correspond to the minimum and maximum in the Alfven speed profile. The CME
heights at GLE particle release are in good agreement with those obtained from
the velocity dispersion analysis (Reames, 2009a,b) including the source
longitude dependence. We also discuss the implications of the delay of GLE
particle release with respect to complex type III bursts and hard X-ray
emission.
-
Maunder Minimum forms an archetype for the Grand minima, and detailed
knowledge of its temporal development has important consequences for the solar
dynamo theory dealing with long-term solar activity evolution. Here we
reconsider the current paradigm of the Grand minimum general scenario by using
newly recovered sunspot observations by G. Marcgraf and revising some earlier
uncertain data for the period 1636--1642, i.e., one solar cycle before the
beginning of the Maunder Minimum. The new and revised data dramatically change
the magnitude of the sunspot cycle just before the Maunder Minimum, from 60--70
down to about 20, implying a possibly gradual onset of the Minimum with reduced
activity started two cycles before it. This revised scenario of the Maunder
Minimum changes, through the paradigm for Grand solar/stellar activity minima,
the observational constraint on the solar/stellar dynamo theories focused on
long-term studies and occurrence of Grand minima.