Aims. The solar activity in the past millennia can only be reconstructed from cosmogenic radionuclide proxy records in terrestrial archives. However, because of the diversity of the proxy archives, ...it is difficult to build a homogeneous reconstruction. All previous studies were based on individual, sometimes statistically averaged, proxy datasets. Here we aim to provide a new consistent multi-proxy reconstruction of the solar activity over the last 9000 yr, using all available long-span datasets of 10Be and 14C in terrestrial archives. Methods. A new method, based on a Bayesian approach, was applied for the first time to solar activity reconstruction. A Monte Carlo search (using the χ2 statistic) for the most probable value of the modulation potential was performed to match data from different datasets for a given time. This provides a straightforward estimate of the related uncertainties. We used six 10Be series of different lengths (from 500–10 000 yr) 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 stability of the long data series was tested. Results. The Greenland Ice-core Project (GRIP) and the Antarctic EDML (EPICA Dronning Maud Land) 10Be series diverge from each other during the second half of the Holocene, while the 14C series lies in between them. A likely reason for the discrepancy is the insufficiently precise beryllium transport and deposition model for Greenland, which leads to an undercorrection of the GRIP series for the geomagnetic shielding effect. A slow 6–7 millennia variability with lows at ca. 5500 BC and 1500 AD in the long-term evolution of solar activity is found. Two components of solar activity can be statistically distinguished: the main component, corresponding to the “normal” moderate level, 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. Conclusions. 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.
Aims. Cosmogenic isotopes provide the only quantitative proxy for analyzing the long-term solar variability over a centennial timescale. While essential progress has been achieved in both ...measurements and modeling of the cosmogenic proxy, uncertainties still remain in the determination of the geomagnetic dipole moment evolution. Here we aim at improving the reconstruction of solar activity over the past nine millennia using a multi-proxy approach. Methods. We used records of the 14C and 10Be cosmogenic isotopes, current numerical models of the isotope production and transport in Earth’s atmosphere, and available geomagnetic field reconstructions, including a new reconstruction relying on an updated archeo- and paleointensity database. The obtained series were analyzed using the singular spectrum analysis (SSA) method to study the millennial-scale trends. Results. A new reconstruction of the geomagnetic dipole field moment, referred to as GMAG.9k, is built for the last nine millennia. New reconstructions of solar activity covering the last nine millennia, quantified in terms of sunspot numbers, are presented and analyzed. A conservative list of grand minima and maxima is also provided. Conclusions. The primary components of the reconstructed solar activity, as determined using the SSA method, are different for the series that are based on 14C and 10Be. This shows that these primary components can only be ascribed to long-term changes in the terrestrial system and not to the Sun. These components have therefore been removed from the reconstructed series. In contrast, the secondary SSA components of the reconstructed solar activity are found to be dominated by a common ≈2400-year quasi-periodicity, the so-called Hallstatt cycle, in both the 14C and 10Be based series. This Hallstatt cycle thus appears to be related to solar activity. Finally, we show that the grand minima and maxima occurred intermittently over the studied period, with clustering near lows and highs of the Hallstatt cycle, respectively.
Direct observations of sunspot numbers are available for the past four centuries, but longer time series are required, for example, for the identification of a possible solar influence on climate and ...for testing models of the solar dynamo. Here we report a reconstruction of the sunspot number covering the past 11,400 years, based on dendrochronologically dated radiocarbon concentrations. We combine physics-based models for each of the processes connecting the radiocarbon concentration with sunspot number. According to our reconstruction, the level of solar activity during the past 70 years is exceptional, and the previous period of equally high activity occurred more than 8,000 years ago. We find that during the past 11,400 years the Sun spent only of the order of 10% of the time at a similarly high level of magnetic activity and almost all of the earlier high-activity periods were shorter than the present episode. Although the rarity of the current episode of high average sunspot numbers may indicate that the Sun has contributed to the unusual climate change during the twentieth century, we point out that solar variability is unlikely to have been the dominant cause of the strong warming during the past three decades.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Solar cycle 23 witnessed the most complete set of observations of coronal mass ejections (CMEs) associated with the Ground Level Enhancement (GLE) events. We present an overview of the observed ...properties of the GLEs and those of the two associated phenomena, viz., flares and CMEs, both being potential sources of particle acceleration. Although we do not find a striking correlation between the GLE intensity and the parameters of flares and CMEs, the solar eruptions are very intense involving X-class flares and extreme CME speeds (average approx. 2000 km/s). An M7.1 flare and a 1200 km/s CME are the weakest events in the list of 16 GLE events. Most (80 %) of the CMEs are full halos with the three non-halos having widths in the range 167 to 212 degrees. 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. For accurate estimation of the CME height at the time of metric type II onset and GLE particle release, we estimated the initial acceleration of the CMEs using flare and CME observations. The initial acceleration of GLE-associated CMEs is much larger (by a factor of 2) than that of ordinary CMEs (2.3 km/sq s vs. 1 km/sq s). We confirmed the initial acceleration for two events for which CME measurements are available in the inner corona. 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 approx 3.09 R(sub s) (median: 3.18 R (sub s) ; range: 1.71 to 4.01 R (sub s) ) for well-connected events (source longitude in the range W20–W90). For poorly connected events, the average CME height at GLE particle release is ∼66 % larger (mean: 5.18 R (sub s) ; median: 4.61 R (sub s) ; range: 2.75–8.49 R (sub s) ). 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. On the other hand, the CME height at metric type II burst onset has no longitudinal dependence because electromagnetic signals do not require magnetic connectivity to the observer. For several events, the GLE particle release is very close to the time of first appearance of the CME in the coronagraphic field of view, so we independently confirmed the CME height at particle release. The CME height at metric type II burst onset is in the narrow range 1.29 to 1.8 R(sub s), with mean and median values of 1.53 and 1.47 R(sub s). 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 increase in CME speed between these two heights suggests an increase in Alfvenic Mach number from 2 to 3. The CME heights at GLE particle release are in good agreement with those obtained from the velocity dispersion analysis, including the source longitude dependence. We also discuss the implications of the delay of GLE particle release with respect to complex type III bursts by approx 18 min (median: 16 in; range: 2 to 44 min) for the flare acceleration mechanism. A similar analysis is also performed on the delay of particle release relative to the hard X-ray emission.
We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and group-sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an ...update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN and GN reconstruction methods while continuing to expand the historical sunspot-number database. Significant progress has been made on the database side while more work is needed to bring the various proposed SN and (primarily) GN reconstruction methods closer to maturity, after which the new reconstructions (or combinations thereof) can be compared with (a) “benchmark” expectations for any normalization scheme (e.g., a general increase in observer normalization factors going back in time), and (b) independent proxy data series such as F10.7 and the daily range of variations of Earth’s undisturbed magnetic field. New versions of the underlying databases for SN and GN will shortly become available for years through 2022 and we anticipate the release of the next versions of these two time series in 2024.
Context. Changes in solar irradiance and in its spectral distribution are among the main natural drivers of the climate on Earth. However, irradiance measurements are only available for less than ...four decades, while assessment of solar influence on Earth requires much longer records. Aims. The aim of this work is to provide the most up-to-date physics-based reconstruction of the solar total and spectral irradiance (TSI/SSI) over the last nine millennia. Methods. The concentrations of the cosmogenic isotopes 14C and 10Be in natural archives have been converted to decadally averaged sunspot numbers through a chain of physics-based models. TSI and SSI are reconstructed with an updated SATIRE model. Reconstructions are carried out for each isotope record separately, as well as for their composite. Results. We present the first ever SSI reconstruction over the last 9000 years from the individual 14C and 10Be records as well as from their newest composite. The reconstruction employs physics-based models to describe the involved processes at each step of the procedure. Conclusions. Irradiance reconstructions based on two different cosmogenic isotope records, those of 14C and 10Be, agree well with each other in their long-term trends despite their different geochemical paths in the atmosphere of Earth. Over the last 9000 years, the reconstructed secular variability in TSI is of the order of 0.11%, or 1.5 W m−2. After the Maunder minimum, the reconstruction from the cosmogenic isotopes is consistent with that from the direct sunspot number observation. Furthermore, over the nineteenth century, the agreement of irradiance reconstructions using isotope records with the reconstruction from the sunspot number by Chatzistergos et al. (2017, A&A, 602, A69) is better than that with the reconstruction from the WDC-SILSO series (Clette et al. 2014, Space Sci. Rev., 186, 35), with a lower χ2-value.
Aims. Miyake et al. (2012, Nature, 486, 240, henceforth M12) recently reported, based on 14C data, an extreme cosmic event in about AD775. Using a simple model, M12 claimed that the event was too ...strong to be caused by a solar flare within the standard theory. This implied a new paradigm of either an impossibly strong solar flare or a very strong cosmic ray event of unknown origin that occurred around AD775. However, as we show, the strength of the event was significantly overestimated by M12. Several subsequent works have attempted to find a possible exotic source for such an event, including a giant cometary impact upon the Sun or a gamma-ray burst, but they are all based on incorrect estimates by M12. We revisit this event with analysis of new datasets and consistent theoretical modelling. Methods. We verified the experimental result for the AD775 cosmic ray event using independent datasets including 10Be series and newly measured 14C annual data. We surveyed available historical chronicles for astronomical observations for the period around the AD770s to identify potential sightings of aurorae borealis and supernovae. We interpreted the 14C measurements using an appropriate carbon cycle model. Results. We show that: (1) The reality of the AD775 event is confirmed by new measurements of 14C in German oak; (2) by using an inappropriate carbon cycle model, M12 strongly overestimated the event’s strength; (3) the revised magnitude of the event (the global 14C production Q = (1.1 − 1.5) × 108 atoms/cm2) is consistent with different independent datasets (14C, 10Be, 36Cl) and can be associated with a strong, but not inexplicably strong, solar energetic particle event (or a sequence of events), and provides the first definite evidence for an event of this magnitude (the fluence >30 MeV was about 4.5 × 1010 cm-2) in multiple datasets; (4) this interpretation is in agreement with increased auroral activity identified in historical chronicles. Conclusions. The results point to the likely solar origin of the event, which is now identified as the greatest solar event on a multi-millennial time scale, placing a strong observational constraint on the theory of explosive energy releases on the Sun and cool stars.
Aims.Using a reconstruction of sunspot numbers stretching over multiple millennia, we analyze the statistics of the occurrence of grand minima and maxima and set new observational constraints on ...long-term solar and stellar dynamo models. Methods.We present an updated reconstruction of sunspot number over multiple millennia, from 14C data by means of a physics-based model, using an updated model of the evolution of the solar open magnetic flux. A list of grand minima and maxima of solar activity is presented for the Holocene (since 9500 BC) and the statistics of both the length of individual events as well as the waiting time between them are analyzed. Results.The occurrence of grand minima/maxima is driven not by long-term cyclic variability, but by a stochastic/chaotic process. The waiting time distribution of the occurrence of grand minima/maxima deviates from an exponential distribution, implying that these events tend to cluster together with long event-free periods between the clusters. Two different types of grand minima are observed: short (30–90 years) minima of Maunder type and long (>110 years) minima of Spörer type, implying that a deterministic behaviour of the dynamo during a grand minimum defines its length. The duration of grand maxima follows an exponential distribution, suggesting that the duration of a grand maximum is determined by a random process. Conclusions.These results set new observational constraints upon the long-term behaviour of the solar dynamo.
The most powerful explosions on the Sun – in the form of bright flares, intense storms of solar energetic particles (SEPs), and fast coronal mass ejections (CMEs) – 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 events 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 coronal 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.
Key Points
Solar, stellar, lunar, and terrestrial data need to be combined
Radionuclide data and ice cores do not tell us much about largest solar flares
We argue for an upper limit for the largest solar flares of about X30
Although sunspot-number series have existed since the mid-nineteenth century, they are still the subject of intense debate, with the largest uncertainty being related to the “calibration” of the ...visual acuity of individual observers in the past. A daisy-chain regression method is usually applied to inter-calibrate the observers, which may lead to significant bias and error accumulation. Here we present a novel method for calibrating the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900 – 1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds
S
S
defined such that the observer is assumed to miss all of the groups with an area smaller than
S
S
and report all the groups larger than
S
S
. Next, using a Monte-Carlo method, we construct a correction matrix for each observer from the reference data set. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot-group records reported by the observers for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series is provided as supplementary material in electronic form and displays secular minima around 1800 (Dalton Minimum) and 1900 (Gleissberg Minimum), as well as the Modern Grand Maximum of activity in the second half of the twentieth century. The uniqueness of the grand maximum is confirmed for the last 250 years. We show that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the eighteenth and nineteenth centuries in some reconstructions.