0%

2019年水星凌日计算

历表:VSOP87
编算:高良超

根据水星凌日规律,2019年11月的水星下合日时将发生水星凌日现象,这一章我们将详细介绍水星凌日的计算。

(一)、求水星与太阳赤经相合时的力学时

从星历表中查出与α日-α水=0的相近的数据前后共四个如下:

11月11日

力学时     α水 α日 α日-α水 △’       △’’

14时 226º33’13.80’’  226º25’03.90’’  -0º08’09’’.90

15时 226º30’06.90’’  226º27’35.85’’  -0º02’31’’.05  +5’38’’.85

16时 226º27’00.15’’  226º30’07.80’’  +0º03’07’’.65  +5’38’’.70  -0’’.15

17时 226º23’53.55’’  226º32’39.90’’  +0º08’46.35’’  +5’38’’.70  +0.00

-f0=+2’31’’.05=+151’’.05

△’½=5’38’’.70=+338’’.7

△’’0+△’’1=-0’’.15

B2=n(n-1)/4

B2

-B2(△’’0+△’’1)

F

n

第一近似

0

0

151.05

0.445970

第二近似

-0.061770

-0.01

151.04

0.445940

第三近似

-0.061769

-0.01

151.04

0.445940

0.445940是时的小数,化为分秒后得26分45.384秒,因而所求水星与太阳赤经相合的时刻为:2019年11月11日15时26分45.384秒(力学时)。

现在以所算出太阳与水星赤经相合的时刻:11月11日15h26m45.4s为引数,从历表查出该时刻水星的黄纬为: 0°01’19.47”,根据凌日判别式肯定有中心凌发生。

如果要计算这次凌日,必先求力学时11月11日12时~19时的太阳、水星要素。

现将各小时的太阳、水星要素分别列成表 1.1和表 1.2

表1.1太阳要素表

力学时

α

δ

s

π

12时

226°20’00.00’’

- 17°24’36.70”

969.26

8.88

13时

226°22’31.95’’

- 17°25’18.08”

969.27

8.88

14时

226°25’03.90’’

- 17°25’59.42”

969.28

8.88

15时

226°27’35.85’’

- 17°26’40.73”

969.29

8.88

16时

226°30’07.80’’

-17°27’22.02’’

969.30

8.88

17时

226°32’39.90

- 17°28’03.27”

969.31

8.88

18时

226°35’11.85’’

- 17°28’44.49”

969.32

8.88

19时

226°37’43.95’’

-17°29’25.68’’

969.33

8.88

表1.2水星要素表

力学时

α

δ

s’

π’

12时

226°39’27.75’’

- 17°31’35.72”

4.97

13.02

13时

226°36’20.70’’

- 17°29’51.32”

4.97

13.02

14时

226°33’13.80’’

- 17°28’06.89”

4.97

13.01

15时

226°30’06.90’’

- 17°26’22.44”

4.97

13.01

16时

226°27’00.15’’

-17°24’37.97’’

4.97

13.01

17时

226°23’53.55’’

- 17°22’53.50”

4.97

13.01

18时

226°20’47.18’’

- 17°21’09.02”

4.97

13.00

19时

226°17’40.80’’

-17°19’24.54’’

4.97

13.00

根据要素使用x=(1-ytanδ)cosδtan(α’-α)

y=δ’-δ+xsinδtanα’-α/2公式,编制较差坐标x与y的表

力学时

x

x’

y

y’

(xx’+yy’)/100

√x’²+y’²

12时

+1113.5’’

-323.1’’

-420.0’’

+146.5’’

-4213.0’’

354.8

13时

+790.4

-323.1

-273.7

+146.3

-2954.2

354.7

14时

+467.3

-323.1

-127.6

+146.1

-1696.3

354.6

15时

+144.1

-323.2

+18.3

+145.9

-439.0

354.6

16时

-179.1

-323.2

+164.0

+145.7

+817.8

354.5

17时

-502.3

-323.3

+309.6

+145.6

+2074.7

354.5

18时

-825.3

-323.3

+455.0

+145.4

+3329.8

354.5

19时

-1148.6

-323.3

+600.1

+145.1

+4584.2

354.4

x、x’,y,y’ 均表为角秒。由线性内插法求得xx’+yy’在

15h.3493=15h20m57s.5

之时为零,即在地心看见的两心距离最短的时刻。于是

x=+31.2

y=69.2   dm=75’’.9=1’15’’.9

再计算在地心所看见的这个现象。因太阳的视半径969.3’’,水星的中心所走的弦的半长度便等于

√969.3²-75.9²=966.3’’

因水星的速度为每时354.5’’,故水星凌日所经历的时间之半为2h.726;因此凌始约在11日12h.62,凌终约在11日18h.08。为了寻求更精确的数字,计算12h.60与12h.65的X与y:

x         y      d

12.60  +919.64  -332.22  977.81

12.65  +903.49  -324.91  960.14

在这个两个时刻,S+S’=974.24’’,S-S’=964.30’’。在表中作d的线性内插,可以求得相切的时刻,并顺便算出方位角Q

x        y      Q

第一次切:12h.6101=12h36m36.4s  +916.4  -330.7  109°51’

第二次切:12h.6382=12h38m17.6s  +907.3  -326.6  109°48’

对于凌终:

x          y       d

18h.05  -841.47  +462.27  960.09

18h.10  -857.63  +469.54  977.75

又S+S’=974.29,S-S’=964.35,由内插法求得

x         y      Q

第三次切:18h.0621=18h03m43.4s  -845.4  +464.0  298°46’

第四次切:18h.0902=18h05m24.7s  -854.5  +468.1  298°43’

1541068962954913.png

2019年11月11日 水星凌日过程图

将地心所见的相切时刻化为地面所见的时刻应加的改正数,介入地心纬度:

ρcocφ’=cosu   ρsinφ’=(1-e)sinu

tanu=(1-e)tanφ (e=1/298.25)

并设A=(y(1-e)(π’-π)cosδ’)/xx’+yy’

Bcosω1=-(y(π’-π)sinδ’)/xx’+yy’

Bsinω1=(x(π’-π))/xx’+yy’

B的符号肯定后便可决定ω1的角,再设

ω=H1’-ω1

H1’是相切时行星在格林尼治的时角。于是地面看到的相切时刻如下表达:

△t=Aρsinφ’+Bρcosφ’cos(ω-λ)

剩下的是计算视差的改正值。在进入日面时可取π’-π=4.14’’,出去时可取4.12’’,由于改正值是对每时秒而定的,故以x’与y’的数值应以3600除之。

切次    A   Bcosω1   Bsinω1     B     ω1     ω

1  +13s.60  +4s.30  -39s.65  -39.88   96.19  -83.62

2  +13s.66  +4s.29  -39s.64  -39.87   96.18  -83.18

3  +19s.21  -6s.02  -36s.79  +37.28  279.29  +184.43

4  +19s.18  -6s.01  -36s.81  +37.30  279.27  +183.99

因此利用下列的公式,可以计算地面看见的相切时刻:

t1  12:36:36.4+13s.60ρsinφ’-39s.88ρcosφ’cos(83°.91-λ*)

t2  12:38:17.6+13s.66ρsinφ’-39s.87ρcosφ’cos(83°.47-λ*)

t3  18:03:43.4+19s.21ρsinφ’+37s.28ρcosφ’cos(184°.14-λ*)

t4  18:05:24.7+19s.18ρsinφ’+37s.30ρcosφ’cos(183°.70-λ*)

接下来求水星凌日时水星在天顶的地点

力学时        恒星时         α水 λ*=恒星时-α水 φ=δ水

t1  12h36m36.4s  239º31’06’’.9  226º37’33’’.6  +12º53’33’’.3  - 17°30’32.03

t2  12h38m17.6s  239º56’29’’.1  226º37’28’’.4  +13º19’00’’.7  - 17°30’29.09”

tm  15h20m57.5s  280°43’09’’.8  226°29’01’’.6  +54°14’08’’.2 - 17°25’45.95”

t3  18h03m43.4s  321°31’18’’.0  226°20’35’’.5  +95°10’42’’.5  - 17°21’02.53”

t4  18h05m24.7s  321°56’41’’.7  226°20’30’’.3  +95°36’11’’.4  - 17°20’59.59”

归纳结果如下:

2019年11月11日 水星凌日(地心视象)

力学时        北京时    方位角         月亮在天顶的地点

11/11     11/11~12                   地理纬度  地理经度

凌始外切 12:36:36.4  20:35:26.9  109°51’     -17°32.5’   -12°36.1’

凌始内切 12:38:17.6  20:37:08.1  109°48’      -17°30.5’  -13°01.6’

凌    甚 15:20:57.5  23:19:48.0              -17°25.8’   -53°56.7’

凌终内切 18:03:43.4  02:02:33.9  298°46’     -17°21.0’   -94°53.3’

凌终外切 18:05:24.7  02:04:15.2  298°43’     -17°21.0’   -95°18.8’

最小角距离=1’15’’.9

至此可以看出这次凌日,在太平洋(除西北部)、北美洲(除西部和北部)、南美洲、大西洋、欧洲(除极北部)、非洲、亚洲极西南部、印度洋西部、南极洲可以看到。

说明:2019年的水星凌日根据VSOP87行星理论结合凌日理论计算而来,给出具体的算例,误差约1秒钟,可供天文爱好者参阅。水星凌日有规律地发生在5月初或11月初,平均每百年出现13~14次水星凌日的现象。