青葉山を流れる小河川の晴天時における水質の経年変化
1114214 木村浩大
1114207 井元 将
1.はじめに面源負荷の一つである森林域からの流出特性を知るために、青葉山植物園を流れる小河川における水質調査を今年も引き続き行った。
過去 5 年間のデータをもとにどのような経年変化があるのかを比較した。
また、2010 年～2014 年までのデータをもとに重回帰分析を行い、それぞれの項目がどのように影響を与えているのかを調べた。
2.調査対象地
調査対象地は仙台市青葉区川内にある、青葉山植物園（面積・約 52 万 m2）である。
図１に示すように、下流を A 地点とし、上流を E 地点とする、計 5 箇所で採水・測定を行った。流下距離と流下時間の測定を行ない表１にまとめた。その結果、流下距離が 520.90m となり、流下時間は、約 2.74hであった。

図1

表１

2014
測定地点
Ａ～B距離（ｍ）
136.00流速（ｃｍ/ｓ）
3.30
流下時間（ｈ）1.14
区間距離（ｍ） 136.00
区間流下時間（ｈ） 1.14

各地点間における流下距離・流下時間

32.00
18.40
0.05

Ｂ～C
69.60
6.10
0.32
169.30
0.80

67.70
4.30
0.44

Ｃ～Ｄ
100.00 46.70
7.90
8.50
0.35
0.15
146.70
0.50

Ｄ～E
68.90
6.60
0.29
68.90
0.29

合計

520.902.74

青葉山植物園３.分析項目および調査内容（1）測定項目
測定項目は、流量、気温、水温、濁度、水素イオン濃度（pH）、電気伝導度（EC）、塩化物イオン濃度（Cl-）、化学的酸素要求量（COD）、総リン（T-P）、総窒素（T-N）、亜硝酸性窒素（NO2-N）、硝酸性窒素（NO3-N）、アンモニア性窒素（NH4-N）、TOC の全 14 種類である。
（2）先行降雨係数 APF(Antecedent Precipitation Factor)
採水前日までに降った雨 5）が、水質に与える影響を知る為に先行降雨係数を用いて検討を行った。下記の（1）式を使用し、2010 年～2014 年までのデータついて計算した。
𝑘APF=∑

𝑡=1

𝑅𝑡( 𝑡 )・・・・・・・・・・・・・・・・・・・・・・（1）

ここで Rt は、対象とする採水日の t 日前の日雨量(mm)を、t は、採水当日からの日数を表す。
（3）重回帰分析
回帰分析は量的変数の予測を目的とした方法である。回帰分析をすると、目的変数を予測する説明変数の関係式が求まる。求めた関係式のことを回帰式と言う。目的変数は予測したい変数、説明変数は予測に使う変数である。回帰分析は、説明変数が 1 つの場合を単回帰分析、複数ある場合を重回帰分析、それぞれで求めた関係式のことを単回帰式、重回帰式と言う 6）。
水温、pH、Cl-、NO3-N 、T-N、T-P、流量、APF（7）（先行降雨係数 7 日間）の 8 項目を、説明変数とし、また、COD を目的変数として重回帰分析によってそれぞれの項目がどのように影響を与えているのかを調べた。
キーワード：面源負荷、青葉山植物園、水質経年変化、先行降雨係数、多変量解析、重回帰分析
No.2－7(中山研究室)-75-

4. 結果と考察―:A地点 ―:B地点 ―:C地点 ―:D地点 ―:E地点

比較検討は、水温・pH・EC・Cl-・COD・TOC・T-N・
NH4-N ・NO2-N・NO3-N・T-P・流量の 12 項目について行う。なお、測定は 6 月～12 月の期間におよそ月 1回のペースで行った。

020

流量(m3/h)(1)各項目の経年変化

水温(℃)

2508

水温：気温と連動して変動している。全地点で同様のpH(-)

傾向が見られる。
流量：毎年、A 地点における変動が大きく、

440

E 地点における変動が小さかった。
EC(mS/cm)

pH：A 地点~D 地点までは、目立った変動は見られない。しかし、Ｅ地点では 2012 年～2013 年にかけて大きく値が減少している。

1515Cl(mg/L)

2010 年の A 地点のみで、減少している。
EC：2010 年においては、測定を行っていない。
2013 年後半に急激に減少する。

77COD(mg/L)

2014 年には A 地点で急な増加が見られた。
Cl-：全地点で値が同様に変動している。A 地点からE 地点にかけて値が増加している

05TOC(mg/L)

COD：水温の変化と同様に増加・減少を繰り返し、2011 年後半から 2012 年後半にかけて A 地点の値が著しく増加し、B 地点～E 地点は減少して

02

TOC：2014 年後半に値が急激に減少 A 地点における増加・減少の変化がある。
2011 年から 2012 年において、急な増加が見ら

T-N(mg/L)いる。

NH4-N(mg/L)

00.3

れ、B 地点～E 地点では減少傾向にあった。
T-N：全地点において同様な増加・減少の傾向が見られるが、2011 年後半から 2012 年半ばにか

00.02

NH4-N：2010 年から 2013 年にかけて増加・減少している。2014 年には、値が 0 になっている。
2014 年の E 地点で微小な増加・減少が見られる。
NO2-N：全地点で 2011 年度に著しく増加している現象が見られたが、2014 年度においては値が 00
1

NO3-N(mg/L)

NO2-N(mg/L)

けて E 地点において値の増加現象が見られる。

NO3-N：2010 年～2011 年にかけて著しい増加・減少が見られる。

T-P(mg/L)

00.04になっている。

0T-P：2011 年以降から減少している6

2014 年に値が著しく増加している。

12 6
12 6
2010
2011

12 62012年

12 6122013
2014

図2. 各項目の経年変化

図2-76-各項目の経年変化

(2)総窒素 T-N の構成による比較
図 3 に総窒素 T-N の構成を、積み上げグラフとして示す。
総窒素 T-N は、T-N＝NH4-N+NO2-N+NO3-N+有機性-N の関係がある。図 3 は、各採水地点における、年ごとの NH4-N・NO2-N・NO3-N・有機性-N の値を積み上げたグラフを示している。全体の高さが T-N濃度となる。
E 地点:2012 年～2013 年に有機性-N の占める割合において減少が見られる。
2012 年～2013 年に有機性-Nの占める割合が少なくなっている。
D 地点:2010 年～2013 年にかけて、有機性-N の占める割合が比較的、小さくなっている。しかし、2014年には、有機性-N の占める割合が大きくなっている。
C 地点:2010 年～2011 年にかけては有機性-N の占める割合に変動が見られるが、2012 年～2013 年にかけては有機性-N の占める割合は小さくなっている。
2014 年に有機性-N の占める割合が大きくなっている。
B 地点:2010 年～2011 年にかけて有機性-N の占める割合が大きくなっている。
2012 年～2013 年には有機性-N の占める割合は少なくなっている。
A 地点:他の測定地点に比べると、有機性-N の占める割合が大きくなっており、NO3-N の占める割合が小さくなっている。
全体として総窒素 T-N は、下流ほど小さくなっており、また、有機性-N の占める割合が増加しており
NO3-N の占める割合が減少している。2014 年には、
NH4-N・NO2-N の値が共にほとんど検出されず、有機性-N と NO3-N のみとなっている。NO2-N においては、値が微小であるためグラフ内で見えにくくなっている。

図 3 総窒素 T-N 構成の積み上げグラフ-77-

(3)重回帰分析の結果
COD を目的変数として標準回帰係数を求めた。その値が大きい説明変数ほどその変化によって目的変数の値が大きく変化する。図 4 に示す、重回帰分析結果から、COD の変化に大きく影響を与えるものは、水温である。また、次に影響を与えるものとして、T-N と Cl-である。pH と NO3-N の値が大きくなると、COD は減少する傾向にある。また、先行降雨係数 APF（7）は大きな影響を与えていることはないと考える。

5.まとめ

図 4 重回帰分析結果

青葉山小河川の水質変化では、経年的に亜硝酸性窒素 (NO2-N)とアンモニア性窒素（NH4-N）の値が減少していることが見られた。
総窒素 T-N における積み上げグラフによる比較においても、亜硝酸性窒素 (NO2-N)とアンモニア性窒素（NH4-N）の値が減少している。上流の E 地点において、NO3-N の占める割合が大きくなっている。しかし、下流の A 地点おいては、有機性-N の占める割合が大きくなっている。
重回帰分析の結果から、水温が最も影響を与えている。また、COD と NO3-N の関係性があり、それは、NO3-Nの値が大きくなると、COD は減少する傾向にある。経年変化のグラフにおいても同様の傾向が見られた。
参考文献
1）千葉 康成、若生 翼：森林域を流れる小河川の水質形成と先行降雨係数の関係：東北工業大学建設システム工学科卒業論文：平成 25 年度
2）永野 知明：森林域を流れる小河川における晴天時の水質形成に関する研究：東北工業大学建設システム工学科卒業論文：平成 24 年度
3）菅原 光平、工藤 洋旗：東北大学植物園内小河川における晴天時の流下に伴う水質変動：東北工業大学建設システム工学科卒業論文：平成 23 年度
4）熊井 淳仁、高橋 恵：東北大学植物園内小河川における季節ごとの水質変化について：東北工業大学建設システム工学科卒業論文：平成 22 年度
5）気象庁：
http://www.data.jma.go.jp/obd/stats/etrn/view/daily_s1.php?prec_no=34&block_no=47590&year=2014&m
onth=12&day=01&view=p1（閲覧日：2014/12/15）
6）内田 治、福島 隆司：例解多変量解析ガイド：東京図書：P24～P25：平成 24 年 3 月 10 日-78-