自然教育園報告（Rept. Inst. Nat. Stu. ）第48号：1−15,2017.

自然教育園の林内気温の特徴①

近藤純正＊・菅原広史＊＊・内藤玄一＊＊
Air temperature in National Park for Nature Study
Junsei Kondo＊, Hirofumi Sugawara＊＊, Genʼichi Naito＊＊

はじめに

森林内の気温については，正しく理解されていない。それには 2 つの理由がある。理由の 1 として気温観測に誤差があること。現在，多方面で使われている通風式気温計や自然通風式（非通風式）気温計は放射影響によって，晴天の日中は高温に，夜間は低温に観測される。
風が弱い晴天時，自然通風式（非通風式）気温計では最大 5℃以上の放射影響の誤差が，ファンモータを使った通風式の比較的高精度といわれている気温計でも 0.3 〜 0.5℃の放射影響による誤差がある（近藤，2014，の表 92.2；近藤，2015a）。
放射影響の誤差により，日当たりのよい林外と日陰の多い林内で観測すると，林外の気温が高く表示され，相対的に林内はかなり低温だという常識が生まれる。
理由の 2 は，木陰の林内では体感温度が低いことによって先入観ができる。また，日中の森林は蒸散によって大気を冷却するという，不正確な常識が生まれている。通常の晴天日中の条件では，葉面は気温より高温となり大気を加熱している。林内の低温は，おもに日陰の作用によるものである。
森林公園内が都市ビル街より高温となる例がある。東京都心部の気温は大手町露場で観測されていたが，2014 年 12 月 2 日から北の丸公園内の露場で観測されるようになった。森林公園内の北の丸露場では，3 月〜 4 月の晴天日中の最高気温はビル街の大手町より平均 1℃ほどの高温である（近藤・菅原・内藤・萩原，2015a）。
北の丸露場は周辺の風通しが悪く，晴天の日中は「日だまり効果」で都心ビル街よりも高温に，夜間は放射冷却が強く低温になる。
前記の不正確な常識を正し，林内の気温について理解を深めることは次の諸問題に生かすことができる。
（1）気象観測露場の環境維持（2）快適な市民公園の整備（3）森林の都市気候に及ぼす影響（4）森林の生態系に及ぼす影響＊東北大学，Tohoku University＊＊防衛大学校，National Defense Academy

自然教育園報告─2─第48号，2017

その目的のために，都心部の森林公園（新宿御苑，明治神宮，代々木公園，北の丸公園），つくば市内や平塚市内の森林公園で気温の観測を行なってきた。また，東京白金台の国立科学博物館附属自然教育園における観測として，2014 年秋から 2015 年 2 月 10 日までの観測（近藤・菅原・内藤・萩原，2015a）と，2015 年 3 月 12 日から 10 月 13 日までの観測（近藤・菅原・内藤・萩原，2016）がある。
これまでの研究結果：
図 1 は，それらのまとめの一部である。縦軸は，森林に隣接する広い芝地広場の気温を基準とした気温差であり，林内気温が芝地広場より高温のときをプラスとする。横軸の木漏れ日率は，林内の直図1

木漏れ日率と気温差の関係，快晴または日射の強い薄曇り日．

（近藤・内藤，2015b，の図 115.3 に緑の楕円と円を加筆）．
黒丸印：見通し良好林，四角印：見通し不良林
上図：4 〜 5 月，赤と黒の塗つぶし印は前日が雨
下図：6 〜 9 月，赤塗つぶし印は大雨後の晴天日，緑と黒塗つぶし印は雨後の晴天日
大きな緑楕円の範囲は，年間の北の丸露場の気温差に相当する．また，大きな緑円の範囲は自然教育園の密な林内の気温差に相当する．

近藤ほか ： 林内気温の特徴

─3─

射光が占める面積の割合である。木漏れ日率は，観測点の周囲約 20 ｍ範囲の目測値であり，雲量の観測に似ている。最下段の横軸は日射量の比「林内日射量÷林外日射量」，同じ横軸の木漏れ日率との関係が分かるように表してある（近藤・内藤，2015a）。
図中のプロットは丸印と四角印で大別し，それぞれ見通し良好・不良を表す。目の高さ（≒気温観測の高さ）で林内を水平にみたときの見通しの良し悪しによって「見通し良好」と「見通し不良」とする。方位の 50％以上の範囲の見通しが，おおむね 30 〜 50 ｍ以下の場合を「見通し不良」，おおよそ 50 ｍ以上まで見える場合を「見通し良好」と定義する。詳細は近藤（2015d）の「林床の木漏れ日率と林内の見通し（詳細）」を参照のこと。
図 1 によれば，気温差は見通し良好林では木漏れ日率＞ 20％の範囲でゼロに近い。しかし，見通し不良林ではプラスとなり季節によって大きく異なる。4 〜 5 月には気温差は＋ 0.8℃前後で大きいのに対し，高温期の 6 〜 9 月には＋ 0.4℃前後で小さい。
いっぽう，木漏れ日率＜ 20％の自然林に近い密な林内では，気温差は 4 〜 5 月には− 0.5℃前後，6 〜 9 月には− 1℃前後である。
ほかに，この図で注目すべきことは，特に密な林内（木漏れ日率＜ 20％）では，降雨の翌日の晴天の日中の気温上昇が小さく気温差は− 2℃前後，通常より 1℃ほど低温である。これは，雨後の林床下の土壌水分が増えて熱容量が大きいために林床面の温度上昇が抑制されることによるものである。
大雨が続いた後に晴天が続いた場合，気温変化が元に復するまでの日数は 4 〜 5 日かかる（近藤・菅原・萩原・内藤，2016 の図 15；近藤・菅原・萩原・内藤，2015）。
図 2 は自然教育園の観測塔で観測された晴天日の気温鉛直分布の例である。日中は樹冠層の上部が日射を受けて高温の極大となり，そこから上層の大気側と下層の林床に向かって顕熱が輸送されていることがわかる。夜間は樹冠層が放射冷却によって低温となり，日中とは逆に上層の大気から樹冠層に向かう下向きの顕熱輸送が生じる。
林内の最下層では，昼夜ともに下ほど気温が低いことから，林床下の土壌層へは顕熱（地中伝導熱）が昼夜ともに入ることがわかる。日中の林床では日射量が林外に比べて 10％以下であり，このわず

図2

自然教育園の観測塔における晴天日の気温の鉛直分布．

（近藤・菅原・内藤・萩原，2016，の図 13 の一部分）．

─4─

自然教育園報告第48号，2017

かの日射量と顕熱が地中へ入り，林床下の地中水分を蒸発させる潜熱のエネルギーとして使われる。
本研究の目的：
本研究の目的は，曇天・雨天も含む年間の条件について林内気温の特徴を調べることである。具体的には，見通し不良林の「開空間」と「密な林内」を対象とする。すなわち林内の「開空間」を代表する北の丸露場の気温と，自然林に近い「密な林内」を代表する観測塔付近の林床上の高度 1.5 ｍの気温について調べる。
図 1 に描かれた大きな緑楕円は，北の丸露場に相当し太陽高度が低い条件も含む範囲を示している。
左方に描いた大きな緑円は自然教育園内の観測塔付近の密な林内に相当する。
気温を比較する際の基準として，自然教育園の樹冠上（高度 19 ｍ，平均の樹高上端から 5 ｍ）の気温を用いる。
この高度 19 ｍの気温は，東京都心部の芝地広場（新宿御苑の芝地広場，明治神宮境内北部の宝物殿前の芝地広場，代々木公園中央広場，北の丸公園の池の隣の芝地広場）の気温とほぼ等しいことから，都心部を代表する気温とみなされる（近藤・内藤，2015c；近藤・菅原・萩原・内藤，2015；近藤・菅原・内藤・萩原，2016）。

観測観測塔
東京白金台（JR 山手線目黒駅の東 500 ｍ）の国立科学博物館附属自然教育園に高さ 20 ｍの観測塔がある。観測塔の位置は図 3 に示してあり，図中の「開空間」，「開空間北」，「塔北 20 ｍ」，「塔南 50 ｍ」
は前報で行った観測地点である。
気温の観測は高度 19.0 ｍ，15.6 ｍ，および 1.5 ｍの 3 高度で行った。高度 1.5 ｍの気温計は，塔の下の観測小屋の影響がないように塔の南 20 ｍの地点に設置した。
観測塔周辺の代表的な樹木は，「塔北 20」の位置では高木スダジイ，亜高木ヤブツバキ，トウネズ
ミモチであり，木漏れ日率は年間を通じて小さい。「塔南 50 ｍ」の付近はコナラ林であり，木漏れ日率の季節変化が大きい。観測塔の周辺は，密な自然林に近いとみなされる。
樹林の着葉が十分になる 5 月上旬以後の季節には，正午前後の林床上の木漏れ日率は 10％以下，すなわち，林床に届く日射量は林外日射量の 8％以下となる（近藤・内藤，2015a）。
観測塔周辺の樹木の平均的な高さは 14 ｍである（直近で 16 ｍの樹木もある）。
気温計
気温の観測では高精度の強制通風式気温計を用いた。この気温計の総合的誤差（通風筒に及ぼす放射影響を含む）は 0.03℃程度である（近藤，2015b）。
センサーは白金抵抗体の Pt1000 オーム，受感部の直径は 2.3 ｍｍを用いている。気温のサンプリ
ングは 10 分間隔である。気温観測用の通風筒のファンモータの電源は 12 ボルトであり，AC100 ボ
ルトの電源に AC アダプターを接続して 12 ボルトとした。ファンモータの電流は，当初は省電力型の 0.06A であったが，長寿命ファンの 0.24A の品に取り換えて使用した。

近藤ほか ： 林内気温の特徴

図3

─5─

自然教育園の観測塔の位置（星印）．

塔北 20 ｍ：観測塔の北 20 ｍの林内（ほとんど常緑樹）
塔南 50 ｍ：同上の南 50 ｍの林内（落葉樹と常緑樹）
開 空 間：旧管理棟跡地の空間（中央部近くに落葉樹）
開空間北 20 ｍ：同上の北 20 ｍの林内（落葉樹と常緑樹）

月平均気温の比較
全体的な傾向を知るために，月ごとの平均気温を表 1 と図 4 にまとめた。参考のために，表 1 には東京の郊外を表すとみなされる府中アメダス（東京農工大学の広い農場）とママ下湧水公園（国立市）の気温も含めてある（近藤・内藤，2017）。
気温変化の特徴を見やすくするために，高度 19 ｍの気温を基準とした気温差も示してある。なお，

自然教育園報告─6─
表1

図4第48号，2017月平均気温の比較（2016 年 5 月の 19 ｍ高度の気温は推定値）．

月平均気温の季節変化．下図は 19 ｍ高度の気温を基準とした気温差の季節変化である．

近藤ほか ： 林内気温の特徴

─7─

気温計のファンモータが故障した 2016 年 5 月の 19 ｍ高度の気温は推定値を示してある。
図 4 によれば，ビル街の大手町の気温（紫プラス印）は 4 〜 10 月には広域を代表する基準値（19ｍの気温）より 0.4 〜 0.6℃ほど高温であり，冬に向かってゼロに近づいている。日照時間が短くなり気温の上昇が抑えられるためと考えられる。
北の丸露場の気温（赤丸印）は，基準値とほとんど同じであるが，10 月からマイナスが大きくなるのは，露場周辺に日陰が多くなることによると考えられる。
密な林内の気温差（緑四角印）は，− 1℃前後で，冬に向かうにしたがってマイナスが大きくなる。
高度 15.6 ｍの気温差（緑破線）は，全期間とも− 0.2℃前後である。図 2 に示した気温鉛直分布の関係が月平均値についても成り立っている。この微小な気温差のマイナスは，大気から森林へ下向きの顕熱輸送が行われているという意味ではない。日中は樹冠から上の大気へ，夜間は逆方向の顕熱輸送が生じる。熱フラックスを詳しく観測した結果によれば 6 年半の平均値は，上向きをプラスとして顕熱輸送量＝＋ 2W/m2 となり，樹冠層から上向きの顕熱輸送がある（近藤・菅原，2016）。

気温日変化の月平均値
図 5 〜 9 は偶数月ごとに示した気温の日変化である。各図の下段は高度 19 ｍの気温を基準にして表した気温差である。各図の下段に注目して見ていこう。
林内の開空間に設置されている北の丸露場の気温差（丸印）について，4 月は正午前後には＋ 0.7℃ほど，深夜から日の出頃にかけて− 0.7℃ほどである。6 月と 8 月は 4 月の傾向とよく似ている。

図5気温の日変化の月平均値（4 月）．

自然教育園報告─8─

第48号，2017

図6

気温の日変化の月平均値（6 月）．

図7

気温の日変化の月平均値（8 月）．

近藤ほか ： 林内気温の特徴

図8

気温の日変化の月平均値（10 月）．

図9

気温の日変化の月平均値（12 月）．

─9─自然教育園報告─ 10 ─第48号，201710 月になると，1 日を通じて気温差は全体としてマイナス側にずれる。12 月には，露場周辺は直射光が少なくなるためか，自然教育園の密な林内（四角印）の気温に近づく。
密な林内（四角印）について，4 月の日中の気温差はゼロに近いが夜間は− 1℃ほどである。6 月と 8 月は 1 日を通して− 1℃前後である。これは落葉樹も完全に着葉した季節であり，林床に日射がわずかしか届かなくなることによると考えられる。大気が乾燥する 10 月には，樹冠上部層の放射冷却で冷えた空気が重力によって林床に下降してくるためか，17 〜 21 時ころの気温差がマイナスで最大となる。12 月になると，この傾向が顕著になり，17 〜 21 時ころ− 2℃前後となる。
高度 15.6 ｍの気温差（破線）は，正午前後にゼロまたは＋ 0.2℃程度のプラスとなるが，その他の時間帯は− 0.2 〜− 0.5℃の範囲内にある。

晴天日の気温日変化
日変化の特徴が顕著になる晴天日について図 10 〜 14 に示した。晴天日とは北の丸公園の科学技術館屋上における日照時間侒 9 時間とする。ただし 1 日の可能日照時間が少なくなる 12 月のみ日照時
間侒 8.8 時間とする。
前節と同じように，自然教育園の高度 19 ｍの気温を基準として気温差の日変化を各図の下段に示した。
北の丸露場の気温差（丸印）は，どの季節も正午頃には「日だまり効果」によって＋ 1℃ほど高温

図 10

晴天日の気温日変化（4 月 12，15，19，20，26，29，30 日の 7 日間平均）．

近藤ほか ： 林内気温の特徴

図 11

晴天日の気温日変化（6 月 2，3，10，11，18，26，27 日の 7 日間平均）．

図 12

晴天日の気温日変化（8 月 5，7，17，25，26，31 日の 6 日間平均）．

─ 11 ─

自然教育園報告─ 12 ─

図 13

図 14第48号，2017

晴天日の気温日変化（10 月 2，15，20，24，26 日の 5 日間平均）．

晴天日の気温日変化（12 月 9，11，17，28，30，31 日の 6 日間平均）．

近藤ほか ： 林内気温の特徴

─ 13 ─になり，夜間は放射冷却によって− 1℃前後の低温となる。10 月と 12 月は露場周辺は直射光が少なくなる関係か，夜間の自然教育園の密な林内（四角印）の気温差に近くなる。
自然教育園の密な林内（四角印）について，落葉樹がまだ完全に着葉していない 4 月には，正午ころ林床の日射量は少ないながらも風が弱いことで気温差はゼロまで上昇するが，夜間は− 1℃程度となる。6 月〜 8 月は 1 日を通して気温差は− 1℃前後で変化している。
10 〜 12 月になると，大気が乾燥し樹冠層の放射冷却が強く，その冷気が林床に降下し，16 時以後は− 1.5 〜− 2.4℃の低温となる。
注意 1（北の丸露場の気温計の誤差）
晴天日中は，気象庁の気温観測用の通風筒は放射影響によって 0.3℃ほど高温に記録される。この誤差を補正すると，正午前後の＋ 1℃は正しくは＋ 0.7℃となり，これまでの開空間で得られている
図 1 の大きい楕円で囲んだ範囲の結果と矛盾しない。なお，図 1 のプロットはすべて高精度の通風式気温計で観測された値である。
注意 2（基準とする気温の選定）
図 1 において，晴天日中の日だまり効果による昇温量の平均値は，季節（気温）によって少し異なり＋ 0.7℃（4 〜 5 月），＋ 0.4℃（6 〜 9 月）であった。本解析の図 10 〜 14 では，丸印で示したように，正午前後の気温差の季節（気温）による違いはほとんど認められない。
その理由の一つは，基準とする気温は図 1 では芝地広場の気温を用い，図 10 〜 14 では樹冠上 5 ｍ（林床上 19 ｍ高度）の気温を用いた。この違いによって気温差にわずかながら差が生じた可能性がある。
高温の夏期は，芝地に比べて森林の蒸発散が盛んで，樹冠上 5 ｍの気温が芝地の 1.5 ｍ高度の気温に比べて相対的に 0.3℃ほど低温であると思われる。その結果，今回の 19 ｍ高度の気温を基準とした気温差が夏期に＋ 0.3℃ほど高温になったと考えられる。
注意 3（夜間の北の丸露場「開空間」と自然教育園「密な林内」の気温差）
晴天日の図 10 〜 14 において，夜間 21 時〜 6 時の気温は北の丸露場「開空間」が「密な林内」に比べて平均的に 0.3℃ほど高温に表されている。北の丸露場では土盛りとなっており，夜間は 0.15℃ほど高温である（近藤，2015c）。この 0.15℃を補正すると，露場下の同じ 1.5 ｍ高度における気温は自然教育園の密な林内の気温よりも平均 0.15℃の高温となる。
要約すると，北の丸公園の開空間のほうが，自然林に近い密な自然教育園の林内に比べて昼夜ともに高温である。前者は，晴天日中には 1℃（12 月〜 4 月）ないし 2℃強（6 〜 10 月）の高温，晴天夜間は平均 0.15℃ほど高温となる。ただし，高度 1.5 ｍの気温の比較である。
北の丸公園の開空間が昼夜ともに高温になる理由は，（1）日中は林床の日射量が相対的に多く林床面が高温になること，（2）その結果，日中に蓄えられた地中熱が夜間の冷却を抑える働きをする。
注意 4（夜間の地表面温度）
上記は高度 1.5 ｍの気温が高温か低温かについて論じてきた。地表面温度と地中温度はこの傾向と異なる。晴天夜間を想定すると，開空間は天空が開けており，密な林内の林床からは天空がほとんど見えない。そのため，地表面の放射冷却は開空間のほうが強く，地表面温度は低温となる。近藤・菅─ 14 ─自然教育園報告第48号，2017

原・内藤・萩原（2015）の図 12 によれば，開空間の地表面温度は林内の地表面温度に比べて，2 〜 3℃も低いことがわかる。

まとめ

東京白金台の国立科学博物館附属の自然教育園において，2016 年 4 月〜 12 月の期間，観測塔の高度 19 ｍで観測した気温を基準として，林床上 1.5 ｍの気温（「密な林内」の気温），および森林公園内に設置されている北の丸露場の気温（「開空間」の気温）を比較した。これらは，いずれも見通し不良林であり，林床の木漏れ日率（日射量）の違いによる気温の特徴が明らかになった。
本論では，「高度 19 ｍの気温」を「基準の気温」として，「気温差＝気温−基準の気温」として定義する。
（1）晴天の日中の気温差は「開空間」の約＋ 1℃に対して「密な林内」では気温上昇は小さく，気温差は− 1℃〜− 1.5℃（6 月〜 10 月）ないしゼロに近い（12 月〜 4 月）である。
（2）夜間の気温差は両者で違いは小さくなるが，開空間のほうが 21 時〜 6 時の平均で 0.15℃ほど高温である。したがって，密な林に比べて開空間は月平均気温（年平均気温）が 0.6℃ほど高く，気温日較差が大きくなる。
北の丸公園の開空間が昼夜ともに高温となるのは，日中の林床に日射量が多く届くことが主な理由である。
ただし，これらは高度 1.5 ｍの気温についての関係であり，地表面温度と地中温度に関しては異なる。
晴天夜間は，開空間の地表面の放射冷却が強く，地表面の温度は密な林内の林床より 2 〜 3℃ほど低温になる。
今回の結果は，これまでに報告してきた前報の主な結果（図 1，図 2）と矛盾しない。本報告は，林内気温のまとめであるので，これまでに得た主な結果も以下にまとめておく。
（3）開空間の北の丸露場では，晴天日中の正午前後の平均気温はビル街の大手町より 0.5℃前後高温となる（近藤・菅原・内藤・萩原，2016；近藤・菅原・内藤・萩原，2015b）。
（4）高度 19 ｍの気温を基準とした林内の気温差は，地域代表の風速（樹冠上の風速）に大きく依存し，昼夜ともに気温差は近似的に風速の逆数に比例する。
（5）見通し良好林では，場所による大きな気温差は生じず，密な林内でも気温差は− 0.5℃前後である。
（6）特に密な林内（木漏れ日率＜ 20％）では，降雨のあった翌日の晴天日中の気温上昇が小さく，気温差は− 2℃前後となり通常より 1℃ほど低温である。
これは，雨後の林床下の土壌水分が増えて熱容量が大きいために林床面の温度上昇が抑制されることによる。大雨が続いた後に晴天が続いた場合，元に復するまでの日数は 4 〜 5 日かかる（近藤・菅原・萩原・内藤，2016 の図 15）。
（7）樹冠上で観測したボーエン比（＝顕熱輸送量／潜熱輸送量）は気温つまり季節によって大きく変化する。晴天日中のボーエン比は，15℃以下の低温期に 2 〜 5，20 〜 25℃で 0.4 〜 2，30℃以上の高温期に 0.1 〜 0.2 となる。つまり，森林から大気への顕熱輸送量は 3 月〜 5 月に最大となるのに対し，蒸散量は 6 〜 8 月に最大となる（近藤・菅原，2016）。

近藤ほか ： 林内気温の特徴

引近 藤 純 正，2014．K92

用文─ 15 ─

献

省 電 力 通 風 筒．http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke92.html（2017.2.20 閲覧）．
近藤純正，2015a．K99

通風筒の放射誤差（気象庁 95 型，農環研 09 Ｓ型）．

http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke99.html（2017.2.20 閲覧）．
近藤純正，2015b．K100．気温観測用の次世代通風筒．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke100.html（2017.2.20 閲覧）．
近藤純正，2015c．K118．北の丸露場における夜間の気温鉛直分布．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke118.html（2017.2.20 閲覧）．
近 藤 純 正，2015d．K121． 空 間 広 さ と 気 温 ─「 日 だ ま り 効 果 」 の ま と め．http://www.asahi-net.
or.jp/~rk7j-kndu/kenkyu/ke121.html（2017.2.20 閲覧）．
近藤純正・内藤玄一，2015a．K113．林内の日射量と木漏れ日率の測定．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke113.html（2017.2.20 閲覧）．
近 藤 純 正・ 内 藤 玄 一，2015b．K115． 新 宿 御 苑 の 気 温 水 平 分 布（2）．http://www.asahi-net.
or.jp/~rk7j-kndu/kenkyu/ke115.html（2017.2.20 閲覧）．
近 藤 純 正・ 内 藤 玄 一，2015c．K116． 東 京 都 心 部 の 代 表 気 温 ─ 大 手 町 露 場 の 代 表 性（ 完 結 報 ）．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke116.html（2017.2.20 閲覧）．
近藤純正・内藤玄一，2017．K140．国立市ママ下湧水公園の気温（8 〜 12 月）．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke140.html（2017.2.20 閲覧）．
近藤純正・菅原広史，2016．K123．東京都心部の森林（自然教育園）における熱収支解析．http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke123.html（2017.2.20 閲覧）．
近藤純正・菅原広史・内藤玄一・萩原信介，2015a．自然教育園内の気温．自然教育園報告，第 46 号，1-15．
近藤純正・菅原広史・内藤玄一・萩原信介，2015b．K107．林内気温の日変化・季節変化，春〜入梅期．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke107.html（2017.2.20 閲覧）．
近藤純正・菅原広史・萩原信介・内藤玄一，2015．K117．自然教育園の林内．気温，6 月〜 10 月．
http://www.asahi-net.or.jp/~rk7j-kndu/kenkyu/ke117.html（2017.2.20 閲覧）．
近藤純正・菅原広史・内藤玄一・萩原信介，2016．自然教育園内の気温．自然教育園報告，第 47 号，1-22．