Next Article in Journal
Constraints on Metastable Dark Energy Decaying into Dark Matter
Previous Article in Journal
Primordial Axion Stars and Galaxy Halo Formation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers

1
Department of Statistics and Data Science, Northwestern University, Evanston, IL 60208, USA
2
Purple Mountain Observatory, Nanjing 210042, China
*
Author to whom correspondence should be addressed.
Universe 2024, 10(9), 370; https://doi.org/10.3390/universe10090370
Submission received: 29 July 2024 / Revised: 7 September 2024 / Accepted: 10 September 2024 / Published: 14 September 2024
(This article belongs to the Section Astroinformatics and Astrostatistics)

Abstract

:
A scientific method that proposes a value Y to estimate a target value ρ is often subject to some level of uncertainty. In the Bayesian framework, the level of uncertainty can be measured by the width of the 68 % interval, which is the range of the middle 68 % of the ranked ρ values sampled from the posterior distribution p ( ρ | Y ) . This paper considers Bayesian knowledge infusion (BKI) to reduce the uncertainty of the posterior distribution p ( ρ | Y ) based on additional knowledge that an event A happens. BKI is achieved by using a conditional prior distribution p ( ρ | A ) in the Bayes theorem, assuming that given the true ρ , its error-contaminated value Y is independent of event A. We use two examples to illustrate how to study whether or not it is possible to reduce uncertainty from 14C reconstruction (Y) of the annual sunspot number (SSN) ( ρ ) by infusing additional information (A) using BKI. Information (A) that SSN is from a year that has a Far Eastern record of naked eye sunspots is found to be not so effective in reducing the uncertainty. In contrast, information that SSN is from a year at a cycle minimum is found to be very effective, producing much narrower 68% intervals. The resulting Bayesian point estimates of SSN (the posterior medians of ρ ) are cross-validated and tested on a subset of telescopically observed SSNs that were unused in the process of Bayes computation.

1. Introduction

We propose a Bayesian knowledge infusion (BKI) method to incorporate additional information to improve a scientific measurement method that is subject to measurement errors. We are interested in a quantity ρ for a target study unit, which is estimated as Y by a scientific method. If the target study unit is known to satisfy an event A, then this additional information may be useful to improve the posterior distribution p ( ρ | Y ) to become a more precise conditional posterior distribution p ( ρ | Y , A ) . The current paper may be regarded as an exercise of applying BKI to improve the estimation of ρ = the annual sunspot numbers (SSN) of a target year in the remote past based on Y = the 14C reconstruction of ρ from tree samples. The additional information is, for example, A = the event that the target year has an ancient record of naked eye sunspot observation.
The evolvement of SSN ( ρ ) reflects the activity of the sun over the history. The sunspot numbers of longer time series are important for studying how solar dynamo works in maintaining solar magnetic field and possible solar influence on climates of the earth; see, e.g., [1]. Recently, [2] recovered the 11-year cycles of SSN since about 1000 A.D. using information of 14C from tree samples. SSN is negatively correlated with the 14C level in the tree samples since, in the lower activity period, the sun generates lower interplanetary magnetic field, which blocks less cosmic rays that increase 14C in the atmosphere and in the plants on the earth (see, e.g., [1]). However, the estimated SSN (Y’s) based on 14C are associated with a large amount of uncertainties, sometimes even being negative. For example, for the year 1566 AD, the reconstruction shows “ S N ” = −18.4 for the “Annual value of the reconstructed sunspot number” and 023.4 for the “1-sigma uncertainty of S N ”. See [3]1. In a Bayesian framework with a flat prior assuming no additional knowledge and Gaussian measurement error, this finding corresponds to estimating the true SSN ρ from a Gaussian posterior distribution with a negative mean 18.4 and a large standard deviation 23.4 , displaying a large amount of uncertainty.
We investigate whether additional knowledge can be used to improve the estimation of the SSN and reduce the uncertainty. For example, the most basic knowledge is that the SSNs have never been negative. They follow a certain distribution that can be derived from the telescopic observations in the later centuries. In the Bayesian statistical perspective, such a distribution on annual sunspot numbers over a collection of years may be regarded as a prior distribution that may be useful to summarize some kind of human knowledge, which can then be infused into the information from the 14C reconstruction, thereby much reducing the level of uncertainty. More details will be given in Section 2. Below, we will consider two kinds of possible additional knowledge and assess how useful they are.
During many of the years that are included in the 14C reconstruction, even before telescopes were invented, there are ancient records of naked-eye observations of the sunspots; see, e.g., a catalogue in [5] covering from 165 BC to AD 1918. We initially investigate whether knowing that a year has a record of naked-eye sunspots (NES) can reduce the level of uncertainty for the 14C reconstruction of the sunspot number of that year. For this purpose, we use a prior distribution of sunspot numbers “conditional on having an NES record”, derived from the observed SSNs with their standard errors during the years 1818–1918 for those years that also have at least one NES record. We found that the knowledge of having a NES record is not very effective for improving the estimation of the SSN. This complements a similar finding on the GSN (group sunspot numbers) by [6], who compared the distributions of the GSN with and without NES observations between 1848 and 1918 and did not find a big difference. The current paper directly shows that the NES information does not further improve the 14C reconstruction of the SSN either.
Next, we want to investigate whether the knowledge of being at a cycle minimum is useful in improving the accuracy of the 14C reconstructed SSN. This is an interesting question since the SSN at the cycle minimum is important for modeling the entire solar cycle, as pointed out by [7,8]2. This is investigated by using a prior distribution that is “conditional on being at a cycle-minimum”, which is constructed from the observed SSNs with their standard errors at the 19 cycle minimums between 1818 and 2023. The information on cycle minima from these years includes the distribution of the observed SSN values, which automatically includes the information that they cannot be negative and also how much higher than 0 they typically are. We notice that this prior distribution, “conditional on being at a cycle-minimum”, does reduce the uncertainty significantly further compared to the “unconditional prior” using all the annual sunspot numbers between 1818 and 2023, ignoring the cycle-minimum information.
To test whether these new estimates from the Bayesian methods are more reliable, we test the new estimates on the observed SSNs at the selected cycle minimums between 1700 and 1899. During this period of time, the 14C reconstruction of the sunspot numbers is also available and displays 12 plausible cycle minimums3, which can be compared with our Bayes estimates, regarding their differences from the cycle minimums of the observed SSNs4.
We find that the Bayes estimates (i.e., the posterior medians) are much closer to the observed SSNs, cutting the root mean square error to less than a half, compared to the original 14C reconstruction truncated above 0. In addition, the new and much narrower interval estimates did not cover the observed SSNs less often than the wider 0-truncated-14C intervals. Here, to be fair, we have used the idea of “cross-validation” in Section 2: We have carefully avoided using the observed SSN data between 1700 and 1899 themselves when obtaining our Bayes estimates and saved this period of data only for comparison at the end.
Here, the closeness between the estimated SSNs and the observed SSNs is measured by the “root mean square error”. The uncertainty is measured by the “average width of the 68% intervals”. These two concepts will be defined in Section 2. In summary, we see that additional information infused through the prior can be either useful or not useful. The Bayesian method we present here is useful in telling which knowledge may be useful and which is not, for reducing the uncertainty in SSN estimation. When the additional information reduces the uncertainty, it is important to use a cross-validation method if possible, to see if the now-narrower 68% intervals will still cover the observed SSNs for at least about 68% times, and also to check if the new point estimates are closer to the observed SSNs overall.
The observed SSN data we use is from “Source: WDC-SILSO, Royal Observatory of Belgium, Brussels”, downloadable from the website [9]. The SSN estimates are available from 1700–2023, with standard errors listed since 1818.

2. Materials and Methods

Bayesian methods found many different applications in astronomy, e.g., on studying galactic nuclei (e.g., [11]), gravitational wave (e.g., [12]), star clusters (e.g., [13]), stellar spectrum(e.g., [14]), and also on sunspots (e.g., [15,16,17]).
Our method adapts Bayesian bootstrap to infuse additional knowledge, improving the recent sunspot number reconstruction from 14C by [2], involving the remote past of more than 1000 years ago. Although there are several versions of Bayesian bootstrap, our method is closest to [18], but we have adapted his method to allow annual sunspot numbers to be not exactly known but subject to different levels of uncertainty, which can be related to fewer number of days of telescopic observation during earlier years.
Let ρ be the true annual sunspot number, say, for year 1566, when there were not yet telescopic records. We have Y, the 14C reconstructed value for ρ , subject to a large amount of uncertainty. The probability density function of Y conditional on ρ , p ( Y | ρ ) , is assumed to be known. We further assume it to be
p ( Y | ρ ) = e ( Y ρ ) 2 / ( 2 σ 2 ) / 2 π σ 2 ,
with a known σ .
Based on the Bayes theorem, the true value of ρ follows a posterior distribution conditional on Y,
p ( ρ | Y ) p ( Y | ρ ) p ( ρ ) ,
where p ( ρ ) is a prior distribution. Additional knowledge about ρ can be incorporated through the choice of the prior distribution. For example, in 1566, there exists a record of NES observation. Let such knowledge be represented by an event A. To incorporate this knowledge, we should compute
p ( ρ | Y , A ) p ( Y | ρ , A ) p ( ρ | A ) = p ( Y | ρ ) p ( ρ | A ) ,
where we assumed
p ( Y | ρ , A ) = p ( Y | ρ ) ,
Which means that given the true ρ , its error-contaminated value Y is independent of event A.
To obtain the “conditional prior” distribution p ( ρ | A ) (conditional on event A), we will use a data set D ( A ) , which satisfies event A and contains information about ρ . Ideally, we could look at years that are similar to 1566, where the true annual sunspot number ρ is known and A also happens (e.g., a NES record also exists). In such an ideal case, the data D ( A ) consists of these sample points of ρ . Then the histogram of these sample points of ρ ’s from these years forms an approximation of p ( ρ | A ) . Using δ x ( ρ ) to denote the point-mass distribution of ρ located at x, then we have
p ( ρ | A ) x D ( A ) δ x ( ρ ) / ( x D ( A ) 1 ) p ^ ,
which could be used to find the estimated quantiles of p ( ρ | A ) . The sample quantities based on p ^ have variations that need to be taken care of. This can be solved by the Bayesian bootstrap that resamples D ( A ) with replacement.
This ideal situation (knowing exactly the ρ ’s for some similar years), however, is not possible due to a lack of telescopic observation back then. Even in the years after 1700, telescopically observed SSN was not so frequent and the true annual value ρ is only estimated to be R, which we will refer to as the “observed SSN” and is subject to an uncertainty that varies with the number of observations in each year. This kind of uncertainty is measured by the standard error, which is needed to apply our method and is available in the dataset we use since 1818. There is also another kind of uncertainty in deciding which subset of years after 1818 are similar to 1566. To accommodate this kind of uncertainty as much as possible in a most conservative approach, we will not infuse any further knowledge on the similarity and allow all the years after 1818 to be considered as possibly similar to the year 1566. Of course, if one wishes, further knowledge on similarity could be added into the event A, sharpening the conditional prior distribution p ( ρ | A ) , at the cost of reducing the data set D ( A ) that is used to estimate p ( ρ | A ) .
In summary, we will use a data set D ( A ) , which consists of R i ’s from all years i S A , where S A consists of the years i { 1818 , , 2023 } for which A i happened (e.g., with NES records). We would like to approximate p ( ρ | A ) by the sample points in D ( A ) . However, these sample points R i are subject measurement errors. Assume that they follow independent normal distributions.
p ( R i | ρ i ) = e ( R i ρ ) 2 / ( 2 s i 2 ) / 2 π s i 2 ,
With a known s i . This same formula will also give the density of p ( ρ i | R i ) approximately, assuming a relatively flat prior of ρ i . This approximation is valid whenever the telescopic-observation uncertainty s i is much smaller than the standard deviation of the true SSN ρ i ’s over different years. To simplify computation, we will also assume similarly that p ( ρ i | R i , A i ) and p ( ρ i | R i , A i , Y i ) follow the righthand side of Equation (6). This allows us to simulate the true SNN ρ i based on the telescopically observed SSN R i alone, if such information is available for year i, and ignore other information. This is usually a good approximation if other information being ignored is much weaker than R i for the purpose of determining ρ i . If not, it may still be regarded as a conservative approach, allowing ρ i to be more uncertain than it really is.
Then we could do the following five steps:
Step I:
Sample the “true SNN”s ρ i from p ( ρ i | R i ) = e ( ρ i R i ) 2 / ( 2 s i 2 ) / 2 π s i 2 independently, for all i S A .
Step II:
Then the ρ i ’s are bootstrapped in a renormalized way to form a possible empirical sample ρ i ’s of p ( ρ | A ) , which is now estimated to be the “predictive” distribution p ( ρ | ρ i , i S A ) for the “future” ρ (in a “future” year 1566), from the simulated “true" SSN data in Step I, using “past” information (from “historical” data ( R i ) i S A ). This is done as follows:
IIa:
{ ρ i , i S A } ’s are bootstrap resampled as { w i , i S A } ’s, and sample standard deviation s w is obtained.
IIb:
{ ρ i , i S A } ’s are bootstrap resampled again as { g i , i S A } ’s, and sample mean g ¯ is obtained.
IIc:
{ ρ i , i S A } ’s are bootstrap resampled again as { f i , i S A } ’s, and we compute
ρ i = ρ ¯ + s ρ ( f i g ¯ ) / s w .
where ρ ¯ is the sample mean, and s ρ is the sample standard deviation for { ρ i , i S A } .
These sub-steps marginalize over uncertainties in the means and variances given the sample of ρ i ’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 ρ ( t ) is sampled from this empirical sample ρ i ’s, weighted by
p ( Y | ρ i ) = e ( Y ρ i ) 2 / ( 2 σ 2 ) / 2 π σ 2 ,
and converted to 0 if negative.
Step IV:
Steps I to III are repeated T times5. We then get ρ ( t ) for t = 1 , , T , which forms an approximate sample of the “posterior distribution” p ( ρ | Y , A , R i , i S A ) .
Step V:
Then we could get the quantiles of ρ from { ρ ( t ) : t = 1 , , T } , with details explained below.
The set { ρ ( t ) : t = 1 , , T } approximates the “posterior distribution”
p ( ρ | Y , A , R i , i S A ) ,
Which combines information Y (14C reconstructed SSN in 1566), A (the event that the year 1566 has a NES record), and information about how SSN relates to the event of NES records based on selected telescopic data of the observed SSNs { R i : i S A } .
Based on this sample { ρ ( t ) : t = 1 , , T } that approximates the “posterior distribution” p ( ρ | Y , A , R i , i S A ) , we can find a triplet of summary statistics
[ a , c , b ] = [ Q ( 0.16 ) , Q ( 0.50 ) , Q ( 0.84 ) ] ,
where Q ( q ) for any fraction q denotes the qth quantile of the sampled values of the true SSN according to the posterior distribution, so that 100q% of these sampled values are below Q ( q ) . For the 14C reconstruction with no additional knowledge infused, [ a , c , b ] = [ Y σ , Y , Y + σ ] . Its zero-truncated version is [ a , c , b ] = [ max { Y σ , 0 } , max { Y , 0 } , max { Y + σ , 0 } ] , which uses the fact that SSN is nonnegative. All these are different methods that describe the probability distribution of the unknown true SSN ρ conditional on the same 14C reconstruction estimate Y, with c being the point estimate of ρ and [ a , b ] being its 68% interval. The [ a , b ] is called a 68% interval since 68% of the ranked sampled values in the middle fall in this interval.

2.1. Bayesian Inference for Many Target Years

The previous procedure was for SSN in one particular target year, e.g., ρ = ρ j , where year j = 1566 . This can be repeated for the ρ j over many target year j’s. For each target year j, we will find a sample { ρ j ( t ) : t = 1 , , T } to approximate the posterior distribution p ( ρ j | Y j , A j , { R i , i S A j } ) . Then for each ρ j we have a triplet of summary statistics [ a j , c j , b j ] = [ Q j ( 0.16 ) , Q j ( 0.50 ) , Q j ( 0.84 ) ] , where Q j ( q ) for any fraction q denotes the qth quantile of the sampled values { ρ j ( t ) : t = 1 , , T } of the true SSN ρ j according to the posterior distribution p ( ρ j | Y j , A j , { R i , i S A j } ) , so that 100q% of these sampled values are below Q j ( q ) . For the 14C reconstruction with no additional knowledge infused, [ a j , c j , b j ] = [ Y j σ j , Y j , Y j + σ j ] . Its zero-truncated version is [ a j , c j , b j ] = [ max { Y j σ j , 0 } , max { Y j , 0 } , max { Y j + σ j , 0 } ] , which uses the fact that SSN is nonnegative. All these are different methods that describe the probability distribution of the unknown true SSN ρ j conditional on the same 14C reconstruction estimate Y j , with c j being the point estimate of ρ j and [ a j , b j ] being its 68% interval.

2.2. Cross Validation

We can compare several different methods for inference about the true annual SSN ρ j over a set of years j S , e.g., S = { 1566 , 1567 , , } . To evaluate the success of the point estimate c j , we would ideally like to see how far it is from the true ρ j , but knowing ρ j would be very difficult in the pre-telescopic years. However, in the post-telescopic years, we can use the “observed SSN” R j , which estimates the value of ρ j from telescopic observations, as its good approximation. We could then compute the “root mean square error”, which is the square root of the average of the ( c j R j ) 2 ’s over all the years j S :
j S ( c j R j ) 2 / j S 1 .
To be objective in prediction, it is preferable for c j (and similarly for [ a j , b j ] ) to be estimated by a “cross-validation” method that w i t h h o l d s the available telescopic information R j for each year j; e.g., computation of c j (and [ a j , b j ] ) uses the 14C reconstruction Y j for year j, but not R j itself.
To evaluate the success of a method for computing the interval estimate [ a j , b j ] , we compute its “empirical coverage” probability of covering R j , i.e., to see how often we find [ a j R j b j ] over all the years j S :
j S I [ a j R j b j ] / j S 1 ,
where I [ e v e n t ] = 1 if e v e n t happens or = 0 if otherwise. We do not wish this empirical coverage probability to be much lower than the intended coverage probability 68 % .
We also care about the “width” of such intervals [ a j , b j ] , defined as | b j a j | . We prefer intervals with “narrower” widths, whose overall behavior is summarized as the “average interval width”, which is the average of | b j a j | ’s over j S :
j S | b j a j | / j S 1 .
The “average interval width” can be used to assess the usefulness of additional knowledge that is infused into the Bayes method in order to improve inference about the true SSN given its 14C reconstruction. Useful additional knowledge infused through the prior distribution will cut down the uncertainty, and make the intervals | b j a j | ’s having a smaller “average interval width”, without hurting the “empirical coverage”.

3. Results

3.1. Using NES Information

Using the year 1566 as an example, when there exists a NES record that year, the 14C reconstruction of the SSN provides a 68% interval [ a , b ] = [ 41.8 , 5.0 ] . (The year 1566 is chosen only for illustration; there are many other years with NES records and 14C reconstructions, which could have been used instead.) We use BKI to combine these two pieces of information.
The telescopically observed SSNs with known standard errors for the 22 years between 1818 and 1918, when NES records exist, are used in a prior distribution to update the 68 % interval obtained from the 14C method, [ a , b ] = [ 41.8 , 5.0 ] , to a Bayesian 68 % interval, [ a , b ] = [ 0.0 , 28.1 ] , according to the posterior distribution. There was a reduction in the width of the interval. However, if the 14C-based interval is zero-truncated to [ a , b ] = [ 0 , 5.0 ] , then combining the telescopic SSN data from NES records actually increases the uncertainty for this particular year. If all the telescopically observed SSNs between 1818 and 2023 are used in the prior distribution, whether or not NES records exists, then the Bayesian posterior 68 % interval actually narrows down slightly, to [ a , b ] = [ 1.6 , 22.4 ] . These imply that NES records provide little additional information. This is consistent with the findings of [6], who show that the distributions of GSN do not differ much between years with or without NES records.
When this is repeated on all the 89 years before 1818 when there is both a 14C reconstruction and an NES record (see Figure 1), we found that the average width of the 68% intervals is { 83.6 , 72.0 , 63.2 , 61.6 } , respectively, for the 14C method, the 14C method with zero-truncation, the Bayesian method using all post-1818 years’ telescopic SSN information, and the Bayesian method using post-1818 telescopic SSN information from only the NES years. There are some small improvements (about 12%) in reducing the uncertainty caused by using the Bayesian methods; however, after using the SSN information from the telescope observations, adding that the year has an NES causes very little further reduction (less than 3%). All these confirm that NES records provide little additional information in estimation of annual SSN, after using the Bayesian method. The combination of using the Bayesian method AND the NES does lead to a moderate decrease in average interval width (about 15%), compared to the 14C method with 0-truncation. This leads us to use cross-validation to check if these slightly narrower intervals cover the observed SSNs less often.

Cross-Validation of the Bayesian Results Using NES Information

The observed SSN (with or without knowing its standard error) is known since 1700. There are 27 NES years since then. We use the later 13 NES years (after 1862) to drive the prior distribution used in the Bayesian method, which incorporates the knowledge that these are NES years. The Bayesian estimates using NES information in later years are then tested on the observed SSNs in the earlier 14 NES years (before 1862). We found that the root mean square error between the Bayesian point estimates (the medians) and the observed SSNs is 50.2, which is only slightly less (less than 8%) than the root mean square error between the 0-truncated 14C estimates and the observed SSNs, which is 54.4. The difference is not statistically significant since out of these 14 NES years, for exactly 7 times the observed SSNs are closer to the Bayesian estimates than to the 0-truncated 14C estimates.
We found that the average interval widths (for the earlier 14 NES years) for the Bayesian method with NES information is 64.4, compared to 75.3 for average interval width for the 0-truncated 14C method. This represents about 14% reduction. The narrower intervals did not reduce the empirical coverage. The Bayesian intervals cover the observed SSN 11/14 times compared to 10/14 times for the empirical coverage of the 0-truncated 14C intervals. However, the reduction of the interval width is not statistically significant, since for 7 times out of these 14 NES years, the Bayesian intervals are wider than the 0-truncated 14C intervals.
In summary, there is some indication that NES information plus the use of the Bayes method may lead to moderately narrower interval estimates compared to the 14C method with 0-truncation, without hurting the ability to cover the underlying SSN, but we did not find the supporting evidence to be statistically significant.

3.2. Using Cycle Minimum Information

3.2.1. Narrowing of 68% Intervals

There are many approximately 11-year cycles of annual SSN since about 1000 years ago, identified by Usoskin and Solanki et al. (2021) using the 14C reconstruction method. The locations of cycle maximums and minimums are only approximate, but are found to be typically within ± 2 years comparing to recent telescopic SSN data. We are interested in whether the information of a cycle minimum can help reduce the large amount of uncertainty in the reconstructed SSN. For example, the year 1363 is a cycle minimum, with 74 ± 58.9 SSN, that is bracketed by two high peaks, one at year 1369 [with SSN in 199.9 ± 80.5 ], and one at year 1357 [with SSN in 206.5 ± 89.3 ], which are 12 years apart. The actual cycle minimum may not be exactly located at 1363, but we would like to improve the estimation of SNN at this cycle minimum around 1363, regardless of in which year exactly it is located6.
The BKI method combines information from the wide 14C interval 74 ± 58.9 with the telescopically observed SSN information from 19 cycle minimums between 1818 and 2023, using their standard errors. The resulting 16%, 50%, and 84% quantiles of the posterior distribution for the SSN are found to be [ a , c , b ] = [ 3.6 , 11.0 , 18.1 ] , respectively. Therefore, the Bayesian interval has a width b a = 18.1 3.6 = 14.5 . In contrast, the 14C method gives an interval width of (74 + 58.9) − (74 − 58.9) = 117.8. We see that the Bayesian interval becomes much narrower (about 1:8) in comparison.
In order to assess whether the information of cycle minimum is influential to the Bayesian method, we incorporated the telescopic information on SSNs for all years without restricting to cycle minimums using their standard errors. This time the Bayesian posterior 0.16, 0.50, and 0.84 quantiles become [ a , c , b ] = [ 16.7 , 60.0 , 109.2 ] , whose range is no longer so much narrower than the 14C interval (about 4:5 in comparison this time). This shows that information about cycle minimum was highly influential in reducing the uncertainty, in contrast to the information of NES records in the previous analysis.
Side-to-side boxplots in Figure 2 confirm that the same phenomenon appears over all the cycle minimum years from the 14C reconstruction before 1818. In addition, the average widths for the 68% intervals are 66.8 for the 14C interval, 30.9 for the 0-truncated 14C interval, 36.5 for the Bayesian interval using telescopic SSN information from all years after 1818, and 12.9 for the Bayesian interval using telescopic SSN information after 1818 only at the cycle minima. Therefore, the cycle minimum information is very useful in reducing the uncertainty in 14C reconstruction of the SSN and reduces the average 68% interval width to less than a half. The usefulness of the cycle minimum information is related to the fact that the SSNs in the cycle minimum years are much smaller than the SSNs in other years, which influence the posterior distribution to focus on small and positive SSN values that allow much less uncertainty. In contrast, in the previous exercise, the SSNs in the years with NES records were not so much different from the SSNs in other years, as found in Usoskin et al. (2015) [6]. That is why the NES information was not so useful.
The results of the Bayesian reconstruction for SSN combining 14C information and cycle minimum information are listed in comparison with the 14C results in Table 1.
In Figure 3, the Bayesian point estimates and 68% intervals for SSNs at plausible cycle minimums are plotted over the years and compared to the 14C reconstruction and to the telescopically observed SSNs in the later years.

3.2.2. Cross-Validation of the Bayesian Results Using Cycle Minima Information

An important question is: will the much narrowed Bayesian intervals miss the observed SSN by a bigger amount, or miss more often?
To test whether the observed SSN in a year is missed by the Bayesian interval, we would like to use the idea of cross-validation and to make sure that this observed SSN to be tested is not used in the prior distribution when computing the Bayesian interval. We use cycle minimum years before 1899 to test if their SSNs are missed since after 1899, the 14C reconstruction data are not available for use in the Bayesian method. This implies that we should only use the data after 1899 to construct the prior distribution for computing the Bayes intervals.
For the SSNs to be tested, we select 11 cycle minimum years for the observed SSNs between 1700 and 1899. For each of these 11 years, a nearest minimum according to the 14C reconstruction is found to be a plausible cycle minimum in the sense of note 3. The SSN cycle-minimum year 1843 is special, since two 14C-plausible-cycle-minimum years are its equally distant nearest neighbors, each being 4 years apart. We take the average for these two nearest neighbor cycle minimums according to the 14C reconstructions.
For the Bayes method, we use the non-overlapping telescopic SSN data from the 12 cycle minimums after 1899 to derive the posterior distribution for each of the 11 cycle minimums before 1899 and compute its 0.16, 0.50, and 0.84 quantiles. The 0.16, 0.50, and 0.84 quantiles are also obtained for the 14C reconstruction method. The details of these quantiles can be seen in Table 2, where the actual telescopic SSN data are also listed for comparison.
In Figure 4, the Bayesian point estimate and 68% intervals for the SSNs at these 11 cycle minimums are plotted over the years before 1899 and compared to the 14C reconstruction and to the unused information on telescopic SSNs. We can see that the Bayesian reconstruction tends to have fewer uncertainties and tends to be closer to the telescopic SSNs.
When we obtain the 14C-based point estimates and 68% intervals, we assume that the uncertainty related to 14C reconstruction is normally distributed, which is alluded to in the discussion of Section 3.4 in [2]. The intervals of the Bayesian method are obtained non-parametrically and may not be symmetric. Now we check the “empirical coverage” of the 68% intervals and the root mean square errors of the corresponding point estimates based on Section 2. For each year j = 1 , , 11 of this set of 11 years, we compute the point estimate and the 68% interval for the Bayesian method and the 14C method with and without 0-truncated, respectively. The 68% intervals based on 14C reconstruction cover the observed SSN 6/11 ≈ 54.5% times, with an average interval width of 71.3 or 41.4 after applying 0-truncation. The 68% intervals from the Bayesian method using the cycle minimum information cover the observed SSN 9/11 ≈ 81.8% times, despite having a much narrower mean interval width of only about 13.4 . The empirical coverage 81.8% is higher than the nominal coverage 68%. This is not surprising, since the empirical coverage, even for a perfect 68% interval over only 11 trials, has a large standard deviation of 0.68 ( 1 0.68 ) / 11 14 % . The higher than-nominal empirical coverage may also be caused by our computation method in Section 2 using some conservation approximations.
If we compare the point estimates to the observed SSNs for these 11 years, the point estimate max { Y j , 0 } ’s based on the 0-truncated 14C reconstruction achieves a root mean square error of 32.4, while the Bayesian estimate (the posterior median Q j ( 0.50 ) ) achieves a root mean square error of about 5.0, only about 1:6 as large.
The absolute differences of these point estimates from the observed SSN R j ’s are compared by taking pairwise differences:
| Q j ( 0.50 ) R j | | max ( 0 , Y j ) R j | , j = 1 , , 11 .
Only 2/11 of such differences are positive. The chance of observing 2 or less such positive differences randomly out of 11 trials is the probability of a Binomial(11, 0.5) random variable being 0, 1, 2, which is only (1 + 11 + 11 * 10/(2 * 1))/211, about 0.033. We have relatively strong statistical evidence that the Bayesian estimates are closer to the telescopically-observed SSN values, compared to the zero-truncated 14C estimates.

4. Conclusions

We apply a Bayesian method to infuse additional information into a “target” measurement result, here being the estimated SSN from the 14C reconstruction. An example of additional information that is not very effective for improving the 14C reconstruction is that a certain year is included in the Far Eastern NES catalog. In contrast, an example of very effective additional information for improving the 14C reconstruction is that a certain year is a cycle minimum. This method may be adapted to other applications for reducing uncertainty with effective additional information.
The BKI method requires choosing a prior distribution properly. In our application, it involves the decision about which years with telescopically observed SSNs are regarded as “similar years” that are similar to a “target year”. The pre-conditions themselves are easy to check, e.g., which of the years in the telescopic era have NES records, or which years are cycle minimums. When the target year satisfies a pre-condition, we can select similar years in the telescopic era that satisfy the same pre-condition. However, any subset of these years could be used, and when we change the subset of “similar years” used in the prior distribution for BKI, the results will differ. For example, the target year 1796 is a cycle minimum according to the 14C reconstruction, and it appears both in Table 1 and Table 2. The prior distribution used for Table 1 uses all 19 telescopically observed cycle minimums between 1818 and 2023. The prior distribution for Table 2 only uses 12 telescopically observed cycle minimums after 1899. This causes the resulting (0.16, 0.50, 0.84) quantiles of the posterior distributions to change from (1.36, 6.18, 13.45) to (1.47, 5.66, 12.77). Although the change is minor in this case, one can imagine that bigger changes could happen when a different subset of years is used in the prior distribution.
The following are some recommendations:
  • 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

W.J. designed the methodology and did the computations. H.J. provided the research direction and suggestions on writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. W.J. is partially supported by a discretionary fund from the Department of Statistics and Data Science, Northwestern University, for covering some related traveling expenses.

Data Availability Statement

Data derived from public domain resources. See related links and references cited in the main text of the paper.

Acknowledgments

We thank the anonymous reviewers for providing comments that are useful for improving this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

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., 76.3 σ ), 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 T = 600 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 ± 2 years of the cycle minimum from the 14C reconstruction.

References

  1. 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]
  2. 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]
  3. 1000-Year Sunspot Series (Usoskin+, 2021). Available online: http://cdsarc.u-strasbg.fr/ftp/J/A+A/649/A141/ (accessed on 11 September 2024).
  4. 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]
  5. 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]
  6. 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]
  7. Hathaway, D.H.; Wilson, R.M.; Reichmann, E.J. The shape of the sunspot cycle. Sol. Phys. 1994, 151, 177–190. [Google Scholar] [CrossRef]
  8. 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]
  9. 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).
  10. 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).
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. Travaglini, G. Bayesian Methods for Reconstructing Sunspot Numbers Before and During the Maunder Minimum. Sol. Phys. 2017, 292, 23. [Google Scholar] [CrossRef]
  17. 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]
  18. 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).
  19. 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).
Figure 1. Bayesian point estimates and 68% intervals (symbolized by red dots with vertical error bars) for the SSNs, computed using post-1818 SSNs at years with NES observations, plotted over the years before 1818, and compared to the 14C reconstruction (solid curve with shade for 68% uncertainty) and to the telescopically observed SSNs (dashed curve).
Figure 1. Bayesian point estimates and 68% intervals (symbolized by red dots with vertical error bars) for the SSNs, computed using post-1818 SSNs at years with NES observations, plotted over the years before 1818, and compared to the 14C reconstruction (solid curve with shade for 68% uncertainty) and to the telescopically observed SSNs (dashed curve).
Universe 10 00370 g001
Figure 2. Information of cycle minimums does further reduce uncertainty of SSN estimated from 14C construction. Each boxplot is from 28 widths of the 68% intervals for the 28 years of plausible cycle minima found from the 14C reconstruction (listed in Table 1). For the boxplot above ‘14C’, the interval width for year j is 2 σ j ; for the boxplot above ‘0-truncated 14C’, the interval width for year j is max { Y j + σ j , 0 } max { Y j σ j , 0 } , where Y j is the 14C reconstruction and σ j is the 1-sigma uncertainty. For both the boxplots above ‘Bayes’ and above ‘Bayes with cycle min’, each interval width is computed as the difference of the 84% quantile and the 16% quantile of the posterior distribution. For ‘Bayes with cycle min’, the prior distribution uses the telescopic observations of annual SSNs, only at years of cycle minima, between years 1818 and 2023. For ‘Bayes’, the prior distribution uses all the telescopic observations of annual SSNs, whether or not at cycle minima, between years 1818 and 2023.
Figure 2. Information of cycle minimums does further reduce uncertainty of SSN estimated from 14C construction. Each boxplot is from 28 widths of the 68% intervals for the 28 years of plausible cycle minima found from the 14C reconstruction (listed in Table 1). For the boxplot above ‘14C’, the interval width for year j is 2 σ j ; for the boxplot above ‘0-truncated 14C’, the interval width for year j is max { Y j + σ j , 0 } max { Y j σ j , 0 } , where Y j is the 14C reconstruction and σ j is the 1-sigma uncertainty. For both the boxplots above ‘Bayes’ and above ‘Bayes with cycle min’, each interval width is computed as the difference of the 84% quantile and the 16% quantile of the posterior distribution. For ‘Bayes with cycle min’, the prior distribution uses the telescopic observations of annual SSNs, only at years of cycle minima, between years 1818 and 2023. For ‘Bayes’, the prior distribution uses all the telescopic observations of annual SSNs, whether or not at cycle minima, between years 1818 and 2023.
Universe 10 00370 g002
Figure 3. Bayesian point estimates and 68% intervals (symbolized by red dots with vertical errorbars) of the SSNs, incorporating telescopic data after 1818, at plausible cycle minimums before 1818, plotted over the years, compared to the 14C reconstruction (solid curve with shade for 68% uncertainty) and to the telescopic SSNs in the later years (dashed curve).
Figure 3. Bayesian point estimates and 68% intervals (symbolized by red dots with vertical errorbars) of the SSNs, incorporating telescopic data after 1818, at plausible cycle minimums before 1818, plotted over the years, compared to the 14C reconstruction (solid curve with shade for 68% uncertainty) and to the telescopic SSNs in the later years (dashed curve).
Universe 10 00370 g003
Figure 4. Bayesian reconstruction for SSNs at 11 cycle minimums before 1899, using post-1899 telescopic data at the cycle minimums, plotted over the years, compared to the 14C reconstruction (solid curve with shade for 68% uncertainty) and to the unused information on telescopic SSNs (dashed curve). The red dots and the vertical error bars display the Bayesian point estimate (posterior median) and the nominal 68% intervals. The horizontal error bars extend to the year of the nearest plausible cycle minimum according to the 14C reconstruction, as defined in note 3, with ties averaged. Therefore, the lengths of the horizontal error bars reflect the errors of placements of the cycle minimums by the 14C reconstruction.
Figure 4. Bayesian reconstruction for SSNs at 11 cycle minimums before 1899, using post-1899 telescopic data at the cycle minimums, plotted over the years, compared to the 14C reconstruction (solid curve with shade for 68% uncertainty) and to the unused information on telescopic SSNs (dashed curve). The red dots and the vertical error bars display the Bayesian point estimate (posterior median) and the nominal 68% intervals. The horizontal error bars extend to the year of the nearest plausible cycle minimum according to the 14C reconstruction, as defined in note 3, with ties averaged. Therefore, the lengths of the horizontal error bars reflect the errors of placements of the cycle minimums by the 14C reconstruction.
Universe 10 00370 g004
Table 1. Table of 16%, 50%, and 84% quantiles of SSNs at years of plausible cycle minimums from 14C reconstruction (last three columns) and from the Bayesian method incorporating the cycle minimum information from later telescope observations (the 2nd to the 4th columns).
Table 1. Table of 16%, 50%, and 84% quantiles of SSNs at years of plausible cycle minimums from 14C reconstruction (last three columns) and from the Bayesian method incorporating the cycle minimum information from later telescope observations (the 2nd to the 4th columns).
Year Q ( 0.16 ) B Q ( 0.50 ) B Q ( 0.84 ) B Q ( 0.16 ) 14 C Q ( 0.50 ) 14 C Q ( 0.84 ) 14 C
17961.366.1813.45−51.5−29.9−8.3
17834.4811.818.4863.4112.7162.0
17712.417.8615.41−43.9−15.512.9
17633.5310.0417.846.761.7116.7
17342.117.1115.05−70.8−40.1−9.4
16093.088.4217.29−29.04.037.0
15842.858.4316.33−33.7−4.325.1
13752.417.7515.43−54.6−24.55.6
13633.5610.9718.0615.174.0132.9
12752.267.014.61−51.7−23.54.7
12502.798.4515.86−42.7−9.523.7
12413.479.6217.4−27.015.257.4
12102.578.4816.18−45.6−11.622.4
11992.057.014.66−59.1−30.9−2.7
11884.139.8517.1−16.228.573.2
11630.925.2512.67−57.5−35.0−12.5
11541.096.013.26−58.5−35.4−12.3
11423.038.5116.78−55.8−13.429.0
11313.6211.0817.8618.364.6110.9
11243.178.8116.57−20.88.036.8
11154.110.6217.635.238.070.8
10932.88.5916.16−39.0−6.326.4
10743.298.8615.98−34.5−3.627.3
10201.095.4911.78−60.1−39.1−18.1
10081.86.8114.85−56.2−29.4−2.6
9971.65.512.43−77.3−51.8−26.3
9881.896.5114.53−77.1−46.0−14.9
9761.075.411.5−86.4−62.0−37.6
Table 2. Table of observed cycle minimum SSN (the 5th column) and its year (the 6th column), 16%, 50%, and 84% quantiles of SSN at years of plausible cycle minimums from 14C reconstruction (last three columns), and from the Bayesian method using the cycle minimum information (the 2nd to the 4th columns). The first column is the year of the nearest-neighbor plausible cycle minimum(s) (defined according to note 3, with ties being averaged), as identified by the 14C method. These 11 cycle minimums are taken from between 1700 and 1899, when both telescope observations of SNN and the 14C reconstructions in [2] are available. The values from the Bayesian method have some minor differences from those on Table 1, since here only non-overlapping cycle minimums after 1899 are incorporated in the prior distribution for Bayesian computation.
Table 2. Table of observed cycle minimum SSN (the 5th column) and its year (the 6th column), 16%, 50%, and 84% quantiles of SSN at years of plausible cycle minimums from 14C reconstruction (last three columns), and from the Bayesian method using the cycle minimum information (the 2nd to the 4th columns). The first column is the year of the nearest-neighbor plausible cycle minimum(s) (defined according to note 3, with ties being averaged), as identified by the 14C method. These 11 cycle minimums are taken from between 1700 and 1899, when both telescope observations of SNN and the 14C reconstructions in [2] are available. The values from the Bayesian method have some minor differences from those on Table 1, since here only non-overlapping cycle minimums after 1899 are incorporated in the prior distribution for Bayesian computation.
Yr14C Q ( 0.16 ) B Q ( 0.50 ) B Q ( 0.84 ) BObsSSNYrObs Q ( 0.16 ) 14C Q ( 0.50 ) 14C Q ( 0.84 ) 14C
18880.75.3312.4110.41889−62.1−39.3−16.5
18781.876.7514.25.71878−48.9−24.10.7
18653.219.4817.4713.91867−9.326.862.9
18562.648.6616.298.21856−27.49.646.6
1847, 18392.658.7517.2718.11843−16.122.360.7
18293.159.0917.0513.41833−33.814.562.8
17961.475.6612.776.81798−51.5−29.9−8.3
17834.1411.6818.6317.0178463.4112.7162.0
17712.237.715.9211.71775−43.9−15.512.9
17633.1710.0818.119.017666.761.7116.7
17342.037.2914.488.31733−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.

Share and Cite

MDPI and ACS Style

Jiang, W.; Ji, H. Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers. Universe 2024, 10, 370. https://doi.org/10.3390/universe10090370

AMA Style

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 Style

Jiang, 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 Style

Jiang, W., & Ji, H. (2024). Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers. Universe, 10(9), 370. https://doi.org/10.3390/universe10090370

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop