日本付近の地磁気経緯度を求める早見表
久 保 木 忠 夫
概 要
最近地球内部の構造を調べるため，臨時の地磁気観測所が多くつくられた。このため観測所の地磁気
経緯度を計算する機会が多し、。この計算式は簡単なものであるが，無味乾燥で，計算に誤りをおかしや
すく，計算に要する時間も短くな
L。
、
このためきわめて簡単に求められる使利な早見表を作成した。この表の使用法はきわめて容易で，地
磁気経度 d・緯度 φ・計算上の偏角
ψ
を求めるの
'
l
こ要する時間はわずか数分であり，求めた値の誤差
は約:f:
O
.
6
'
できわめて正確である。
早見表は
1
5
組の表で構成されてし、る。!原点の北地磁気極の経緯度を
2
9
1
.
O
O
E
，
7
8
.
5
0N
として，地理的
経緯度 1
。毎の
A.
φ
・
ψ
の主値に.'の補正値を 2個づっ加算するだけで暗算でも求められる。もし原
点をずれた値に Lたいときは，それにさらに 2倒の補正値を加えればよ
L、。このほかに度分換算表，比
例撹入表がある。
この早見表は三角関数を用いる計算式とは異って，単なる加算だけで求めるものであるから誤り計算
をすることが全くなし、。
また現在使用されている日本の各観測所の地磁気経緯度の値を示し，その
2
，
3
の問題点を述べた。
また誤差に関する計算についても少しふれた。
1 . は し が き
日本付近には割合に多くの地磁気観測所がある
Oと く に 地 球 内 部 開 発
CUpperm
a
n
t
l
e
p
r
o
j
e
c
t
)
の研究の一部として，
1
9
6
4
年から
1
9
6
5
年の
IQSY
の観測期間と期を同じくして，
日本では数多くの臨時地磁気観測所が建設され，観測が行なわれた
o
このとき著者は資料の
整理の一部を分担して，各観測所の地磁気経緯度の計算を行なった。
また柿聞などの地磁気経緯度の他は，ごくわず.かで、あるが，それぞれ公表されている値が
異なっている。この原因を山口又新が調査した結果，それは北地磁気極の経緯度の違いと計
算の誤りであることが分ったく
1)0地磁気経緯度の計算式はよく知られているぬように，球面三角法で容易に求まる簡単な式
である
Oしかも一般には
0
.
1
0まで求めれば十分であり，計算は容易なはずである。しかし
実際に求めてみると，角度の存在範聞のきめ方に手数がかかり誤りをおかし易く，意外に多
くのi
時間を要する。慣れないと 2時間もかかることがあり，労力も馬鹿にならなし、。
快界の任意の地点、の地磁気経緯度を求めるため ，
A
.
G
.
M
c
N
i
s
h
は 北 地 磁 気 極 を
φ
。
=78.5
0N
，Ao
=29
1
.
0
0E
として求めるノモグラムを作成した向。これはある地点の地理的
経緯度 φ，λ を与えたときに地磁気経緯度 0
，
A
との差 φ-0，λ
-A
および地理的子午線
と地磁気子午線のなす角
ψ
を，それぞれ
φ
，λ の関数として与えたものである
oそ れ ぞ れ
の{直は:t
0
.1
0_
_
0
.
2
0まで求められる。
現在まで多くの観測所は，このノモグラムを利用して地磁気経綜度 0
，
A
および角
ψ
を
求めて，それを使用してきた。
地磁気経緯度の使用目的から言えば:t
0
.
1
0の誤差があっても十分であるが，報告した
f
i
t
t
がたとえ
0
.
1
0の差であっても，ずれていると都合の悪いことが多く，印刷誤りなどととられ
る向きもある。また北地磁気極も最近の資料によるとゆ。
=
7
8
.
5
0N
，
λ
。
=29
1
.
0
0E
とは少し
ず、れた位置になっている。
このため著者は日本付近の観測所の地磁気経純度 0
，
A
および V を求める早見表を作っ
た。これは φ。 =78.5
0N ， λ。 =291.0
0E として日本付近の経緯度 1
0毎の地磁気経綜度~0 ，
A
および角
ψ
を求め，この主値に地理的経緯度の'の値に相当する分を補正として加算して
求める
oしたがって
2
個の補正をすればよいし，北地磁気極が変ったときは，それに相当す
る補正をさらに
2
伺加えればよ
L
、。この早見表は使用法が簡単で，便利で， しかも計算は正
確で迅速にできる。
数百個に及ぶ地点の値を求めるならば，電子計算機を利用する方が賢明であるが，わず‘か
2
，
3
個の計算をする場合はその準備作業にのみ労力をとられて，非常に不利である。しか
も観測所はそれほどつぎつぎと建設されるのでないから，このような早見表の力がはるかに
すぐれている。そして将来のことも考えて，早見表は:t
1
'
の精度が得られるようにした。
また最近の観測所の地磁気経緯度についても 2， 3ふれてみたいと思う。とくに若干の差
のある値について出典を明らかにしておき，混乱を避けるようにした。
近頃では補正磁気経緯度が利用され，羽倉によって多くの観測所の値が求められている∞。
これは実際の地球磁場の分布から求めた経緯度へ極付近の現象を取り扱うのに好都合であ
る。しかしこのために従来の地磁気経純度が不必要となる訳ではないし，とくに地磁気観測
所の表示要素には必ず必要なので，著者の作った早見表も今後利用されることが多いと考え
る
。
2
.
地磁気極の経緯度
地磁気の極の位置は多くの観測結果を用いて，伏角が
-90
0または
9
0
0，あるいは水平分
力が零になる南磁極と北磁極が求められて
L
、る。一方地球上の多くの資料を用いてこれを球
関数に展開して，地球の中心においた双極子による磁場におきかえて，それの示す雨，北地
磁気極の位置が求められる
o
最近のように，観測値とくに極地方の備が多くなると，地磁気の極の経総度もどんどん新
しいものに更新されていくことになる。
ここでは一般に知られている主な値について第 1表に示した山， (的。この表で
Chapman
，
B
a
r
t
e
l
s
の
1
9
4
0
年の値は
B
a
u
e
r
の値の四捨五入したものであるが，この値が一般に地磁気
経緯度を求めるときの北地磁気極の経緯度として用いられている。
No.7
，
8
，
9
はほとんど同
じ資料を用いて計算されたもので，総度は同じであるが経度に若干の差がある。このほかに
もいくつか公表された値がある
o日本付近の地磁気経緯度を求める早見表
65
Table 1 The p
o
s
i
t
i
o
n
o
f
t
h
e
geomagnetic p
o
l
e
North p
o
l
e
50uth p
o
l
e
No
Author
!Epoch;Latitude
JLongitude
L
a
t
i
t
u
d
e
l
Longitude
i
Note
I Geo
・ ， _
!Geo- I
n'_I
Geo
・
tn，
_
I
Geo-I
'Dip
1 .::~:::-I
Dip i :
:
:
:
:
-
I
Dip
1 -:..::::I
Dip
1 -:..:-:-i
l
i9
0
0
l
mq-
190
。
i
ma
・
写
190
。
lmq-! 90
。
l
f￨
-
- Inetic I
-
- I
n
e
t
i
c
I
-
- I
n
e
t
i
c
I
-
- I
n
e
t
i
c
I
1 5hmidt
1
8
8
5
-
7
80
32'N -
2
9
10
3
0
'
E
ー ￨ ー ￨ ーAdmiralty c
(
E
n
g
l
i
s
h
)
h
a
r
t
1922710N
.
.
.
，
.
.
.
.
.
.
.
.
.
.
.
ー2
6
40
E
ー7
30
5I - I
1
5
60
E
I
-..~....
.
.
.
.
.
.
.
，
I 1
.
.
.
，
.
.
"
.
.
.
.
1
3 Bauer
1
9
2
2
-
7
80
-32' -
2
9
0
0
-
5
2
'
ーI
4 Chapman
，
Bart
出(
1
倒0
)
ー7
8
.
50
-
2
9
1
。
5 Meyer
194
ら。--
2 6 0 J -
ーにご!ー
iI
三
M m r
1952 十竺ざ二一竺~:~~.7CJ
-_
I
叩
。
L
二
1-I
l
e
x
p
a
n
d
I
7 Finch-Leaton
附 ー7
8
.
3
10 -
291ω
。
- 1
ーi
ーI- I
u
n
t
i
l
I I I
Pg
I
8 V
す
ま
i
-
忠
f
A
)
師 一 泊4
5
0
2
9
0
・8
9
0
-
1 - I - I - 1P
!
i
i
7
噌
1
5
2
2
L
)
向 て
M
。 一 遍 五
-l-l-l-lP11!
今後新しい観測値が増えたとしても，一般に用いられている
NO.4 の値から大幅にずれ
た値になるとは考えられないから，実用上はこの値をそのまま使用しでもそれ程の問題は起
こらなし、。しかし地磁気経緯度を使用するときは，それが求められたときの北地磁気極の経
緯度がどの値であったかは知っておく方がよし、。そうすれば新しい北地磁気極の経緯度の値
に換算するとき都合がよいし，値のずれによって混乱の起こることがな L、
。
3
.
地磁気経緯度と子午線聞の交角(計算上の偏角〉の求め方
3
.
1
第
1
の方法
第
1図の説明図の球面三角形で考えてみる由。
N:
地理的北極
B:
。北地磁気極(北磁極ではない)
P:
観測所
B
点の経度 λ0，
緯度
φ
。
P
点の経度 λ，
緯度
φ
P
点の地磁気経度
A
，
地磁気緯度
e
計算から求まる磁気子午線と地理的子午線と
のなす角(ここでは
H計算上の偏角"という〉
を
ψ
とする。
ここで経度は東経，緯度は北緯，計算上の偏
角
ψ
は東偏をそれぞれ正とする
c図では
ψ
は
1
1
.
_0
0 ー 今E
F
i
g
.
1 I
l
l
u
s
t
r
a
t
i
n
g
t
h
e
c
a
l
c
u
l
a
t
i
o
n
o
f
t
h
e
geomagnetic c
o
o
r
d
i
n
a
t
e
s
(
P
a
r
t
1
)
西偏であるからー
ψ
と表示した。
B
，
P
点の経緯度が分れば，球間三角法で容易に
A
，t
J
，
[
ψ
は求まる
O対数計算が可能な
ように
，
P
点より
NB
に震直線
PD
を引き補助角
ND=K
を用いる
O第 1
図で
NB=8
1，
NP=8
，
BP=
θ
，
NBP=180
0
ー伊!
，
!
NPB=
ψ
，
BNP=
伊1，
DBP=
伊2，
ND=K
とすれば
81
=900
ーゆ
0
，
8=900
ー
φ
，
θ=900-0
伊
1
=
λ
一
λ
0
，伊
2=
/
1
となる
oChapman
の求めた式
(
2
)
と若干記号が異なるが本質的には変らな
L
、。求める式は
次のように表わされる。
または
t
a
n
K
=
t
a
n
8
c
o
s
伊1
…
・
…
.
.
.
.
.
・
H・
.
.
…
・
・
・
・
・
・
・
・
・
一
…
-
・
・
・
・
・
・
・
・
…
…
…
.
.
.
.
.
.
…
・
・
・
・
(
1
)
刷 物=tan
例豆 旦 豆 一 …
H・
H・
-
…
H・
H・-………
H・
H・
-
…
H・
H・
-
・
・
…
・
・
(2)
i
n
(K-8
I)c
o
t
θ=cos
伊2c
o
t
(K-8
1)………・・・……・…………・・………・・……・
(3)
sinψ=-sin 81
ま
号
(4)
LIsm
伊
。
=
-sin 81
盃
五
6
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・ (5)
これは別に面倒な計算式でないが，角度のとり方に誤りをおかし易く，意外に計算に手間
どるものである。
角度の存在範聞は
McNish
に詳しく述べられている仰が，簡単に日本付近の例で述べる。
(1)
式で伊
1は負号であるから
tanK<O
ゆえに
900~K<1800または
2
7
00壬
K<3600
しかし
K
の定義より
OO~K孟 1800ゆえに
K
の存在範囲は次のように限定される。
9
0
0 ;;;玉K<1800
………・…・・……・・・………・…・・……・・・…・・……・・・・・…・・・・・
(6)
この
K
，
λ，
φ
0
，
A
.
oを考慮すると公式(2)から
t
a
n
伊
2>0
ゆえに
0<
伊
2
:
:
三
9
00
または
1
8
00<
伊
2
豆
2
7
00
(3)
式から
0< 伊2~900ならば
9
00<
θ
孟
1
8
00
または
2
7
00<
θ
孟
3
6
00
180
0< 伊2~2700
ならば
00<
θ
孟
9
00
または
1
8
00<
θ
孟
2
7
00
しかし
θ
の値は
00
孟
θ
孟
1
8
00
したがって θ と伊
2の聞には次の制限がつく。
1800< 伊2~玉 2700と
00<
θ
孟
9
00
・
・
・
・
・
・
…
…
・
…
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
…
・
・
・
・
・
・
・
・
・
(7)
または
0< 伊2~五 900と
9
00<
θ
孟
1
8
00
・………・・…
(8)
一方(4)と(5)式から
s
i
n
(
)
s
i
n
伊I=sinθsin
伊2
・
H・
H・
.
.
…
…
.
.
.
.
・
H・....………・・・・…・・
(9)
(9)
式で
伊
1
くO したがって
s
i
n(
)
s
i
n
伊
1<0
ゆえに
s
i
n
θsin
伊
2<0
ゆえに
00
三
五
θ孟1800
と
1800三伊
2~三3600... …...………・…・・………(10)日本付近の地磁気経緯度を求める早見表
67
または
1800~θ 三3600と
伊豆伊2
三
~1800... ・ H ・..…・・・・・・…...・...……… (11)したがって伊
2，
θ
の聞で共通に存在するのは， (7)，
(8)，
(
1
0
)，
(
1
1
)式・のうちで(7)と
(
1
0
)
式だけであり，その関係から制限の条件として次式をうる
O0
0
<θ
三
900
で
180
0
<
伊
2
;
'
;
三
2700
………・・・・・・……・・…・…・・・・・…・・・…・
(12)
したがって地磁気経度
A
は次式になる。
1
8
0
0<A
孟
2700
・
・
H・
H・
.
.
.
.
.
・
H・
.
.
.
.
.
・
H・
.
.
.
・
H・-………・・・・・………...…………
(13)
以上の (6)と
(12)
式が計算のとき必要になる。
次に計算の一例を柿同について示寸。計算は最終的な値が
0.1
0単位でよいから ，l
'単位
の三角関数表を使用した。この計算例はだ足の感じがするが，後に述べる著者の早見表が便
Table 2 Example o
f
c
a
l
c
u
l
a
t
i
o
n
o
f
t
h
e
magnetic c
o
o
r
d
i
n
a
t
e
a
t
Kakioka (
P
a
r
t
1
)
1 尚
L_F
…
r
J
.
~_~
Value. _
_
_
_
，
_
一 円
(
1
)
Ao2
9
10
IE
a
s
t
l
o
n
g
i
t
u
d
e
(2)
仇
7
8
.0
5
I
North l
a
t
i
t
u
d
e
l ( 3 ) 1 A
l 附1
1
r
i
E
制l
吋
t
山(4)
ゆ
360 1
4
'
North l
a
t
i
t
u
d
e
(5)
。
5
30 4
6
'
9
0
。
ー
(4)
(6 )
l
o
g
t
a
n
(
)
0
.
1
3
5
0
2
5
l
o
g
t
a
n
(5)
(7)
Cll-1500 4
9
'
(3)
ー(
1
)
(8)
l
o
g
c
o
s
C
P
l
9
.
9
4
1
倒6
(ー)
l
o
g
c
o
s
(
7
)
(9)
l
o
g
t
a
n
K
0
.
0
7
6
0
7
1
(ー)
(6)+(8)
(
1
0
)
K
1
3
00 0
0
'
(9 )
→a
n
t
i
l
o
g
.
.
e
q
u
a
t
i
o
n
(6)
(
1
1
)
。
l1
1
0 3
0
'
90
。
一
(2)
(
1
2
)
K
.
-
(
}
.
1
1
80 3
0
'
(
1
0
)
ー
(
1
1
)
(
1
3
)
l
o
g
t
a
n
判9
.
7
4
7
0
2
3
l
o
g
t
a
n
(
7)
(
1
4
)
l
o
g
s
i
n
K
.
9
.
8
8
4
2
5
4
l
o
g
s
i
n
(
1
0
)
(
1
5
)
l
o
g
sin(K-(}.)
9
.
9
4
3
8
9
9
l
o
g
s
i
n
(
1
2
)
(
1
6
)
l
o
g
t
a
n
C
P
2
9
.
6
8
7
3
7
8
(
1
3
)
+ (
1
4
)
ー
(
1
5
)
(
1
7
)
C
P
2
2
0
50 5
7
'
(
1
6
)→a
n
t
i
l
o
g
.
.
e
q
u
a
t
i
o
n
(
1
2
)
.
(
1
3
)
(
1
8
)
A
2
0
50 5
7
'
East
(
1
9
)
l
o
g
c
o
s
C
P
2
9
.
9
5
3
8
4
5
(ー)
l
o
g
c
o
s
(
1
8
)
(
2
0
)
l
o
g
c
o
t
(K
-
(
)
，
)
9
.
7
3
4
7
6
4
(ー)
l
o
g
c
o
t
(
1
2
)
(
2
1
)
l
o
g
c
o
t
θ
9
.
6
8
8
6
ω
(
1
9
)
+ (
2
0
)
(
2
2
)
e
6
30 5
9
'
(
2
1
)→a
n
t
i
l
o
g
.
.
e
q
u
a
t
i
o
n
(
1
2
)
(
2
3
)
φ
2
60 0
1
'
90
。
一
(
2
2
) North
(
2
4
)
l
o
g
s
i
n
，
)
(
9
.
2
9
9
6
5
5
l
o
g
s
i
n
(
l
l
)
(
2
5
)
l
o
g
s
i
n
C
P
2
9
.
6
4
1
侃4
(ー)
l
o
g
s
i
n
(
1
7
)
(
2
6
)
l
o
g
s
i
n
(
)
9
.
9
凶6
6
7
l
o
g
s
i
n
(
5)
(
2
7
)
l
o
g
s
i
n
1
J
!
9
.
0
3
4
0
5
2
(ー)
(
2
4
)
+ (
2
5
)
ー
(
2
6
)
(
2
8
)
-
w
_60 1
3
'
(
2
7
)→a
n
t
i
l
o
g
.
，
e
q
u
a
t
i
o
n
(
4
)
(
2
9
)
v
60 1
3
'
(
E
a
s
t
)
利なものであることを実証するため，とくに比較し易いようにのせた
oこの計算で
K
の備と符号，仰の値と符号が最も誤り易し、。大体の
A
，θ
，
ψ
はわかって
いるからそれを参考にすればよいが
，
K
を誤ると予想債にならなくなり混乱する。
この計算はなれれば験算まで入れて
3
0
分
.
-
4
0
分でできるが，まごつくと
1
時間半もかかり，
手元に適当な参考文献とくに
McNish
の論文加でもないと
2
時間くらいかかることがある。
従来の値が古い北地磁気極の経緋度であるので，新しいものに改めるために再計算する場
合などは，つまらぬ労力を用いることになる。しかも値がたかだか
0
.
1
0単位でよいから余
計につまらぬ労働であると感じる。
3
.
2
第
2
の方法
本質的には全く変らないが，角度の計算
の誤りを少なくするため，第
2
図のように
考える。この方法は日本付近ではとくに好
都合である。
第 1図と同じく
N
点;地理的北極
B
点;北地磁気極，経度 λ0，緯度
φ
o
P
点;観測所
経度
λ
，純度
φ
-P
点の地磁気
経度
A
，
純度
e
とする
oG
e
o
g
J
;
t
1
P
h
i
c
，
a
l
n
o
n
h
凶
e
ー - TE
ここで経度は東経，緯度は北緯，
H
算上
の偏角
ψ
は東偏をそれぞれ正とする。図
は第 1図と同じく西偏に書いてあるから
F
i
g
.
2
g
e
o
m
I
a
I
g
l
u
n
s
e
t
t
r
i
a
c
t
c
i
n
o
g
o
t
r
d
h
i
e
n
c
a
t
e
a
s
l
(
c
u
l
P
a
a
t
r
i
t
o
2
n
o
)
f
t
h
e
一
ψ
と表示しておく。
第
2
図で
とすれば
NB=8
1
，
NP=
θ
，
BP=()
NBP=
伊1NPB=
ψ
， BNP=180
0ー 伊ZDNP=
伊2，
BD=K
θ
1
=90-(to
，
(
;
}
=
9
0
0 -(/)，
θ=90
0ーφ
伊1=A-180
0(
A
は東経)
伊2=180
0ー
λ
(
。
-λ)
となる。したがって (1
)
-
-
(
5)に相当する式は次のようになる。
c
o
t
θ
cot(K
ー
θ
1
)
=
一
一
一
一
-
…
.
.
.
・
H・
.
.
.
.
.
・
H・
.
.
.
.
.
・
H・
H・
H・
.
.
.
.
.
・
H・
.
.
…
.
.
.
・
H・
.
.
…
…
(
1
4
)
¥..u:.CP2s
i
n
(K-8
1)t
a
n
伊l
=
t
a
n
伊2…
.
.
.
.
.
・
H・
'
"
・
H・
.
.
.
.
・
H・
.
.
.
・
H・
.
.
.
.
.
・
H・
.
.
.
.
.
・
H・
.
.
(
15
)
sinK
tanK
t
a
n
θ=
二 ご7
・
.
.
.
.
.
.
…
・
・
…
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
ο
6
)
¥..u:'CPtまたは
日本付近の地磁気経緯度を求める早見表
sinψ=sin 8
1ま
τ
(
1
7
)
~ Sln CD2=
s
i
n
(Jl
一子三三…...・
H・
"
…
.
.
.
・
H・..………'"・
H・
.
.
…
…
.
.
.
・
H・
.
.
(
1
8
)
s
m
σ
次に日本付近の角度の存在範囲を考える。
(
1
4
)
式で
0<θ<900
，
0<
伊
2<900
であるから
0<K-8
1<90
0……… .
.
.
・
H・
.
.
.
.
.
・
H・
.
.
…
.
.
.
・
H・
.
.
.
.
.
・
H・
-
…
・
…
.
.
.
.
.
.
・
H・
.
.
(
1
9
)
したがって
0<K<90
0.
.
.
・
H・
.
.
.
.
.
・
H・
-
・
…
…
.
.
.
・
H・..……・・………...・
H・
.
.
.
.
.
.
.
.
…
・
(
2
0
)
(
1
5
)
式より
0<
伊
1<900
または
2
7
00<
引く
3
ω
。
…
.
.
.
・
H・
.
.
…
.
.
.
・
H・
.
.
…
.
.
.
・
H・
.
.
…
…
(
2
1
)
(
1
6
)
式より
t
a
n
8
・
c
o
s
伊
l=tanK
・
・
H・
H・
.
.
.
.
・
H・
-
…
・
・
…
H・
H・
.
.
.
.
.
・
H・
.
.
.
.
・
H・
.
.
.
・
H・
.
.
.
.
・
H・
(
2
2
)
(
2
0
)
式より
θ
，伊
lに次の関係がある
o-900<
伊
，
<900
ならば
O
話。
<900
または
1
8
00
三
8<2700
または
900<
伊
1<2700
ならば
9
00<8
孟
1
8
00
または
2700<
θ
三
3600
しかし
8
の値は
伊豆
8~三180
0したがって
O
と伊の聞には次の制限がつく。
-90
0<伊1<90
0と
0~8<900 …………'"・ H ・...・ H ・..……...・ H ・ (23)
または
900<
伊
1<2700
と
9
00<8
三
1
8
00
.
.
.
・
H・
.
.
…
.
.
.
・
H・"………(
2
4
)
(
1
7
)
と
(
1
8
)
式から
したがって
ゆえに
s
i
n
8
s
i
n
伊
l
=
s
i
n
θ
s
i
n
伊
z
…
…
.
.
.
・
H・..………'"・
H・
.
.
.
.
.
…
・
・
…
・
…
(
2
5
)
s
i
n
8
s
i
n
伊1>0
0
;
;
;
;
三
O
三
1
8
00
と
0
孟納豆
1
8
00
………・… .
.
.
.
・
H・
.
.
.
.
.
.
・
H・
.
.
…
(
2
6
)
1800~三θ 孟 3600と
1800~伊1 三3600・………
H・
H・
..，…....・ H ・ '..(27)6
9
(
2
3
)，
(
2
4
)
，
(
2
6
)
，
(
2
7
)
式で
8
，伊
1
の共通に存在するのは
(
2
3
)と
(
2
6
)
式および(
2
4
)と
(
2
6
)
式の
2
つの組合せであるから，
(
2
6
)
式のみとなる。しかるに
(
2
1
)
式 で 引 に 制 限 が つ く の で
次の条件が求めるものとなる。
0<8<900
と
0<
伊
，
<900
.
.
.
・
H・
.
.
.
.
.
・
H・
.
.
…
…
.
.
.
・
H・"…・・…・………・・・
(
2
8
)
以上の(1
9
)
，
(
2
0
)
，
(
2
8
)
式から
K-(}I
，
K
，
θ
，
伊
，
~ìl 、ずれも第一象限内にのみ存在し
ているから，符号は一切考えないでよ
L
。
、
したがって(1
4
)
-
-(
1
8
)
式は角度の符号のとり違いが無くなり，
計算誤りがなくなり，
3
.
1
で、述べた第
1
の方法よりはるかにすぐれてし、る。しかしこの方法も場所が大幅に変ると符号
が変るからまた若干の変形の必要がある。
前節と同様に計算例を示す。
この計算法は角度がすべて第一象限内であるから，計算誤りが少なし、。しかし計算の験算
まで入れると
3
0
.
.
.
.
.
.
6
0
分の時間を必要とし，決して少ない時間ではない。
3.3
地理的経緯度，北地磁気極の経緯度の変化に対する地磁気経緯度の変化
観測所の地理的経緯度が少し変ったり，北地磁気極の経緯度が変ったりした場合に，地磁
Table 3 Example o
f
c
a
l
c
u
l
a
t
i
o
n
o
f
t
h
e
m
ag
n
e
t
i
c
c
o
o
r
d
i
n
a
t
e
a
t
Kakioka (
P
a
r
t
2
)
N
o
.
F
a
c
t
o
r
(
1
)
Ao(2)
や
q(3)
A
(4)
や
(5)
θ
(6)
l
o
g
c
o
t
5
9
(
2
9
0
0
-
1h
(7)
(8)
l
o
g
c
o
s
C
P
2
(9)
l
o
g
c
o
t
(
K
-0
1
)
(
1
0
)
K-0
1
(
1
1
)
0
1
(
1
2
)
K
(
1
3
)
l
o
g
t
a
n
C
P
2
(
1
4
)
l
o
g
sin(K-0
1)(
1
5
)
l
o
g
s
i
n
K
(
1
6
)
l
o
g
t
a
n
'1'1
(
1
7
)
C
P
l
(
1
8
)
d
(
1
9
)
l
o
g
tanK
(
2
0
)
l
o
g
C
O
S
C
P
l
(
2
1
)
l
o
g
t
a
n
0
(
2
2
)
。
(
2
3
)
φ
(
2
4
)
l
o
g
s
i
n
O
l
(
2
5
)
l
o
g
s
i
n
C
P
l
(
2
6
)
l
o
g
s
i
n
θ
(
2
7
)
l
o
g
s
i
n
1
J
r
(
2
8
)
型TValue
2
9
1
。
7
8
.
5
。
1
4
00 1
1
'
3
60 1
4
'
5
30 4
6
'
9
.
8
6
4
9
7
6
2
90 1
1
'
9
.
9
4
1
0
4
6
9
.
9
2
3
9
3
0
5
00 0
0
'
1
10 3
0
'
6
10 3
0
'
9
.
7
4
7
0
2
3
9
.
8
8
4
2
5
4
9
.
9
4
3
8
9
9
9
.
6
8
7
3
7
8
2
50 5
7
'
2
0
50 5
7
'
0
.
2
6
5
2
3
6
9
.
9
5
3
8
4
5
0
.
3
1
1
3
9
1
6
3
。
ー
5
9
'
2
6
。
ー
0
1
'
9
.
2
9
9
6
5
5
9
.
6
4
1
0
6
4
9
.
9
0
6
6
6
7
9
.
0
3
4
0
5
2
60-13'
Note
s
t
l
o
n
g
i
t
u
d
e
r
t
h
l
a
t
i
t
u
d
e
5
t
l
o
n
g
i
t
u
d
e
r
t
h
l
a
t
i
t
u
d
e
-(4)
ogcot(5)
80
0ー{(
1
)
ー
(3)
}
og c
o
s
(
7
)
)
ー (8)
)→
a
n
t
i
l
o
g
.
ー(
2
)
0
)
+
(
1
1
)
ogtan(7)
og 5
i
n
(
1
0
)
og 5
i
n
(
1
2
)
3
)
+
(
1
4
)
ー(
1
5
)
6
)
→a
n
t
i
l
o
g
.
7
)
+
1
8
00
E
a
s
t
g
t
a
n
(
1
2
)
og c
o
s
(
1
7
)
9
)
ー(
2
0
)
1
)
→a
n
t
i
l
o
g
.
。
一
(
2
2
) North
g
s
i
n
(
1
1
)
og s
i
n
(
1
7
)
og s
i
n
(
5)
(
2
4
)
+
(
2
5
)
ー2
6
)
(27) ー~antilog.E
a
s
t
気経緯度がどの程度変るかについては， (1)--(5)式を微分すれば求まる
oこれについては
山口の計算がある
ω
が，若干誤りがあるので訂正した。
(
t
a
n
θ+cot (
)
)
8
8
ー
t
a
n
伊1
O
伊l dK=-\V-~' ~，
~""'V:: ~~_
v
:
-
:
.
'Y，...'Y'.
.
.
・
H
・
.
.
.
.
.
・
H
・
.
.
…
…
…
.
.
.
.
.
.
.
.
・
H
・
H
・
H
・
-
…
.
.
.
・
H
・
(29)
K+cotK
c
otK-cot (
K
-
(
)
1)
t
ー←，
a
n
(
)
+
c
_
.
.
o
.
.
.
v
t
(
)
~O
(
)
伊2
t
a
n
K+cot K •
t
a
n
伊2+Cot
伊2
+{~m 伊1+ ∞t 伊1
似K
ー 似(
K
-
(
)
1
)
t
a
n
伊I
l
早l
t
a
n
伊2+Cot
伊2
t
a
n
K
+cot
K
t
a
n
伊2+Cot
伊/0
伊tc
o
t
(
K
-
(
)
1
)
δ
(
)
1
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
・
・
・
・
・
・
・
・
・
(
3
0
)
t
a
n
伊2+C
O
t
伊2
日本付近の地磁気経緯度を求める早見表
dθ=t_t
竺空
E
，.c
:
o
tK
ー
cot(K
ー
θ
1
L
.
ぬ
n8+ωt8
{tanθ+cotθtan
K
+cot
K
t
a
n
伊2
+
C
o
t
伊2
+
_
!
a
n
8+cotθtan
.
(K-8
1)
+
ー
c
o
t
(K-81
)
1
- rt
a
n
K+cot K
t
a
n
θ+cotθJ
+[Uan
(
/
)
1
+∞t~-型企∞t(K-8
1
2. t~n ~~~t=J
・ 旦 位
l
t
a
n
伊
2+
'
c
o
t
伊
2
t
a
n
K
+cot K
t
a
n
伊
2+C
o
t
伊
2
' t
a
n
θ+cotθ
t
a
n
(K-8
1)+cot
(K
-(1
)
t
a
n
伊1
1
.
.
一
一
旬nθ+C
瓦吾一一
.
t
a
n
玄示。
tKJO
伊1
( 州
'K-8
1L._t
旬a
n
軌一一一
t
旬a
叫
n
(K
一
θ
川
o
t(K
一仇
-
一
「
δ0
仇
1.
…
.
リ
.
，
…
.
口
"
…
.
υ
.
.
日.
.
…
(
σ
3
1
υ
)
t
a
n
伊2+C
∞
o
悦t
伊2
句nθ+c
∞
o
叫t θ t a nθ+c
∞
ο
瓜
t θ j
I
. L _ ， _ . _ _ ' -C
o
t
K-cot (K-8
1)t
a
n
8+cot8
'
-
_
_
.
_
_
.
L
I
1
)
8ψ=itanψ
l
~_...."・
c
_
_
o
V
t
伊T "2
・
t
a
n
K
+cot
K
・
t
a
n
伊2+Cot
一 泊nψcot8tS8
伊
![Table 8 L i s t o f Geomagnetic c o o r d i n a t e s o f magnetic o b s e r v a t o r i e s i n ]apan ( c o t i n u e d )](https://thumb-ap.123doks.com/thumbv2/123deta/5825500.541649/19.730.85.649.203.973/table-l-geomagnetic-magnetic-v-e-apan-n.webp)
