大阪平野における地下温暖化の実態○有本弘孝 1・谷口真人 2・濱元栄起 3・岸本安弘 4・水間健二 5・小林 晃 61(株)地域地盤環境研究所・2 総合地球環境学研究所・3 埼玉県環境科学国際センター

4

国土交通省近畿地方整備局企画部・5 大阪府環境農林水産部・6 関西大学環境都市工学部

1. はじめに都市域における都市化の進行やヒートアイランド現象は，地表面温度を上昇させ，地中では主に熱拡散によって表層から深さ方向に地下温度を徐々に上昇させている。すなわち，“地下の温暖化”である。地下温度が上昇すると地層から化学成分の溶出や汚染物質の離脱が促進され，土壌水や地下水の水質だけでなく地下の微生物活動にも影響が生じる可能性がある。これら地下温暖化に伴う地下環境への影響は，隠れた熱汚染問題であり，現在，各方面で研究が進められている。一方で地下の温暖化には，地下温度上昇分の地中熱が蓄積されているため，その蓄熱分も活用できるという利点もある。地
下温暖化の将来にわたる変化を予測し，その対策や活用を講じるためにも地下温度の正確な客観的事実データを得ておくことは重要な課題であるといえる。
大阪平野における地下温暖化の調査研究は，1998 年，2003 年，2011 年，谷口ほか 1) 2) により主に大阪湾沿岸部と郊外を対象に実施された（以降，一斉調査と呼ぶ）。その一方で 2009 年から，都市化が古くから始まった大阪都心部の中之島周辺を対象とした調査研究が，筆者らが所属した「地下水地盤環境に関する研究協議会」の前地下水水質研究委員会(2013年まで活動)によって実施され，その成果の一部を北岡ほか 3) 5)，有本ほか 4) 6) 7)，濱元ほか 8) により報告済みであり，現在は同協議会に設置された地中熱研究委員会(委員長：関西大学 小林晃教授)に引き継がれている。今後，地下温暖化による地下環境の変化がもたらす地中構造物等への影響を調査研究する前段階として，また，大阪平野における地中熱利用ポテンシャルの調査研究を行うためにも，大阪平野における地下温暖化の実態を正確で詳細に調査研究することは極めて有意義であると考えられる。
本稿は，地中熱研究委員会の研究成果の一部として，大阪平野の 2011 年一斉調査で得られた 30 地点の地下温度（地下水位観測井の井戸管内の水温鉛直分布）の観測結果について，過年度一斉調査からの経年変化，大阪平野の地下温暖化の実態（地下温度の高さ，地下温度の上昇速度等）について報告するものである。
2. 地下の温暖化を示す地温鉛直分布
地上の温暖化に伴う地温形成のメカニズムは，熱的性質が一様な半無限の媒体（地層）で表面（地表）の温度がステップ状に上昇したとして，媒質内温度分布の変化を熱伝導理論を用いて計算すると理解しやすい。図-1 は，地表面温度 15.0 ℃，地温勾配0.030 ℃/m の直線温度分布を温暖化開始前の定常的な地温分布として，地表面温度がステップ状に 3.0 ℃上昇した時の地温鉛直分布の変化 5) を示したものである。この図より，地下温暖化を示す地温分布の特徴は，時間が大きく経過しない範囲
では，ある深さで極小となる地温分布ができ，その極小となる深さは時間とともに深くなりつつ全体的に温度が上昇することである。この温度の上昇部分(増分) が地下温暖化で蓄積された地中熱(蓄熱) であり，地表面温度が一旦上昇した後，温 度 （℃）
0

深 度 （G.L.-m）

2040

161y
2510201718

気温変化，人工排熱等

1920寒冷化

蓄熱 （地下温暖化）

ΔT

地下温暖化
（地温の逆転現象，蓄熱）

T050ｙｔT0 = 15.0 ℃
Δ T = 3.0 ℃
地盤：一層モデル
地温勾配G ：0.030 ℃/m
熱拡散率κ ：5.0×10-7 m2/sec

60
200

80500y

地表面温度の変動

Gの異常

100

100T0

温暖化

T（℃）

温度 Ｔ

年周（季節）変動
GL-10～-20m
攪乱域
都市毎に異なるT（z）＝T0 + G z
G : 地温勾配（℃/100m）

深度 ｚ

15100m G2～4 ℃

非攪乱域

地殻熱流量
（= 地温勾配×熱伝導率）

1000y

図-2 都市域における地温鉛直分布の一般的な特徴(模式図)の変化 5)

図-1 地表面温度がステップ状に上昇した時の地温鉛直分布の変化 5) 一部加筆

Actual States of Subsurface Warming in Osaka Plains, JapanHirotaka Arimoto1, Makoto Taniguchi2, Hideki Hamamoto3, Yasuhiro Kishimoto4, Kenji Mizuma5, Akira Kobayashi6 ( 1Geo-ResearchInstitute, 2Research Institute for Humanity and Nature, 3Center for Environmental Science in Saitama, 4Kinki Regional Development
Bureau, Ministry of Land, Infrastructure, Transport and Tourism, 5Osaka Prefectural Government, Department of Environment,Agriculture, Forestry and Fisheries, 6Kansai University )KEY WORDS: Urbanization, Osaka Plains, Observation well, Temperature
71地表面温度の上昇が無くても，地温は地層の熱伝導により上昇する。

B1321A53 46

図-2 は，地盤に地下水流動がない仮定で，都市域における地温鉛

117121415 1691718C192021

直分布の一般的な特徴を示した図である。地温鉛直分布は，地下10深部からの熱供給(地殻熱流量)と地表面温度 T0 により決定されるが，地表面温度変動のうち，年周変動以外の長周期的変動がより地下深部に伝導していき，地温勾配 G の異常(地温の逆転現象)をもたらす。この地温分布の攪乱域が地下温暖化(蓄熱)の範囲であ22り，その深さは都市毎に異なる。大阪で約 100m，東京で約 150mという解析結果 2)も報告されている。

FＦ

242526D2829 30

3. 地下温度の観測(1) 大阪平野で地温観測した井戸をそれぞれ図-3 と図-4 に示す。前者の井戸(●印)は，河川流域水3132

2011 年大阪平野一斉調査と大阪都心部の地温観測井の位置図

E34Observation well

文観測のための単管井戸 10 本，広域地盤沈下観測のための単管井

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

戸 8 本および二重管井戸 12 本の合計 30 本であり，図中には気象庁の測候所を▲印(A～F)で示している。後者の井戸は，地下水協JR大阪駅

議会が所有する単管井戸 6 本（N1～N6）と，大阪市が所有する単管井戸の馬場町(Ⅱ)と二重管井戸の中之島 B である。二重管井戸

N1N3N2とは，地下水位観測用の外管と地盤沈下観測用の内管(ロッド)の中之島Ｂ二重管構造(写真-2(a)(b)参照)を意味する。
(2) 地温の観測方法

N6N5N4御堂筋

N2

従来からの地温観測は，観測井管内の水温鉛直分布の測定が行

18. 馬場町(Ⅱ)
大阪管区気象台N3N1

われている。井戸管内の水の温度は，井戸管の材質や井戸管と地0.5 km
Googleマップから引用

層との密着度などの影響を受けるが，時間変動が小さい場合は，周囲の地層の温度にほぼ等しく，近似的にその深さの地温を表し

図-4 大阪都心部における観測井の位置図ているものと考えられる。通常の地温勾配では，大口径の管内でも水の対流は生じない。
単管井戸における地温観測を写真-1 に示す。観測は，温度と圧力を同時に計測できる測定器(メモリーを内蔵)を手動の巻き上げ機を用いて井戸管中を約 2 cm/sec の速度(記録間隔 1 sec)で降下させることによって行った。降下速度は経験上，2cm/sec 程度がベストであり，降下速度が 5 cm/sec 程度よりも速いと，温度の応答に若干のずれが生じたり，また，水に乱れを生じさせる可能性がある。測定器は JFE アドバンテック社製の Compact-TD 計で，水深 500 m まで計測可能である。
一方，二重管井戸における地温観測を写真-2 に示す。二重管井戸では，井戸深部での外管径の段落ちと外管と内管の間にあるセントライザーの存在(写真-2(b)参照)により，外管と内管の隙間が狭いことから，単管井戸の測定器よりも直径の小さい温度測定器を採用した。観測は，温度センサーケーブルに 1m 間隔の目印を付けた温度センサーを井戸内に手動で降下させながら，温度を深度 1m 間隔に高精度デジタル温度計で読み取る方法で行った。水深 400 m まで計測可能である。
ワイヤ

自記水位計
（ﾌﾛｰﾄ式）

観測器具

202mm

(b)井戸管頭部

温度測定精度
分解能
水圧測定精度
分解能±0.05 ℃
0.001 ℃
±0.3 F.S.
0.76 cm

φ40mm
温度センサー部（サーミスタ）
圧力感知部(c)地温観測状況(a)単管井戸(d)温度測定器

写真-1 単管井戸における地温観測
自記水位計（ﾌﾛｰﾄ式）

内管地盤沈下
自記記録計

φ24mm

温度測定精度分解能
水深測定精度観測器具
外管ｾﾝﾄﾗｲｻﾞｰ

高精度
ﾃﾞｼﾞﾀﾙ
温度計(b)井戸管頭部(a)二重管井戸(c)地温観測状況

写真-2 二重管井戸における地温観測
72

温度センサー（サーミスタ）(d)温度測定器±0.1 ℃
0.01 ℃
1cm 程度

4. 大阪平野における地温観測結果

温 度 （℃）
15(1) 大阪平野の地温状態

17192123250

2011 年春，連名者谷口を研究代表とする共同研究グループにより，大阪平野の 30 地点で地温鉛直分布の一斉調査が行われた。
この一斉調査は 1998
年と 2003 年の一斉調査を引き継いだもので，同じ井戸での観測を行って地5013

下温暖化の進行状態を把握することが目的である。2011 年一斉調査の地温
深 度 （G.L.-m）

観測結果をまとめて図-5 に示す。G.L.-50m または-100m 地点の地温が低い
観測井 13，31 および 32 を除いて，深層(＞ 50 m) における地温は，深度100に対して，2 本の平行線（地温勾配：0.033 ℃/m，間隔：3.5 ℃）の間に収まる。一方，浅層(＜ 50 m) では，深層の地温勾配から外れ，高温側に湾曲する。これは明らかに，図-1 で示した地下温暖化の特徴を示している。

0.033℃/ｍ

3132

150

3 回の一斉調査(一部の観測井では 1998 年は未調査) で得られた地温分

0.033℃/ｍ

布の経年変化を図-6 と図-7 に示す。観測地点それぞれで土地利用や地表面被覆形態，さらには地層状態に違いがあるにも係わらず地温は上昇して200いる。
このような地温上昇は 34.泉南を除く他の観測井全てで認められたこ

3.5 ℃とから，大阪平野において都市化の進行している地域では，地下温暖化が図-5 2011 年大阪平野一斉調査による着実に進行していると考えられる。

地温鉛直分布（観測値）7) 一部加筆

なお，34.泉南の地温鉛直分布は，先述した一般的な特徴（図-2 参照）とは異なる特異な S 字状の分布を示していることから，地温データの扱いには注意する必要がある。

19230

19 1505050501501000

23 1705050100150

温 度 （℃）
1921100250

18.馬場町250

19.港(Ⅱ)B

10015014.庭窪 1-3200100

19

19150

10015022.堺 5-3100

9.十三

1921

2325

100

15020017.鴻池 22123

温 度 （℃）
27 17025

1950

150

100

15020024.堺 2-3

図-6 地温鉛直分布の繰り返し観測結果(1)7350

温 度 （℃）
1910025019

温 度 （℃）

19 17016.南郷200

23 170

深 度 G.L.-（m）

深 度 G.L.-（m）

深 度 G.L.-（m）17175015.新森小路

217.吹田

21 150

温 度 （℃）

21 150200

200

250

10050100

21.加美東

50

50

100

200

13.交野

温 度 （℃）

19 170171719150

20.長瀬
温 度 （℃）
15100

温 度 （℃）

21 17050

10050

6.豊中

50

50

0： 1998年 観測値：200
2003年 観測値
： 2011年 観測値

200

25 150150

凡 例19

温 度 （℃）

23

深 度 G.L.-（m）21

深 度 （G.L.-m）

深 度 G.L.-（m）

19

12.点野

200

温 度 （℃）
17100

11.鳥飼西200

1750150

150

10.柴島

19 15017

深 度 G.L.-（m）

200

17

深 度 G.L.-（m）

100

15

深 度 G.L.-（m）

深 度 G.L.-（m）

深 度 G.L.-（m）0100

温 度 （℃）

温 度 （℃）

温 度 （℃）

215.曽根19

深 度 G.L.-（m）17100

4.口酒井第 150

温 度 （℃）

温 度 （℃）

19 17017

深 度 G.L.-（m）

15

100

50

19 150

深 度 G.L.-（m）

温 度 （℃）

3.野間

50

17

深 度 G.L.-（m）

100

2.石橋

50

19 15017

温 度 （℃）

温 度 （℃）

温 度 （℃）

19 150

深 度 G.L.-（m）50

17

深 度 G.L.-（m）100

19 15017

深 度 G.L.-（m）50

温 度 （℃）

温 度 （℃）

21 150

深 度 G.L.-（m）

1.荒牧

19

深 度 G.L.-（m）50

深 度 G.L.-（m）

深 度 （G.L.-m）

0100

温 度 （℃）

19 17017

深 度 G.L.-（m）

温 度 （℃）
15250

25.堺 A-3

21

23

50

100

150

25 17
023

21 15
0

19150

100150

温 度 （℃）

19 15017100150

19 17017150

凡 例192123

25

50

50

100

温 度 （℃）

温 度 （℃）

19 15
050

50

50

100

17

深 度 G.L.-（m）50

21

深 度 G.L.-（m）019

温 度 （℃）

温 度 （℃）

温 度 （℃）
25 170深 度 G.L.-（m）23

深 度 G.L.-（m）21

深 度 G.L.-（m）

深 度 G.L.-（m）19

深 度 G.L.-（m）

温 度 （℃）
17100

150

100150： 1998年 観測値
200200200
200： 2003年 観測値

200200200： 2011年 観測値
250

26.堺 B-3250

28.岸和田 3250

29.岸和田 2250

30.岸和田 1250

31.貝塚-2250

32.泉佐野250

34.泉南

図-7 地温鉛直分布の繰り返し観測結果(2)300(2) 地温の低い観測地点の考察
図-5 において，観測井 13，31，32 の地温が他の観測井よりも相対的に低いことについて，一つの考察を試みる。地下温暖化前後の地温鉛直分布の形状は，地下水流動の有無で，図-8 に示す 4 つの模式図に大別される。すなわち，(a)：一定の地温勾配で表される地下温暖化前の定常的な地温分布(b)：(a)の状態から地中鉛直方向に地下水涵養・流入(Recharge)や地下水流出(Discharge)といった地下水流動(groundwater flow)の影響を受ける場合の地温分布（地下温暖化前の定常的な地温分布と言ってもよい）(c)：(a)の状態から地表温暖化(surface warming)が生じた場合の地温分布(d)：(b)の状態から(c)の地表温暖化が生じた場合の地温分布
標高
O.P.(m)13.交野
半径100m
田畑

観測井

神社

緑地

田畑

(a)定常状態小学校

0グラウンド

50m

13.交野

31.貝塚-2半径100m河川

観測井(b)g.f.の影響

河川

水田
河川

050m

31.貝塚-2

32.泉佐野大阪湾(c)s.w.の影響

更地
佐野漁港

観測井
半径100m

更地

野球グランド阪神高速道路湾岸線(d)s.w.と g.f. 両方の影響
0

図-8 地下温度に及ぼす地表温暖化(s.w.)と地下水流動(g.f.)の影響 9)

50m

更地

32.泉佐野図-9 観測井 13，31，32 周辺の空中写真と地層断面図
74である。図-9 は，地温の低い 13.交野，31.貝塚-2 および 32.泉佐B
A

野の井戸周辺の空中写真と地層断面図を示したものである。13.
交野では，観測井の周囲に田畑や緑地，小学校敷地周りには水路が存在するうえ，砂質土層が主体の地層構成となっている。31.
貝塚-2 では，2 つの中小河川に挟まれ，砂質土層優勢の地層構成Cにある。32.泉佐野では，大阪湾の佐野漁港に近接し，観測井の周囲はグランドや更地の状態であるが，層厚 10m 程度の粘性土層FＦが水平に堆積している。前二者は，雨水浸透しやすい地表面被覆

Dと水平に堆積した厚い粘性土層が存在しない地層構成から，地下水 涵養が生 じてい る可能性 もあり ， 地下温 暖化の 弱い(d)のRecharge タイプに近い。一方で後一者は，水平で厚い粘性土層温度（℃）の凡例の存在から地下水涵養の可能性は低く，それよりも地温勾配が小17.0以上 17.5以上 18.0以上 18.5以上 19.0以上
17.0未満 17.5未満 18.0未満 18.5未満 19.0未満

さい状況下での(c)タイプに近いと推測される。

E

上述したような地形，地質に基づく全観測井における地温鉛直

Observation well

分布のタイプ区分については，今後の研究課題としたい。
13.交野(1) G.L.-40m 地温の分布と19その上昇速度
地下温暖化を評価する当面の指標として，G.L.-40m 地温とその上昇速度を採用する 7)。G.L.-40 m 地温の採用は，2003 年大阪平野一斉調査の結温 度 （℃）

地下温暖化の実態

20

さは GL-約 80 m

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

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

0.06
0.04

180.02170.00

16-0.02

15-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)： GL-40m地温

大阪都心部

10.柴島

温度上昇速度（℃/年）

5. 大阪平野における

図-10 G.L.-40m 地温の平面分布図 7) ｶﾗｰ図化まであり，その中

19 ,20 ,21 ,22 ,24 ,25 ,26 ,28 ,29 ,30 ,31 ,32 ,34

観測井No.

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

間深さをとったものである。

表-1 G.L.-40m 地温とその上昇速度 7) 一部加筆

図-10 は，各々観測井の G.L.-40 m 地温を大阪平野の白地図へプロットした結果である。プロットしたデータには，図-3 に記載のない郊外

測候所
No.の観測井（大阪市大正区の 1 本，柏原市の 1 本）で観測したデータも含まれている。また，気象庁測候所と観測井の位置関係から，観測井を測候所毎に目安としてのグループ化も行っている。図によると，地温は郊外ほど低く，反対に沿岸部や都心部では高い傾向を示すことがA
（豊中）

B（枚方）

分かる。
G.L.-40 m地温とその上昇速度についてまとめた結果を表-1 と図-11に示す。表では，図-3 および図-10 で示した気象庁測候所 A～F 毎にG.L.-40 m 地温の上昇速度の平均値を参考値として付記した。図表から，
C以下のことが分かる。

（大阪管区）① G.L.-40 m 地温とその上昇速度は，観測井毎に異なる。
② G.L.-40 m 地温の最も高い観測井は，大阪都心部の N3 西天満の大阪都心部

19.27 ℃で，次いで同じ都心部の N6 道修町である。
③ G.L.-40 m 地温の上昇速度は，特異な地温鉛直分布を示す 34.泉南を除いて，0.001～0.054 ℃/年の範囲にある。

F（八尾）④ G.L.-40 m 地温の上昇速度が最も大きい観測井は，大阪都心部のN6 道修町の 0.054 ℃/年で，次いで同じ都心部の 18.馬場町と N3D（堺）

西天満である。これら井戸の地温上昇速度は 0.050 ℃/年以上の速さで，大阪管区気象台の過去 130 年間の気温平均上昇速度0.020 ℃/年や過去 63 年間の 0.028 ℃/年の約 2 倍程度大きい。
⑤ 大阪管区気象台の気温平均上昇速度 0.020 ℃/年以上となる観測井の数は，計 23 本で，全 34 本の約 2/3 にも及ぶ。
⑥ 地温上昇速度が最も小さい観測井は，大阪都心部の N5 西船場で，次いで 4.口酒井，N1 福島である。いずれも 0.01 ℃/年未満である。

75E（熊取）観測GL-40m
観測
井戸
地温
所名
No.
（℃）
荒牧
117.30石橋
218.48*野間
317.25
口酒井第1417.15
曽根
517.03豊中
617.26
交野
1315.99吹田
718.77
十三917.59柴島
1017.89
鳥飼西
1116.67
点野
12
***
庭窪1-3
1417.42新森小路
1518.85
南郷
1617.44鴻池21717.93
18
N6
N3
N2
N1NK
N5
N4

192021222425
26
28
2930
313234

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

G.L.-40m地温の上昇速度 (℃/年）

2003～2011*
0.017

平均値

***

0.035
0.005
0.035
0.029
0.021
0.038
0.018
0.047
0.020

0.024

0.021

***

0.037
0.041
0.033
0.029

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.0300.032（都心部を除く）

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年のデータ
泉南：特異な地温鉛直分布であり，地温データの扱いには要注意

G.L.-40 m 地温の上昇速度が気温平均上昇速度を上回っている井戸は，一つの解釈と14.庭窪1-3

して，都市化の進行に伴う地表面温度の上昇速度が気温の上昇速度よりも大きい可能性が考えられる。

17.鴻池2中之島B

18.馬場町(Ⅱ)(2) 大阪平野の地下温暖化

19.港(Ⅱ)B

大阪平野において，図-12 に示す観測井で地下温暖化の深さを分析した。そのうち，

24.堺2-3

今回分析した結果を図-13 に示す。図より，各井戸の破線を地下温暖化が始まる前の定常的な地温分布を表すものとすると，地下温暖化の深さは井戸毎に異なるが，概ねG.L.-70～-90 m 程度と推測され，既に報告した大阪都心部の 18.馬場町(大阪管区気象台の近傍)の深さ 80 m

7)28.岸和田3と同程度である。このことは，大阪都心部以外の地点でも，（18.馬場町と中之島 B の分析は，文献 7 にて報告済み）ローカルな地下温暖化が早くから起こっていた可能性がある。

図-12 地下温暖化を分析した
観測井の位置図

6. まとめ大阪平野および大阪都心部において，地下の温暖化は着実温 度 （℃）に進行中である。
地下温暖化の進行性の指標として，G.L.-40m

15

17井は，大阪都心部の N6 道修町の 0.054 ℃/年で，次いで同じ都心部の 18.馬場町である。これら観測井は大阪管区気象台の気温平均上昇速度の約 2 倍程度の大きさである。また，観測

深 度 （G.L.-m）

観測井は，大阪都心部の N3 西天満の 19.27 ℃で，次いで同G.L.-40m 地温の上昇速度でみると，速度の最も大きい観測

2123250

地温とその上昇速度で評価すると，G.L.-40m 地温の最も高い
じ都心部の N6 道修町である。

19地下温暖化の深さ
概ね70～90m

50

24.堺2-3
17.鴻池2

100
14.庭窪1-3

19.港(Ⅱ)B

28.岸和田3150

井 34 本中 23 本において，G.L.-40m 地温の上昇速度は，過去130 年間の大阪管区気象台の気温平均上昇速度よりも大きい。
地下温暖化の深さから評価すると，大阪都心部以外でも深200いところで概ね G.L.-90 m 程度と推測され，大阪都心部以外

図-13 大阪平野における地下温暖化の状況

でもローカルな地下温暖化が早くから起こっていた可能性がある。
謝辞：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., and Jago-on, K. : Combined effects of urbanization and global warming on subsurface temperature infour 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.
7）有本弘孝，北岡豪一，谷口真人，濱元栄起：大阪都心部における地下温暖化の実態，Kansai Geo-Symposium 2013 地下
水地盤環境・防災・計測技術に関するシンポジウム論文集，地盤工学会関西支部，pp.53-58，2013.
8）濱元栄起，有本弘孝，北岡豪一，谷口真人：大阪都心部における地下温暖化履歴の推定，Kansai Geo-Symposium 2013 地
下水地盤環境・防災・計測技術に関するシンポジウム論文集，地盤工学会関西支部，pp.59-64, 2013.
9）Taniguchi, M. Shimada, J. Tanaka, T. Kayane, I. Sakura, Y. Shimano, Y. Dapaah-Siakwan, S. and Kawashima, S. : Disturbances oftemperature-depth profiles due to surface climate change and subsurface water flow, 1. An effect of linear increase in surface
temperature caused by global warming and urbanization in the Tokyo metropolitan area, Japan, Water Resources Research, Vol.35,No.5, pp.1507-1517, 1999.

76