Silviase/jgov_v1.4_full_ann huggingface dataset:
データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{light} $ の境界線とした（式17）。

$$ f_{light}=\frac{0.003\times PPFD+1-\left[(0.003\times PPFD+1)^{2}-4\times0.003\times PPFD\times0.9\right]^{1/2}}{2\times0.9} $$

図(1)-3-Bに、ミズナラの1st flush薬における $ g_{s} $ 実測値の $ g_{\max} $ に対する相対值 $ (g_{s}/g_{\max}) $ と $ g_{s} $ 測定時の飽差（VPD, kPa）との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{VPD} $ の境界線とした（式18）。

$$ f_{VPD}=\min(1,-0.81\times VPD+3.3) $$

図(1)-3-Cに、ミズナラの1st flush薬における$g_{\mathrm{s}}$ 実測値の$g_{\mathrm{max}}$ に対する相対値$(g_{\mathrm{s}}/g_{\mathrm{max}})$ と$g_{\mathrm{s}}$ 測定時の気温（temp,℃）との関係を示した。この$g_{\mathrm{s}}/g_{\mathrm{max}}$ の分布をもとに、データ値集団の境界に位置する$g_{\mathrm{s}}/g_{\mathrm{max}}$ の近似式を$f_{\mathrm{temp}}$ の境界線とした（式19）。

$$ \begin{aligned}&f_{temp}=\left(\frac{temp-T_{min}}{T_{opt}-T_{min}}\right)\times\left(\frac{T_{max}-temp}{T_{max}-T_{opt}}\right)^{t}\quad t=\frac{T_{max}-T_{opt}}{T_{opt}-T_{min}}\\ &\\ ここで、T_{min}=15\ ^\circ C\\ &T_{opt}=32\ ^\circ C\\ &T_{max}=40\ ^\circ C\\ \end{aligned} $$

$ T_{min}, T_{opt} $ および $ T_{max} $ は、それぞれ気孔が反応する気温域における最低気温、 $ g_{s} $ が $ g_{max} $ となる気温、気孔が反応する気温域における最高気温を示している。

図(1)-3-Dに、ミズナラの1st flush薬における $ g_{s} $ 実測値の $ g_{\max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と出薬日から $ g_{s} $ 測定時までの有効積算温度（EATL，°C・day）との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{phen} $ の境界線とした（式20）。

$$ \begin{aligned}&EATL<800\ ^\circ C\cdot day 的とき、\\&\quad f_{phen}=0.0013\times EATL\\&800\leq EATL\leq1700\ ^\circ C\cdot day 的とき、\\&\quad f_{phen}=1\\&EATL>1700\ ^\circ C\cdot day 的とき、\\&\quad f_{phen}=-0.00052\times EATL+1.9\end{aligned} $$

EATLは出薬からの有効積算温度であり、出薬日からの薬の成長過程の指標として用いた。

図(1)-3-Eに、ミズナラの1st flush菓における$g_{s}$実測値の$g_{\max}$に対する相対値$\left(g_{s}/g_{\max}\right)$と$g_{s}$測定時の育成土壌における体積含水率（SWC，%）との関係を示した。この$g_{s}/g_{\max}$の分布をもとに、データ値集団の境界に位置する$g_{s}/g_{\max}$の近似式を$f_{\mathrm{SWC}}$の境界線とした（式21）。
$$ f_{SWC}=\left[1+\left(\frac{SWC}{13}\right)^{-0.4}\right]^{-1} $$

図(1)-3-Fに、ミズナラの1st flush薬における $ g_{s} $ 実測値の $ g_{\max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と大気オブン濃度 $ \left(\mathrm{O}_{3},\mathrm{ppb}\right) $ との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{O3conc} $ の境界線とした（式22）。

$$ f_{O3conc}=\left[1+\left(\frac{O_{3}}{210}\right)^{1.1}\right]^{-1} $$

## 5 ）シラカンバの春葉の水蒸気気孔コンダクタンスと環境要因の関係

図(1)-4-Aに、シラカンバの春葉における $ g_{s} $ 実測値の $ g_{max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と $ g_{s} $ 測定時の光合成有効放射束密度（PPFD, $ \upmumol\ m^{-2}\ s^{-1} $ ）との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{light} $ の境界線とした（式23）。

$$ f_{light}=\frac{0.005\times PPFD+1-\left[\left(0.005\times PPFD+1\right)^{2}-4\times0.005\times PPFD\times0.9\right]^{1/2}}{2\times0.9} $$

図(1)-4-Bに、シラカンバの春葉における $ g_{s} $ 実測値の $ g_{\operatorname*{max}} $ に対する相対値 $ (g_{\mathrm{s}}/g_{\operatorname*{max}}) $ と $ g_{s} $ 測定時の飽差（VPD, kPa）との関係を示した。この $ g_{s}/g_{\operatorname*{max}} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\operatorname*{max}} $ の近似式を $ f_{VPD} $ の境界線とした（式24）。

$$ f_{VPD}=\min(1,-1.1\times VPD+3.86) $$

図(1)-4-Cに、シラカンバの春葉における $ g_{s} $ 実測値の $ g_{\max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と $ g_{s} $ 測定時の気温（temp,℃）との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{temp} $ の境界線とした（式25）。

$$ f_{temp}=\left(\frac{temp-T_{min}}{T_{opt}-T_{min}}\right)\times\left(\frac{T_{max}-temp}{T_{max}-T_{opt}}\right)^{t}\quad t=\frac{T_{max}-T_{opt}}{T_{opt}-T_{min}} $$

$$ \begin{aligned} ここで、T_{min}&=15\ ^\circ C\\ T_{opt}&=25\ ^\circ C\\ T_{max}&=42\ ^\circ C\end{aligned} $$

$ T_{min}, T_{opt} $ および $ T_{max} $ は、それぞれ気孔が反応する気温域における最低気温、 $ g_{s} $ が $ g_{max} $ となる気温、気孔が反応する気温域における最高気温を示している。
図(1)-4-Dに、シラカンバの春葉における$g_{\mathrm{s}}$実測値の$g_{\mathrm{max}}$に対する相対値$(g_{\mathrm{s}}/g_{\mathrm{max}})$と出葉日から$g_{\mathrm{s}}$測定時までの有効積算温度（EATL, °C·day）との関係を示した。この$g_{\mathrm{s}}/g_{\mathrm{max}}$の分布をもとに、データ値集団の境界に位置する$g_{\mathrm{s}}/g_{\mathrm{max}}$の近似式を$f_{\mathrm{phen}}$の境界線とした（式26）。

$$ \begin{aligned}EATL<150\ ^\circ C\cdot day 的とき、\\f_{phen}=0.004\times EATL+0.4\\150\leqq EATL\leqq450\ ^\circ C\cdot day 的とき、\\f_{phen}=1\\ EATL>450\ ^\circ C\cdot day 的とき、\\f_{phen}=-0.00038\times EATL+1.2\end{aligned} $$

EATLは出薬からの有効積算温度であり、出薬日からの薬の成長過程の指標として用いた。

図(1)-4-Eに、シラカンバの春葉における $ g_{s} $ 実測値の $ g_{max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と $ g_{s} $ 測定時の育成土壤における体積含水率（SWC，%）との関係を示した。この $ g_{s}/g_{max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{max} $ の近似式を $ f_{SWC} $ の境界線とした（式27）

$$ f_{SWC}=\left[1+\left(\frac{SWC}{14}\right)^{-0.4}\right]^{-1} $$

図(1)-4-Fに、シラカンバの春葉における $ g_{s} $ 実測値の $ g_{\max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と大気オブン濃度 $ \left(\mathrm{O}_{3},\mathrm{ppb}\right) $ との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{O3conc} $ の境界線とした（式28）。

$$ f_{\mathrm{O3c o n c}}=\left[1+\left(\frac{\mathrm{O}_{3}}{90}\right)^{6}\right]^{-1} $$

## 6 ）シラカンバの夏葉の水蒸気気孔コンダクタンスと環境要因の関係

図(1)-5-Aに、シラカンバの夏葉における $ g_{s} $ 実測値の $ g_{max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と $ g_{s} $ 測定時の光合成有効放射束密度（PPFD, $ \upmumol\ m^{-2}\ s^{-1} $ ）との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{light} $ の境界線とした（式29）。

$$ f_{light}=\frac{0.005\times PPFD+1-\left[\left(0.005\times PPFD+1\right)^{2}-4\times0.005\times PPFD\times0.7\right]^{1/2}}{2\times0.7} $$

図(1)-5-Bに、シラカンバの夏葉における$g_{\mathrm{s}}$実測値の$g_{\mathrm{max}}$に対する相対値$(g_{\mathrm{s}}/g_{\mathrm{max}})$と$g_{\mathrm{s}}$測定時の飽差（VPD, kPa）との関係を示した。この$g_{\mathrm{s}}/g_{\mathrm{max}}$の分布をもとに、データ値集団の境界に位置する$g_{\mathrm{s}}/g_{\mathrm{max}}$の近似式を$f_{\mathrm{VPD}}$の境界線とした（式30）。

$$ f_{VPD}=\min(1,-0.4\times VPD+1.77) $$
図(1)-5-Cに、シラカンバの夏葉における $ g_{s} $ 実測値の $ g_{max} $ に対する相対値 $ (g_{s}/g_{\max}) $ と $ g_{s} $ 測定時の気温（temp,℃）との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{temp} $ の境界線とした（式31）。

$$ f_{temp}=\left(\frac{temp-T_{min}}{T_{opt}-T_{min}}\right)\times\left(\frac{T_{max}-temp}{T_{max}-T_{opt}}\right)^{t}\quad t=\frac{T_{max}-T_{opt}}{T_{opt}-T_{min}} $$

$$ \begin{aligned} ここで、T_{min}&=13\ ^\circ C\\ T_{opt}&=28\ ^\circ C\\ T_{max}&=41\ ^\circ C\end{aligned} $$

$ T_{min}, T_{opt} $ および $ T_{max} $ は、それぞれ気孔が反応する気温域における最低気温、 $ g_{s} $ が $ g_{max} $ となる気温、気孔が反応する気温域における最高気温を示している。

図(1)-5-Dに、シラカンバの夏葉における$g_{\mathrm{s}}$実測値の$g_{\mathrm{max}}$に対する相対値$(g_{\mathrm{s}}/g_{\mathrm{max}})$と出葉日から$g_{\mathrm{s}}$測定時までの有効積算温度（EATL, °C・day）との関係を示した。この$g_{\mathrm{s}}/g_{\mathrm{max}}$の分布をもとに、データ値集団の境界に位置する$g_{\mathrm{s}}/g_{\mathrm{max}}$の近似式を$f_{\mathrm{phen}}$の境界線とした（式32）。

$$ \begin{aligned}EATL&<150\ ^{\circ}C\cdot day 的とき、\\&f_{phen}=0.0067\times EATL\\150&\leqq EATL\leqq400\ ^{\circ}C\cdot day 的とき、\\&f_{phen}=1\\ EATL&>400\ ^{\circ}C\cdot day 的とき、\\&f_{phen}=-0.00043\times EATL+1.2\end{aligned} $$

EATLは出菓からの有効積算温度であり、出菓日からの葉の成長過程の指標として用いた。

図(1)-5-Eに、シラカンバの夏葉における $ g_{s} $ 実測値の $ g_{\operatorname*{max}} $ に対する相対値 $ (g_{s}/g_{\operatorname*{max}}) $ と $ g_{s} $ 測定時の育成土壤における体積含水率（SWC，%）との関係を示した。この $ g_{s}/g_{\operatorname*{max}} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\operatorname*{max}} $ の近似式を $ f_{SWC} $ の境界線とした（式33）。

$$ f_{SWC}=\left[1+\left(\frac{SWC}{13}\right)^{-0.6}\right]^{-1} $$

図(1)-5-Fに、シラカンバの夏葉における $ g_{s} $ 実測値の $ g_{\operatorname*{max}} $ に対する相対値 $ (g_{\mathrm{s}}/g_{\mathrm{max}}) $ と大気オブン濃度 $ \left(\mathrm{O}_{3},\mathrm{ppb}\right) $ との関係を示した。この $ g_{s}/g_{\max} $ の分布をもとに、データ値集団の境界に位置する $ g_{s}/g_{\max} $ の近似式を $ f_{03conc} $ の境界線とした（式34）。

$$ f_{\mathrm{O3c o n c}}=\left[1+\left(\frac{\mathbf{O}_{3}}{100}\right)^{7}\right]^{-1} $$
## 7 ）水蒸気気孔コンダクタンスとオゾン吸収量の推定式の構築

Stomatal conductance multiple modelは、気象条件等の環境要因が $ g_{s} $ に相加的に作用することを前提としたモデルである。このモデルを用い、各環境要因に対する $ g_{s} $ の応答変数を $ g_{max} $ に掛け合わせることで $ g_{s} $ の推定値 $ (g_{\mathrm{sto}}) $ を算出する推定モデル $ ^{[10], [11], [12]} $ を構築した（式35）。

$$ g_{sto}=g_{max}\times f_{phen}\times f_{light}\times max\left\{f_{min},\left(f_{temp}\times f_{VPD}\times f_{SWC}\times f_{O3conc}\right)\right\} $$

なお、 $ f_{\min} $ は $ g_{\min} $ と $ g_{\max} $ の比 $ (g_{\min}/g_{\max}) $ である。

図(1)-6に、本研究で開発した各樹種の薬の $ g_{s} $ 推定モデルから算出した $ g_{s} $ の推定値と $ g_{s} $ の実測値との関係を示した。ブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬、シラカンバの春葉および夏葉における $ g_{s} $ の実測値と推定値との間に有意な正の相関関係が認められ、その回帰直線の傾きはそれぞれ1.02, 0.83, 0.87, 0.82および0.80であり、切片はそれぞれ53, 89, 112, 166および194であった。これらの結果から、各樹種の薬の $ g_{s} $ を気象条件等の常時測定値から推定することが可能になった。また、 $ g_{s} $ の推定値を式4に代入し、式2と式3および風速とオゾン濃度の常時測定値を用いて $ F_{st} $ を推定することで、ブナ、コナラ、ミズナラおよびシラカンバの薬のオゾン吸収量を推定することが可能になった。
大気オゾン濃度 (ppb)

図(1)-1 ブナの1st flush薬の水蒸気気孔コンダクタンス $ (g_{\mathrm{s}}) $ の実測値の最大水蒸気気孔コンダクタンス $ (g_{\mathrm{max}}) $ に対する相対値 $ (g_{\mathrm{s}}/g_{\mathrm{max}}) $ と光合成有効放射束密度(A)、飽差(B)、気温(C)、出薬からの有効積算温度(D)、土壤体積含水率(E)、大気オゾン濃度(F)との関係.
大気オゾン濃度 (ppb)

図(1)-2 コナラの1st flush薬の水蒸気気孔コンダクタンス $ (g_{\mathrm{s}}) $ の実測値の最大水蒸気気孔コンダクタンス $ (g_{\mathrm{max}}) $ に対する相対値 $ (g_{\mathrm{s}}/g_{\mathrm{max}}) $ と光合成有効放射束密度(A)、飽差(B)、気温(C)、出薬からの有効積算温度(D)、土壤体積含水率(E)、大気オゾン濃度(F)との関係.
大気オゾン濃度 (ppb)

図(1)-3 ミズナラの1st flush薬の水蒸気気孔コンダクタンス $ (g_{s}) $ の実測値の最大水蒸気気孔コンダクタンス $ (g_{\mathrm{max}}) $ に対する相対値 $ (g_{\mathrm{s}}/g_{\mathrm{max}}) $ と光合成有効放射束密度(A)、飽差(B)、気温(C)、出薬からの有効積算温度(D)、土壤体積含水率(E)、大気オゾン濃度(F)との関係.
光合成有効放射束密度（ $ \mu $ mol m $ ^{-2} $ s $ ^{-1} $ ）

出薬からの有効積算温度（℃・day）

氖温 (℃)

大気オゾン濃度 (ppb)

図(1)-4 シラカンバ春葉の水蒸気気孔コンダクタンス $ (g_{s}) $ の実測値の最大水蒸気気孔コンダクタンス $ (g_{\mathrm{max}}) $ に対する相対値 $ (g_{s}/g_{\mathrm{max}}) $ と光合成有効放射束密度(A)、飽差(B)、気温(C)、出薬からの有効積算温度(D)、土壌体積含水率(E)、大気オゾン濃度(F)との関係.
氖温(℃)

大気オゾン濃度 (ppb)

図(1)-5 シラカンバの夏葉の水蒸気気孔コンダクタンス $ (g_{s}) $ の実測値の最大水蒸気気孔コンダクタンス $ (g_{\mathrm{max}}) $ に対する相対値 $ (g_{s}/g_{\mathrm{max}}) $ と光合成有効放射束密度(A)、飽差(B)、気温(C)、出薬からの有効積算温度(D)、土壌体積含水率(E)、大気オゾン濃度(F)との関係.
図(1)-6 ブナの1st flush薬(A)、コナラの1st flush薬(B)、ミズナラの1st flush薬(C)、シラカンバの春薬(D)および夏葉(E)の水蒸気気孔コンダクタンス( $ g_{s} $ )の実測値と推定値( $ g_{sto} $ )との関係．図中の実線は回帰直線であり、破線は1:1の直線を示す。回帰直線の式とその決定係数を図中に示した。Pearson's correlation test: ***p<0.001.
## 8 ）環境要因による葉の水蒸気気孔コンダクタンスの律速程度

図(1)-7に、2013年における1.5倍O_{3}区で育成したブナ、コナラおよびミズナラの1st flush薬とシラカンバの春葉および夏葉の出薬日からの単位葉面積あたりの積算オゾン吸収量（COU, mmol m^{-2}）とオゾン濃度の日平均値の経時変化を示した。ブナの1st flush薬のCOUは8月に最も高く、コナラの1st flush薬およびミズナラの1st flush薬のそれは6月に最も高かった。また、シラカンバの春葉および夏葉のCOUは5月に最も高かった。2013年におけるブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬、シラカンバの春葉および夏葉のCOUは、それぞれ38.5, 25.4, 35.0, 37.2および35.1 mmol m^{-2}であった。

1.5倍O₃区で育成したブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬、シラカンバの春葉および夏葉の日中の大気オゾン濃度の1時間値の積算値（Daylight AOT0, ppm h）は、5月に最も高かった。同様に、1.5倍O₃区で育成したシラカンバの春葉および夏葉のCOUは、5月に最も多かった。これに対して、1.5倍O₃区で育成したブナの1st flush薬のCOUは8月に最も高く、コナラとミズナラの1st flush薬のCOUは6月に最も高かった。このように、ブナ、コナラおよびミズナラの1st flush薬において、オゾンのDaylight AOT0とCOUが最も高くなる月は異なっていた。地中海地域の夏季においても、このような相違が生じることが報告されている $ ^{13)} $ 。地中海地域の夏季においては、強い日射、高い気温および大気の乾燥によって光化学反応が促進され、オゾン濃度が高くなる $ ^{14)} $ 。しかしながら、夏季における大気の乾燥は気孔を閉鎖させるため、気孔を介した葉のオゾン吸収量は比較的少なかった $ ^{13)} $ 。このようなオゾンのAOT0とCOUが最も高くなる時期の相違は、樹木に対するオゾンの影響を評価する際に、AOT0ではなく、COUを用いるべきであることを示している。水蒸気気孔コンダクタンスの最大値（g_{max}）は、シラカンバの夏葉＞ミズナラの1st flush薬＞ブナの1st flush薬＞シラカンバの春葉＞コナラの1st flush薬の順に高かった。一方、2013年の1.5倍オゾン区の出薬日からの単位葉面積あたりの積算オゾン吸収量（COU, mmol m^{-2}）は、ブナの1st flush薬＞シラカンバの春葉＞シラカンバの夏葉＞ミズナラの1st flush薬＞コナラの1st flush薬の順に高かった（図(1)-7）。これらの結果は、g_{max}の高い葉が必ずしもCOUが高い訳ではなく、環境要因に対するg_{s}の応答がCOUの決定要因として重要であることを示している。

図(1)-7 2013年における1.5倍O_{3}区で育成したブナ、コナラおよびミズナラの1st flush薬とシラカンバの春葉および夏葉の出薬日からの単位薬面積あたりの積算オゾン吸收量(COU, mmol m^{-2})と大気オゾン濃度の日平均値の経時変化.
図(1)-8に、1.5倍O_{3}区で育成したブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬、シラカンバの春葉および夏葉の2013年の春季（4～5月）、夏季（6～8月）および秋季（9～10月）におけるLimiting degreeを示した。Limiting degree（D_{factor}:D_{phen},D_{light},D_{temp},D_{VPD},D_{SWC},D_{O3conc}）とは各環境要因による水蒸気気孔コンダクタンス（g_{s},mmol H_{2}O m^{-2}s^{-1}）の律速程度であり、各環境要因に対するg_{s}の応答変数（f_{factor}:f_{phen},f_{light},f_{temp},f_{VPD},f_{SWC},f_{O3conc}）を用いて以下のように定義した（式36）。

$$ D_{factor}=1-f_{factor} $$

春季においては、ブナ、コナラおよびミズナラの1st flush薬では $ D_{phen} $ が最も高く、それぞれ0.72, 0.63および0.71であった。シラカンバ春葉および夏菓の $ D_{phen} $ はそれぞれ0.09および0.19であり、ブナ、コナラおよびミズナラの1st flush薬のそれと比較して低かった。また、シラカンバの春葉および夏菓では $ D_{temp} $ が最も高く、それぞれ0.32および0.26であった。夏季においては、ブナおよびミズナラの1st flush薬では $ D_{light} $ が最も高く、それぞれ0.19および0.25であった。コナラの1st flush薬では $ D_{temp} $ が最も高く、0.019であった。シラカンバの春葉および夏菓では $ D_{phen} $ が最も高く、それぞれ0.40および0.40であった。秋季においては、ブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬、シラカンバの春葉および夏菓において $ D_{phen} $ が最も高く、それぞれ0.27, 0.53, 0.52, 0.90および0.95であった。

日本のオゾン濃度は、春季（3～5月）に比較的高い $ ^{15),16} $ 。したがって、樹木に対するオゾンの影響評価において、春季における葉のオゾン吸収量は重要である。1.5倍O $ _{3} $ 区で育成したブナ、コナラおよびミズナラの1st flush葉における春季のCOU（5.6,4.1および5.2 mmol m $ ^{-2} $ ）が一成長期間のCOU（38.5,25.4および35.0 mmol m $ ^{-2} $ ）に占める割合は比較的低く、それぞれ14.5,16.0および15.0%であった。一方、1.5倍O $ _{3} $ 区で育成したシラカンバの春葉およびシラカンバの夏葉における春季のCOU（19.3および14.7 mmol m $ ^{-2} $ ）が一成長期間のCOU（37.2および35.1 mmol m $ ^{-2} $ ）に占める割合は高く、それぞれ51.8%および42.0%であった。このような春季におけるCOUの樹種間差異の原因を明らかにするために、1.5倍O $ _{3} $ 区で育成したブナ、コナラおよびミズナラの1st flush葉とシラカンバの春葉および夏葉の春季におけるLimiting degree（D $ _{factor} $ :D $ _{phen} $ ,D $ _{light} $ ,D $ _{temp} $ ,D $ _{VPD} $ ,D $ _{SWC} $ ,D $ _{O3conc} $ ）を計算した。その結果、ブナ、コナラおよびミズナラの1st flush葉ではD $ _{phen} $ が最も高く、それぞれ0.72,0.63および0.71であった（図(1)-8）。一方、シラカンバの春葉および夏葉のD $ _{phen} $ は、ブナ、コナラおよびミズナラの1st flush葉のそれと比較して著しく低く、それぞれ0.09および0.19であった。その他のLimiting degree（D $ _{light} $ ,D $ _{temp} $ ,D $ _{VPD} $ ,D $ _{SWC} $ ,D $ _{O3conc} $ ）に樹種による著しい違いは認められなかった。f $ _{phen} $ のパラメータの一つであるphen_bは、葉の完全展開や成熟に伴って気孔が最も開くようになる時（f $ _{phen} $ =1）の有効積算温度（EATL,℃day）を示している。ブナの1st flush葉、コナラの1st flush葉、ミズナラの1st flush葉、シラカンバの春葉および夏葉のphen_bは、それぞれ850,600,800,150および150℃・dayであった。この結果は、シラカンバの春葉および夏葉は、ブナ、コナラおよびミズナラの1st flush葉に比べて、出葉後に短期間でg $ _{s} $ が高くなることを示している。以上のことから、ブナ、コナラおよびミズナラの1st flush葉は成熟が比較的遅いため、春季において大気オゾン濃度は比較的高かったが、気孔を介した葉のオゾン吸収量は低いことが明らかになった。これに対して、シラカンバの春葉および夏葉は成熟が早く、出葉後に短期間でg $ _{s} $ が高くなるため、春季のオゾン吸収量が高いことが明らかになった。
図(1)-8 1.5倍O_{3}区で育成したブナの1st flush葉(緑)、コナラの1st flush葉(黄)、ミズナラの1st flush葉(青)、シラカンバの春葉(黒)および夏葉(赤)の2013年の春季(4～5月)、夏季(6～8月)および秋季(9～10月)におけるLimiting degree. Limiting degree (D_{factor}: D_{phen}, D_{light}, D_{temp}, D_{VPD}, D_{SWC}, D_{O3conc})とは各環境要因による水蒸気気孔コンダクタンス(g_{s})の律速程度であり、各環境要因に対するg_{s}の応答変数(f_{factor}: f_{phen}, f_{light}, f_{temp}, f_{VPD}, f_{SWC}, f_{O3conc})を用いて定義した。D_{factor}=1-f_{factor}.
## （２）日本の落葉広葉樹の積算純光合成量に対するオゾンの影響の現状評価

## 1 ）積算純光合成量に対するオゾンの影響

図(1)-9に、サブテーマ（1）とサブテーマ（2）のオゾン暴露実験から得られたブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬およびシラカンバの春薬の積算純光合成量の相対値（%）と出薬日からの単位薬面積あたりの積算オゾン吸收量（COU, mmol m⁻²）との関係を示した。積算純光合成量の相対値とは、浄化空気区（CF区）の積算純光合成量に対する各オゾン処理区の積算純光合成量の相対値である。ブナの1st flush薬の積算純光合成量の相対値（Y, %）と単位薬面積あたりの積算オゾン吸收量（X, mmol m⁻²）との回帰曲線（ウェイブル関数，Y=exp[-(X/ω)^]、ω、λは定数）は、Y=100×exp[-(X/83.2)^3.08]であった。コナラの1st flush薬の積算純光合成量の相対値（Y, %）と単位薬面積あたりの積算オゾン吸收量（X, mmol m⁻²）との回帰曲線は、Y=100×exp[-(X/69.0)^2.59]であった。ミズナラの1st flush薬の積算純光合成量の相対値（Y, %）と単位薬面積あたりの積算オゾン吸收量（X, mmol m⁻²）との回帰曲線は、Y=100×exp[-(X/66.5)^3.51]であった。シラカンバの春薬の積算純光合成量の相対値（Y, %）と単位薬面積あたりの積算オゾン吸收量（X, mmol m⁻²）との回帰曲線は、Y=100×exp[-(X/48.3)^4.41]であった。いずれの樹種においても、COUの増加に伴って積算純光合成量の相対値は低下した。オゾンによる積算純光合成量の相対値の低下程度には樹種間差異があり、シラカンバ>コナラ>ミズナラ>シラカンバの順に高かった。この結果は、気孔から薬内に同量のオゾンを吸收しても、オゾンによる積算純光合成量の低下程度が樹種によって異なることを示している。

○ 赤城(2012年) △ 赤城(2013年) □ 八王子(2013年)

図(1)-9 赤城(〇:2012年，△:2013年)および八王子(□:2013年)におけるブナの1st flush薬、コナラの1st flush薬、ミズナラの1st flush薬およびシラカンバの春薬の積算純光合成量の相対値(%)と出薬日からの単位薬面積あたりの積算オゾン吸收量(COU, mmol m⁻²)との関係.積算純光合成量の相対値とは、浄化空気区(CF区)の積算純光合成量に対する各オゾン処理区の積算純光合成量の相対値である。 ブナ: Y=100×exp[-（X/83.2）³.⁰⁸] コナラ: Y=100×exp[-（X/69.0）².⁵⁹] ミズナラ: Y=100×exp[-（X/66.5）³.⁵¹] シラカンバ: Y=100×exp[-（X/48.3）⁴.⁴¹]
## 2 ）現状のオゾンによる日本に分布する落葉広葉樹の積算純光合成量低下率の評価

図(1)-10に、2009～2011年における日本のブナ分布域の日中（日射量：>50 W m^{-2}）の大気オブン濃度の1時間値の積算値（Daylight AOT0）の年平均値と1st flush薬の単位薬面積あたりの積算オブン吸収量（COU）の年平均値を示した。ブナ分布域におけるオブンのDaylight AOT0および1st flush薬のCOUの全国平均値は、それぞれ80.0 ppm hおよび43.0 mmol m^{-2}であった。

図(1)-11に、2009～2011年における日本のコナラ分布域のオブンのDaylight AOT0の年平均値と1st flush薬のCOUの年平均値を示した。コナラ分布域におけるオブンのDaylight AOT0および1st flush薬のCOUの全国平均値は、それぞれ82.3 ppm hおよび27.9 mmol m^{-2}であった。

図(1)-12に、2009～2011年における日本のミズナラ分布域のオゾンのDaylight AOT0の年平均値と1st flush薬のCOUの年平均値を示した。ミズナラ分布域におけるオゾンのDaylight AOT0および1st flush薬のCOUの全国平均値は、それぞれ76.0 ppm hおよび35.6 mmol m^{-2}であった。

図(1)-13に、2009～2011年における日本のシラカンバ分布域のオゾンのDaylight AOT0の年平均値と春葉のCOUの年平均値を示した。シラカンバ分布域におけるオゾンのDaylight AOT0および春葉のCOUの全国平均値は、それぞれ58.6 ppm hおよび31.7 mmol m^{-2}であった。
図(1)-10 2009～2011年における日本のブナ分布域のオゾンのDaylight AOT0の年平均値と1st flush葉の単位葉面積あたりの積算オゾン吸収量(COU)の年平均値.
図(1)-11 2009～2011年における日本のコナラ分布域のオゾンのDaylight AOT0の年平均値と1st flush葉の単位葉面積あたりの積算オゾン吸収量(COU)の年平均値.
図(1)-12 2009～2011年における日本のミズナラ分布域のオゾンのDaylight AOT0の年平均値と1st flush葉の単位葉面積あたりの積算オゾン吸收量(COU)の年平均値.
図(1)-13 2009～2011年における日本のシラカンバ分布域のオゾンのDaylight AOT0の年平均値と春葉の単位葉面積あたりの積算オゾン吸収量(COU)の年平均値.
日本に分布するブナ、コナラ、ミズナラおよびシラカンバの2009～2011年におけるオブンによる積算純光合成量の低下率（%）を積算純光合成量の相対値（%）と出薬日からの単位薬面積あたりの積算オブン吸収量（COU, mmol m⁻²）との関係（図(1)-9）から推定した。

図(1)-14に、積算オゾン吸收量（COU）に基づいて算出したブナの1st flush葉の積算純光合成量の低下率を示した。積算純光合成量の低下率は1～32%であり、その全国平均値は12%であった。

図(1)-15に、積算オゾン吸収量（COU）に基づいて算出したコナラの1st flush葉の積算純光合成量の低下率（%）を示した。積算純光合成量の低下率は2～17%であり、その全国平均値は10%であった。

図(1)-16に、積算オゾン吸収量（COU）に基づいて算出したミズナラの1st flush葉の積算純光合成量の低下率（%）を示した。積算純光合成量の低下率は1〜32%であり、その全国平均値は12%であった。

図(1)-17に、積算オゾン吸収量（COU）に基づいて算出したシラカンバの春葉の積算純光合成量の低下率（%）を示した。積算純光合成量の低下率は9〜58%であり、その全国平均値は16%であった。

オゾンによる積算純光合成量の低下率(%)

図(1)-14 出葉日からの単位葉面積あたりの積算オゾン吸収量(COU, mmol m^{-2})に基づいて算出した2009～2011年におけるブナの1st flush葉の積算純光合成量の年平均低下率(%).
図(1)-15 出葉日からの単位葉面積あたりの積算オゾン吸収量（COU, mmol m^{-2}）に基づいて算出した2009～2011年におけるコナラの1st flush葉の積算純光合成量の年平均低下率(%)．

オゾンによる積算純光合成量の低下率(%)

図(1)-16 出葉日からの単位葉面積あたりの積算オゾン吸収量(COU, mmol m^{-2})に基づいて算出した2009～2011年におけるミズナラの1st flush葉の積算純光合成量の年平均低下率(%)。
図(1)-17 出葉日からの単位葉面積あたりの積算オゾン吸収量(COU, mmol m^{-2})に基づいて算出した2009～2011年におけるシラカンバの春葉の積算純光合成量の年平均低下率(%)。

## 3 ）日本における森林衰退に対するオゾンの関与

本研究の結果を基に、日本における森林衰退に対するオゾンの関与を検討する。日本において、コナラやミズナラの森林衰退が観察されているが、主な原因としてはナラ枯れが指摘されており $ ^{17,18} $ 、それらの森林衰退に対するオゾンの関与の判断は困難である。一方、ブナ林の衰退の主な原因として、オゾンが指摘されている $ ^{19,20,21,22,23} $ 。例えば、丹沢山地においては、現状濃度レベルのオゾンによるブナ苗の葉のクロロフィル含量（SPAD値）の低下、落葉の早期化および個体乾物成長の低下などが認められている $ ^{21} $ 。

図(1)-18に、日本における森林衰退や葉の異常が観察されているブナの分布地域 $ ^{20} $ ， $ ^{22} $ ， $ ^{23} $ と2009～2011年におけるブナの単位葉面積あたりの積算オゾン吸収量（COU）に基づいて算出した積算純光合成量の低下率を示した。なお、葉の異常の判断指標は、その大きさ、色、形、空間的密度および枯葉数である $ ^{23} $ 。丹沢山地（神奈川県）、英彦山（福岡県）、脊振山（福岡県）、函南原生林（静岡県）、富士山（静岡県）および立山（富山県）におけるオゾンによる積算純光合成量の低下率は15～20%であり、その全国平均値（12%）に比べて高かった。さらに、臥竜山（広島県）におけるオゾンによる積算純光合成量の低下率は20～25%であり、その全国平均値に比べて極めて高かった。したがって、日本に分布するブナ林の衰退や葉の異常にオゾンが関与している可能性が十分にある。
オゾンによる積算純光合成量の低下率(%)

図(1)-18 日本における森林衰退や葉の異常が観察されているブナの分布地域と2009～2011年におけるブナの単位葉面積あたりの積算オゾン吸收量(COU)に基づいて算出した積算純光合成量の年平均低下率(%).

## 5 ．本研究により得られた成果

## （1） 科学的意義

サブテーマ（1）では、光合成有効放射束密度、気温、飽差、有効積算温度、土壌体積含水率および大気オゾン濃度に基づいて、日本の代表的な落葉広葉樹であるプナ、コナラ、ミズナラおよびシラカンバの葉の水蒸気気孔コンダクタンスの推定式を構築し、気孔を介した葉のオゾン吸収量を算出することが可能になった。その結果、比較的高濃度のオゾンが観測された春季においては、プナ、コナラおよびミズナラの1st flush葉は未成熟であったため、気孔開度が比較的低く、オゾン吸収量が低いことが明らかになった。一方、シラカンバの春葉および夏葉は葉の成熟が早く、春季においても葉のオゾン吸収量が比較的高いことが明らかになった。

薬の積算オブン吸収量に基づいて、日本に分布する落葉広葉樹4種の積算純光合成量に対するオブンの影響評価を行った結果、ブナ、コナラおよびミズナラの1st flush葉とシラカンバの春葉の積算純光合成量のオブンによる年平均低下率はそれぞれ12%, 10%, 12%および16%と推定された。したがって、アジア大陸からの越境大気汚染対策や日本国内のオブンの前駆物質の排出源対策を行うことで、日本の森林の保護だけでなく、その地球温暖化防止機能を向上させることができることを科学的に示した。
## （２）環境政策への貢献

## ＜行政が既に活用した成果＞

環境省の越境大気汚染・酸性雨対策検討会の生態影響分科会の下にあるオゾン等の植物影響評価ワーキンググループ（座長：伊豆田猛）において、本研究で得られた研究成果などを提示し、平成26年3月に公表された環境省の越境大気汚染・酸性雨長期モニタリング（平成20～24年度）最終報告書の内容（119～123ページ）の検討過程で多大な貢献をした。同ワーキンググループで実施しているオゾンによる植物影響のパイロットモニタリングの地点選定過程で本研究で得られた研究成果を提示し、福岡県（英彦山）、新潟県（八海山）および北海道（摩周湖外輪山）におけるオゾン濃度の測定と植物被害の評価が開始され、その結果が平成26年3月に公表された環境省の越境大気汚染・酸性雨長期モニタリング（平成20～24年度）最終報告書に掲載された（123～129ページ）。

## ＜行政が活用することが見込まれる成果＞

本研究では、気孔を介した葉のオゾン吸収量に基づいた日本の落葉広葉樹のCO₂吸収・固定能力に対するオゾンの影響の評価手法を確立したため、日本の現状レベルのオゾンによる落葉広葉樹の年間CO₂吸収量の低下程度が評価できるようになった。したがって、大気汚染対策と温暖化対策のコベネフィットアプローチに対して科学的知見を提供できる。

本研究によって、日本におけるオゾンによる森林衰退の危険地域（ホットスポット）の選択が可能になったため、越境大気汚染による生態系影響の発現予測やその対策に貢献できる。

本研究による森林のCO₂固定能へのオゾンの影響評価は、気候変動に関する政府間パネル（Intergovernmental Panel on Climate Change, IPCC）の報告書などの国際的な科学的知見の集積に直接的に貢献し、温暖化防止を目的として東アジア各国の連携の下で国際的大気汚染対策を推進するための科学的根拠を提供できる。

本研究は、現在、オゾンなどのガス状大気汚染物質の測定が検討されている東アジア酸性雨ネットワーク（EANET）におけるオゾン測定地点の設定やオゾンの植物影響評価地点の選定などに対して科学的根拠が提供できるため、多大な貢献が期待できる。

## 6 ．国際共同研究等の状況

特に記載すべき事項はない。
## 7 ．研究成果の発表状況

## （1） 誌上発表

## <論文（査読あり）>

1) M. YAMAGUCHI, M. WATANABE, H. MATSUMURA, Y. KOHNO and T. IZUTA: Asian J. Atmos. Environ., 5, 65-78 (2011) “Experimental Studies on the Effects of Ozone on Growth and Photosynthetic Activity of Japanese Forest Tree Species”

2) M. WATANABE, M. YAMAGUCHI, H. MATSUMURA, Y. KOHNO, T. KOIKE and T. IZUTA: Asian J. Atmos. Environ., 5, 205-215 (2011) “A Case Study of Risk Assessment of Ozone Impact on Forest Tree Species in Japan”

3）伊豆田 猛：大気環境学会誌, 47, 1, A12-A15（2012）「葉のオゾン吸収量に基づいた植物に対するオゾンの影響評価」

4) M. WATANABE, M. YAMAGUCHI, H. MATSUMURA, Y. KOHNO and T. IZUTA: Eur. J. Forest Res., 131, 475-484 (2012) "Risk Assessment of Ozone Impact on Fagus crenata in Japan: Consideration of Atmospheric Nitrogen Deposition"

5) F. AZUCHI, Y. KINOSE, T. MATSUMURA, T. TANOMATA, Y. UEHARA, A. KOBAYASHI, M. YAMAGUCHI and T. IZUTA: Environ. Pollut. 184: 481-487 (2014) "Modeling Stomatal Conductance and Ozone Uptake of Fagus crenata Grown under Different Nitrogen Loads"

6) T. KOIKE, M. WATANABE, Y. HOSHIKA, M. KITAO, H. MATSUMURA, R. FUNADA and T. IZUTA: Climate Change, Air Pollution and Global Challenges: Understanding and Solutions from Forest Research (R. Matyssek et al., eds.), Elsevier. (2013) “Effects of Ozone on Forest Ecosystems in East and Southeast Asia”

## くその他誌上発表（査読なし）>

1）公害防止の技術と法規 編集委員会 編：新・公害防止の技術と法規 2012 大気編 I、（社）産業環境管理協会、168-179（2012）「植物に対する大気汚染物質の影響（執筆担当：伊豆田猛）」

2）伊豆田 猛、船田 良：北方林業，65, 23-24（2013）「我が国における対流圏オゾンの植物影響評価に関する試み」

3）公害防止の技術と法規 編集委員会 編：新・公害防止の技術と法規 2013 大気編 I、（社）産業環境管理協会、175-186（2013）「植物に対する大気汚染物質の影響（執筆担当：伊豆田猛）」

4）公害防止の技術と法規 編集委員会 編：新・公害防止の技術と法規 2014 大気編 I、（社）産業環境管理協会、116-128（2014）「植物に対する大気汚染物質の影響（執筆担当：伊豆田猛）」

## （2） 口頭発表（学会等）

1）伊豆田 猛：大気環境学会関東支部講演会（2011）「葉のオゾン吸収量に基づいた植物に対するオゾンの影響評価」
2）伊豆田 猛:第28回エアロゾル科学・技術研究討論会（日本エアロゾル学会）特別講演（2011）「植物に対する越境大気汚染物質の影響」

3）安土文鹿、母袋朗、松村友絵、伊藤伸一郎、伊豆田猛：第52回大気環境学会（2011）「ブナ苗の成長、光合成および炭素固定量に対するオゾンの影響」

4）渡辺 誠，山口真弘，松村秀幸，河野吉久，小池孝良，伊豆田 猛：第123回日本森林学会（2012）「日本の森林樹種に対するオゾンのリスク評価，窒素沈着の考慮」

5）山口真弘、渡辺 誠、松村秀幸、伊豆田 猛：第123回日本森林学会（2012）「ブナ苗の光合成活性と窒素代謝に対するオゾンと土壌への窒素負荷の複合影響」

6）伊豆田 猛：大気環境学会 植物分科会 全国講演会（2012）「大気汚染と植物」

7）渡辺 誠、山口真弘、松村秀幸、河野吉久、伊豆田 猛：大気環境学会 植物分科会 全国講演会（2012）「窒素沈着量を考慮したブナに対するオゾンのリスク評価」

8）黄瀬佳之、安土文鹿、松村友絵、上原唯、鹿又友彰、小林亜由美、山口真弘、伊豆田猛：第53回大気環境学会年会（2012）「ブナの葉のオゾン吸収量推定モデルの構築」

9）安土文鹿、松村友絵、上原唯、鹿又友彰、黄瀬佳之、小林亜由美、山口真弘、伊豆田猛：第53回大気環境学会年会（2012）「葉のオゾン吸収量に基づくブナ苗に対するオゾンと土壌への窒素負荷の複合影響の評価」

10) T. IZUTA : Workshop on Atmospheric Deposition in East Asia 2012, Bangkok, Thailand, (2012) "Effects of ozone on plants"

11) 山口真弘、渡辺 誠、松村秀幸、河野吉久、伊豆田 猛：第123回日本森林学会大会テーマ別シンポジウム（2012）「葉のオゾン吸収量に基づくブナ苗に対するオゾンと土壌窒素負荷の複合影響の評価」

12) 松村秀幸、伊豆田 猛：第53回大気環境学会年会（2012）「同化箱を用いて測定した数種落葉広葉樹の葉のオゾン吸収と純光合成との関係」

13) 山口真弘、安土文鹿、松村友絵、上原唯、鹿又友彰、黄瀬佳之、小林亜由美、松村秀幸、伊豆田猛：第124回日本森林学会大会テーマ別シンポジウム（2013）「葉のオゾン吸収量に基づくブナ苗に対するオゾンと土壌窒素負荷の複合影響の評価」

14) 山口真弘、黄瀬佳之、安土文鹿、上原唯、小林亜由美、鹿又友彰、松村秀幸、伊豆田猛：第54回大気環境学会（2013）「ブナ、コナラ、ミズナラおよびシラカンバの葉の気孔を介したオゾン吸収量の推定」

15) 松村秀幸、山口真弘、黄瀬佳之、米倉哲志、伊豆田猛：第54回大気環境学会（2013）「落葉広葉樹4種の葉の純光合成速度におよぼす1成長期間のオゾン暴露の影響とオゾン吸收量との関係」

16) 安土文鹿、上原 唯、黄瀬佳之、鹿又友彰、小林亜由美、山口真弘、伊豆田 猛：第54回 大気環境学会（2013）「葉のオゾン吸収量に基づくブナ、ミズナラおよびシラカンバの純光合成速度におけるオゾン感受性の評価」

17) M. YAMAGUCHI and T. IZUTA: International Symposium on Agricultural Meteorology 2014, Sapporo, Japan, (2014) “Effects of air pollutants on photosynthetic activity of Asian forest tree species”
18) 山口真弘、安土文鹿、黄瀬佳之、小林亜由美、鹿又友彰、上原唯、松村秀幸、伊豆田猛：第125回日本森林学会大会テーマ別シンポジウム（2014）「葉のオゾン吸収量に基づくブナ、コナラ、ミズナラおよびシラカンバの葉の純光合成速度に対するオゾンの影響評価」

## （3） 出願特許

特に記載すべき事項はない

（4）シンポジウム、セミナー等の開催（主催のもの）

特に記載すべき事項はない

## （５）マスコミ等への公表・報道等

1）東京新聞（2013年7月28日、サンデー版、「光化学スモッグ」）

## (6) その他

「国民との科学・技術対話」（平成22年6月19日）を推進するために、平成23～25年度に合計17回の講演等を行い、本研究プロジェクトで得られた成果を一般市民や学生などに紹介した。

1）伊豆田猛：府中市エコサマースクール（2011）「府中市の大気はだいじょうぶ？」

2）伊豆田 猛：オゾン等による植物影響シンポジウム（北海道弟子屈町）（2011）「植物に対する越境大気汚染物質の影響」

3）伊豆田猛：東京大学森林科学セミナー（2012）「樹木に対する越境大気汚染物質の影響」

4) Izuta, T.: JICA Capacity Development (2012) "Effects of transboundary air pollutants on plants"

5）伊豆田猛：津軽半島環境サミット（2012）「植物に対する越境大気汚染物質の影響」

6）伊豆田猛：北海道大学講演会（2012）「対流圏オゾンが農林作物の成長に及ぼす影響の研究の動向」

7）伊豆田 猛：北海道立総合研究機構 環境科学研究センター 対流圏オブン等による植物影響講演会（2012）「越境大気汚染の影響，北海道の農作物と樹木への影響は？」

8) 伊豆田 猛：北海道立総合研究機構 環境科学研究センター 対流圏オゾン等による植物影響

検討会（2012）「オゾンの植物影響評価には何が必要か，暴露試験の変遷と装置の特徴」

9）伊豆田猛：横浜市立フロンティアサイエンス高等学校（2013）「大気汚染と植物」

10) 伊豆田 猛：府中市講演会（2013）「大気汚染と植物」

11) 伊豆田 猛：東京都立科学技術高等学校（2013）「大気汚染と植物」

12) 伊豆田 猛：慶熙大学校講演會（2013）“Effects of ozone on plants”

13) 伊豆田 猛：東京都立武蔵野北高等学校（2013）「大気汚染と植物」

14) 伊豆田 猛：神奈川県農業技術センター講演会（2013）「ブナに対するオゾンの影響」

15) 伊豆田 猛：北海道大学講演会（2013）「越境大気汚染と植物」

16) 伊豆田 猛：東京農工大学科学技術展（2013）「樹木に対するオゾンの影響に関する研究」

17) 伊豆田 猛：京都大学生存圏研究所国際シンポジウム（2013）“Effects of transboundary air pollutants on plants”