Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers
Abstract
:1. Introduction
2. Materials and Methods
- Step I:
- Sample the “true SNN”s from independently, for all .
- Step II:
- Then the ’s are bootstrapped in a renormalized way to form a possible empirical sample ’s of , which is now estimated to be the “predictive” distribution for the “future” (in a “future” year 1566), from the simulated “true" SSN data in Step I, using “past” information (from “historical” data ). This is done as follows:
- IIa:
- ’s are bootstrap resampled as ’s, and sample standard deviation is obtained.
- IIb:
- ’s are bootstrap resampled again as ’s, and sample mean is obtained.
- IIc:
- ’s are bootstrap resampled again as ’s, and we compute
These sub-steps marginalize over uncertainties in the means and variances given the sample of ’s, which would lead to a correct predictive t-distribution (e.g., [19], first formula on p16) in a special case involving normal distributions. - Step III:
- Then one is sampled from this empirical sample ’s, weighted by
- Step IV:
- Steps I to III are repeated T times5. We then get for , which forms an approximate sample of the “posterior distribution” .
- Step V:
- Then we could get the quantiles of from , with details explained below.
2.1. Bayesian Inference for Many Target Years
2.2. Cross Validation
3. Results
3.1. Using NES Information
Cross-Validation of the Bayesian Results Using NES Information
3.2. Using Cycle Minimum Information
3.2.1. Narrowing of 68% Intervals
3.2.2. Cross-Validation of the Bayesian Results Using Cycle Minima Information
4. Conclusions
- When we are not sure, we propose to use all years available that satisfy the pre-condition to account for the uncertainty in our knowledge as much as possible. This is the approach we take in the current paper.
- It may be safer to report a result together with the corresponding assumption, e.g., “Assuming that the target year 1796 is similar to the cycle minimum years after 1899, then the posterior quantiles for the SSN in year 1796 will be (1.47, 5.66, 12.77)”.
- Subject matter knowledge may help to tell if a subset of years are similar to the target year, or none of the years in the telescopic era is similar to the target year, and a scaling up or long-term de-trending (say) of the SSN values has to be done first before applying BKI. (So we may regard BKI as a mathematical framework that can be made useful with additional scientific knowledge).
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
1 | It is noted that this issue with negative sunspots does not affect the series by [1] or [4]. However, [1] reconstructed the solar activity decadely without resolving the 11-year cycles. [4] did not directly reconstruct the sunspot numbers but instead estimated the solar modulation potential, which, as [2] pointed out, characterizes the flux intensity of galactic cosmic rays and is not straightforward to convert into quantities useful for Sun-Earth relations. |
2 | Hathaway et al. (1994) found that in a parametric model of the SSN cycle involving a cubic polynomial rising phase followed by a Gaussian tail decline, two parameters are most important: the starting time and the amplitude. The amplitude parameter can be estimated well using only about 2–3 years of data after the start of the cycle, near the cycle minimum. This is why the placement and size of the cycle minimum are very important for fitting the whole cycle. We study the size of the cycle minimum in this paper but not the placement, since the latter is already determined quite well by the 14C reconstruction, according to Section 4.1.3 in [2]. |
3 | The “plausible cycle minimums” in this paper are minimums found from the 14C reconstruction of the SNN series, such that each minimum is bracketed by two adjacent maximums, one on each side, with neither one being too low after accounting for its standard error (i.e., ), and located a reasonable number of years apart (between 7–16). According to all the telescopically observed SSN data since 1700, the value 76.3 is the lowest of all cycle maximums, and all adjacent cycle maximums are located 7–16 years apart. |
4 | The cycle minimums of the observed SSNs are found from the local minimums in the observed SSN dataset from [9]. A flat local minimum at years 1711.5 and 1712.5 are averaged. Some minor local minimums that are at least 1 year away from any cycle minimum listed in the following website (before the current cycle) are omitted: [10]. |
5 | We used for all our computations, and np.random.seed(101) and random.seed(10) in Python codes. |
6 | We only examine the values of the sunspot minimums, since their times are reconstructed by the 14C method quite successfully already, according to Section 4.1.3 in [2]. They found by comparison with the direct sunspot series in later years, that the true year of the cycle minimum is usually located within years of the cycle minimum from the 14C reconstruction. |
References
- Solanki, S.K.; Usoskin, I.G.; Kromer, B.; Schüssler, M.; Beer, J. Unusual activity of the Sun during recent decades compared to the previous 11,000 years. Nature 2004, 431, 1084–1087. [Google Scholar] [CrossRef] [PubMed]
- Usoskin, I.G.; Solanki, S.K.; Krivova, N.A.; Hofer, B.; Kovaltsov, G.A.; Wacker, L.; Brehm, N.; Kromer, B. Solar cyclic activity over the last millennium reconstructed from annual 14C data. Astron. Astrophys. 2021, 649, A141. [Google Scholar] [CrossRef]
- 1000-Year Sunspot Series (Usoskin+, 2021). Available online: http://cdsarc.u-strasbg.fr/ftp/J/A+A/649/A141/ (accessed on 11 September 2024).
- Brehm, N.; Bayliss, A.; Christl, M.; Synal, H.A.; Adolphi, F.; Beer, J.; Kromer, B.; Muscheler, R.; Solanki, S.K.; Usoskin, I.; et al. Eleven-year solar cycles over the last millennium revealed by radiocarbon in tree rings. Nat. Geosci. 2021, 14, 10–15. [Google Scholar] [CrossRef]
- Yau, K.K.C.; Stephenson, F.R. A revised catalogue of Far Eastern observations of sunspots (165 BC to AD 1918). Q. J. R. Astron. Soc. 1988, 29, 175–197. [Google Scholar]
- Usoskin, I.G.; Arlt, R.; Asvestari, E.; Hawkins, E.; Käpylä, M.; Kovaltsov, G.A.; Krivova, N.; Lockwood, M.; Mursula, K.; O’Reilly, J.; et al. The Maunder minimum (1645–1715) was indeed a grand minimum: A reassessment of multiple datasets. Astron. Astrophys. 2015, 581, A95. [Google Scholar] [CrossRef]
- Hathaway, D.H.; Wilson, R.M.; Reichmann, E.J. The shape of the sunspot cycle. Sol. Phys. 1994, 151, 177–190. [Google Scholar] [CrossRef]
- Wilson, R.M.; Hathaway, D.H.; Reichmann, E.J. On the importance of cycle minimum in sunspot cycle prediction. NASA Tech. Publ. TP-3648 1996, 16, 1–11. [Google Scholar]
- Sunspot Data (WDC-SILSO, Royal Observatory of Belgium, Brussels). Available online: https://www.sidc.be/SILSO/DATA/SN_y_tot_V2.0.txt (accessed on 11 September 2024).
- Minima and Maxima of Sunspot Number Cycles (NOAA). Available online: https://www.ngdc.noaa.gov/stp/space-weather/solar-data/solar-indices/sunspot-numbers/cycle-data/table_cycle-dates_maximum-minimum.txt (accessed on 11 September 2024).
- Li, Y.R.; Wang, J.M.; Ho, L.C.; Du, P.; Bai, J.M. A Bayesian approach to estimate the size and structure of the broad-line region in active galactic nuclei using reverberation mapping data. Astrophys. J. 2013, 779, 110. [Google Scholar] [CrossRef]
- Fan, X.; Messenger, C.; Heng, I.S. A Bayesian approach to multi-messenger astronomy: Identification of gravitational-wave host galaxies. Astrophys. J. 2014, 795, 43. [Google Scholar] [CrossRef]
- Shao, Z.; Xie, X.; Chen, L.; Zhong, J.; Hou, J.; Lin, C.C. Bayesian Inference of Kinematics and Mass Segregation of Open Cluster. Int. Astron. Union. Proc. Int. Astron. Union 2015, 12, 265–266. [Google Scholar] [CrossRef]
- Kang, X.; He, S.Y.; Zhang, Y.X. A novel stellar spectrum denoising method based on deep Bayesian modeling. Res. Astron. Astrophys. 2021, 21, 169. [Google Scholar] [CrossRef]
- Yu, Y.; van Dyk, D.A.; Kashyap, V.L.; Young, C.A. A Bayesian analysis of the correlations among sunspot cycles. Sol. Phys. 2012, 281, 847–862. [Google Scholar] [CrossRef]
- Travaglini, G. Bayesian Methods for Reconstructing Sunspot Numbers Before and During the Maunder Minimum. Sol. Phys. 2017, 292, 23. [Google Scholar] [CrossRef]
- Velasco Herrera, V.M.; Soon, W.; Hoyt, D.V.; Muraközy, J. Group sunspot numbers: A new reconstruction of sunspot activity variations from historical sunspot records using algorithms from machine learning. Sol. Phys. 2022, 297, 8. [Google Scholar] [CrossRef]
- Hjort, N.L. Bayesian and Empirical Bayesian Bootstrapping. Preprint Series. Statistical Research Report 1991. Available online: https://www.duo.uio.no/bitstream/handle/10852/47760/1/1991-9.pdf (accessed on 11 September 2024).
- Murphy, K.P. Conjugate Bayesian Analysis of the Gaussian Distribution. Def 2007, 1, 16. Available online: https://www.cs.ubc.ca/~murphyk/Papers/bayesGauss.pdf (accessed on 11 September 2024).
Year | B | B | B | C | C | C |
---|---|---|---|---|---|---|
1796 | 1.36 | 6.18 | 13.45 | −51.5 | −29.9 | −8.3 |
1783 | 4.48 | 11.8 | 18.48 | 63.4 | 112.7 | 162.0 |
1771 | 2.41 | 7.86 | 15.41 | −43.9 | −15.5 | 12.9 |
1763 | 3.53 | 10.04 | 17.84 | 6.7 | 61.7 | 116.7 |
1734 | 2.11 | 7.11 | 15.05 | −70.8 | −40.1 | −9.4 |
1609 | 3.08 | 8.42 | 17.29 | −29.0 | 4.0 | 37.0 |
1584 | 2.85 | 8.43 | 16.33 | −33.7 | −4.3 | 25.1 |
1375 | 2.41 | 7.75 | 15.43 | −54.6 | −24.5 | 5.6 |
1363 | 3.56 | 10.97 | 18.06 | 15.1 | 74.0 | 132.9 |
1275 | 2.26 | 7.0 | 14.61 | −51.7 | −23.5 | 4.7 |
1250 | 2.79 | 8.45 | 15.86 | −42.7 | −9.5 | 23.7 |
1241 | 3.47 | 9.62 | 17.4 | −27.0 | 15.2 | 57.4 |
1210 | 2.57 | 8.48 | 16.18 | −45.6 | −11.6 | 22.4 |
1199 | 2.05 | 7.0 | 14.66 | −59.1 | −30.9 | −2.7 |
1188 | 4.13 | 9.85 | 17.1 | −16.2 | 28.5 | 73.2 |
1163 | 0.92 | 5.25 | 12.67 | −57.5 | −35.0 | −12.5 |
1154 | 1.09 | 6.0 | 13.26 | −58.5 | −35.4 | −12.3 |
1142 | 3.03 | 8.51 | 16.78 | −55.8 | −13.4 | 29.0 |
1131 | 3.62 | 11.08 | 17.86 | 18.3 | 64.6 | 110.9 |
1124 | 3.17 | 8.81 | 16.57 | −20.8 | 8.0 | 36.8 |
1115 | 4.1 | 10.62 | 17.63 | 5.2 | 38.0 | 70.8 |
1093 | 2.8 | 8.59 | 16.16 | −39.0 | −6.3 | 26.4 |
1074 | 3.29 | 8.86 | 15.98 | −34.5 | −3.6 | 27.3 |
1020 | 1.09 | 5.49 | 11.78 | −60.1 | −39.1 | −18.1 |
1008 | 1.8 | 6.81 | 14.85 | −56.2 | −29.4 | −2.6 |
997 | 1.6 | 5.5 | 12.43 | −77.3 | −51.8 | −26.3 |
988 | 1.89 | 6.51 | 14.53 | −77.1 | −46.0 | −14.9 |
976 | 1.07 | 5.4 | 11.5 | −86.4 | −62.0 | −37.6 |
Yr14C | B | B | B | ObsSSN | YrObs | 14C | 14C | 14C |
---|---|---|---|---|---|---|---|---|
1888 | 0.7 | 5.33 | 12.41 | 10.4 | 1889 | −62.1 | −39.3 | −16.5 |
1878 | 1.87 | 6.75 | 14.2 | 5.7 | 1878 | −48.9 | −24.1 | 0.7 |
1865 | 3.21 | 9.48 | 17.47 | 13.9 | 1867 | −9.3 | 26.8 | 62.9 |
1856 | 2.64 | 8.66 | 16.29 | 8.2 | 1856 | −27.4 | 9.6 | 46.6 |
1847, 1839 | 2.65 | 8.75 | 17.27 | 18.1 | 1843 | −16.1 | 22.3 | 60.7 |
1829 | 3.15 | 9.09 | 17.05 | 13.4 | 1833 | −33.8 | 14.5 | 62.8 |
1796 | 1.47 | 5.66 | 12.77 | 6.8 | 1798 | −51.5 | −29.9 | −8.3 |
1783 | 4.14 | 11.68 | 18.63 | 17.0 | 1784 | 63.4 | 112.7 | 162.0 |
1771 | 2.23 | 7.7 | 15.92 | 11.7 | 1775 | −43.9 | −15.5 | 12.9 |
1763 | 3.17 | 10.08 | 18.1 | 19.0 | 1766 | 6.7 | 61.7 | 116.7 |
1734 | 2.03 | 7.29 | 14.48 | 8.3 | 1733 | −70.8 | −40.1 | −9.4 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Jiang, W.; Ji, H. Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers. Universe 2024, 10, 370. https://doi.org/10.3390/universe10090370
Jiang W, Ji H. Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers. Universe. 2024; 10(9):370. https://doi.org/10.3390/universe10090370
Chicago/Turabian StyleJiang, Wenxin, and Haisheng Ji. 2024. "Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers" Universe 10, no. 9: 370. https://doi.org/10.3390/universe10090370
APA StyleJiang, W., & Ji, H. (2024). Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers. Universe, 10(9), 370. https://doi.org/10.3390/universe10090370