大阪都心部における地下温暖化の実態○有本弘孝 1・北岡豪一 2・谷口真人 3・濱元栄起 41(株)地域地盤環境研究所・2 岡山理科大学理学部・3 総合地球環境学研究所・4 埼玉県環境科学国際センター

1. はじめに都市域における都市化の進行やヒートアイランド現象は，地表面温度を上昇させ，地中では熱拡散によって表層から深さ方向に地下温度を徐々に上昇させている。すなわち，“地下の温暖化”である。地下温度が上昇すると地層から化学成分の溶出や汚染物質の離脱が促進され，土壌水や地下水の水質だけでなく地下の微生物活動にも影響が生じる可能性がある。これら地下温暖化に伴う地下環境への影響は，隠れた熱汚染問題であり，現在，各方面で研究が進められている。
一方で地下の温暖化には，その温暖化以前の地中熱を対象とした再生可能エネルギー以上に地中熱が蓄積されているため，その蓄熱分を活用できるという利点もある。地下温暖化の将来にわたる変化を予測し，その対策や活用を講じるためにも地下温度の正確な客観的事実データを得ておくことは重要な課題であるといえる。
大阪平野における地下温暖化の調査研究は，1998 年，2003 年，2011 年，谷口ほか 1) 2) により主に大阪湾沿岸部と郊外を対象に実施された。その一方で 2009 年から，都市化が古くから始まった大阪都心部の中之島周辺を対象とした調査研究が，筆者らが所属する「地下水地盤環境に関する研究協議会」の地下水水質研究委員会(委員長：大阪工業大学 青木一男教授)によって実施され，その成果の一部を北岡ほか 3) 5) と有本ほか 4) 6) により報告済みである。今後，地下温暖化による地下環境の変化がもたらす地中構造物等への影響を調査研究する前段階として，今，都市化の最も進んだ都市中心部における地下温暖化の実態を正確で詳細に調査研究することは極めて有意義であると考えられる。なお，本報告において，大阪都心部とは JR 大阪環状線内（ただし西側は新なにわ筋まで）のエリアとする。
本稿は，大阪都心部の８つの観測井で観測された地下温度（地下水位観測井の井戸管内の水温鉛直分布）について，その温度分布を大阪平野の約 30 地点のそれと比較し，大阪平野から見た大阪都心部の地下温暖化の実態（地下温度の高さ，地下温度の上昇速度）について報告するものである。

温 度 （℃）

2. 地下の温暖化を示す地温鉛直分布15

地上の温暖化に伴う地温形成のメカニズムは，熱的性質が一様な半無限の媒0体（地層）で表面（地表）の温度がステップ状に上昇したとして，媒質内温度分布の変化を熱伝導理論を用いて計算すると理解しやすい。図-1 は，地表面温度温分布として，地表面温度がステップ状に 3.0 ℃上昇した時の地温鉛直分布の変化5)を示したものである。この図より，地下温暖化を示す地温分布の特徴は，時間が大きく経過しない範囲では，ある深さで極小となる地温分布ができ，その極小となる深さは時間とともに深くなりつつ全体的に温度が上昇することである。この温度の上昇部分（増分）が地下温暖化で蓄積される地中深 度 （G.L.-m）

15.0 ℃，地温勾配 0.030 ℃/m の直線温度分布を温暖化開始前の定常的な地205102017181920

蓄熱
50ｙ
1006020080

θ s 0 = 15.0 °C
地温勾配： 0.030 °C/m ，
熱拡散率 κ : 5.0 × 10 −7 m 2 /sec

100

図-1
3. 地下温度の観測

1y240

熱（蓄熱）であり，地表面温度が一旦上昇した後，地表面温度の上昇が無くても，地温は地層の熱伝導により上昇する。

16Δθ = 3.0 °C

500y
1000y

地表面温度がステップ状に上昇した時の地温分布の変化 5)(1) 大阪都心部で地温観測した井戸
地温の観測井戸は，地下水協議会が所有する広域地下水位観測井のうちの 6 井戸（N1～N6）と，大阪市が所有する中之島 B および馬場町(Ⅱ)観測井で，それらの位置を図-2 に示す。なお，N3 観測井は上町断層の西方約 150 m 離れ，馬場町観測井は大阪管区気象台の観測位置に近接している(離隔約 30 m)。N1～N6 観測井の孔底深度は，N3 井戸では第 1洪積砂礫層(Tg1)の被圧地下水位を観測するために G.L.-約 40 m，その他の井戸では第 2 洪積砂礫層(Tg2)を対象とするために G.L.-60 m 前後である。一方，中之島 B および馬場町観測井の孔底深度は，大阪広域地盤沈下を観測するために前者では G.L.-約 180 m，後者では G.L.-約 150 m と深い。
(2) 地温観測日
表-1 に地温観測日を示す。2011 年 11 月～2012 年 8 月に実施した N1～N6 観測井の地温観測は，四季における地温分
Actual states of subsurface warming in the central area of Osaka City, Japan
Hirotaka Arimoto1, Koichi Kitaoka2, Makoto Taniguchi3, Hideki Hamamoto4 ( 1Geo-Research Institute, 2Okayama University ofScience, 3Research Institute for Humanity and Nature, 4Center for Environmental Science in Saitama )KEY WORDS: urbanization, Osaka area, observation well, temperature

布変化の把握を目的としたものである。

JR大阪駅(3) 地温の観測方法

N3N2

井戸管内の水の温度は，井戸管の材質や井戸管と地層N1との密着度などの影響を受けるが，時間変動が小さい場

中之島Ｂ合は，周囲の地層の温度にほぼ等しく，近似的にその深御堂筋

N5さの地温を表しているものと考えられる。通常の地温勾

N4配では，大口径の管内でも水の対流は生じない。

N6N3N2観測は，温度と圧力を同時に計測できる測定器(メモリ馬場町(Ⅱ)
大阪管区

N1ーを内蔵)を手動の巻き上げ機を用いて井戸管中を約 20.5 kmcm/sec の速度（記録間隔 1 sec）で降下させることによっGoogleマップから引用て行った。測定器は JFE アドバンテック社製の

図-2

大阪都心部における観測井の位置

Compact-TD 計で，温度の測定精度は±0.05 ℃，分解能は表-10.001 ℃，水深の精度は 0.3%F.S.，分解能は 0.76 cm，水深 500 m まで計測可能である。
4. 中之島周辺の N1～N6 観測井の地温観測結果
図-3（次ページ参照）は，N1～N6 観測井全てについて，観測された地温鉛直分布，年平均気温(大阪管区気象台)，ならびに土質柱状図，井戸スクリーン位置を示したものである。各々観測井における地温データの始まりは，井戸管中の水位面を表している。

観測井No.，名称
N1（福島公園）
N2（西梅田公園）
N3（西天満公園）
N4（中之島西公園）
N5（西船場公園）
N6（市道道修町線）
NK（中之島B）
18.馬場町(Ⅱ)

地温観測日

2009年

2012年

2011年

12/15
4/28

11/62/4

5/198/119/29

12/14－12/15
2005/10/4
1998/12/4

2003/10/15

2011/6/14（注）地温データの所有者
N1～N6観測井，中之島B観測井：地下水水質研究委員会の地温WG（主査：北岡教授）
馬場町(Ⅱ)観測井：総合地球環境学研究所の研究プロジェクト（研究代表者：谷口教授）

図-3 の観測結果には，以下のような特徴が見られる。

19.4① 地温そのものと地温分布のパターンが観測井によってかなり異なっている。

N3 0.053 ℃/年（G.L.-38.5m 地点）② 地表面から深度 13～20 m 6) までの地温には季節変動が生じている。
③ 季節変動のない深度 20 m 以深に着目すると，地温が全般的に低い N1 と N5
観測井では，極小を示すような分布ではなく深度に対して直線的な分布となっており，その時間変化は小さい。ところが，N3，N4 および N6 観測井では時間経過とともに地温がほぼ一様に増加している。
④ N2 観測井では深度約 30 m を中心に上凸，N6 観測井でも深度約 40 m を中心に上凸の温度分布のパターンを示す部分が見られる。この上凸の中心深度の土層は，N2 では砂質土層（As(L)），N6 では砂礫層（Tg1）と粘性土層（Ma12）の層境に対応している。

G.L.-40m 地点の地温 （℃）

19.2

N6 0.054 ℃/年

19.018.8
18.6

N2 0.036 ℃/年

18.4

N4 0.018 ℃/年

18.2
18.0
N1 0.007 ℃/年

17.8

上述③の地温の一様増加について検討する。図-4 は，地温の季節変動深さの下端(深度 20 m)から井戸孔底(深度 60 m)の中間深度である G.L.-40 m 地点について，地温の経時変化を示したもので，地温の変化を直線近似した地温上昇速度も明示している。図によると，地温が高い井戸ほど地温の上昇速度が大きくなっている。地温の17.6

N5 0.001 ℃/年

20092011

2013

時 刻
図-4 G.L.-40m 地点の地温経時変化

高い N3 および N6 観測井の地温上昇速度は 0.05 ℃/年で，こ年～2012 年）の気温の平均上昇速度 0.020 ℃/年や過去 63 年間（1950 年～2012 年）の平均上昇速度 0.028 ℃/年よりも約2 倍程度大きい。一方，地温の最も低い N5 観測井では，地温はほとんど上昇していない。
5. 大阪都心部の地下温暖化の実態

19

年平均気温 （℃）の速度は，図-5 に示す大阪管区気象台の過去 130 年間（18831950 - 2012（63年間）
y = 0.028 x - 38.999

18171615

1883 - 2012（130年間）
y = 0.020x - 23.341

14131880

2011 年春，連名者谷口を研究代表とする共同研究グループ

1900

1920

1940

1960

1980

2000

2020

観測年

5.1 大阪平野の地温状態
図-5大阪管区気象台における気温の経年変化により，大阪平野の約 30 地点で地温鉛直分布の一斉調査が行われた。この一斉調査は 1998 年と 2003 年の一斉調査を引き継いだもので，同じ井戸での観測を行って地下温暖化の進行状態を把握することが目的である。2011 年一斉調査時の観測井位置図を図-6 に，地温分布の全観測結果を図-7 に示す。図-7 によれば，深層（＞ 50 m）におけるほとんどの地温は，深度に対して，2 本の平行線（地温勾配：0.033 ℃/m，間隔：3.5 ℃）の間に収まる。一方，浅層（＜ 50 m）では，深層の地温勾配から外れ，高温側に湾曲する。これは，

温 度 （℃）
15161718

温 度 （℃）
192016-101718192021-102012 2011 2009 20102012 2011 2009 2010土質柱状図

00As(U)As(U)1020As(L)30090428

Tg1091215

40111106
12020450120519

深 度 （G.L.-m）

深 度 （G.L.-m）

10Tg2

20As(L)

30090428
11110612020450120519

12081112092960

年平均気温（大阪管区）
井戸スクリ－ン70

年平均気温（大阪管区）
井戸スクリ－ン

70(a) N1（福島公園）(b) N2（西梅田公園）

温 度 （℃）

温 度 （℃）
1718

1920

21

1722

18

1920

21

22-10-10

2012 2011 2009 20102012 2011 2009 2010

＋0 ＋＋ ＋

10

100As

20As

30090428

Tg1091215

40111106
12020450120519
120811

深 度 （G.L.-m）

深 度 （G.L.-m）

Tg2

12081112092960

Tg1

091215403040120519
120929年平均気温（大阪管区）-1018(d) N4（中之島西公園）

温 度 （℃）
1920

1718

1920

21

22-10

2012 2011 2009 20102012 2011 2009 2010

＋0 ＋

0As(U)

As(U)

1020

40As(L)
090428Tg1

09121411110612020450120519
1208116070

深 度 （G.L.-m）

深 度 （G.L.-m）

1030

Tg2

井戸スクリ－ン

温 度 （℃）
17

Tg1

1202045070(c) N3（西天満公園）

16As(L)091214120811

井戸スクリ－ン

15090428
11110660

年平均気温（大阪管区）

70

20Tg2120929

60

As(U)

As(U)

20As(L)

30091215

12020412051950

120811Tg2

12092912092960

年平均気温（大阪管区）

Tg111110640

年平均気温（大阪管区）

Tg2

井戸スクリ－ン

井戸スクリ－ン

70(e) N5（西船場公園）(f) N6（市道道修町線）
土質柱状図の凡例：砂礫図-3

大阪中之島周辺における地温鉛直分布(観測値)と土質柱状図：砂：粘性土温 度 （℃）

B
132

1A
5
3 4
67

15

19

21

23

25

01112
14
15 16
91718
C
1920
21

3.5 ℃

1710

深 度 （G.L.-m）

2250

FＦ

2425
26
D

100

28
29 30
320.033℃/ｍ

150

310.033℃/ｍ

E342003.5 ℃

図-62011 年大阪平野一斉調査における観測井戸の位置図

図-7図-1 で示した地下温暖化の特徴を明らかに示し

全観測井戸における地温鉛直分布（観測値）

1921

温 度 （℃）

温 度 （℃）

温 度 （℃）
17ている。

231725192115

230005050501719

３回の一斉調査（一部の観測井では 1998 年調査

これら土地利用形態に違いがあるにも係わらず地温は上昇している。このような地温上昇は一つの観測井（後出する表-2 の No. 34 泉南）を除く他の100150観測井全てで認められたことから，大阪平野全体において地下温暖化が着実に進行していると考え

200

られる。

観測日
左側の分布から順に
1998/12/04
2003/10/15
2011/06/01

250

100

100

150

150

200

観測日
1998/12/04
2003/10/15
2011/06/14

19.港Ⅱ(B)

観測日

20018.馬場町(Ⅱ)
250

図-85.2 地下温暖化の進行性を評価する指標

深 度 （G.L.-m）

共公益施設用地，No.19 港は市街地の総合公園近傍，No.13 交野は郊外の小学校用地に位置する。

深 度 （G.L.-m）の一例を図-8 に示す。No.18 馬場町は都心部の公深 度 （G.L.-m）

はされていない）で得られた地温分布の経年変化

250

2003/11/11
2011/05/25

13.交野

地温鉛直分布の繰り返し観測結果の一例

地下温暖化の進行性を評価する指標として，

温 度 （℃）

図-1 の地温分布変化を考慮すると，以下の３つが考えられる。

16① 地温鉛直分布の極小値の深さとその地温

171819200② 地下が温暖化する前の定常地温分布からの乖離開始深さ（地下N1 N510温暖化の深さ）③ GL-40 m 地点の地温とその上昇速度

N2中之島B
N3N4N620①については，観測井全てで地温分布に必ずしも明瞭な極小値が認

30

浅く，地下が温暖化する前の定常的な地温分布が定かでない。③について，G.L.-40 m 地点は，2003 年大阪平野一斉調査の結果から，大阪平野の地下温暖化の深さは GL-約 80 m2)まであり，その中間深さをとったものである。この③については，観測井の孔底深度が 40 m 以上あれば評価することができ，①や②の指標よりも指標値の算定に客観

深 度 （G.L.-m）

められない。②については，ほとんどの観測井の深さは 100 m 未満と405060性があることから，当面の指標として③を採用することとする。以後，

80

GL-40 m 地点の地温を GL-40 m 地温と呼ぶ。

905.3 大阪平野から見た大阪都心部の地温状態(1) 地温鉛直分布
大阪都心部の地温状態について，大阪平野全体の地温分布から見た

18.馬場町70100

図-9

大阪平野から見た大阪都心部の地温鉛直分布（観測値）

21

場合の検討を行う。地温鉛直分布について，大阪平野から見たB大阪都心部の状況を図-9 に示す。N1～N6 観測井のデータは，A2011 年一斉調査と同じ年の観測データ(2011 年 11 月)である。
季節変動の影響を受けない深度 20 m から 50 m 程度までの浅層の地温に着目すると，N3 西天満の地温が最も高く，次いで N6C道修町が高い。そして，No.18 馬場町，N2 西梅田，N4 中之島西と続く。これらの地温分布は，大阪平野全体から見れば地温FＦの高い範囲に位置している。一方，中之島 B や N5 西船場，N1福島の地温は，大阪平野全体から見れば中程度の高さと言えよDう。
(2) G.L.-40m 地温の分布とその上昇速度
各々観測井の地温鉛直分布から G.L.-40 m 地温を読み取り，温度（℃）の凡例
17.0以上 17.5以上 18.0以上 18.5以上 19.0以上
17.0未満 17.5未満 18.0未満 18.5未満 19.0未満

その温度を大阪平野の白地図へプロットした結果を図-10 に示Eす。プロットしたデータには，表-1 に記載のない郊外の観測井（大阪市大正区の 1 井戸，柏原市の 1 井戸）で観測したデータも含まれている。また，気象庁測候所と観測井の位置関係から，

図-10

G.L.-40m 地温の平面分布図

観測井を測候所毎に目安としての13.交野高い傾向を示す。このことは，郊外よりも沿岸部，都心部ほど都市化の進行に伴う地下温暖化が進んでいることを意味している。

温 度 （℃）

図-10 よると，地温は郊外ほど低く，反対に沿岸部や都心部では： GL-40m地温

大阪都心部

10.柴島

： GL-40ｍ地温の上昇速度

20

過去63年間の気温平均上昇速度
過去130年間の気温平均上昇速度

190.060.04

180.02170.0016-0.0215G.L.-40 m 地温とその上昇速度

温度上昇速度（℃/年）

グループ化も行っている。

-0.04
18 N6 N3 N2 N1 NK N5 N4

1, 2, 3, 4, 5, 6, 13, 7, 9,10,11 ,12,14,15,16,17

についてまとめた結果を表-2 と19 ,20 ,21 ,22 ,24 ,25 ,26 ,28 ,29 ,30 ,31 ,32 ,34

観測井No.

図-11

G.L.-40m 地温とその上昇速度の分布

図-11 に示す。表では，測候所毎に G.L.-40 m 地温の上昇速度の平均値を参考値として付記した。
表-2 と図-11 によると，以下のことが分かる。
① No.34 泉南を除いて，G.L.-40 m 地温の上昇速度は観測地点に

表-2
測候所
No.

よって異なり，0.001～0.054 ℃/年の範囲にある。特に大阪都心部ではその差異が顕著である。

A（豊中）② 地温上昇速度が最も大きい地点は，大阪都心部の N6 道修町で，次いで，No.18 馬場町，N3 西天満である。これら地点

B（枚方）の地温上昇速度は 0.050 ℃/年以上の速さで，大阪管区気象台の過去 130 年間の気温平均上昇速度 0.020 ℃/年や過去 63年間の気温平均上昇速度 0.028 ℃/年の約 2 倍程度大きい。
③ 大阪管区気象台の気温平均上昇速度 0.020 ℃/年以上となる観測地点の数は，計 23 地点で，全 34 地点の約 2/3 にも及ぶ。
④ 一方で，地温上昇速度が最も小さい地点は，大阪都心部のC（大阪管区）

大阪都心部

N5 西船場で，次いで No.4 口酒井，N1 福島である。いずれも 0.01 ℃/年未満である。
G.L.-40 m 地温の上昇速度が気温平均上昇速度を上回っている
地点は，一つの解釈として，都市化の進行に伴う地表面温度の上昇速度が気温の上昇速度よりも大きい可能性が考えられる。これ

F（八尾）
D（堺）

らの解析や，土地利用状態の変化と地表面被覆形態の違いが地温に及ぼす影響等については，今後の研究課題としたい。
5.4 大阪都心部の地下温暖化
大阪都心部に位置する N1～N6 観測井，中之島 B および No.18
馬場町の地温鉛直分布を重ねて図-12 に示す。図より，N1～N6E（熊取）

G.L.-40m 地温とその上昇速度

観測
GL-40m
G.L.-40m地温の観測
井戸
地温上昇速度 (℃/年）
所名No.
（℃） 2003～2011年* 平均値
荒牧
17.30
0.017
1
石橋
***
18.48*
2
野間
3
17.25
0.035
0.024
口酒井
17.150.005
4
曽根
17.030.035
5豊中
17.260.029
6
交野
0.021
1315.990.021
吹田
718.770.038
十三
17.59
0.018
9
柴島
10
17.89
0.047
鳥飼西
11
16.67
0.020
点野
***
***
12
庭窪
14
17.42
0.037
0.041
15 新森小路 18.85
南郷
1617.440.033
鴻池
1717.93
0.029
0.032
18N6
N3
N2
N1
NK
N5
N4

192021
2224
25
26
28
2930
313234

馬場町
道修町
西天満
西梅田
福島
中之島B
西船場
中之島西港
長瀬
加美東
堺5
堺2
堺A
堺B
岸和田3
岸和田2
岸和田1貝塚
泉佐野
泉南

18.63
18.97
19.27
18.41
17.86
17.86*
17.74
18.31

0.053
0.054
0.053*
0.036
0.007－0.001
0.018

18.52
17.55
17.49
17.91
18.84
18.83
19.01
18.62
18.40
17.19
17.26
17.77
18.67

0.024
0.025
0.029
0.013
0.030
0.023
0.016
0.036
0.018
0.035
0.030
0.017
-0.030（都心部を除く）

0.027
0.021

0.034

（注） *** ：G.L.-40m±5m内に地温データなし－ ：地温上昇速度のデータなし
18.48* ：G.L.-44.0m地点での地温17.86* ：2005/10/04の地温0.053* ：G.L.-38.5m地点での地温上昇速度
大阪都心部の地温上昇速度：2008年～2012年のデータ

温 度 （℃）

観測井と中之島 B では，中之島 B の地温勾配で表され13る温度分布 a 線を，No.18 馬場町では馬場町の地温勾0始まる前の定常的な地温分布を表すものとすると，大込まれる。
6. まとめ大阪平野および大阪都心部において，地下の温暖化は着実に進行中である。観測 34 地点中 23 地点において，G.L.-40m 地温の上昇速度は，過去 130 年間の大阪管区気象台の気温平均上昇速度よりも大きい。
地下温暖化の進行性の指標として，G.L.-40m 地温と深 度 （G.L.-m）

G.L.-20m で約 3℃，G.L.-40m で約 1.5℃の地温上昇が見

1921232527

年平均気温20

阪都心部における地下温暖化の深さは地点毎に異なる
化 の 進 行 し た N3 西 天 満 で は , 温 暖 化 開 始 前 か ら

17-20配で表される温度分布 b 線を，それぞれ地下温暖化がが，概ね G.L.-50～80 m 程度と推測される。最も温暖

154060N1N5160180

地下温暖化の範囲

概ね50～80m程度

中之島B120140

a線，b線を地下温暖化
前の定常状態を表すものとした場合の

1 8 . 馬場町(Ⅱ)

80100

N2N3N6N4

a線 （中之島B）
深度 90～120 m間
地温勾配 0.0220 ℃/m
地表切片 16.68 ℃

b線

a線

b線 (馬場町(Ⅱ)）
深度 110～150 m間
地温勾配 0.0409 ℃/m
地表切片 15.48 ℃（注） N1～N6の地温：2012/02/04（冬季）と2012/08/11（夏季）の平均値
中之島Bの地温：2005/10/04 の値
18.馬場町(Ⅱ)の地温：2011/06/14 の値
年平均気温：大阪管区，2011/08/12～2012/08/11 の平均値

図-12

大阪都心部における地下温暖化の状況

その上昇速度で評価すると，都心部の N3 西天満のG.L.-40m 地温は全観測地点の中で最も高く（19.27 ℃），その上昇速度も大きい（0.05 ℃/年）。これは大阪管区気象台の気温平均上昇速度の約 2 倍程度の大きさで，ローカルな地下温暖化が早くから起こっていた可能性がある。
中之島 B および馬場町の地温勾配で表される温度分布を，地下温暖化が始まる前の定常的な地温分布を表すものとすると，大阪都心部における地下温暖化の深さは観測地点毎に異なるが，概ね G.L.-50～80 m 程度と推測される。
謝辞：
本報告は，地下水地盤環境に関する研究協議会に設置された「地下水水質研究委員会」の研究成果として報告するものである。中之島周辺の地温観測に際しては，中之島高速鉄道（株），京阪電鉄（株），地下水協議会から数々の便宜供与をいただいた。また，2011 年の大阪平野一斉調査は，住友財団研究助成金（研究代表者：谷口真人）により実施されたもので，地温観測に際しては，国土交通省近畿地方整備局淀川河川事務所ならびに猪名川河川事務所，大阪府環境農林水産部，大阪市環境局から数々の便宜供与をいただいた。あわせてここに記して謝意を表します。
参考文献
1）Taniguchi, M. and T. Uemura: Effects of urbanization and groundwater flow on the subsurface temperature in Osaka, Japan,Physics Earth Planetary Inter., 152, pp.305-313, 2005.
2）Taniguchi, M. Uemura, T., Jago-on, K. : Combined effects of urbanization and global warming on subsurface temperature in four
Asian cities. Vadose Zone Jour. 6(3), 591-596, 2007.
3）北岡豪一，吉岡龍馬，有本弘孝：地下温度の鉛直分布に記録されている気候温暖化履歴の解析，地下水地盤環境に関するシンポジウム 2006 発表論文集，地下水地盤環境に関する研究協議会，pp.115-124，2006.
4）有本弘孝，北岡豪一，谷口真人，吉岡龍馬，上村剛史：大阪地盤における地下温度鉛直分布の地域性，地下水地盤環境に関するシンポジウム 2007 発表論文集，地下水地盤環境に関する研究協議会，pp.51-60，2007.
5）北岡豪一，有本弘孝，吉岡龍馬：温暖化に伴う大阪平野の地温モニタリング，地下水地盤環境に関するシンポジウム
2010 発表論文集，地下水地盤環境に関する研究協議会，pp.69-80，2010.
6）有本弘孝，北岡豪一，谷口真人，濱元栄起：大阪中之島周辺における地下温度鉛直分布の特徴，地下水地盤環境に関するシンポジウム 2012 発表論文集，地下水地盤環境に関する研究協議会，pp.27-38，2012.