(大学院) 圏間相互作用特論 II
このウェブページは http://macroscope.world.coocan.jp/ja/edu/ris/kenkan/ です。
第1回は、研究科授業一覧どおりの日時 (10月6日 2限)、場所 (1号館 1306号室) でおこなう予定です。
ネットワーク接続可能なノートパソコンを用意してください。
実習の回は、研究室 (3号館 3410) でおこなう予定です。
各回の内容
- 10月 6日 (1306) 序論、進めかたの相談
- (10月 13日 休日) (3410) パソコンを端末とするためのソフトウェアのインストール
- 10月 20日 (1306) 暖かさの指数について文献を読む
- 10月 27日 (3410) [暖かさの指数を計算してみる] ← このリンクさきページで説明します。
- 11月 4日 (火、月曜授業日) (3410)
- 11月 10日 同じ時間に生物多様性のミニシンポジウム (8号館 8201) があるので、
この授業は 2限にはやらないことにします。
9日に材料を用意したので、各自作業して、17日に報告してください。
- 改訂されたプログラムで乾湿指数を計算しなおす。
- [暖かさの指数を計算してみる] への補足。
乾湿指数を計算する際に、降水量を月ごとに頭打ちにするところの、
適切と思われるやりかたがわかりました。(下にリンクするブログ記事で説明しました。)
- それにしたがってプログラム hia_to_csv.py をつくったので、実行してみてください。
結果として hia.csv ができます。
(cd crucl ; cp /home/masudako/ex/crucl/hia_to_csv.py ; python hia_to_csv.py )
- それを地図に表示するプログラムは 同じ crucl ディレクトリにある map_hia_eastasia.py ですが、
これは前回の map_hix_eastasia.py の入出力ファイル名と表示する文字列をちょっと変えただけです。
- 横軸に乾湿指数、縦軸に暖かさの指数をとった散布図を、植生分類ごとにつくる。
- [気候と植生の関係を見るための可視化 (1)] のウェブページをつくり、
前回作業した年降水量と年平均気温の散布図をふくめました。
- そこに、横軸に乾湿指数、縦軸に暖かさの指数をとった散布図を
植生分類ごとにつくるプログラム scat_wi_hia_by_veg.py を 用意しました。
上記の hia.csv をつくったあと、これをやってみてください。
(cd clim_veg_w1 ; cp /home/masudako/ex/clim_veg_w1/scat_wi_hia_by_veg.py ; python scat_wi_hia_by_veg.py )
- 選択されたもの以外の植生の地点についても薄い色で散布図に点をうつと、
選択された植生がそのうちでどこにあるかがわかりやすくなりますが、
計算時間がかかります。
「backgr = True」としたばあいにかぎってこれをやることにしました。
- 軸の向きのとりかたはいろいろ考えられるので、
上下は bottomwarm に True, False のいずれをあたえるか、
左右は righthumid に True, False のいずれをあたえるかで
切りかえられるようにしました。
- cf. 今西・吉良 (1953) の 第84図: 湿潤→, 温暖↓
- cf. Budyko [ブディコ] のダイアグラム: 乾燥→, エネルギー大↑
- ぜひやってほしいというわけではありませんが、
もし、ちがうバージョンの乾湿指数をつかった作図をしたければ、
プログラム中で読むファイル名
(「dfwi = pd.read_csv」でつかわれる filenamein にあたえる文字列) の hia を hix や hiik にかえてください。
そのときは、横軸につける文字列 (ax.set_xlabel で指定) や、
図を保存する画像ファイル名 (imagefilename で指定) もかえてください。
- Ramankutty の植生分類のうち「Mixed forest」は、いろいろな混合をふくんでいるので、
散布図のなかの分布範囲がしぼられないのは当然でしょう。
- 11月 17日 (3410) [暖かさの指数などと植生を対比してみる] (つづき)
- 植生の分布をみる
- Ramankutty の植生分類の地理的分布をみるため、
/home/masudako/ex/clim_veg_w1 のディレクトリに、
(地図上の表示と基本的に同じなのですが) 横軸に経度、縦軸に緯度をとった
散布図の形で表示するプログラム scat_lon_lat_by_veg.py を追加したので、
動かして作図してみてください。
- 植生分布と気候要因 (吉良による指数) との関係をみる
- 植生分類ごとに、地理的分布の図と、「暖かさの指数・乾湿指数」の
2軸での散布図を見くらべて、分布の特徴を見ていきます。
- さらに条件をしぼった点の分布を作図しながら考えようと思います。
- 寒さの指数についても見てみようと思います。
ただし、寒さについては理屈の面からは最低気温のほうがよさそうなので、
吉良の寒さの指数にはそれほど時間をかけないことにします。
- 標高との関連もみる
- 気候は、緯度・経度のほかに標高によっても変わります。
- 気温については、大まかに対流圏の気温減率を反映して、標高の高いところほど地上気温が低く、
そのことは CRU のデータセット作成でも反映されています。
- 降水量と地形との関係は、風上斜面で降水が多いなどの傾向がありますが、
方位や数量は地域ごとにまちまちで、
CRU データセットでは地形依存性は計算にはいっておらず、
観測点のない標高帯の降水量はあまり正確ではないでしょう。
- ともかく、植生の地理的分布をみるうえで、
地面の標高を考慮したほうがよいことはあるでしょう。
- 標高のデータセットとして、
アメリカ合衆国 NOAA National Geophysical Data Center (現 National Centers for Environmental Information) で編集された
GLOBE という 緯度経度 0.5分 (約 1 km) の格子データを
緯度経度 0.5 度で、ひとまず平均したものをつくりました。
- データセットの紹介は、これから書く予定です。
([別ページ]でつかった「ETOPO1」は海底の水深がはいっていましたが、
GLOBE は海の地点には値がありません。)
- 平均後の標高データは /raid/globe_dem/hdegav/elev_hdeg.r4l 、
それを読んで地図上に表示するプログラム例は /raid/globe_dem/py/map_elevhdeg.py にあります。
- [2025-11-18] 11月17日の授業にまにあわなかった、散布図を標高によって色わけするプログラムを、
/home/masudako/ex/clim_veg_w1/ に追加しました。
- scat_wi_hia_by_topo_and_veg.py ... 植生型ごとに、縦軸に暖かさの指数、横軸に乾湿指数をとって散布図をつくる。点を標高の階級別に色わけをする。
- scat_lon_lat_by_topo_and_veg.py ... 植生型ごとに、横軸に経度、縦軸に緯度をとって点をプロットする。点を標高の階級別に色わけをする。
- いずれも、標高の階級の境界を thershs、色を colors というリストで指定するようにしました。
色のところに 'none' を指定すると、その階級については点を打ちません。
- 標高の利用としては、(経度を省略して) 緯度と標高の2次元で見ることもやってみようと思います。
- 11月 24日 (1306) 文献を読む。
- 新井 正, 2005: 気候区分における水文要素 -- その歴史的展開。
地球環境研究 (立正大学), 7, 1-7.
(『地球環境研究』の論文は、地球環境学部ウェブサイトの「刊行物」のところにあるが、この存続は不確か。
このたび、永井先生 と 図書館職員のかたがたの努力により、
立正大学リポジトリ https://rissho.repo.nii.ac.jp/ に置かれ、こちらが存続する見こみ。)
- 近藤 純正、徐 健青, 1997:
ポテンシャル蒸発量の定義と気候湿潤度。
天気 (日本気象学会), 44: 875-883.
(『天気』の論文は、日本気象学会のサイトの『天気』のページ https://www.metsoc.jp/tenki/ から検索で PDF 取得可能。)
- (そのあとの回で、Budyko (または 内嶋)、酒井、Woodward)
- ...
Linux 上のエディタについて [2025-10-27, 2025-11-03]
気象のセミナーでは、エディタとしては emacs をつかっている。
研究室のマシン hamana と niigata では、
X Window の機能をつかわない emacs がインストールしてある。
他のサイトで X Window の機能をつかわないで実行するためには
オプション「-nw」が必要になることがある。
emacs よりも直観的に操作できるエディタがいくつかある。
多人数が同時にログインするときにはすすめられないが、
少人数のときはさしつかえなさそうだ。
これまでためしたうちでは featherpad がよさそうなので、
fpad という別名でもつかえるようにした。
他のものもふくめた情報はつぎの記事参照。
Python の matplotlib での色の指定 [2025-11-03]
色は文字列で指定する ('red' のように)。
文字列は、英語の色の名まえがいろいろ予約されているほか、
色を RGB (赤、緑、青) に分解してそれぞれ16進数2けたで表現した6文字の あたまに「#」をつけた形で指定することもできる
('#FF0000' のように)。
今学期の構想 (10月6日現在のもの。今後書きかえます。)
2025年度は、増田が単独で担当します。
2023年度、岩崎先生と半分ずつを分担したときと同様な内容をふくらませて、
そのとき見送った世界の植生と気候との関係をディジタルデータで見る実習をしたいと思っています。
構想は、つぎの文章に書いたことからはじまります。
この構想は、1学期の [熱収支・水収支特論]
とつながっています。
その科目を履修ずみである必要はありませんが、必要に応じてそちらの内容 (さしあたりつぎのもの) を自習してください。
これからやっていくことは、だいたい次のように予定しています。
- 植生の分布を制約する気候条件について述べたいくつかの文献を読む。
- ケッペン (Köppen)
- 吉良 竜夫 の「暖かさの指数」=「温量指数」 (および、「寒さの指数」、「乾湿指数」)
- 吉良 (1949) 日本の森林帯。吉良 (1971) 『生態学からみた日本の自然』 [読書メモ] に収録。
- 川喜田 二郎 (1999) 『環境と人間と文明と』 [読書メモ] ... 1940年代の着想を回顧している。
- 今西 錦司、吉良 竜夫 (1953) 生物地理。福井 英一郎 編 『自然地理 II』 [読書メモ] 第5章。
乾湿指数の説明をふくむ。[スキャンした画像 (臨時リンク)]
- (以下、文献検討中)
- Holdridge の biotemperature
- Budyko (ブディコ)、内嶋 善兵衛 ... 地表面の正味放射と降水量にもとづく「放射乾燥度」
- Thornthwaite (ソーンスウェイト) の可能蒸発散量、それをつかった Lieth (リート) と Box の光合成生産量評価
- 酒井 昭 ... 凍結について
- Woodward (内嶋訳 1993) 『植生分布と気候変化』
- 生態学用語、とくにバイオーム (biome、生物群系) について学ぶ。
- 地表面のエネルギー収支 (熱収支) と水収支について
- 植生分類と気候要素の両方の地理的分布のデータをつかって関連を解析する。
- たとえば、Budyko のダイアグラムは、
縦軸にエネルギーの指標 (正味放射)、
横軸にそれと水分の指標 (降水量) との比をとった 2次元平面に、
おそらく各地点の植生コードをプロットして、それぞれの植生コードのひろがりをまとめたものである。
つぎのページ (アクセス制限あり) の下の図を参照。
- 同様なことを、現代 (1961-2020年 のうちのいずれかの30年ぐらいの期間) のデータでやってみたい。
格子点 (あるいは、マス目) 型のデータをつかえば、計算機処理はできる。
- ためしにやった例: Masuda 2000 都立大学 地理 英文紀要 [text]
- 緯度経度 1度格子。
- 植生データは Matthews (1984), 気候モデルのために専門地図帳類から編集されたものだが、品質が不ぞろい。
- 気候データは New et al. (1999), East Anglia 大学 CRU で編集された気温と降水量の平年値。
- 散布図の配置は Budyko 型だが、エネルギーの指標は Holdridge による average biotemperature である。
- しかし、つぎのような問題点がさけられない。
- 格子点の値は、格子点で観測・観察されたものではなく、空間内挿、あるいは衛星で観測可能な量からの推定がはいっている。
- 別々につくられたデータの同じマス目の値は、マス目の内にある同じところを見た結果とはかぎらない。
- とくに、同じマス目のなかで標高の不均一性が大きいばあい、見ているものがちがう可能性が大きくなる。
- 標高の情報や観測点密度の情報も参照して、比較的信頼できる部分について論じることにしたい。
- 対象は世界と日本があり、データのくわしさがことなる。
今年度は世界をやろうと考えているが、受講生と相談して決める。
- 作業は、研究室の Linux マシン「hamana」に、パソコンを端末として login し、
Python のプログラムを動かす形でおこなう。
この形の授業を経験していない受講生には、準備作業の時間を別にとる。
参考文献
- 松山 洋、増田 耕一 編, 2021: 『大気と水の循環』。朝倉書店, 136 pp. ISBN 978-4-254-16076-5.
[読書メモ]
... 第1～6章、8章を増田が書いた。
過去の教材ページ
(リンクさきがなくなっているばあいや、別のものにおきかわってしまうばあいもあります。)
関連の教材ページ
(更新) 2025-11-18
増田 耕一 (MASUDA Kooiti)