Journal of Astronomy and Space Sciences
The Korean Space Science Society
Research Paper

Analysis of Solar and Lunar Motions in the Seonmyeong Calendar

Ki-Won Lee†http://orcid.org/0000-0002-4097-1827
Daegu Catholic University, 13-13 Hayang-ro, Hayang-eup, Gyeongsan 3410, S. Korea
Corresponding Author Tel: +82-53-850-2573, Fax: +82-53-850-3800 E-mail: leekw@cu.ac.kr

© The Korean Space Science Society. All rights reserved. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: May 8, 2019; Revised: May 31, 2019; Accepted: Jun 1, 2019

Abstract

We investigate solar and lunar motions in the Seonmyeong (SM) calendar that was compiled by Xu, Ang of the Tang dynasty (A.D. 618-907) in China and used for 71 years from 822 to 892. This calendar was also used in Korea during the Goryeo dynasty (A.D. 918-1392) and in Japan for 823 years from 862 to 1684, the longest time among the three countries. Referring to historical documents of China, Korea, and Japan, we analyze the calendrical methods of calculating the daily apparent movements of the Sun and Moon in the SM calendar, which were considered their unequal motions, and compare the movements with the results of modern calculations for three periods in the Goryeo dynasty: 919, 1155, and 1392 years (i.e., the beginning, middle, and ending of the dynasty, respectively). We find that a quadratic equation was employed to obtain the daily movement of the Sun using physical quantities on the instant of each solar term, which was tabulated in its calendar book such as the Goryeosa (History of the Goryeo Dynasty). For quantitative analysis, we compute the mean absolute difference (MAD) of the daily apparent movement between the SM calendar and modern calculations and obtain 0.33, 0.30, and 0.31 arcmin for the periods of 919, 1155, and 1392 years, respectively. Meanwhile, we find relatively large MAD values in the daily movement of the Moon: 0.217, 0.284, and 0.240 degrees for each corresponding year. An interesting point is that the MAD value in the lunar motion shows the maximum in 1155 years, and is the minimum in the solar motion. In conclusion, we believe that this study will facilitate in the understanding of the SM calendar further, particularly in the calendrical methods of calculating sunrise, sunset, and eclipse times.

Keywords: historical astronomy; Seonmyeong calendar; ephemeris; Sun; Moon

1. INTRODUCTION

The Xuanming (Manifest Enlightenment) calendar was compiled by Xu, Ang (birth and death years are unknown) of the Tang dynasty (A.D. 618–907) in A.D. 822 and has been used for 71 years in China (Lee 1985). This calendar was also used in Korea and Japan, which was called the Seonmyeong and Senmyo calendars, respectively. According to the preface of the astronomical chapter of the Goryeosa (History of the Goryeo Dynasty), the Goryeo court (A.D. 918-1392) of Korea used the Seonmyeong calendar (hereinafter SM calendar) from the beginning of the dynasty by succeeding it from the Silla kingdom (B.C. 57–A.D. 935) to the end of the dynasty, particularly for calculating solar and lunar eclipses although the Shoushi (Granting the Seasons) calendar of the Yuan dynasty (A.D. 1271–1368) in China was enforced in 1281 (Seo 2012, Chang 2014, Jun 2014). Therefore, the SM calendar was utilized for 475 years in Korea, at the least. Meanwhile, this calendar was introduced to Japan in A.D. 862 through the Balhae kingdom (A.D. 698–926) of Korea and used for 823 years until 1684 (Kim 2012).

In this study, we analyzed the calendrical methods of solar and lunar motions in the SM calendar referring to the Xintangshu (New History of Tang Dynasty) of China, Goryeosa of Korea, and Seonmyoreki (Seonmyo Calendar) and Chokei-Senmyoreki-Sanpo (Calendrical Method for the Seonmyo Calendar of the Chokei, shortly Sanpo) of Japan. With regard to Japanese literature, the former is preserved in the National Library of Korea and is known to be made of a Japanese woodblock in 1644 but the author is unknown. Meanwhile, the latter was compiled by Ando, Yueki in 1654 (Soma 2014) and the Chokei (Changqing in Chinese) is the Chinese reign-title used for four years from 821 to 824.

Together with the Shoushi calendar, the SM calendar is a well-studied Chinese calendar, particularly by Japanese scholars. Sato (1968) reviewed the calendar in the study on the calendar history of Japan. Uchida (1975) utilized the Sanpo for identifying Japanese lunisolar calendar dates with Julian calendar dates. However, Soma (2014) found some mistakes in the daily movement of the Sun, which is tabulated in the Sanpo. According to our examination, the same mistakes exist in the Seonmyoreki preserved in Korea. ISOBDU (1998) compared the contents of the SM calendar recorded on various versions of the Goryeosa with the Sanpo. Qu (2008) studied Chinese calendars including the SM calendar in terms of mathematical techniques (see also Qian 1981, Martzloff 2016). In Korea, the SM calendar was first studied by the scholars of North Korea in 1964 when translating the Goryeosa written in Chinese characters into Korean (EDS 1991). Recently, SADU (2016) retranslated the astronomical chapter of the Goryeosa including detailed annotations and some calculation examples.

The SM calendar is known as a good Chinese calendar, particularly for the calculations of eclipses (Lee 1985). According to the work of Hahn et al. (2016), a quadratic equation of type x(x-a) was employed in the SM calendar to consider the horizontal parallax of the Moon and was utilized in subsequent Chinese calendars such as the Shoushi calendar. However, the accuracy of the SM calendar has not been analyzed through comparison with modern calculations. The purpose of this study is to investigate the accuracy of the solar and lunar motions by the SM calendar. The remainder of this paper is structured as follows. In Section 2, we introduce the units and astronomical constants in the SM calendar, which are used in this study. We analyze the numerical tables relating to solar and lunar motions listed in calendar books and compare with the results of modern calculations in Section 3. Finally, we summarize our findings in Section 4.

2. UNITS AND ASTRONOMICAL CONSTANTS

In the SM calendar, Il (Day), Bun (Minute) and Cho (Second) were used as the units of an hour (hereinafter, dH, minH, and secH, respectively) and Do (Degree), Bun (Arcminute), and Cho (Arcsecond) as those of an angle (hereinafter, dA, minA, and secA, respectively): 1 dH = 84 or 8400 minH = 672 or 67200 secH (i.e., 1 minH = 8 secH) and 1 dA = 84 or 8400 minA = 25200 or 2520000 secA (i.e., 1 minA = 300 secA). Because two scales were used in units of minH and minA (i.e., 84 and 8400 called the Gakbeob and Seonmyeong-Tongbeob, respectively), we distinguished them using the superscripts G (i.e., Gakbeob) and T (i.e., Seonmyeong-Tongbeob) such as minHG or minAT.

As in other Chinese calendars, the beginning of the calendrical calculation for a given year, Y, was to determine the instant of the winter solstice in its previous year, Y-1, called the Cheonjeong-Dongji (hereinafter, CJDJ) of year Y. In addition, although the calendrical epoch of the SM calendar was the instant of the CJDJ of 822 (i.e., instant of the winter solstice in 821), that of the CJDJ of a given year was obtained by counting from the Super Epoch (i.e., an ancient epoch where the beginnings of the solar, lunar, and sexagenary cycles coincided). The length between the Super Epoch and calendrical epoch of the SM calendar, Sangwon-Jeoknyeon, was 7070138 years. The abolition of the Super Epoch is a remarkable characteristic in the Shoushi calendar compared to its previous calendars (Sivin 2008).

Regarding the solar motion, the length of a tropical year called the Jangse was 3068055 minHT or 365 dH 2055 minHT, and a year was divided into 24 sub-seasons (hereinafter, 24 solar terms) in Chinese calendars. A method of determining the instants of 24 solar terms is to equally divide the length of the tropical year into 24 intervals, as mentioned in Section 3.1. In this case, the length of a solar term is 15 dH 1835 minHT 5 secH and this value was called the Jungjeol. The celestial sphere was divided into 28 lunar loges (refer to Mihn et al. 2017), and the length of the celestial perimeter called the Jucheon is 365 dA 2153 minAT in the Goryeosa. However, the value of the Tang dynasty, i.e., 365 dA 2153 minAT 299 secA, is also presented in cases of the Jucheon as well as some other constants. When converting the Tang value into units of secA, it becomes 920446199 that is recorded as the Sangsu in the Goryeosa. Therefore, it appears that the Goryeo astronomer used up to the unit of minA for the value of the Jucheon presumably to simplify the calculations. From the lengths of the tropical year and celestial perimeter, we could deduce that the mean motion of the Sun in the ecliptic was nearly 1 dA per day although this value was not explicitly defined in historical literature. In addition, the lengths of a synodic month called the Hapcheok and of an anomalistic month called the Yeokjuil was 29 dH 4457 minHT and 27 dH 4658 minHT 19 secH, respectively. Similarly, the 24 solar terms were utilized as the cardinal points in relation to the solar motion, and the periods of 7 dH 7465 minHT and 14 dH 6529 minHT derived from an anomalistic month were utilized as the reference points in relation to the lunar motion as mentioned below. The values of 7465 and 6529 were called the Chosu of the 7th and 14th day, respectively, and the values subtracted the Chosu from 8400, i.e., 935 and 1871, were called the Malsu of the 7th and 14th days, respectively. Finally, although the mean motion of the Moon among the stars was also not mentioned, it could be deduced as 13 dA 31 minAG per day, as estimated in Section 3.2. In Table 1, we summarize the values of astronomical constants in the SM calendar that are used in this study.

Table 1. Summary of astronomical constants in the SM calendar used in this study.
jass-36-2-87_T1
Download Excel Table

3. SOLAR AND LUNAR MOTIONS

3.1 Solar motion

The inequality of the solar motion has been known as early as A.D. 570 in China (Needham 1959). In the Goryeosa, the apparent movements of the Sun considering its unequal motion are tabulated with respect to the solar term, as in the case of the Chongxiu-Daming calendar (refer to Choi et al. 2018a). A year was divided into 24 solar terms, as mentioned above, and two methods were utilized for determining the instants of the solar terms: mean and corrected (literal translation) methods. The former method is to equally divide the length of the tropical year into 24 intervals while the latter method is to equally divide the ecliptic into 24 intervals conceptually. Currently, the corrected method is utilized for determining the dates of 24 solar terms (KASI 2018). In both methods, the starting point is the instant of a winter solstice of the preceding year. With regard to calculating the daily movement of the Sun from each solar term, only the calculation procedure is explained in the Goryeosa. In the Japanese literature, meanwhile, the amounts of daily Sun movement in each solar term are tabulated.

We present the ready reckoner relating to the solar motion in Table 2, which is tabulated in the Goryeosa. In column 1 of the table, i represents the solar term notated with an Arabic number. For reference, i= 1, 7, 13, and 19 represent the winter solstice, spring equinox, summer solstice, and autumnal equinox, respectively. In column 2, ∆USG is the unequal amount of the solar motion in units of minHG during the length of each solar term by the mean solar term method (i.e., Jungjeol) compared to its mean motion. In column 3, JG is the length of the ith solar term by the corrected solar term method given in units of dH, minHT, and secH. Because the Sun moves faster and slower around the winter and summer solstices, respectively, than its mean motion, the values of JG could be obtained by subtracting or adding the value of ΔUSG×100 to the Jungjeol (the units of ΔUSG is minHG not minHT). As an example, the value of JG in the second solar term, JG2, is 14 dH 5235 minHT 5 secH by subtracting +50×100 minHT to the Jungjeol (i.e., 15 dH 1835 minHT 5 secH). In this sense, the value of 2335 minHT for JG18 in the Goryeosa is the typographical error of 2435 minHT, which is the same that for JG7. In all tables, the values marked with the asterisk symbol are the corrected values in this study. In column 4, USG is the accumulation of ΔUSG until the preceding solar term; hence, the value of USG in the first term, US,1G, is zero. In column 5, ΔUSM is the value dividing ΔUSG×100 by the mean motion of the Moon per day, u¯MG. Although the value of u¯MG is not defined, it was 13 dA 31 minAG (or 13.3689 dA), according to our estimation in the next section. Adopting this value, we found that the values of ΔUSM were the results rounded half up. In addition, the value of ΔUSMin the 14th solar term, ΔUS,14M, should be +374 recorded as the value of the Tang dynasty in the Goryeosa. The last column, USM, is the accumulation of the ΔUSM until the preceding solar term; hence, the value of USM in the first term, US,1M, is zero. The value of US,17M is the typographical error in the Goryeosa. Although the Goryeosa states that the value of the Tang dynasty is 1246, it is recorded as 1346 in the version of the Xintangsh that we refer to in this study.

Table 2. Ready reckoner in relation to the solar motion.
jass-36-2-87_T2
Download Excel Table

As shown in Table 2, the solar motion is symmetrical around two solstices and equinoxes. The typographical errors discussed above can be verified from this perspective. In reality, however, the Earth revolves around the Sun and its motion is fast around the perihelion and slow around the aphelion. To estimate the accuracy of the solar motion in the SM calendar, we calculated the daily apparent movement of the Sun following the calculation procedures described in the Goryeosa. In Appendix A, we present the procedures using mathematical symbols and formulas to facilitate the understanding of contemporary concepts. As shown in the appendix, a quadratic equation was used to consider the inequality of the Sun, as in the Shoush calendar (refer to Choi et al. 2018b).

3.2 Lunar motion

The length of the anomalistic month, 27 dH 4658 minHT (19 secH), was divided into two ranges referring to the apogee (A) and perigee (P). Hereinafter, Sectors I and II stand for the A-P and P-A ranges, respectively. To calculate the position of the Moon for a given time, the SM calendar used the daily apparent movement of the Moon in each sector, which was given in a tabular form. In Table 3, we present the table listed in the Goryeosa in relation to the lunar motion. In the table, column 1, denoted as mA or mP (integer), represents the day number since the instants of the apogee and perigee passages, respectively (hereinafter, m is the common designation for mA and mP). Because the day number was counted from the first day (i.e., m=1), the last day number in each sector was the 14th day. In column 2, uMG is the apparent movement of the Moon in units of minAG on its orbital during the day (i.e., 8400 minHT) on the mth day. In column 3, ΔuMG is the difference of two successive values in uMG for m≤13, i.e., ΔuMGm=uMGm+1uMGm, and 0 for m = 14. In column 4, uMD is the same quantity with uMG but converted in units of dA and minAG, which is listed in a separate table in the Goryeosa. For instance, because the value of uMG on the third day in Sector II is 1203 minAG (= 14×Gakbeob + 27 minAG), hence that of uMD on the day is 14 dA 27 minAG. Suppose that the u¯MG is the mean motion of the Moon per day and δuMG is the difference between uMG and u¯MG. In columns 6 and 7, ΔuMM and UMM are the physical quantities dividing δuMG by u¯MG and accumulating ΔuMM until the preceding day, respectively; hence, UMM is zero for m=1. Strictly speaking, the value of ΔuMM for the mth day, ΔuMM (m), is

Δ u M M m = α × δ u M G m u ¯ M G m / 8400 min H T ,
(1)

Table 3. Ready reckoner in relation to the lunar motion.
jass-36-2-87_T3
Download Excel Table

where α is 1 except for m = 7 and 14. As shown in Table 3, two values are listed in ΔuMM for m = 7. These values can be obtained from Eq. (1) using the values of the Chosu and Malsuv on the 7th day. That is, two values in each sector can be obtained using α of 7465/8400 and 935/8400. Meanwhile, the value of ΔuMM for m = 14 is the result using the value of Chosu on the 14th day, i.e., α = 6529/8400. However, the value of the daily mean motion of the Moon is not stated in historical literature, as mentioned above. In addition, the daily apparent movement of the Moon in Sector I is not symmetrical to the movement in Sector II, as shown in column 1. Hence, a simple method of estimating the value of u¯MG is using the average of the minimum and maximum apparent movements of the Moon, i.e., 1012 and 1234 minAG, respectively. Subsequently, we can obtain the u¯MG of 1123 minAG (= 13 dA 31 minAG) or 13.1766 degrees. This value is similar to that defined in the Chongxiu-Daming calendar, i.e., 13.1763 degrees (Choi et al. 2018a) or as known in modern astronomy, i.e., 13.176358 degrees (Cox 2001). Applying the u¯MG value of 1123 minAG to Eq. (1), we obtained a good agreement with the values of ΔuMM except for several cases indicating differences of approximately 1 minHT. In our opinion, the discrepancies might be caused by the significant figures in the mean motion of the Moon estimated in this study.

3.3 Accuracy estimation

We estimated the accuracy of the solar and lunar motions by the SM calendar when the calendar was used in Korea by the Goryeo dynasty, and not in China by the Tang dynasty. Hence, we selected the periods of 919, 1155, and 1392 (i.e., the beginning, middle, and ending periods of the Goryeo dynasty, respectively) and compared the amounts of daily Sun and Moon movements listed in the Goryeosa with the results of modern calculations for each period. As mentioned earlier, the first step of the calendrical calculation in the SM calendar is to determine the instant of the CJDJ of a given year. However, in this study, we will not discuss the procedures of calculating the instant of the CJDJ in detail because they are well explained in the work of SADU (2016) where examples are provided. In principle, the instants of the CJDJ and those of the apogee passage are determined by accumulating the lengths of a tropical year and anomalistic month since the Super Epoch. However, it is noteworthy that the day of the winter solstice (also CJDJ) should be placed in the lunar November in Chinese calendars. In other words, a lunar November is the month containing the winter solstice day. In addition, it is noteworthy that the calendar days are obtained as a sexagenary circle instead of Arabic numbers. Hence, we referred to the work of Han (2001) to convert lunisolar days obtained as a sexagenery circle into Julian calendar days. Because the hour system of the SM calendar was based on the apparent solar time (AST), we also converted the AST at the Gaeseong (i.e., a capital of the Goryeo dynasty) into the mean solar time in Korean standard time (KST, i.e., at the meridian of 135° E) by correcting the longitudinal difference, equation of time (ET), and ∆T (i.e., time difference between the terrestrial time and universal time) (see also Lee et al. 2011). We assumed the longitude of the Gaeseong as 126° 33ʹ E, i.e., the longitudinal difference is -33.8 min (Lee 2017), and referred to the work of Morris and Stephenson (2004) for the value of ∆T (e.g., 34.5 min in 918). Finally, we used the astronomical algorithms of Meeus (1998) and the ephemeris of Standish et al. (1997) in modern calculations.

To estimate the accuracy of the solar motion, we first reproduced the instants of the CJDJ of the years used as the examples in the work of SADU (2016) and obtained the same results. Subsequently, we determined the instants of CJDJ of 919, 1155, and 1392 and summarized the results in Table 4 including the perihelion passage times in those days from modern calculations. In the table, YY, MM, DD, and HH represent the year, month, day, and hour, respectively. In the lunisolar date, the day is expressed with the notation utilized by Lee et al. (2012) and the unit of hour is minHT as that in the AST. For instance, #55 November 1154 in the lunisolar calendar corresponds to 15 December 1154 in the Julian calendar, and 1905 minHT in AST at the Gaeseong to 4.55 h in KST by correcting the longitudinal difference of -33.8 min, equation of time of -0.72 min, and ∆T of 14.7 min. According to our study, the difference between the CJDJ instant of 1155 and the perihelion passage time at that time was the smallest, i.e., -0.14 days, among the three periods. After determining the CJDJ instant, we computed the dates of the 24 solar terms by accumulating the Jungjeol and calculated the daily Sun movement by the SM calendar since the instant of each solar term throughout the procedures described in Appendix A. Subsequently, we compared the movements for 15 days in each solar term (i.e., for a total of 360 days) with the results of modern calculations. Hence, the interval between the last day in a solar term and the first day in the next solar term (i.e., next solar term day itself) is 1 dH 1835 minHT 5 secH, and not 8400 minHT (i.e., a day).

Table 4. Summary of the CJDJ instants of 919, 1155, and 1392 by the SM calendar and the perihelion passage times in those days by modern calculations.
jass-36-2-87_T4
Download Excel Table

Figure 1 shows the differences in the daily apparent movement of the Sun between the SM calendar, SSM, and modern calculations, SMC, for the periods of 919, 1155, and 1392. The horizontal axis is the elapsed date since the CJDJ instant of each year and the vertical axis is SiSMnSiMCn in units of (modern) arcmin, where i is the solar term from 1 to 24 and n is the day number since each solar term from 0 to 14. In the figure, the red solid, green dot, and blue short-dash lines represent the results for the periods of 919, 1155, and 1392, respectively; the vertical dotted lines represent the elapsed dates of the spring equinox, summer solstice, and autumn equinox, from left to right, since the CJDJ instant of each year. For a quantitative comparison, we also estimated the MAD value defined as follows:

MAD = 1 N i = 1 i = N X S M X M C i ,
(2)

jass-36-2-87_F1
Fig. 1. Differences in daily solar motion between the SM calendar, SSM, and modern calculations SMC: red solid, green dot, and blue short-dash lines represent the results for the periods of 919, 1155, and 1392, respectively. The horizontal axis is the elapsed dates since the CJDJ instant of each year and the vertical axis is SiSMnSiMCn in units of (modern) arcmin. The vertical dotted lines are the elapsed dates of the spring equinox, summer solstice, and autumn equinox, from left to right, since the CJDJ instant of each year.
Download Original Figure

where XSM and XMC are physical quantities from the SM calendar and modern calculations, respectively, and N is the number of data used in the calculations. We found that the MAD values are 0.33, 0.30, and 0.31 arcmin for the periods of 919, 1155, and 1392, respectively. In addition, the pattern in Figure 2 is similar to that in the result of the Chongxiu-Daming calendar of the Jin dynasty for the solar motion: the maximum differences are indicated around the winter and summer solstices, and the minimum around the spring and autumn equinoxes (refer to Choi et al. 2018a).

jass-36-2-87_F2
Fig. 2. Differences in daily lunar motion between the SM calendar, MSM, and modern calculations, MMC: red solid, green dot, and blue short-dash lines represent the results for the periods of 919, 1155, and 1392, respectively. The horizontal axis is the day number since the instant of the apogee passage and the vertical axis is MSM (m)-MMC (m) in units of (modern) degrees, which are averaged over the 12 cycles of anomalistic month length for each day number.
Download Original Figure

In the lunar motion section of the Goryeosa, the methods of determining the elapsed dates of four phases of the Moon (i.e., new, first quarter, full, and last quarter moons) since the instant of an apogee passage are described. In addition, the instant of the apogee passage to calculate the elapsed date of the new moon in lunar November was computed using the length of the Yeonyeo (i.e., the difference in instants between the CJDJ and apogee passage around that time). For this study, however, we obtained the instant of the apogee passage immediately before the CJDJ instant (hereinafter, the CJAP instant) by accumulating the length of the anomalistic month from the Super Epoch. In Table 5, we present the instants of the CJAP referring to the instants of the CJDJ of 919, 1155, and 1392 together with the results of modern calculations.

Table 5. Summary of the instants of CJAP of 919, 1155, and 1392 by the SM calendar and the differences with apogee passage times in those days by modern calculations.
jass-36-2-87_T5
Download Excel Table

To estimate the accuracy of the lunar motion in approximately one tropical year, we first determined 12 successive instants of the apogee passage in each year by accumulating the anomalistic month length to the CJAP instant. Subsequently, we calculated the daily Moon movement during each anomalistic month length using the DE405 ephemeris. As indicated in Section 3.2, however, the lunar movement listed in Table 3, uMG, is the daily amount (i.e., for 8400 minHT) for 28 days since the instant of an apogee passage although the length of the anomalistic month is 27 dH 4658 minHT 19 secH. Therefore, in modern calculations for the daily lunar movement, we determined the beginning point of each day by accumulating 8400 minHT from the instant of an apogee passage but reset the point on the days of mA= 8 and mP= 1 and 8 using the Chosu of the 7th and 14th days. For example, the beginning point on mA= 8 is 7 dH 7465 minHT and not 8 dH 0 minHT if the instant of each apogee passage is 1 dH 0 minHT. From the comparison with the movements listed in the Goryeosa, we found that the MAD values are 0.217, 0.284, and 0.240 (modern) degrees for the periods of 919, 1155, and 1392, respectively. In Figure 2, we present the differences in daily Moon movement between the SM calendar, MSM, and modern calculations, MMC. The horizontal axis is the day number since the instant of an apogee passage, and the vertical axis is MSMm¯MMCm¯ in units of (modern) degrees, which are averaged over the 12 cycles of anomalistic month length for each day number. In the figure, the legend is the same as that in Figure 1: red solid, green dot, and blue short-dash lines represent the results for the daily Moon movement for the periods of 919, 1155, and 1392, respectively.

4. SUMMARY

We investigated the solar and lunar motions in the SM calendar of the Tang dynasty in China. According to the preface of the Goryeosa, this calendar was used in Korea during the Goryeo dynasty succeeding the Silla Kingdom. Although it was known that the Goryeo dynasty introduced the Shoushi calendar in 1281, the SM calendar had been utilized until the end of the dynasty for the prediction of solar and lunar eclipse times. Primarily referring to the Goryeosa, we analyzed the calendrical methods of the SM calendar to calculate the solar and lunar motions. Subsequently, we compared the motions with modern calculations using the DE405 ephemeris for three periods: 919, 1155, and 1392, i.e., the beginning, middle, and ending of the Goryeosa dynasty, respectively. Regarding the solar motion, we found that a quadratic equation was employed to calculate the daily Sun movement. Furthermore, we found that the MAD values were approximately 0.31 arcmin for the three periods, although the minimum was indicated in 1155, where the CJDJ instant of 1155 was the closest to the perihelion passage time among the three periods. Meanwhile, we found that relatively large MAD values for the lunar motion: 0.247 degrees on average. In addition, we found that the MAD value indicated the maximum in 1155, unlike the cases of the solar motion. This might be the reason for the largest difference in apogee passage time between the SM calendar and modern calculations in 1155 among the three periods analyzed in this study.

In conclusion, we believe that this study will facilitate in the understanding of calendrical methods in the SM calendar, such as the sunrise and sunset times, eclipse times, and so forth.

APPENDIX A.

CALCULATION METHODS OF THE DAILY SUN MOVEMENT

According to the work of Qu (2008), the method of calculating the daily Sun movement was based on the Dayan calendar of the Tang dynasty. In the following, we formulated the calculation procedures explained in the Goryeosa. First, calculate the physical quantities of the Jungryul (JR), Hapcha (HC), and Jungcha (JC):

J R i = Δ U s , i G × α / G i
(A1)

H C i = J R i J R i + 1
(A2)

J C i = H C i × J G i / J G i + J G i + 1 ,
(A3)

where i represents the ith solar term; ΔUS,iG and JGi are physical quantities in the ith solar term, which are given in Table 2, and α is the Gakbeob (i.e., 84). In addition, the value of JRi+1 in Eq. (A.2) should be substituted as JRi-1 in the calculations of HCi for i = 6, 8, 12, and 18 (i.e., solar terms incurring the sign change in the value of ΔUSG compared to each subsequent solar term). Subsequently, compute the following quantities: Ilcha (IC), Choryul (CR), Malryul (MR), Chojeongryul (CJR), and Maljeongryul (MJR).

I C i = J C i × 2 × β / J G i
(A4)

C R i = J C i + J R i , M R i = J C i J R i
(A5)

C J R i = C R i I C i / 2 , M J R i = M R i + I C i / 2 ,
(A6)

where β is 100. In the SM calendar, the addition or subtraction operation differs with respect to the solar term in calculating the values of CRi, MRi, CJRi, and MJRi. However, these values can be obtained from Eqs. (A.5) and (A.6) if adopting the signing system given in Table 2. Finally, the daily inequality of the Sun (or equation of center of the Sun) in the ith solar term is obtained by cumulatively subtracting ICi from CJRi. That is, the unequal amount of the solar motion on the nth day in the ith solar term, ΔUS,iG (n), is

Δ U S , i G n = C J R i j = 0 j = n I C i for 0 n 14 Δ U S , i G n = M J R i × γ for n = 15 ,
(A7)

where n is the day number (integer) from 0 (i.e., day of each solar term) to 15, and γ is the decimal fraction when converting the length of the Jungjeol into units of d_H (i.e., 1835+5/88400 ).

According to our examination, however, the values of α and β should be a unity (see also Martzloff 2016). This can be verified by summing ΔUS,iG(n) for n from 0 to 15. That is, the value of n=0n=15ΔUS,iG(n) for a given solar term should equal to that of ΔUS,iG In addition, we reproduced the amounts of ΔUSM on the nth day in the ith solar term, ΔUS,iM(n), using the same method to obtain ΔUS,iG(n). Except for ΔUS,2M(n), in which their errors were already discussed by Soma (2014), we found an exact agreement with the values listed in the Seonmyoreki, a Japanese SM calendar preserved in the National Library of Korea, and the Sanpo. Finally, the unequal amount of solar motion until the nth day in the ith solar term, US,iG(n), is calculated as follows:

U S , i G n = U S , i G + j = 1 j = n Δ U S , i G n = U S , i G + j = 1 j = n ( C J R i k = 0 k = j I C i ) .
(A8)

It can be rewritten as

U S , i G n = U S , i G + n × C J R i n n 1 2 × I C i .
(A9)

This equation indicates that a quadratic equation was utilized to consider the inequality of the Sun.

ACKNOWLEDGEMENTS

This work was supported by research grants from Daegu Catholic University in 2017.

References

1.

Chang DI, The First Day in the Lunar Calendar in Koryeosa, History Education Review 52, 183-232 (2014).

2.

Choi GE, Lee KW, Mihn BH, Models of solar and lunar motions in the Chinese Chongxiu-Daming calendar, Advances in Space Research 62, 3237-3248 (2018a).

3.

Choi GE, Lee KW, Mihn BH, Li L, Ryu YH, Ahn YS, Investigating Chinese mathematical techniques to calculate sunrise and sunset times in Datongli, Astronomical Notes 339, 520-532 (2018b).

4.

CoxN, Allen’s Astrophysical Quantities. 4th Eds (Springer, New York, 2001).

5.

Editorial Department of Sinseowon (EDS), North Korean Translation of the Goryeosa: Volume 5 (Sinseowon, Seoul, 1991).

6.

HahnYH, LeeEH, Kang MJ, Chijeongsan-Naepyeon 1 & 2 (Institute for the Translation of Korean Classics, Seoul, 2016).

7.

HanBS, Arrangement of Chronological Tables on Korea: Volume II (Yeungnam University Press, Gyeongsan, 2001).

8.

Institute for Oriental Studies of Daito Bunka University (IOSDBU), Study on the Senmyoreki of the Calendar Chapter in the Goryeosa (Institute for Oriental Studies of Daito Bunka University, Tokyo, 1998).

9.

Jun YH, Calendrical Calculation System and Annual calendar of the Goryeo Dynasty, The Journal of Korean Medieval History 39, 193-257 (2014).

10.

Kim IG, A Study of Koryo Calendar System by Translation of the Koryo-sa Yokji, Journal of Seokdang Academy 4, 99-143 (2012).

11.

Korea Astronomy and Space Science (KASI), Korean Astronomical Almanac 2019 (SMBook, Seoul, 2018).

12.

LeeES, The Principle Analysis of the Calendar (Jeongeumsa, Seoul, 1985).

13.

Lee KW, Daylight Observations of Venus with Naked Eye in the Goryeosa, Journal of Astronomy and Space Sciences 34, 67-73 (2017)

14.

Lee KW, Ahn YS, Mihn BH, Verification of Calendar Days of the Joseon Dynasty, Journal of the Korean Astronomical Society 45, 85-19 (2012).

15.

Lee KW, Ahn, YS, Yang, HJ, Study on the system of night hours for decoding Korean astronomical records of 1625-1787, Advances in Space Research 48, 592-600 (2011).

16.

Martzloff JS, Astronomy and Calendars – The Other Chinese Mathematics: 104 BC – AD 1644 (Springer, Berlin, 2016).

17.

MeeusJ, Astronomical Algorithms (Willmann-Bell Inc., Virginia, 1998).

18.

Mihn BH, Choi G, Lee YS, Astronomical Instruments with Two Scales Drawn on Their Common Circumference of Rings in the Joseon Dynasty, Journal of Astronomy and Space Sciences 34, 45-54 (2017)

19.

Morrison LV, Stephenson FR, Historical values of the Earth’s clock error and calculation of eclipses, Journal for the History of Astronomy 35, 327-336 (2004).

20.

NeedhamJ, Science and Civilisation in China: Volume 3 (Cambridge University Press, Cambridge, 1959).

21.

SatoM, History of Japanese Calendar Study (Surugadai Press, Tokyo, 1968).

22.

Seo KS, A Review on the preface of Yeokju, Goryeosa through the development of Goryeo Calendrial Science, Chonam Historical Review 47, 29-70 (2012).

23.

Seokdang Academy of Dong-A University (SADU), Korean Translation of the Goryeosa: 14 (Kyungin Publishing, Seoul, 2016).

24.

Sivin N, Granting the Seasons: The Chinese Astronomical Reform of 1280, With a Study of Its Many Dimensions and an Annotated Translation of Its Records (Springer, New York, 2008).

25.

SomaM, Mistakes in the Xuanming Calendar, in The Eighth International Conference on Oriental Astronomy, Hefei, China, 26-28 March, 2014.

26.

StandishEM, NewhallXX, Williams, JG, Folkner, WF, JPL Planetary and Lunar Ephemeris (CD-ROM) (Willmann-Bell Inc., Virginia, 1997).

27.

QianB, Zhongguo shuxue shi (Kexue chubanshe, Beijing, 1981).

28.

QuA, Zhongguo shuli tianwenxue (Kexue chubanshe, Beijing, 2008).

29.

UchidaM, Nihon Rekijitsu Gentn (Yuzankaku Inc., Tokyo, 1975).