土壌の物理性 第 139 号 平成 30 年 7 月 20 日発行（年 3 回発行） 昭和 45 年 7 月 31 日 学術刊行物承認

ISSN 0387-6012

土壌の物理性Journal of the Japanese Society of Soil Physics

第 139 号2018 年 7 月

土壌物理学会Japanese Society of Soil Physics

土壌の物理性
第 139 号目2018 年 7 月

次

巻頭言

竹内晴信

...

1

亀山幸司 · 岩田幸良 · 宮本輝仁 · 佐々木康一

...

3

...

15

...

23

会務報告

........................................................................

25

編集後記

........................................................................

27

畑地灌漑特集論

文

HYDRUS-1D を用いた砂丘未熟土畑における灌水操作方法の検討

研究ノート
Handmade Becquerel meter using commercial scintillation survey meter
S. YOSHIDA, S. SHIOZAWA, T. ZAITSU, H. YAMANO
and K. NISHIDA

土粒子
土壌物理をめぐる研究環境

井上久義

表紙写真の説明
ロール状の鉛板を円筒状に加工して簡易な遮蔽容器を作成し，その内部に市販のサーベイメータ用のシンチレーションプロ−ブをセットして試料の放射能を測定できるようにした装置．試料をセットしたときと，しないときの γ 線カウント数の差から，試料の放射能を求めることができる．今号掲載の「Handmade
Becquerel meter using commercial scintillation survey meter」をご参照下さい．

第 16 回（ 2018 年度）
土壌物理学会賞（論文賞）候補の推薦（公募）について
土壌物理学会では，下記の要領で学会賞候補（推薦）を公募いたします．
学会賞種類：論文賞
対 象 論 文：2017（平成 29）年度に「土壌の物理性」（第 136，137，138 号）に掲載された論文（original paper）
推 薦 期 限：2018（平成 30）年 8 月 15 日（水）必着
推薦書に必要事項をご記入いただき，学会事務局（事務局長）までお送り下さい．推薦書様式は，
学会ホームページ

http://js-soilphysics.com/prz

下部の “論文賞推薦書” をダウンロードしてご記入願います．
表

彰：2018（平成 30）年 10 月 27 日（土）2018 年度大会にて

2018 年度土壌物理学会大会のご案内

日時：10 月 27 日（土） 9:00 ∼ 18:00

場

所：北海道大学農学部

参加費：3,000 円 （要旨集代として．ただし，学生会員は無料）

シンポジウムのプログラム，情報交換会，昼食弁当の申し込み方法等は，8 月 10 日までに土壌物理学会ホームページ（https://js-soilphysics.com/conf）上に掲載します．また，大会当日に，学生会員の入会手続を行った学生は，参加費が無料 となります．

1. 開会 · 事務連絡 農学部 4 階大講堂

9:00 ∼ 9:15

2. ポスター · セッション 農学部 N11，N21，N31 教室 9:30 ∼ 11:00
ポスターは，9:30 までに所定の位置に掲示してください．
発表者は，:::::::::::
学会員のみとします．土壌物理学会発表要領に基づいて作成した発表要旨（A4，2 ページ（200 字程度の研究紹介を含む），pdf 形式）を，学会ホームページ（https://js-soilphysics.com/conf）
へアップロードして下さい．8::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
月 20 日（2 ヶ月前）受け付け開始，9 月 27 日（1 ヶ月前）締切です．

3. 土壌物理学会総会

農学部 4 階大講堂

11:15 ∼ 12:15

総会では，土壌物理学会賞（論文賞）の授与を行います．

4．第 60 回シンポジウム

農学部 4 階大講堂

13:15 ∼ 18:00

テーマ「これからの持続的農業を考える；土壌中の物質循環 · 微生物 · 共生からの視点」
総合司会

柏木淳一

事務局長

趣旨説明

石黒宗秀

会長

北海道大学農学研究院

北海道大学農学研究院（1） 土壌生産力の支配因子を求めて水野直治

元酪農学園大学農獣医学研究科教授 · 東京農業大学 · 北海道立農業試験場（2） 近年における水田土壌の変化と持続的水稲生産に向けた対応
金田吉弘

秋田県立大学生物資源科学部（3） 直接分析から土壌中の微生物 – 元素 – 鉱物相互作用を調べる
光延 聖

愛媛大学農学部

（4） 植物共生科学から考える農耕地生態系の物質循環と持続的農業池田成志

農研機構北海道農業研究センター大規模畑作研究領域（5） 農薬や肥料に依存した現代稲作への警鐘 ∼ 江戸時代に開発された水田の多数回中耕除草法が意味するもの
粕渕辰昭（6） 総合討論

5. 企業展示

山形大学農学部名誉教授

司会

石黒宗秀

会長

北海道大学農学研究

竹内晴信

副会長

北海道立総合研究機構十勝農業試験場

農学部 N11, N21, N31 講義室

6. 情報交換会

18:15 ∼ 20:00

エンレイソウ（北大構内）

7. 大会会場への交通手段
JR 札幌駅から，徒歩で約 15 分．詳細は，北大農学部 HP のアクセスマップをご覧ください．
「https://www.agr.hokudai.ac.jp/i/access」
問い合わせ先:土壌物理学会事務局事務局長〒060-8589

柏木淳一

札幌市北区北 9 条西 9 丁目

北海道大学農学部

電話 011-706-3641 E-mail: kashi@env.agr.hokudai.ac.jp

8. 農業農村工学会土壌物理研究部会（土壌物理学会 共催）
日時：10 月 26 日（金）13:30 ∼ 17:00

場

所：北海道大学農学部 4 階大講堂

テーマ「土壌環境と気候変動（仮）」「気候変動が農地土壌の物理環境に及ぼす影響の予測（仮）」

弘前大学

加藤千尋

「大気 CO2 濃度の増加と温暖化が農耕地からの温室効果ガス排出におよぼす影響（仮）」
農研機構

常田岳士

北海道大学大学院工学研究院

山田朋人「気候変動の影響を考慮した今後の治水に関する取り組み（仮）」問合わせ先：
農業農村工学会土壌物理研究部会事務局〒840-8502 佐賀県佐賀市本庄町 1 佐賀大学農学部生物環境科学科 徳本家康
電話：0952-28-8755 E-mail：yasu@cc.saga-u.ac.jp http://www.jsidre.or.jp/dojou/

評議員選挙の実施について（予告）
土壌物理学会選挙管理委員
本年は，評議員の選挙の年です．現在の評議員は，平成 31 年 3 月 31 日に任期が満了しますので，会則第 6 条（2）
「評議員（イ）15 名正会員の中から選挙によって選出される．」および役員選出規定（1970 年 11 月 18 日改定）に基づいて，選挙を実施します．
有権者のみなさまには，9 月上旬に投票用紙等選挙関係書類をお送りしますので，ご投票下さい．ご協力をお願い申し上げます．

会費請求について
昨年度，事務局の手違いから皆様に会費請求を怠っており，多くの方々にご迷惑をお掛けしました．不手際がありましたこと，深くお詫び申し上げます．
未納の方におかれましては，本年度を含めて 2 年度分の会費を請求させて頂きますのでどうぞご容赦下さい．近日中に皆様のお手元に届きますよう，会費請求の準備を進めておりますので，請求書が届きましたら納付下さいますようお願い致します．
なお，郵便振込以外のお支払い方法を希望されます場合，また納付状況に対するお問合せにつきましては，お手数ですが，spkaikei@ml.affrc.go.jp までお申し出下さるようお願い致します．

J. Jpn. Soc. Soil Phys.
土壌の物理性No. 139, p.1 ∼ p.2 (2018)

農業生産の現場から考える土壌物理研究竹内晴信 1

大それたことを仰せつかってしまった
一昨年，学会事務局が北海道に移管された際，新学会長に就かれることとなった石黒先生より副会長の打診を戴いた．従前より，「学会は真理の探求を極め，またそれを活かした技術を実現し，それらのエッセンスを後世に伝える場である」と認識していたのであるが，いざその場に入って舵取り役の補佐を仰せつかるとは思ってもみなかった．
本学会は，「土壌の物理的現象に関連する科学と技術の研究領域を対象にする」ことを謳い，様々な専門分野の研究者，技術者達が学際領域の発展を夢見て設立されたもので，60 年に及ぶ歴史を持っている．日々，眼前の技術的サポートに振り回されてろくに論文も書いていない自分がいったい何をやって貢献すれば良いのやら，誠に困った事態になってしまった．
さて，学際領域ではあるけれども，土壌の物理的な性質や現象に関わる科学的知見を活用する場面として真っ先に考えられるのは，土木工学と農学 · 環境学の分野であろう．本学会を盛り立てておられる会員諸氏も多くが後者を軸足においておられる．なので，筆者の少ない経験の中から，農業生産の視点から見た土壌物理性に関わる技術を展望
してみることとしたい．
作物生産上の問題
古いデータで恐縮だが，1960 ∼ 70 年代に国主導で行われた地力保全基本調査事業により，国内の農耕地土壌
288.7 万 ha の生産力阻害要因の抽出，集計が行われた．これによると，相当大きな阻害要因を持つ不良土壌（第 III，
IV 等級）の面積は水田の 38.9 %，畑地では 69.2 % に達している．この等級区分は，以下に示す要因項目，すなわち，表土の厚さ，有効土層の深さ，表土の礫含量，耕耘の難易，※湛水透水性，※酸化還元性，＃土地の乾湿，自然肥沃度，養分の豊否，障害性，災害性，＃傾斜，＃侵蝕の各項目毎に I ∼ IV の評価を行い，最も低いカテゴリをもって当該土壌の生産力分級としたものである（※水田のみ，＃畑のみ）．
北海道の畑地についてみると，第 III，IV 等級となる要因項目として，自然肥沃度，養分の豊否といった化学的要因で不良を示しているものが 30 % を超えて最も多かったことから，物理性不良が突出していたわけではない．しかし，近年の資材多投型農業生産体系の下では，こうした養分含量の不足は大幅に改善され，むしろ過剰蓄積の弊害が顕在化している．
一方，1979 年から続く耕地土壌の実態を北海道内の広域モニタリング調査によりとりまとめた成果によると，水田，普通畑，草地とも近年（2010 年頃）の作土深はその 20 ∼ 30 年前と比較として差が認められない．心土のち密度，仮比重とも横ばいで推移しており，巷で言われる「近年土壌が薄く硬くなってきた」はデータの裏付けがない．
逆に，この 30 年間の作土の全炭素含量は，水田で横這い，普通畑で漸減（−1.7 %），草地で大幅増（+2.6 %）となっている．また，有効態リン，交換性カリウムなどの養分は調査した時代毎に大きな変化を示す例が見られ，基本的に増加傾向にある．易分解性の窒素は全炭素と同様に減少傾向であった．これらのことから，土壌が本来持つ物理的特性は長年維持されたままである一方で，炭素や NPK 等の養分量は営農活動の影響を大きく受けていると考えられる．
したがって，養分供給能の改善が図られた段階では，もともと不良な物理的要因の改善要望が表面化していくことが理解される．
土壌物理性の診断
それでは，土壌物理性の何をどの程度まで（どうやって）改善すれば良いのだろうか．
地力増進法（昭和 59 年法律第 34 号）では土壌の性質の改善目標が定められている．すなわち，普通畑では，作土の厚さ 25 cm 以上，主要根群域の最大ち密度が山中式硬度計で 22 mm 以下，粗孔隙量 10 vol. % 以上，0 – 40 cm の易
有効水分保持能 20 mm 以上，腐植含有量 3 % 以上，であり，これらの達成に向け基本的な改善方策も示されている．
これを受け，各都道府県においては地域の実情を加味して各々独自の土壌診断基準値を設定し，改善を図っている．
例えば北海道では，（水田で）すき床層の貫入抵抗値 0.5 ∼ 1.5 MPa，収穫期地耐力 0.25 MPa 以上，（畑地で）作土の1 地方独立行政法人北海道立総合研究機構

十勝農業試験場

土壌の物理性

2

第 139 号(2018)

砕土率 70 %，飽和透水係数（cm s−1 ) −3 ∼ −4 オーダ，耕盤層の判定指針 1.5 MPa 以上，などが独自の基準となっている．
蛇足ながら，ユニークな技術開発事例として，「有機農業を行うための診断基準値」を紹介しておきたい．道立農試（現 · 北海道立総合研究機構）では，有機物重点利用で露地野菜を栽培する際のほ場適性を，窒素肥沃度水準と土壌物
理性の両面から評価し，土壌の粘土含量，心土のち密度，腐植含量によって 5 段階に区分した．これによると，特別栽培でキャベツを栽培した場合に水準 I のほ場では，慣行栽培比で 90 以上の収量が期待され，水準 V のほ場では同
比 30 の収量しか見込めない，というものである．水準 V に相当するのは，例えば腐植 3 % 以下では土性 CL より細粒質でち密度 24 以上の条件である．
土壌物理性の改善に向けた取り組み
このように改善目標はある程度定量的に示されているのだが，化学性と違ってこれらの基準値を簡単に達成することは困難なことが多い．大規模な土木的工事を行わない限り，あるいは土壌の構成要素そのものを変えない限り，その抜本的な改良が行えないという問題もある．そうではあるが，政策的に公共事業が続けられてきたことにより，ほ場の地下水位は低下し，作土厚が確保されるなど，劣悪な条件の改善は進んだ．
営農対策としては，近年は大型トラクタの利用が一般化してきている中で，ほ場管理作業の中に心土破砕が組み込まれるようになった．またレーザーレベラーの普及で，ほ場均平化も進んでいる．さらに，適切な輪作体系や有機物の投入が行われていれば，土壌炭素の消耗や踏圧による土壌の圧密固化への影響は小さいものと考えている．
このように，外的，内的努力により，耕地土壌の物理性も改善が進んでいると理解されるが，広域的に見るとそれはまだ十分な水準ではない．また，不適条件（水分）でのトラクタ走行や耕起砕土によって土壌を練り返すことが微細孔隙の閉塞を招き透水性低下をもたらすなど，引き続き解決を要する問題点も未だ多く残されている．改めてそれらを集約するなら，①根域土層深の確保（拡大），②硬さ（堅さ，固さを含む）の改善，③透排水性の向上，④侵食対策，といったところになろう．
私たちはこれから何をすべきか（自問自答）
農地の大規模化，粗放化（= 低コスト管理）が推奨され，土木的工事だけが突出して進められるのなら，決して「良い土」の創出にはならないであろう．化学性や生物性とのバランスが必要なことは言うまでも無いが，環境への影響や土壌劣化に対して足を引っ張らないような改善対策を行わなければいけない．
今後の土壌改良の視点として，例えば一度行った改良のアクションがどの程度の持続性を持つものなのか検証した例は少ない．心土破砕で形成した破砕孔の変化や，練り返した土壌の物理的特性の時間的変化などを知りたい．もと
もと団粒化や亀裂の形成が見られない土壌で，粗大孔隙を低コストで短期間に形成する画期的技術はできないものか．
リモートセンシングにより土壌物理性を広域的に把握する取り組みは，これまでにも行われてきたが，近年の急速に進歩した ICT の活用で，その精度や簡便性は大幅に高まっている．しかし，何が作物の生育を規制しているのか，
単純ではないその関係を面的にひもとくことが重要と思っている．
近年勢いを増す気象災害による土壌侵食への対応も頭に浮かぶが，土壌の改善だけでは対処できることではないであろう．
ところで，現場においては，決して生産性向上のみを追い求めているわけではない．環境保全の見地からは，硝酸性窒素をはじめとした栄養塩類の土壌からの溶脱 · 水系への流出の定量評価が一通り行われた．生産管理技術としての硝酸性窒素流亡抑制は，主に窒素施用の適正化と堆肥等有機物の施用上限量の設定で対処され，降水量と土壌の保水性よりカテゴライズした土地の窒素環境容量も設定されている．地域の土地利用と水系の水質との関連性は多く検討されているが，それをどのように社会や生産技術に還元するかの提案が少ないように感じる．
近年は温室効果ガスフラックスの評価と排出抑制対策に関する研究も多く，土壌への炭素隔離の研究も進みつつある．この部分では，そもそも農地への有機物投入量の絶対的不足があって，農業生産の系内だけ見ても炭素の循環がうまくいっていないことの証左であり，今後は研究面というよりも環境保全型の生産技術の実践に対するインセン
ティブを高めていく等の対策が必要であろう．
例えば透排水性の向上と栄養塩類の流亡，土壌微生物の活性化と蓄積炭素の分解消耗に代表されるように，生産性と環境保全機能にトレードオフが生じた場合，何を優先するかの判断をサポートするような評価指標が必要になっているのかもしれない．
以上，思いつくままに並べて何やら総説みたいな文章になってしまったが，農業生産を念頭に置いた土壌物理の分野で，会員諸兄の研究ネタ探しのきっかけとなってくれればもっけの幸いである．

J. Jpn. Soc. Soil Phys.
土壌の物理性No. 139, p.3 ∼ 12 (2018)

HYDRUS-1D を用いた砂丘未熟土畑における灌水操作方法の検討
亀山幸司 1 · 岩田幸良 1 · 宮本輝仁 1 · 佐々木康一 2

Optimization of irrigation regime in a sand dune soil field using HYDRUS-1D
Koji KAMEYAMA1 , Yukiyoshi IWATA1 , Teruhito MIYAMOTO1 and Koichi SASAKI2

Abstract: In this study, we investigated the optimum irrigation regime for a sand dune soil field using a numericalsimulation model. For this purpose, we used HYDRUS1D to simulate the soil water regime in the root zone of anirrigated sand dune soil and evaluated the effects of irrigation scenario (matric potential threshold for irrigation anddaily irrigation depth) on the cumulative irrigation amount,cumulative water discharge from the root zone, ratio of cumulative irrigation amount over potential evapotranspiration, and ratio of transpiration reduction using the triggered
irrigation module in HYDRUS-1D. When a soil water
characteristic parameter determined from undisturbed soil
samples in the laboratory was used, volumetric water contents simulated by HYDRUS-1D differed from field measurements (RMSE = 0.031 − 0.039). Estimation accuracy
was improved using a soil water characteristic parameter
derived from an inverse analysis (RMSE = 0.007 − 0.009).
Based on the results of scenario simulation with different matric potential thresholds for irrigation and daily irrigation depths, the matric potential threshold set at −75
cm and the daily irrigation depth of 3.2 mm are recommended to decrease transpiration reduction and water discharge from the root zone.
Key Words : irrigation regime, numerical simulation, soil
water movement, triggered irrigation

1. はじめに三里浜砂丘地のような砂丘未熟土の圃場に灌漑施設が整備されると，野菜やメロン，スイカ等の灌水効果が高
い新規作物の導入が推奨される．これらの作物は一般に高い土壌水分状態で栽培される．例えば，砂丘地圃場でのスイカ栽培では，圧力水頭が −25 ∼ −32 cm の範囲で維持されていたとの報告例がある（橋本ら, 2010）．また，井上 · 竹内（2000）は砂地圃場でホウレンソウを栽培するときの最適な圧力水頭の範囲は −20 ∼ −30 cm であったと報告した．
高い土壌水分量を維持するためには頻繁に灌水を実施する必要があるが，これを人力で行おうとすると多大なる労力を要する．そのため，これまでに畑地灌漑の自動化が検討されてきた（例えば，井上 · 竹内, 2000）．灌水を自動的に行うためには，圃場で体積含水率や圧力水頭をモニタリングし，測定値がある閾値を超えたときに必要な量の灌水を行う必要がある．このような畑地灌漑の自動化は，安価な土壌水分計の普及により実用化されてきている（Blonquist et al., 2006; Vellidis et al., 2008）．
灌水を開始する際の土壌水分の閾値や灌水量，灌水時
間などの灌水パラメータは栽培作物や栽培環境によって異なるため，灌水の自動化を行うためには個々の圃場条件に応じた灌水パラメータを設定する必要がある．しかし，圃場試験による灌水パラメータの取得には労力と時間を要する．そこで近年，数値解析を応用して灌水パラメータを決定する試みが始まった．これまでに土福井県北部に位置する三里浜砂丘地では，畑地灌漑用
水として地下水を利用している．しかし，海に近いため，地下水が塩水化することで塩害や灌漑用水の不足を引き起こしている（勝田, 2007）．この対策として，国営九頭竜川下流農業水利事業の中で，農業用水の再編により新たに生み出された用水が三里浜砂丘地に供給される予定である．この水源転換による将来の安定した灌漑用水供給を想定し，三里浜砂丘地では新規導入作物の検討や灌壌水分 · 溶質移動汎用プログラムが整備され，畝間灌漑
や点滴灌漑，地中灌漑等の異なる灌水方法について，土壌水分や溶質の動態が把握された（Abbasi et al., 2004;

Bastiaanssen et al., 2007; Subbaiah, 2013; El-Nesr et al.,2014）．日本においても，畑地灌漑の用水諸元の一つである日消費水量を精度良く求めるために HYDRUS が使用されている（中村ら, 2016; 弓削 · 阿南, 2016）．その結果，様々な灌水方法に汎用プログラムを適用する際の注

水効果の検証，肥料の溶脱状況の検証などの営農実証試

意点や，数値解析を行うときの境界条件の設定方法が明験が行われている（河原, 2013）．

らかになってきた．さらに，HYDRUS には，土層内のあ1 Institute for Rural Engineering, NARO, 2-1-6 Kannondai, Tsukuba,

る深さの圧力水頭が閾値に達した時に予め設定された水Ibaraki 305-8609, Japan. Corresponding author: 亀山幸司, 農研機構

量と時間によって間断灌漑を行うオプション（triggered

農村工学研究部門．

irrigation）がある（Dabach et al., 2013）．Marković et al.
（2015）は HYDRUS に組み込まれた間断灌漑オプション

2 Regional Agriculture Division, Fukui Prefectural Government, 3-17-1

Ote, Fukui city, Fukui 910-8580, Japan.
2017 年 6 月 13 日受稿 2018 年 4 月 2 日受理を使用して様々な栽培条件でシナリオ解析を行い，最適

土壌の物理性

4

第 139 号(2018)な灌水パラメータの選定が可能となることを示した．ま土であることや地下水面がないことを確認した．このた，Müller et al.（2016）は半乾燥地でナスを栽培する際，際，地表面から 7.5 cm，22.5 cm，45 cm の深さにおい

点滴灌漑を行うための最適な灌水パラメータを圃場実験て 100 cm3 の円筒試料を 3 点ずつ採取し，これらの試と数値解析により求めた．しかし，三里浜砂丘地のよう料の飽和透水係数および水分特性曲線を測定した．飽和に保水性の極端に低い砂丘未熟土において，HYDRUS を

透水係数は定水位法により測定した．また，水分特性曲

用いて灌水パラメータを決定した事例はほとんどない．

線は，−31.6 cm（pF 1.5）以上では砂柱法，−31.6 cm ∼砂丘未熟土の圃場において比較的高水分状態を維持す

灌水することにより根群域下端からの浸透損失を抑え−15, 800 cm（pF 1.5 ∼ 4.2）では加圧板法により測定した．各土層の水分特性曲線の測定結果に van GenuchtenMualem モデルを適合させるため，非線形回帰プログ
ラムの RETC（PC-Progress s.r.o.; van Genuchten et al.,1991）を用いた．

られることは知られている（Hillel, 1998）．ただし，砂室内試験により得られた水分特性曲線を Fig. 1 に示丘未熟土畑における少量頻繁灌水について，下方浸透し，水分特性曲線から得られた van Genuchten-Mualem

を極力抑えるための灌水パラメータを検討した事例はモデルパラメータを Table 1 に示した．また，室内試験

るように灌水を実施する場合，下方浸透が卓越すること
による灌水の損失や肥料の溶脱が懸念される（亀山ら,2016）．作物によって消費された水分を少量ずつ頻繁に少ない．一方，近年の研究では HYDRUS を用いて根群により測定された飽和透水係数も Table 1 に示した．水

域下端の水フラックスを計算し，畑地灌漑計画に役立て

分特性曲線については深さによる大きな変動は見られなることも検討されている（成岡ら, 2015; 岩田ら, 2016）．

かったが，飽和透水係数については深さによる変動が見

HYDRUS を用いて様々な灌水条件のもとで下方浸透量

られた．

を計算することで，砂丘未熟土圃場における灌水の浸透損失や肥料の溶脱を抑制するような灌水パラメータの決定が可能になると考えられる．

2.3 土壌水分量 · 灌水量 · 気象要素の測定
栽培期間中の土層内の土壌水分量，ハウス内の灌水量，日射量および温湿度を測定した．播種後の 2014 年

そのため，まず，現地の砂丘未熟土畑で生産者によっ10 月 21 日に，深さ 7.5 cm，15 cm，22.5 cm に土壌水て慣行的に行われている土壌水分を比較的高めに維持す

分センサー（Decagon 5TE）を設置した．センサーの設る少量頻繁灌水の下での土壌水分動態を HYDRUS によ

置深さは，コカブの貯蔵根長（15 cm 程度）を考慮してり再現する．次に，同様のモデルにより灌水パラメータ決定した．なお，試坑を掘った後，試坑の鉛直断面に対を変化させたときのシナリオ解析を行い，作物の消費水

して水平方向に各センサーの回路部を含む全体を挿入し量を満たしながら下方への浸透損失をできるだけ抑えるた．また，土壌水分センサーは 3 本のロッドが鉛直方向

灌水操作方法を検討する．

に並ぶように断面に挿入した．土壌水分センサーの設置
場所は，ハウスの側面から 1 m 程度離れたコカブの列間

2. 調査方法の中央である．土壌水分量の測定は，ハウス内の 2 地点で実施した．土壌水分センサーの出力値は，試験圃場か

2.1 試験圃場と栽培状況ら採取した土壌を用いて実験室で求められた校正式によ試験圃場は，三里浜砂丘地内（福井県坂井市三国町）にり体積含水率に変換した（亀山ら, 2016）．灌水量は，ハ位置するビニールハウス温室（幅 5.5 m × 長さ 85 m（室ウスの側面から 1 m 程度離れた位置に雨量計（Decagon

内面積，約 5 a）
）である．土壌は砂丘未熟土に分類され，ECRN-100）を設置して測定した．コカブが成熟した際砂分を 90 % 以上含む砂土である（亀山ら, 2016）．このにも葉が雨量計の測定を邪魔しないように，雨量計の開温室において，新規導入作物の検討のための実証試験が口部が地表面から 30 cm 高い位置になるように設置し実施され，前々作にメロン，前作にコカブが栽培された．

た．また，全天日射計（Decagon PYR）および温湿度セ

前作終了後，作土をロータリーにより耕耘した．2014

ンサー（Decagon VP-3）は，ハウス内において側面から

年 10 月 17 日に基肥として化成肥料（日の本化成 3 号）
を約 0.1 kg m−2 散布した後，コカブを播種した．コカブの列間は約 20 cm であり，株間は約 10 cm であった．栽
培期間中，ビニールハウスの左右に固定された片側噴霧式の灌水ホースで散水灌漑を行った．なお，その時の地表面の水分状態に基づいて栽培農家が適宜判断し，灌水を行うタイミングや灌水量を決定した．播種から 50 日
目以降，出荷可能な大きさ（直径が 5 cm 程度）に育った
株から順に収穫が行われた．コカブの栽培期間は, 2014年 10 月 19 日 ∼ 2015 年 1 月 12 日であった．

2.2 土壌の物理性の測定
試験に先立って同ハウス内に深さ 100 cm 程度の試坑を掘り，土壌断面調査を実施し，深さ 100 cm までは砂Table 1 室内試験により測定された水分特性曲線から得られた van Genuchten-Mualem モデル（（2）式）パラメータ（θr , θs , αVG , n）及び飽和透水係数（Ks ）測定値．
van Genuchten-Mualem model (equation (2)) parameter (θr ,θs , αVG ，n) determined from soil water characteristic curves
and saturated hydraulic conductivity (Ks ) measured fromundisturbed soil samples in the laboratory.

θr
θs
αVG
3
−3
3
−3
（cm）（m m ）（m m ）（cm−1 ）
Depth

Ks

n

（cm d−1 ）

7.5

0.07

0.44

0.06

3.57

2,012

22.5

0.07

0.44

0.07

3.50

1,524

45

0.06

0.43

0.07

3.39

2,586

論文：HYDRUS-1D を用いた砂丘未熟土畑における灌水操作方法の検討

5

ここで，θ ：体積含水率（m3 m−3 ），t ：時間（d），z：鉛
直深さ（cm），h：圧力水頭（cm），K(h)：不飽和透水係数（cm d−1 ），S(h)：根の吸水速度（d−1 ）である．
不飽和透水係数および水分特性曲線は van Genuchten-

Mualem モデル（van Genuchten，1980）により近似した．

θ = θr +θs − θr[1 + |αVG h|n ]m

)m ]2
[
(
1
K(h) = Ks Se l 1 − 1 − Sem

Se =Fig. 1 室内試験と逆解析から得られた水分特性曲線．
Soil water retention curves obtained from laboratory experiments and inverse analysis.

θ − θrθs − θr

(2)

(3)

(4)

ここで，θr ：残存体積含水率（m3 m−3 ），θs ：飽和体積

含水率（m3 m−3 ），αVG ，n，m = 1 − 1/n：曲線形状パラ
メータ，Ks：飽和透水係数（cm d−1 ），Se：有効水分飽和度（–），l ：間隙結合係数（0.5）である．

3.2 根の吸水速度と可能蒸発散速度の設定（1）式中の根の吸水速度 S(h) は次式により計算される．

S(h) = α (h)β (z)Tp(5)

ここで，Tp は可能蒸散速度（cm d−1 ），α (h) は吸水速度の制限因子となる圧力水頭による水ストレス応答関数，

Fig. 2 圧力水頭と根の吸水抑制の関係を表す S-shape モ
デル（
（6）式）
（図中の h50 は，根の吸水速度が可能蒸散速度の 50 % に抑制される時の圧力水頭を示す）．
Plot of the S-shaped function（equation（6）） for describing

β (z) は正規化された根密度分布（cm−1 ）であり，根群域を深さ方向に積分すると 1 になる．
HYDRUS-1D では，（5）式中の水ストレス応答関数
α (h) として Feddes モデル（Feddes et al., 1978）と Sshape モデル（van Genuchten, 1987）を選択することができる．ここでは，少ないパラメータ数で α (h) が計算できる S-shape モデルを用いた（Fig. 2）．

root water uptake reduction with different matric potential.
h50 represents a matric potential at which the root water uptake rate is reduced by 50 % from the potential rate.

1( )pα (h) =
1+

h
h50

(6)

0.5 m 程度離し，高さ 1 m の位置に設置した．全ての測定データを，30 分間隔でロガー（Decagon Em50b）に記制される時の圧力水頭（cm）であり， p は曲線形状パラ録した．

メータである．

ここで，h50 は根の吸水速度が可能蒸散速度の 50 % に抑（5）式中の Tp は可能蒸発散速度（ETp ）をもとに与え3. 数値解析方法られる．そこで，ハウス内の日射量，日最高気温，日最低気温，平均相対湿度の測定値を用い，Penman-Monteith

3.1 土壌水分移動に関する支配方程式

式（Allen et al., 1998）により ETp を算定した．その際試験圃場では，畝などの顕著な凹凸はなく，散水灌漑に必要となるハウス内の日平均風速は，農林水産省農村により土壌面への水分補給が空間的に一様に行なわれ振興局（2016）に倣い，屋外の日平均風速の 15 % と仮ていたため，土壌水分動態の数値解析を HYDRUS-1D

定した．なお，屋外の風速データとして圃場から直線距

ver. 4.16.0110（Šimůnek et al., 2015）を用いて行った．
HYDRUS-1D では土層内の鉛直 1 次元の土壌水分移動

離で約 5 km の位置にある三国観測所の観測値を援用しが次式により計算される．

式により可能蒸発速度（Ep ，cm d−1 ）と Tp に配分した．

(
)
∂θ
∂
∂h
=
K(h) + K(h) − S(h)
∂t
∂z
∂z

(1)

た．算定された ETp を，葉面積指数（LAI ）をもとに次Ep = ETp exp(−k × LAI)

(7)

土壌の物理性

6

第 139 号(2018)

均一に分布していると仮定し，根の吸水割合を表す（5）Tp = ETp [1 − exp(−k × LAI)](8)

式の β (z) を地表面から深さ 15 cm で同じ値（1/15）に設定した．

ここで，k は作物群落による放射の減衰に関わる係数で既存の研究では，根の給水に係る S-shape モデル（
（6）

あり，デフォルト値である 0.463（Šimůnek et al., 2015）

土壌断面調査により，100 cm 深さまでは砂土である式）のパラメータの h50 は −100 ∼ −5, 000 cm， p は 2∼ 3 が用いられている（Dabach et al., 2013; Dabach etal., 2015; Müller et al., 2016; Segal et al., 2006; Skaggs
et al., 2006）．砂土の水分保持は毛管水領域で急激に低ことが確認されたため，計算領域として深さ 100 cm の下し，それ以下の領域では殆どゼロとなる．このため，土層を設定した．土壌断面調査では，水分特性曲線につSkaggs et al.（2006）が指摘しているように，砂土の h50

いては深さによる大きな変動は見られなかったが（Fig.

は一般に大きくなることが考えられる．現地土壌の水分1），飽和透水係数については深さによる変動が見られた（Table 1）．このため，3 層の成層を仮定し，地表面から深さ 15 cm までを第 1 層，深さ 15 cm から 30 cm までを
第 2 層，30 cm 以深を第 3 層とした．また，計算期間は，
土壌水分センサーを設置した翌日の 2014 年 10 月 22 日
から栽培期間最終日の 2015 年 1 月 12 日までとした．

特性を見たとき，各層共に −10 cm から水分が急激に低計算対象期間の最初にあたる 2014 年 10 月 22 日 0 時る（Segal et al., 2006）．これらのパラメータを用いる場

を用いた．

3.3 計算領域および初期 · 境界条件

下し，−100 cm 以下では殆ど水分が保持されていないことが確認され（Fig. 1），−100 cm 以下では根の吸水がかなり抑制されると考えられる．そこで，今回の計算では

h50 = −100 cm，p = 3 を用いた．なお，同様のパラメー
タが砂質土における観葉植物栽培に対して適用されていの時点における体積含水率の測定値は，深さ 7.5 cm で合，−20 cm 以下から吸水抑制が始まり，−160 cm 以下

0.154 m3 m−3 ，深さ 15 cm で 0.147 m3 m−3 ，深さ 22.5
cm で 0.146 m3 m−3 であり（いずれも 2 地点の測定値のでは 80 % 以上吸水が抑制される計算となる（Fig. 2）．

3.6 蒸発散量推定に係るパラメータの設定平均），深さ方向に大きな違いは見られなかった．そこ

可能蒸発散量を推定するため，コカブの葉面積指数

で，深さ 100 cm までの全層の初期体積含水率を一定と（LAI ）は以下のように設定した．収穫期（播種から 50

し，初期条件として 0.15 m3 m−3 を与えた．上部境界条日目）における 1 株あたりのコカブの葉面積は大和ら

件は大気境界条件とした．また，竹内 · 長江（1990）によ（1992）を参考に 0.05 m2 と仮定した．植栽密度（列間る三里浜砂丘地での測定事例によれば，地下水位は 4.5m 程度であったと報告されている．砂土では一般に毛管20 cm，株間 10 cm）から 1 m2 あたりの株数を 50 本と仮定すると，収穫期における 1 m2 あたりの LAI は 2.5上昇が生じにくく，毛管水帯が浅い．このため，地下水となる．それ以降は，出荷可能な大きさになった株から面が深さ 100 cm 付近の水移動に与える影響は極めて小順に収穫されるため，LAI は減少する．このため，播種さいと考えられたため，下部境界条件は自由排水条件と日の LAI は 0，播種後 50 日目の LAI は 2.5，栽培期間終した．

了日（2015 年 1 月 12 日）の LAI は 0 と仮定し，その間

3.4 水分特性曲線パラメータの設定の期間については線形補間した LAI を使用して，蒸発速

各土層の水分特性曲線の測定結果に RETC を適用し度と蒸散速度の割合を計算した．

て推定した van Genuchten-Mualem モデルのパラメータ（Table 1）を用いて数値解析を行った．計算結果と体積

4. 現地調査結果と土壌水分動態の

含水率の実測値を比較し，大きな乖離が見られる場合は，再現性の検討

逆解析による van Genuchten-Mualem モデルのパラメー
タの最適化を先のパラメータを初期値として実施した．
逆解析を行う際，体積含水率の実測値と計算値の Root

4.1 現地調査結果
栽培期間中の日灌水量，日射量，最高 · 最低温度，相対mean square error（RMSE ）ができる限り小さくなるよう

湿度，可能蒸発量，可能蒸散量，可能蒸発散量を Fig. 3に最適化を行った．

に示す．栽培初期から 11 月までは頻繁に灌水が行われ√たが，地上部の生長が止まる 12 月以降はあまり灌水が行
2(Xi − xi )n

われなかった（Fig. 3（a））
．頻繁に灌水が行われた 10 月
22 日から 11 月 30 日までの積算灌水量は 41.8 mm であり，日灌水量は 0.2 ∼ 5.4 mm d−1（平均 1.7 mm d−1 ）とここで，Xi は実測値，xi は計算値，n はデータ数を示す． 日によって大きく異なった．また，灌水の間隔は 0 ∼ 3
なお，逆解析のための実測値として，深さ 7.5，15，22.5
日と不規則であった．日射量や気温の変動（Figs. 3（b），
cm の体積含水率の日平均値（2 地点の測定値の平均）を （c））によって農家が灌水のタイミングを変えたことが，RMSE =

∑i=1
n

(9)

与えた．

3.5 根の吸水に係るパラメータの設定
現地圃場ではコカブの貯蔵根および側根は深さ約 15
cm 以内に観察された．そこで，根は深さ 15 cm までに日灌水量や灌水間隔に影響を与えたと考えられる．

10 月と 11 月の平均日射量は 16.6 MJ m−2 d−1 と値が大きく，ハウス内の日最高気温も平均 26.4 ◦ C と高い傾向があったが，12 月以降は平均日射量が 8.5 MJ m−2 d−1

論文：HYDRUS-1D を用いた砂丘未熟土畑における灌水操作方法の検討

7

Fig. 3 試験圃場の（a）日灌水量，（b）日射量，（c）最高 · 最低温度，（d）相対湿度の測定値および（e）可能蒸発量，可能蒸散量，可能蒸発散量の推定値．

(a) Daily irrigation depth, (b) solar radiation, (c) maximum and minimum temperatures, (d) relative humidity and(e) estimated potential evaporation, transpiration and evapotranspiration in an experimental field.

と減少し，ハウス内の日最高気温も平均 14.9 ◦ C と低下した（Figs. 3（b），（c））．

（b），（c））．
4.2 土壌水分動態の再現性

11 月までの ETP は 0.7 ∼ 2.7 mm d−1 （平均 1.6 mm−1d ）であったが（Fig. 3（e）），11 月までは頻繁に灌水がを用いて計算した体積含水率と実測値との比較を Fig.

行われていたため（Fig. 3（a）），土壌の体積含水率は全4 に示した．室内実験から得られた水分特性曲線パラての深さで 0.12 ∼ 0.17 m3 m−3 で推移していた（Figs. 4

メータを用いた場合，計算値と実測値は乖離しており室内実験により決定された水分特性曲線パラメータ

（a）
，
（b）
，
（c）
）
．水分特性曲線（Fig. 1）から圧力水頭は（RMSE = 0.031 ∼ 0.039），土壌水分量の変化幅も大きく−100 cm（pF 2.0）よりも湿潤側に相当する水分状態で異なった（Figs. 4（a）
，
（b）
，
（c）
）
．特に灌水により体積あったと判断されることから，圃場容水量よりも多い水

含水率が高い状況にあった栽培初期から 11 月にかけて分状態であったと推察される．12 月以降は灌水が殆ど行

計算値が低くなった．

われず，ETP は 0.3 ∼ 2.6 mm d−1（平均 1.1 mm d−1 ）と原因を探索するため，蒸発散量のない条件で計算した蒸発散が無視できなかった（Fig. 3（e））．そのため，全ての深さで体積含水率は緩やかに減少した（Figs. 4（a），

が，計算結果は大きく変わらず，RMSE は殆ど変わらなかった．このため，蒸発散量の推定誤差が 11 月までの

土壌の物理性

8

第 139 号(2018)

Fig. 4 室内パラメータと逆解析パラメータを用いた場合の体積含水率の計算値ならびに 2 地点の測定値の経時変化（a）深さ 7.5 cm，（b）深さ 15 cm，（b）深さ 22.5 cm．
Time series of simulated volumetric water contents and measured volumetric water contents at two measurement positions.
Simulated values were calculated using HYDRUS-1D with soil water characteristic parameter obtained from laboratory
data (solid line) and inverse analysis (dashed line). (a) depth 7.5 cm, (b) depth 15 cm and (c) depth 22.5 cm.

実測値と計算値の違いの原因ではないと考えられた．さらに，初期条件の影響についても検討したが，初期体積2.14（15 – 30 cm），3.63（30 cm 以深）であった．その結果を Fig. 1 に示すと，表層ほど推定値は実測値から外含水率を変化させても同様に計算結果への影響は殆どれ，推定された水分特性曲線は実測した曲線よりも勾配見られなかった．そのため，室内試験による水分特性曲が緩やかであることがわかる．現地の土壌水分動態を再線パラメータが実際の試験圃場の水分特性曲線を十分に現するためには，表層における重力排水後の土壌水分量

反映していないことが，体積含水率の計算結果が測定結が 0.15 m3 m−3 程度となるよう水分特性曲線を修正する

果と乖離する主要因であると判断した．室内実験による必要があった．この要因として現地土壌の不均一性に加水分特性曲線パラメータは迅速かつ正確に得ることが可

え，採取試料と現地土壌で土壌構造や水分分布，封入空

能である一方，圃場条件における水分特性曲線の代表値

気量等が異なったこと（中村ら, 2016）が考えられる．

とならない場合もあることが指摘されている（例えば，Sonnleitner et al., 2003; 中村ら, 2016）．

5. シナリオ解析による灌水操作方法の検討

逆解析では，推定するパラメータをできる限り減らす
ことが推奨されている（Rassam et al., 2004）．これまでに

5.1 灌水パラメータの設定

報告されている砂土の水分特性曲線（例えば，Nakamura
et al., 2004; 坂井 · 取出, 2007; 藤巻ら, 2012）を比較した

初期の 10 月 22 日から 11 月 30 日までを対象に，作物

ところ，曲線勾配は幅を持つことが明らかとなった．その消費水量を満たしながら下方への浸透損失をできるだ

こで，van Genuchten-Mualem モデルのパラメータのう

け抑える灌水操作方法を検討した．HYDRUS-1D に組

ち曲線勾配に関わる n を推定パラメータとして，体積含み込まれた間断灌漑オプションを使用すると，灌水を開水率の実測値と計算値の RMSE ができる限り小さくな

始する際の圧力水頭の閾値とそれを設定する深さ，およるように逆解析を行った．

び 1 回あたりの灌水量（灌水強度と 1 回あたりの灌水時

コカブの栽培期間のうち，頻繁に灌水が行われた栽培

推定パラメータを用いることにより RMSE は 0.007 ∼間）を設定して土壌水分動態の数値解析ができる．本研

0.009 となり，計算結果は大きく改善された（Figs.4（a），

究では，間断灌漑オプションを用いたシナリオ解析によ（b），（c））．逆解析から得られた n は 1.99（0 – 15 cm），

り，灌水を開始する際の圧力水頭の閾値や日灌水量を変論文：HYDRUS-1D を用いた砂丘未熟土畑における灌水操作方法の検討

9化させた場合，計算を実施した期間である 40 日間の積の大部分が深部に浸透しやすい土層条件が形成されてし算灌水量や深さ 100 cm における積算下方浸透量がどのまうことが要因として考えられる．

ように変化するかを評価した．

灌水を開始する際の圧力水頭の閾値と積算流出量の関灌水を開始する際の圧力水頭の閾値は −250，−200，

係をみると，閾値が −100 cm 以下では，日灌水量に依−150，−100，−75，−50，−25 cm とした．これらの値は，らずほぼ一定の積算流出量となった．これに対して，閾

黒ボク土圃場のコカブ栽培では圧力水頭 −100 ∼ −1, 000

値が −100 cm 以上になると下端からの浸透損失量が増cm（pF 2 ∼ 3）を目安に灌水を行っているとする報告（高えはじめ，−50 cm 以上で積算流出量もほぼ頭打ち状態

野 · 猪野, 2008）や，砂土圃場では比較的高水分で灌水管となった（Fig. 5（b））．しかし，日灌水量が 4.8 mm の理されていることが多いとする報告（井上 · 竹内, 2000;

ときは閾値を高くすると積算流出量も増加した．Fig. 5

橋本ら, 2010）を参考に設定した．圧力水頭の閾値を設（a）と Fig. 5（b）を合わせてみると，閾値が −100 cm

定する深さは根群域（15 cm）の中心の 7.5 cm とした．

以上のときには，灌水量が増加すると積算流出量も灌水

現地の日灌水量と ETp （Fig. 3）を参考に，シナリオ量の増加分とほぼ同じ量だけ増加することから，灌水の解析で用いる日灌水量を設定した．現地調査結果から，ほとんどは作物に有効利用されないと考えられる．この10 月 22 日から 11 月 30 日までの日灌水量は 0.2 ∼ 5.4
mm d−1 （平均 1.7 mm d−1 ）であり，ETp は 0.7 ∼ 2.7
mm d−1 （平均 1.6 mm d−1 ）であった．このため，シナリオ解析で用いる日灌水量は 1.6（1 × ETp ），2.4（1.5 ×
ETp ），3.2（2 × ETp ），4.8（3 × ETp ）mm とした．また，その際の灌水強度は，現地データを参考に 10 cm d−1 とした．現場では灌水間隔が 1 日以上だったことから，灌水開始から 24 時間経過までは灌水を実施せず，これ以降で深さ 7.5 cm の圧力水頭が上記の閾値よりも低下しように土層を高水分状態に管理しようとすると，深部へた時点で灌水を実施する設定にした．灌水量は，灌水時間を変化させることで調整した．
初期条件は，現地の土壌水分動態を再現する数値解析と同様に，深さ 100 cm までの全層の体積含水率を 0.15
m3 m−3 とし，上部境界条件は大気境界条件，下部境界
条件は自由排水条件とした．
数値解析結果をもとに，積算 ETP に対する積算灌水量の割合（I · ETP−1 ）と，積算 TP に対する積算 TP と根によの浸透水量も増加する．そのため，作物の吸水を阻害しない程度に灌水を開始する際の圧力水頭の閾値を下げて灌水を実施することが，下方への浸透損失を抑制するためには望ましいと考えられる．
灌水を開始する際の圧力水頭の閾値および日灌水量とI · ETP−1 の関係から（Fig. 5（c）），閾値が −50 cm 以上か
つ日灌水量が 2.4 mm 以上の場合，I · ETP−1 は 1.5 以上の過剰灌水となり（Fig. 5 （c））
，浸透損失量が増加する可能性が示された（Fig. 5（b））．灌水を開始する際の圧力
水頭の閾値が −100 cm 以下の場合には，I · ETP−1 は 1 以
下で灌水不足となる（Fig. 5（c））．I · ETP−1 が 1 となる
圧力水頭の閾値をグラフから求めると，日灌水量が 1.6mm，2.4 mm，3.2 mm，4.8 mm のとき，それぞれ −50
cm，−65 cm，−75 cm，−85 cm であった．灌水を開始する際の圧力水頭の閾値を下げて日灌水量を増やすこと
により，I と ETP のバランスがとられることがわかる．

る積算吸水量の差の割合（蒸散抑制率）を算定し，灌水操灌水を開始する際の圧力水頭の閾値および日灌水量と作の結果の評価指標とした．I · ETP−1 は可能蒸発散量に蒸散抑制率の関係（Fig. 5（d））から，閾値が −50 cm 以

対する灌水量の割合であることから，作物の消費水量に上の場合は蒸散抑制率が低く，−75 cm 以下の場合は閾

対する灌水量の妥当性を評価できる．一般に，塩類集積値の低下に伴って蒸散抑制率が上昇する結果となった．

を発現させないためにも灌水量は蒸発散量よりも僅かにまた，閾値が −50 cm 以上の場合は日灌水量が多いほど多いことが望ましいとされている（Dabach et al., 2013）．

蒸散抑制率が低下する傾向が見られた．閾値が −75 cm

蒸散抑制率は，土壌水分低下によって生じる根の吸水抑以下では日灌水量が 1.6 ∼ 3.2 mm の場合は同様の蒸散制割合を表しており，蒸散抑制率が大きくなると収量低抑制率を示し，日灌水量が 4.8 mm では他の日灌水量の下が懸念される（例えば，Segal et al., 2006）．

場合よりも蒸散抑制率が 2 ∼ 7 % 低くなる傾向が見られ

5.2 試験圃場における灌水操作の最適化

た．水ストレス応答関数（Fig. 2）を反映して，閾値が灌水を開始する際の圧力水頭の閾値や日灌水量を変化−100 cm の時の蒸散抑制率が約 0.5，−160 cm の時の蒸散抑制率が約 0.8 と計算された．今回の計算では，水分特性曲線を考慮して（Fig. 1），水ストレス応答関数が設させた場合の積算灌水量，下端（深さ 100 cm）からの積
算流出量，I · ETP−1 ，蒸散抑制率を Fig. 5 に示す．
灌水を開始する際の圧力水頭の閾値が高くなるにつれ定されている．しかし，灌水パラメータと蒸散抑制率のて，積算灌水量も増加する（Fig. 5（a））．また，日灌水

関係を今後より詳細に検討するためには，圧力水頭と作量が多いほど積算灌水量も多くなることもわかる．灌水

物根による吸水抑制との関係を実験的に確認することもを開始する際の圧力水頭の閾値が −50 cm 以上では，日重要と考えられる．

灌水量によらず，積算灌水量がほぼ頭打ち状態となった．

砂丘未熟土においては，灌水を開始する際の圧力水頭

これは灌水を行っても設定した圧力水頭まで回復しないの閾値を −50 cm 以上に設定した場合，土壌水分量が高まま，設定した量の灌水が連日行われたためと考えられく（Fig. 1），灌水の大部分が深部に浸透しやすくなるこる．砂丘未熟土では灌水を開始する際の圧力水頭の閾値とが示唆された（Fig. 5（b））
．一方，閾値を −100 cm 以を高く設定した場合，不飽和透水係数も高くなり，灌水下に設定した場合，土壌水分量が極端に小さくなり（Fig.

土壌の物理性

10

第 139 号(2018)(a) Cumulative irrigation amount (I)(b) Cumulative discharge from the bottom(c) Ratio of cumulative irrigation amount over(d) Transpiration reduction ratio

potential evapotranspiration (I × ETP −1 )Fig. 5 日灌水量と灌水を開始する際の圧力水頭の閾値を変化させた場合の（a）積算灌水量（I ），（b）下端（100cm）からの積算流出量，（c）可能蒸発散量に対する積算灌水量の割合（I × ETP−1 ），（d）蒸散抑制率の計算結果．
(a) Cumulative irrigation amount (I), (b) Cumulative discharge from the bottom, (c) Ratio of cumulative irrigation
amount over potential evapotranspiration (I × ETP−1 ), and (d) Transpiration reduction ratio, simulated by HYDRUS-1D
with different daily irrigation depths and thresholds of matric potential for irrigation.

1），作物根による水分吸収が抑制されるようになると推定された（Fig. 5（d））．この結果は，砂丘未熟土の畑に

これまで砂丘畑で報告されているものより低い値となっおいては，粗孔隙が多いことに由来する特殊な水分特性

べて少ないことがあげられる．Dabach et al.（2013）も曲線により，灌水を開始する際の圧力水頭の閾値の最適

ETP の増加に伴い灌水を開始する際の圧力水頭の閾値が範囲が狭くなる可能性を示唆している．

高くなることを数値解析結果から示している．

た原因として，今回の栽培期間における ETP が夏場に比Fig. 5（c）から，I · ETP−1 = 1 を満たすような日灌水量

6. おわりにと圧力水頭の閾値組み合わせは 1.6 mm と −50 cm，2.4

mm と −65 cm ，3.2 mm と −75 cm，4.8 mm と −85 cmと判断された．Fig. 5（d）に示したように，これらの組

砂丘未熟土畑において作物の消費水量を満たしながらみ合わせではいずれも蒸散が抑制されるが，閾値が −75
cm の場合，日灌水量が 3.2 mm 以下では灌水量の増加分が積算流出量の増加に寄与する割合が小さい（Fig. 5
（a）
，
（b）
）
．そのため，10 月から 11 月の期間のコカブ栽
培では，深さ 7.5 cm における灌水を開始する際の圧力
水頭の閾値を −75 cm，日灌水量 3.2 mm に設定することにより，消費水量と灌水量のバランスが良く（Fig. 5（c）），作物根の吸水抑制が少なく（Fig. 5（d）），浸透損失をできるだけ抑えた（Fig. 5（b）
）灌水を行うことがで下方への浸透損失をできるだけ抑える灌水操作方法をきると判断した．灌水を開始する際の圧力水頭の閾値が反映していないことが明らかになった．そこで，逆解析検討するため，HYDRUS-1D により灌水制御下の根群域の土壌水分動態予測を行うと共に，灌水を開始する際の圧力水頭の閾値と日灌水量を変えてシナリオ解析を行った．
室内実験から得られた水分特性曲線パラメータを用いて HYDRUS-1D により計算した体積含水率は，現地圃
場の体積含水率の実測値と乖離が見られ，室内試験による水分特性曲線パラメータが実際の圃場の特性を十分に論文：HYDRUS-1D を用いた砂丘未熟土畑における灌水操作方法の検討

11により水分特性曲線パラメータを最適化し，これを用い

Dabach, S., Shani, U. and Lazarovitch, N. (2015): Optimal ten-て HYDRUS-1D により計算した結果，体積含水率の推

siometer placement for high-frequency subsurface drip irriga-

定精度は大幅に改善された．

tion management in heterogeneous soils. Agricultural Water

灌水を開始する際の圧力水頭の閾値と日灌水量を変化

Management, 152: 91–98.

させた場合の計算結果から，砂丘未熟土においては，灌El-Nesr, M.N., Alazba, A.A. and Šimůnek, J. (2014): HYDRUS

水を開始する際の圧力水頭の閾値を −50 cm 以上に設定

simulations of the effects of dual-drip subsurface irrigation

した場合，灌水の大部分が深部に浸透しやすいことが示

and a physical barrier on water movement and solute transport

唆された．一方，同一の条件で閾値を −100 cm 以下に設in soils. Irrigation Science, 32: 111–125.

定した場合では，作物根による水分吸収が抑制されるよ

Feddes, R.A., Kowalik, P.J. and Zarandny (1978): Simulation of

うになると推定された．この結果は，砂丘未熟土の畑で

field water use and crop yield. John Wily & Sons, New York,は，粗孔隙が多いことに由来する特殊な水分特性曲線により，灌水を開始する際の圧力水頭の閾値の最適範囲が狭くなることを示唆している．計算結果をまとめると，

10 月から 1 月にかけて行われるコカブのハウス栽培では，灌水を開始する際の圧力水頭の閾値を −75 cm，日
灌水量 3.2 mm に設定することにより，I · ETP−1 は約 1となり，作物根の吸水抑制や下方への浸透損失も少なくできると判断された．このように，具体的な営農条件を設定してシナリオ解析を行うことによって灌水パラメー
タが定量化されることで，節水かつ安定した作物収量を確保するための水管理計画の策定が可能となると考えられる．

NY.
藤巻晴行, 布和, 天野貴久, 井上光弘 (2012): 砂丘砂における地中
点滴灌漑下の水分溶質移動の数値解析. 日本砂丘学会誌, 59:

1–10.
橋本岩夫, 千家正昭, 伊藤健吾, 丸山利輔 (2010): 砂丘地露地畑のスイカ栽培における栽培管理用水の実態. 農業農村工学会論文集, 78: 289–303.

Hillel, D. (1998): Environmental Soil Physics. Academic Press,
San Diego, CA.
井上光弘, 竹内芳親 (2000): 土壌水分制御によるハウス内の砂
ベッドおよび砂地野菜栽培 — 感圧センサーの開発とその適用 —. 農業および園芸, 75: 600–606.
岩田幸良, 成岡道男, 宮本輝仁, 中村俊治, 松宮正和, 亀山幸司

謝辞(2016): 水フラックスの計算による有効土層深決定方法の検討: — HYDRUS を用いた数値シミュレーションによる新しい手法の提案 —. 農業農村工学会論文集, 84: I 175–I 183.

福井県坂井市農家の川合芳彦氏には，実証試験のため

亀山幸司, 岩田幸良, 佐々木康一, 成岡道男, 宮本輝仁 (2016): 樹の圃場をお借りした上，圃場管理等でご尽力いただきま

皮由来バイオ炭の砂丘地圃場への施用による土壌の保水性 ·した．また，実証試験の実施に当たっては，北陸農政局

保肥性改善効果. 農業農村工学会論文集, 84: I 65–I 74.

九頭竜川下流農業水利事業所にご協力 · ご助言をいただ

勝田英郎 (2007): 福井県三里浜砂丘砂の開発と農業. 日本砂丘学

きました．ここに記して感謝申し上げます．

会誌, 54: 79–83.
河原あゆみ (2013): 九頭竜川下流地区における畑作農業振興の

引用文献

取組. 畑地農業, 659: 2–8.

Marković, M., Filipović, V., Legović, T., Josipović, M. and Tadić,V. (2015): Evaluation of different soil water potential by field

Abbasi, F., Feyen, J. and van Genuchten, M.Th. (2004): Twodimensional simulation of water flow and solute transport below furrows: Model calibration and validation. Journal of Hydrology, 290: 63–79.
Allen, R.G., Pereira, L.S., Raes, D. and Smith, M. (1998): Cropevapotranspiration — guidelines for computing crop water requirements, FAO Irrigation and drainage paper 56. Food andcapacity threshold in combination with a triggered irrigation
module. Soil and Water Research, 10: 164–171.
Müller, T., Ranquet Bouleau, C. and Perona, P. (2016): Optimizing drip irrigation for eggplant crops in semi-arid zones
using evolving thresholds. Agricultural Water Management,
177: 54–65.
Nakamura, K., Harter, T., Hirono, Y., Horino, H. and Mitsuno, T.

Agriculture Organization of the United Nations, Rome.

(2004): Assessment of root zone nitrogen leaching as affected

Bastiaanssen, W.G.M., Allen, R.G., Droogers, P., D’Urso, G. andby irrigation and nutrient management practices. Vadose Zone

Steduto, P. (2007): Twenty-five years modeling irrigated anddrained soils: State of the art. Agricultural Water Management, 92: 111–125.
Blonquist Jr., J.M., Jones, S.B. and Robinson, D.A. (2006): Precise irrigation scheduling for turfgrass using a subsurface electromagnetic soil moisture sensor. Agricultural Water Management, 84: 153–165.
Dabach, S., Lazarovitch, N., Šimůnek, J. and Shani, U. (2013):
Numerical investigation of irrigation scheduling based on soil
water status. Irrigation Science, 31: 27–36.

Journal, 3: 1353–1366.
中村公人, 大串祥子, 池浦康広, 田中宣多 (2016): 畑地用水計画のための HYDRUS-1D を用いた土壌水分移動解析における土壌水分特性パラメータの推定例 (畑地灌漑特集). 土壌の物理性, 134: 25–40.
成 岡 道 男, 宮本 輝 仁, 岩 田 幸 良, 亀 山 幸 司, 中 村 俊 治 (2015):

HYDRUS-1D を使ったユリ畑の上向き補給水量の算定. 農業
農村工学会誌, 83: 565–570.
農林水産省農村振興局 (2016): 土地改良事業計画設計基準
画「農業用水 (畑)」. 農業農村工学会.

計土壌の物理性

12

Rassam, D., Šimůnek, J. and van Genuchten, M.Th. (2004):
HYDRUS-2D による土中の不飽和流れの計算 (取出伸夫,
井上光弘監訳). 農業土木学会土壌物理研究部会 HYDRUS グ
ループ.

第 139 号(2018)

高野幸成, 猪野 誠 (2008): コカブ横縞症に対する複数の発生要因の相互関係. 千葉農総研研報, 7: 1–8.
竹内睦雄, 長江亮二 (1990): 電気探査による地下水流動モニター法の研究. 応用地質, 31: 12–18.

坂井 勝, 取出伸夫 (2007): 砂質土と団粒土に対する水分移動特性関数の検討 (モデル特集論文). 107: 63-77.

van Genuchten, M.Th. (1980): A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil

Segal, E., Ben-gal, A. and Shani, U. (2006): Root water uptake
efficiency under ultra-high irrigation frequency. Plant and Soil,
282: 333–341.

Science Society of America Journal, 44: 892–898.
van Genuchten M.Th. (1987): A numerical model for water andsolute movement in and below the root zone. U.S. Salinity

Šimůnek, J., Šejna, M., Saito, H., Sakai, M. and van Genuchten,

Laboratory, USDA, ARS.

M.Th. (2015): The HYDRUS-1D software package for simu-

van Genuchten, M.Th., Leij, F.J. and Yates, S.R. (1991): The

lating the one-dimensional movement of water, heat, and mul-

RETC Code for quantifying the hydraulic functions of un-tiple solutes in variably-saturated media, version 4.17, HY-saturated soils, version 1.0. EPA Report 600/2-91/065, U.S.

DRUS software series 3. Department of Environmental Sci-

Salinity Laboratory, USDA, ARS.

ences, University of California Riverside, Riverside, CA.

Vellidis, G., Tucker, M., Perry, C., Kvien, C. and Bednarz, C.

Skaggs, T.H., Shouse, P.J. and Poss, J.A. (2006): Irrigating for-(2008): A real-time wireless smart sensor array for schedulingage crops with saline waters: 2. Modeling root uptake andirrigation. Computers and Electronics in Agriculture, 61: 44–

drainage. Vadose Zone Journal, 5: 824–837.

50.

Sonnleitner, M.A., Abbaspour, K.C. and Schulin, R. (2003): Hy-

大和陽一, 前川 進, 稲垣 昇, 寺分元一 (1992): カブの生育に及ぼ

draulic and transport properties of the plant-soil system esti-

す構造物の被陰の影響とそのモデル化について. 農業気象,

mated by inverse modelling. European Journal of Soil Science,

48: 19–28.
弓削こずえ, 阿南光政 (2016): 畑地における作物根の分布と土壌54: 127–138.
Subbaiah, R. (2013): A review of models for predicting soil water dynamics during trickle irrigation. Irrigation Science, 31:

面蒸発の変動を考慮した土壌水分動態解析と消費水量の定量化 (畑地灌漑特集). 土壌の物理性, 134: 17–24.

225–258.

要旨

本研究では，砂丘未熟土畑において，作物の消費水量を満たしながら下方への浸透損失をできるだけ抑える灌水操作方法を検討するため，砂丘未熟土畑で行われている土壌水分を比較的高めに維持する少量
頻繁灌水の下での土壌水分動態を再現すると共に，灌水パラメータを変化させたときのシナリオ解析を行った．室内実験から得られた水分特性曲線パラメータを用いて土壌水分量を計算した結果，現地圃場
の実測値との乖離が見られた．そこで，逆解析から得られた水分特性曲線パラメータを用いて灌水を開始する際の圧力水頭の閾値や日灌水量を変化させた場合のシナリオ解析を実施したところ，灌水を開始する際の圧力水頭の閾値を −75 cm，日灌水量 3.2 mm に設定することにより，根の吸水を抑制することなく，下方への浸透損失も少なくできると判断された．
キーワード：灌水制御，数値解析，土壌水分動態，間断灌漑True TDR, TDT 土壌水分センサー

True TDR（時間領域反射法）方式

ハンディモニタ・ロガー

TDT（時間領域透過法）方式

CACC-TDR-315L

C-M1001

CACC-SEN-SDI

定価：44,000 円（税抜）

定価(本体のみ)：100,000 円（税抜）

定価：34,000 円（税抜）【仕様】
CACC-TDR-315L、CACC-SEN-SDI
要素

体積含水率

導電率

温度

測定範囲

0～100%

0～5.0 dS/m

TDR:-４0℃～+６
0℃
TDT:1℃～+50℃

精度

±2%

TDR：
0.1dS/m(0.1mS/cm)±2.5%
TDT：±0.2 dS/m

TDR:±0.2(-20-50℃)
TDT:±1℃(1-50℃)±1% FS

分解能

0.1%

-

-

-

±1%
TDR:±2%
EC 安定性 0-5ds/m
TDT：±1%
導電率= S/cm = 1/(Ω・cm)
1 [S/m] = 10 [dS/m] = 10 [mS/cm] = 10,000[uS/cm]-

温度安定性 1-50℃

誘電率-

True TDR, TDT センサーは共に反射(TDR)および透過(TDT)信号波形を取得することができます。波形が出力できるので粘土質土壌、穀物の水分量測定など、多様な物質の水分量測定可能性があります。
【TDT 取得波形の例】 キャンベル社製ロガーCR800 を使用

Handy 土壌水分モニター/ロガー

C-M1001 仕様

電源：アルカリ単 3 電池 X3 本
外部電源端子：アルカリ単 3 電池 X3 本 または マイクロ USB typeB
対応センサ：TDT(Acclima 社製)、TDR(Acclima 社製)、CS650(Campbell 社製)
記録メディア：μSD カード 16GB まで FAT32
寸法：80W x 140L x 40Dmm
オプション ：拡張ターミナル(5ch、10ch 対応）
分岐コネクター（2ch）
外部電源（マイクロ USB 端子付）
キャリングケース

クリマテック株式会社
本社：東京都豊島区池袋 4-2-11 ＣＴビル 6F
TEL 03-3988-6616 FAX 03-3988-6613

札幌営業所：札幌市東区北２２条東８丁目４－５
TEL 011-711-9921 FAX 011-711-9922

型式︓WDR-1-3WET 型

主な特徴● 現場で瞬時に⼟壌⽔分/EC/温度を測定可能● ハウス栽培、露地栽培・植物⼯場等、あらゆる環境に適します● 従来では難しかった体積⽔分率 40%以上の測定も可能● ⼤型カラー液晶搭載による⾼い視認性と、タッチパネルによる直感的な操作性● 計測データはメモリーカードに記録され、1 測定当たり最⼤ 65,000 件の記録が可能

体積含⽔率(VWC)、電気伝導度(EC)、温度(℃)の同時計測を⼿軽に⾏うことが出来る携帯型データロガーです。乾電池で動作させる事が出来、電源の確保が困難なフィールドでの計測に最適です。
■仕様
測定項目 ︓体積⽔分率、電気伝導度(EC)、温度
測定範囲 ︓⽔分 0〜100%、EC 0〜7mS/cm、温度-10〜50℃
測定精度 ︓⽔分±3%、EC±5%、温度±1℃
測定間隔 ︓1/10/30 秒、1/10/30 分、1/3/12 時間
ディスプレイ︓3.5 インチ タッチパネル⽅式
記憶容量 ︓1 測定あたり 65,000 件
サイズ重量︓100×130×35mm、290g(本体)

開発・製造元 株式会社Ａ・Ｒ・Ｐ
販売元

本社 〒114-0024 東京都北区⻄ケ原 1-46-16

営業部－ 弊社ホームページ www.fujiwara-sc.co.jp－

TEL 03-3918-8111 FAX 03-3918-8119

千葉 営業所 TEL 0438-64-0800 FAX 0438-64-0820
つくば営業所

TEL 0298-39-4500 FAX 0298-39-4512

J. Jpn. Soc. Soil Phys.
土壌の物理性No. 139, p.15 ∼ 20 (2018)

Handmade Becquerel meterusing commercial scintillation survey meter
Shuichiro YOSHIDA1 , Sho SHIOZAWA1 , Takuya ZAITSU2 , Hiroaki YAMANO3
and Kazuhiro NISHIDA1

Abstract: After devastating nuclear disasters, the Becquerel meter should be very important for monitoring thecontamination of our surrounding environment and judging the risk of exposure. However, the commercial Becquerel meters are too expensive for the citizens to purchase. Even the availability for the administrative agency
is not ensured in emergency. Considering these risks, aself-produced Becquerel meter using a commercially available scintillation survey meter is proposed. The sample formeasurement should be set in a cylindrical bottle having
a diameter of 3 cm and the height of 5 cm. Although theshielding of the system with lead is not perfectly given,
the subtraction of the background dose enables the measurement of the radiocesium concentration of the contaminated soil samples with sufficient accuracy. The theoretical error analysis revealed the relative error depends on the
total counts and the contribution of the background radiation. The counts required for the measurement with 10% of error is 4500 if the contribution of the backgroundradiation is 80 % of the total counts, while the required
counts reduced to 600 if the contribution of backgroundis 50 %. The correction methods to apply this system tothe sample with variable thickness in the bottle is also presented. Because the relationship between the counts rate
and the radiocesium concentration depends on the proportion of 134 Cs to 137 Cs, the calibrations are generally required at each measurement. However, no significant difference was revealed for the calibration lines measured atJan. 2013 and Aug. 2013 during which the proportion of
134 Cs highly decreased. Thus, any frequent calibrations
are not necessary for the practical use of the system. Whenthe aforementioned conditions are satisfied, the proposed
Becquerel meter meets the practical demand for knowing
the radioactivity of soil samples.
Key Words : radiocesium concentration, radioactive contamination, simple measurement, soil contamination, contamination of agricultural products
1 Graduate School of Agricultural and Life Sciences, The University ofTokyo, 1-1-1, Yayoi, Bunkyo, Tokyo, 113-8657, JAPAN. Corresponding
author: Shuichiro Yoshida, 東京大学大学院農学生命科学研究科.
2 Ministry of Agriculture, Forestry and Fisheries, 1-2-1 Kasumigaseki,
Chiyoda-ku, Tokyo 100-8950, Japan.
3 Kajima Corp. 3-1, Motoakasaka 1-chome, Minato-ku, Tokyo 107-8388,
Japan.
2017 年 12 月 22 日受稿 2018 年 5 月 23 日受理

1. Introduction
On unforeseen contamination by radionuclides in a nation scale such as an accident of nuclear power station, notonly researchers or governmental officers who should prepare apparatus to monitor the situation, but citizens also eager to know the real situation. At the first step, they wouldpurchase dose meter or survey meter, to monitor the surrounding environment. But it is impossible to measure theconcentration of radionuclide in samples of soil or foods
without having an expensive professional instrument. For
the official purposes, the methods to measure the concentration of radionuclides contained in foods or soil is defined in the guidelines published by the government. For
example, Ministry of Health Labor and Welfare (2012) established a guideline on the analytical method of radionuclides in foods in which the analysts must use a germanium
semiconductor spectrometer (Ge spectrometer) having relative efficiency larger than 15 % and surrounded by leadshield with a thickness of 10 to 15 cm. NaI (TI) spectrometer which is much cheaper than Ge spectrometer havebeen also used for screening or voluntary analysis of foods.
Although the lower limit of the measurable concentrationby the system is as dilute as 1 Bq kg−1 , it is yet too expensive and bulky for citizens to purchase. The availability is also questioned after devastating accidents. In fact,
the shortage of the instruments was reported just after theaccident of Fukushima Daiichi Nuclear Power Station in
2011, leading researchers on the front-line to seek alternative methods to measure radioactivity of soil samples using
commercially available scintillation survey meter (Nemoto
et al., 2013). They proposed to measure the radiation fromthe surface of the 1 kg of soil sample packed and filled
in a 2 L beaker. The background radiation is subtracted
from the total counts to determine the net radiation fromthe sample. Although the principle to measure the radioactivity of samples is common for all methods, the mass ofthe sample, the shielding capacity, and the uniformity of16

土壌の物理性第 139 号(2018)

Fig. 1 Schematic of the proposed apparatus.

the geometry are the highly substantial factors to establish
useful methods. The method proposed by Nemoto et al.
(2013) does not require any apparatus except the scintillation survey meter. However, the large mass of soil sample
is necessary to assure the consistent accuracy.
In-situ measurement of vertical profile of radiocesium
concentration in soil inexpensively was also studied.
Suzuki et. al. (2013) arranged photodiodes along a rod
to detect the vertical radiation profile from which the concentration profile is analyzed. Due to the low sensitivity ofthe photodiode sensors, long counting time is required for
accurate measurements. But the required time can be reduced by replacing the sensors to more sensitive ones, withadditional cost.
The present report aims to provide an information about
making a cheap but reliable Becquerel meter from a commercially available gamma ray survey meter and rolled
lead sheet. The accuracy of the apparatus is thoroughly
discussed to present the applicability of the method.

2. Design of the apparatus
Measurement of radiocesium concentration of contaminated soil samples by arranging a commercial scintillation probe (Model 5530) and survey meter (Model 5000CYPHER) was tested. The gamma sensitivity of the scintillation probe is approximately 30,000 cpm µ Sv−1 h−1 for
137 Cs. according to the manufacturer. The survey meter
including the probe approximately costs 4,000 US dollars.
The GM (Geiger-Muller) tube survey meters are cheaper
than scintillation survey meters, though the measuring time
would have been longer due to the lower sensitivity. Fig.

1 shows the design of the instrument. The instrument consists of a shielded sample holder and shielded scintillation
probe.
A cylindrical shielding cap (a) and shielding frame (b)were made of lead sheets which is soft enough to shape.
The cost for the 6.5 kg of rolled lead sheet used for theshielding cap and frame was about 200 US dollars. We
wrapped multiple thin lead sheets having thickness of 3
mm around a cylindrical mold having a diameter of 35mm, each of which was glued with adhesive agent (Super
Glue). The resultant thickness of the shielding wall was 20mm, which attenuates direct gamma beam to 7 %. Because
the lead wall alone is not enough to shield the backgrounddose completely, subtraction of the background radiation isnecessary. The procedure will be discussed later. The topof the cylindrical cap (a) was filled with a top lead cap (h)having thickness of 20 mm. The probe (c) must be set at thepredefined position in the cylindrical shielding to be fixedby the screws when used. Thus, the probe can be takenout from the frame easily to be used as a survey meter tomonitor the dose rate in environment. The NaI scintillator(d) having a diameter of 25.4 mm and a length of 25.4 mm(1 inch) is installed at the top of the probe, where gamma
beam from the sample and the background is converted tolight and then the light is converted to electronic pulse signal by the photomultiplier. The probe is connected to thesurvey meter (g) (analyzer) which displays count of photon
in a predefined period. Because the system we used cannot detect the energy spectrum, the response depends onthe energy. The major energies of the gamma beam fromradiocesium are 661 keV (137 Cs), 605 keV (134 Cs) and 796
keV (134 Cs) (e.g. Malonda and Carles, 2012). Due to the

研究ノート：Handmade Becquerel meter using commercial scintillation survey meter

17

Fig. 2 Blank counts per 10 seconds; (a) Fluctuation during 600 seconds; (b) Correlation incounts per 10 seconds between a couple of the apparatus put next to each other.

dependency of the sensitivity on energy, the relationshipbetween the counts and the dose can change when the proportion of 134 Cs and 137 Cs changes due to time. This drawback will be discussed later. A couple of the measurement
sets: probes, caps, frames and the controller, were constructed to compare the individual specificity.

tors influencing the accuracy, duration of the measurement,
the background dose rate, and the thickness of the sample
in the bottle were examined. All the measurements were
performed in an unshielded laboratory in the University ofTokyo located in Bunkyo-ku, Tokyo. Only the Set A wasused for the measurement except for the comparison between the probes.

3. Methods of calibration
and accuracy test

4. Results and discussion

The contaminated soil samples were taken from a profile
of a paddy field located in Iitate Village, Fukushima, Japanand the bottom of a reservoir in Motomiya City, Fukushima
Japan. The samples were subject to drying in an oven at
105 ◦ C. The highly contaminated sediment of the agricultural reservoir was mixed with non-contaminated volcanic
ash soil to make a specimen with variety of predefined radiocesium concentration, which was 100 %, 80 %, 60 %,40 %, 20 % and 10 % of the original sediment by weight.
The volcanic ash soil was taken from the experimental fieldof The University of Tokyo located in Nishi-Tokyo city
before the accident of Fukushima Daiichi Nuclear Power
Station. The radiocesium concentration of these samples
were measured by the Ge semiconductor detector run byRadioisotope (RI) Center of the Graduate School of Agricultural and Life Sciences, The University of Tokyo. The
samples were also subject to the counting of photons by
the proposed system. The contribution of the sample tothe detected count of photons were evaluated by subtracting the background from the total. The background canbe determined by measuring the radiation without setting
the samples, to be called blank. Among the expected fac-

The variation of the background radiation was examinedusing a couple of the proposed apparatus which were thesame model of the survey meter but the sensibility to the
gamma ray fluence are different due to the shielding capability. Fig. 2 shows the temporal variation of the background detected by the probes (Set A and Set B) in 10 seconds. Observed numbers of photons in space at predefinedtime interval is widely known to follow Poisson distribution (e.g. Malonda and Carles, 2012) represented by P(λ ).
When the probability variable X follows the Poisson distribution and the mean value is λ , the variance of X is also
λ . The observed mean (Mean) and variance (σb 2 ) shownin Fig. 2a are consistent with this theoretical knowledge.
Because the mean count is proportional to the duration of
the measurement, the theory tells that the standard deviation increases proportionally with the root square of duration of measurement. This means that the longer we spent
for the measurement, the less the relative error becomes.
This principle is important in measuring dilute samples.
Because the count due to the sample alone (Ns ) is evaluated by subtraction of the count due to the background(Nb ) from the measured count (N):

土壌の物理性

18

第 139 号(2018)

dition because Nb /N cannot be reduced. For instance, thetotal counts required for the concentration measurementwith 10 % of error (σs /Ns = 0.1) is 4500 if the contribution of the background (Nb /N) is 80 %, while the required
counts reduced to 600 if the Nb /N is 50 %. In general, therelative contribution of Nb : Nb /N is not known in advance.
Thus, the rate of counts should be measured for the sample with the minimum radioactivity among the samples forsome minutes to roughly estimate Nb /N.

Fig. 3 Theoretical estimate of relative error in the countdue to sample Ns as a function of the proportion of the background radiation Nb to the total count N.

Ns = N − Nb(1)

the range of error in the count due to the sample can be
estimated by eq. (2) based on the principle of propagationof error:

σs2 = σ 2 + σb2(2)where σs 2 , σ 2 and σb 2 are variance of Ns , N and Nb , respectively. Since N = σ 2 and Nb = σb 2 ,σs2 = N + Nb

(3)

Then the relative error in Ns can be written as√σs
N + Nb=Ns
N − Nb(4)

The diagram of this relationship is provided as Fig. 3,
which is useful to estimate the minimum total count to obtain the count due to the sample with practical accuracy.
We should usually wait until the total count meets the con-

Although two probes were set on the same table, thecounts in 10 seconds measured by the both probes were
never synchronized to each other. Fig. 2b examines thecorrelation between the counts by both probes. It is obvious no correlation was found in the temporal trends ofthe counts in 10 seconds detected at the neighboring location. This basic characteristic show real-time subtractionof the background does not work to increase the accuracyof measurements.
Table 1 shows the variation in the averaged blank counts
per minute (CPM) due to the measurement dates. The measurement was kept for an hour and the counts per every
minute were averaged. The variation due to the dates canbe negligible compared to the fluctuation observed in every
minute. Therefore, frequent measurement of blank CPM is
not required unless the environmental radiation is known
to be considerably change. However, needless to write, themeasurements of samples and blanks must be done in thesame place because the background radiation highly depends on the place.
Fig. 4 shows the relationship between total radiocesium concentration of mixture of the contaminated sediment and andosol measured by Ge semiconductor detector
and Counts in 400 seconds by the proposed apparatus. The
temporal decay of the concentration is corrected based on
the date of the measurements. The relationship was so linear that the CPM obtained by the proposed apparatus canbe converted to concentration. However, we should notethat the relationship could change depending on height,
density and 134 Cs/137 Cs ratio of samples.
The thickness (height) of the sample perpendicular tothe scintillator affects the self-attenuation of the gamma

Table 1 Variation in averaged blank CPM due to the measurement dates.

Set A Meanσ
Set B Meanσ

2016/3/11

2016/3/14

2016/3/14

2016/3/15

2016/3/21

PM

AM

PM

AM

PM

490

491

490

488

492

21.7

23.3

25.5

24.3

21.7

311

309

312

310

307

18.0

17.1

18.1

17.0

17.6

研究ノート：Handmade Becquerel meter using commercial scintillation survey meter

19

beam. Fig. 5 shows the relative decrease in the CPM perunit weight of soil with an increase in the thickness of thesample in the bottle. Thickness of 33 mm was applied as
the reference. The 10 % of decrease in thickness caused
5 % increase in CPM for paddy soil; the 30 % decreasecaused 15 % of increase in CPM for the sediment. The dependency of CPM on the thickness had a slight difference
between the samples having different bulk density. The
dependency of CPM on the thickness for the denser sample (paddy soil) was less than that for the looser sample(sediment). If the thickness of the sample is adjusted tothe standard thickness, the difference in the soil density ispractically less important.
Fig. 4 Relationship between total radiocesium concentration measured by Ge semiconductor detector and counts per
400 seconds by the proposed apparatus (Set A). The sample
was the mixture of sediment in an agricultural reservoir andnon-contaminated upland soil. The temporal decay of theconcentration is corrected based on the date of the measurements.

Because the measurement of numerous samples usually
took so much time that the radiocesium concentration of
the samples were subject to the effect of natural decay.
Thus, the measured concentrations should be converted
into the ones at the reference date by considering the halflife period and existing ratio of 134 Cs and 137 Cs. However,since the survey meter we used cannot detect the energy
spectrum, the response depends on the component of theenergy of the gamma ray. Fig. 6 compares the differencein the relationship between the radiocesium concentration
and CPM due to the measurement date: 23 Jan. and 29
Aug. 2013, during which the proportion of 134 Cs to thetotal Cs should decreased from 0.34 to 0.28. Because thecovariance analysis reveals no significant difference in the
regression coefficients, the calibration is not frequently required unless the measurement spans many years.

5. Conclusion
Fig. 5 Dependency of CPM on thickness of sample inbottle. Normalized to the CPM for the sample with thickness of 3.3 cm. Set A is used for these measurements.

Self-produced Becquerel meter using commerciallyavailable scintillation detector is proposed. Sufficient accuracy can be ensured if the measurement is conducted foradequate time determined from the background radiationand the concentration of the sample. Shieling by lead sheet
helps reduce the effect of background radiation and shorten
the time required for the accurate measurement. Although
any types of detector can be principally used, the NaI scintillation probe should be the best regarding both the cost
and the sensitivity. Detectors with a function of counting
photons in any time interval is highly recommended to produce the system.

References
Fig. 6 Temporal change in relationship between CPM by
the proposed apparatus (Set A) and total radiocesium concentration. The samples were taken from a profile of a contaminated paddy field. The concentrations were correctedto the day when the CPM measurement was performed.

Nemoto, T., Sato, M., Ono, Y., Fujimura, K., Nohei, N. andSuzuki, Y. (2013): A simple method for the measurementof radiocesium in soil samples using a NaI(Tl) scintillation
counter. Research Report of Fukushima Agricultural Technol-

土壌の物理性

20

第 139 号(2018)

ogy Center, Special Issue on countermeasures for radionu-

for Radioactive Substances in Food, Notice No. 0315 Article

clides, 15–18. (in Japanese)
Malonda, A.G. and Carles, A.G (2012): Radioactivity counting

Suzuki, S., Iwase, H., Noborio, K., Mizoguchi, M., Kobayashi, D.

statistics, in L’Annunziata, F., Michael (eds.) Handbook of Ra-

and Ito, T. (2013): Development of a device for measuring the

dioactivity Analysis 3rd Edition, pp.163–189, Elsevier Inc.

vertical distribution of radioactivityin soil using photodiode.

Ministry of Health, Labor and Welfare (2012): Testing Methods

要

4 of the Department of Food Safety, March 15, 2012.

Journal of Japanese Society of Soil Physics, 125: 29–33.

旨

原子力事故の後には，周辺環境の汚染モニタリングや被爆リスクの判断のためにベクレルメータが不可欠である．しかし，市販のベクレルメータは一般市民には高価であり，しかも非常時には行政機関でさえもすぐに入手することはできない．これらのリスクを考慮し，比較的入手が容易な市販のシンチレー
ションサーベイメータを用いてベクレルメータを自作し，その精度を検証した．作成したベクレルメー
タは，シンチレーションプローブの周囲に鉛板を巻き付けることで，簡易な遮蔽を行うとともに，ブランク測定によるバックグラウンドを差し引くことで，汚染試料の測定において実用上十分な精度を得ることができた．
キーワード：放射性セシウム濃度，放射能汚染，簡易計測，土壌汚染，農産物汚染

● デカゴン社(Decagon Devices Inc.)とUMS社(UMS AG)は2017年1⽉に合併
し、デカゴン社とUMS社の社名を[METER Group, Inc.]に変更いたしました。

ATOMOS-41

複合型気象計測ユニット

EM60
ソーラーデータロガー

1つチャンネルで14の気象観測項⽬を同時測定しますEM60ソラーデータロガーは、すべてのMETER社製センサーとの接続が可能です[技術仕様］◎ 測定項目・雷の回数
・日射量
・雷からの距離・降雨量
・傾き
・最大降雨量・コンパス方向
・気温
・気圧
◎ 寸法
・水蒸気圧
直径: 10 cm 高さ: 34 cm
・相対湿度
・風速
◎対応データロガー
・風向EM60, EM60G
・瞬間最大風速[特徴］
・最大6種類のセンサーを接続可能
・ソーラー充電機能を装置、長期間の計測が可能
・すべてのMETER社製センサーの接続が可能
・温度センサーと気圧センサーを内蔵
・8MBメモリ(40,000〜80,000データを保存)
・GPS搭載、設置場所の管理が可能
・自動時刻同期を搭載

アイネクス株式会社(⽇本総代理店)〒144-0035 東京都⼤⽥区南蒲⽥2-16-1
テクノポートカマタセンタービル
TEL：03-5713-0388
FAX：03-5713-1388
Web site：www.ai-nex.co.jp

バイアル瓶脱気装置セット（アクリル製）
EN-1003★特徴・小型容器（バイアル瓶等）を素早く容易に減圧状態にする事ができます・15ml のバイアル瓶では、同時に約 35 本の減圧が可能です
・組立式ですので部品の交換が容易です★セット内容・バイアル瓶脱気装置本体（真空計・バルブ付）
・真空ポンプセット（ホース・オイル・オイルミストトラップ付）★定価（税別）脱気装置本体：￥250,000
真空ポンプセット：￥130,000★製造・販売元（有）エンドウ理化〒001-0910
札幌市北区新琴似 10 条 7 丁目 3 番 16 号
Tel：（011）763-1088
Fax：（011）763-1667
Mail：endo-rika@mb.infosnow.ne.jp※本仕様は予告なく変更する場合がございますJ. Jpn. Soc. Soil Phys.
土壌の物理性No. 139, p.23 ∼ 24 (2018)

土壌物理をめぐる研究環境
井上久義 1

出だしから私事で大変恐縮ですが，私は 2013 年 3 月

方が来られる場合が多く，私たち専門外の人間の通り一に定年退職し，その後，再任用職員として同じ職場で働遍の説明ですますことはできません．どうしても専門のき，この 2018 年 3 月に任期満了で二度目の退職をしま研究チームに対応をお願いすることになります．そしてした．この 5 年間は広報普及室で勤務し，問い合わせ ·これが，多いときには週に 2〜3 回重なるときがあり，見学の対応，一般公開やサイエンスカフェ等イベントのチームでの研究の進行を相当邪魔しているだろうことは実施，所内ゼミの開催等々の仕事をしましたが，研究室想像に難くありません．

では味わえない面白い体験ができ結構楽しませてもらいました．

加えて感じたことは，研究に対する管理 · 運営の優先の傾向です．サイエンスカフェ自体は，研究における基問い合わせでは，多くは専門チームに回せば良いので

礎学問は特別なものではなく，現実生活に深く密接に関すが，中には「除草剤をまいたらコ̇ブ̇のようなものがでわっていることをわかりやすく地域の皆さんに理解してきて困る」のような問いあわせがあり，瘤のようなもの

もらうために面白いイベントだと私は思っています．加ができるのは聞いたことがなくメーカーもしくは普及所

えて，その準備段階では，自分の知識を再整理する上でに問い合わせてもらうようにお願いしました．後でゆっ非常に役に立ちます．しかし，実際に携わってみると，くり考えると，コ̇ブ̇は瘤ではなく昆布のことかと思い至年 2 回することに決まっているからとにかくそれを卒なりました．大きな雑草類が抑制されるとコケなどの地衣くこなすこと自体が目的となっているのではと思われる類が優勢になりますが，それが昆布のことなのかと一人節が多々見られました．わずか 4 研究チーム，研究者が苦笑していました．見学では，毎年 5 月の連休の頃，近30 人を割るような職場で，これを毎年繰り返せば研究側

くの小学校の新 5 年生が見学にやってきてくれ，そこでの負担がかなり大きくなります．せっかくのいいイベン

研究センターでの仕事の紹介を行いました．3，4 年生のトなので回数を減らしてでももっと大切に行えばいいの

社会 · 理科の指導要領を取り寄せ，どこまで話すか，どにと思っていました．

のような言葉が使えるかなどについて勉強し話をしまし

この他，所内ゼミの担当も行いました．ゼミはお互いたが，子供たちの笑顔や一生懸命聞いてくれている顔々の行っている研究の情報交換，ブラッシュアップというが記憶に残っています．また，サイエンスカフェでは土意味から重要なものと考えています．特に，いろいろな壌物理の話を 2018 年 2 月にさせてもらい，いい思い出分野の研究者が集まる地域の研究センターでは大切ではとなっています．

ないかと思っています．これを進める中で，「ゼミは所広報の仕事は個人的には楽しむことはできましたが，の規定にあるのだから，すべきものとして順番にやらせ

研究面から考えるとこれでいいのかなと思うところがればいい」という意見を幾度となく伺いました．どうも多々あり憂鬱な気分にさせられることが度々ありまし

こうした考え方になじめず，その都度めぼしい人に打診た．

して可能な人に話していただくという形で進めていきま一番感じたのは，あまりにもイベントが多すぎること．

した．

これは，「効果的 · 効率的な研究開発を総合的に推進し，こうしたことは広報普及内での小さな経緯にしか過ぎ

それらを広く社会に移転する」という機構の中期目標にませんが，もうひとつ大事な点は，研究という場を離れ，もかかわることでしょうが，各種ビジネスフェアへの展

広報普及という此岸から研究という彼岸を眺めることが示，マッチングフェスタ，一般公開，サイエンスカフェ

できた点にあります．

そして所内見学等々があります．所内見学では生産者の
1 元 農業 · 食品産業技術総合研究機構

究拠点2018 年 5 月 14 日受稿

西日本農業研究センター四国研2018 年 5 月 27 日受理

現在，研究チームの置かれている環境は非常に厳しいものがあります．チームにおける研究の基本的な活動資金というべき交付金が年々減らされ，研究活動はもとよ24

土壌の物理性第 139 号(2018)

り，重要なスタッフであるパートさんの賃金もそれのみ

告の作成，次の職場の探索等で忙殺されてしまいます．

では賄いきれない状況にあります．このため，チーム運大学でも同様で，昨年の新聞では，40 歳以下の若手研営は競争資金に頼らざるを得なくなります．プロジェク

究者の 6 割以上が任期付任用になっていると報道されてトに参画すると，年 2 回の推進会議と現地検討会の開催

います．また，このような厳しい現状の中，博士課程へが必要となります．さらに，課題評価，中間評価，会計の進学者が大幅に減少しているとも言われています．こ処理のための膨大な資料作成が求められます．プロジェ

うした状況を見ていくと果たして基本的な問題についてクトの期間は長くて 5 年，下手をすると 3 年となります．

3 年などの場合，初年度は試験地の選定，測定機器の設どこが研究するのか心配になります．上のような Natureの報告が出るのも宜（むべ）なるかなと思われます．

置等々で忙殺され，中 1 年置いて 3 年目には成果のとり

現在，地域の研究センターでは土壌のことがわかる人

まとめ，そして，次のプロジェクトへの応募の準備とながほとんどいなくて困っていると聞いています．作物栽

ります．それに加えて，広報普及から雨あられと降って培や育種の現場で新たな品種や技術開発を進めている

くる各種イベントや見学の対応が重なってきます．各研中で，排水や灌漑を通じて土壌と水の関係や土壌特性の究チーム，特にチーム長のご苦労はいかばかりかと察し把握などの問題が避けて通れない重要な課題となっていて余るものがあります．

ます．一方では，土壌物理においては最近，pF（中村,チーム長の方々は 40 歳台位の人が多いのではないか2018）や TDR による水分計測（吉本ら, 2017）の特に作と思います．ちょうど研究に脂の乗った頃合いで，研究土における不安定性についての問題提起がなされていま者として非常にポテンシャルの高い方々ですが，とてもす．このことは栽培技術の開発にとって基本的かつ fatal

自分の研究を進めていく余裕などないのではないかと察な課題となります．これらの問題に対して，地域研究セ

せられます．チーム員も資料作りなどで忙殺される中，ンターにおいては，たまたま在籍する土壌物理研究者のプロジェクトで 3 ∼ 5 年でなにがしかのアウトプット

個人的な熱意や才覚，好意に依存しているのが実態ではを出すことが求められます．現場の問題を解決するためないかと思います．栽培技術等の開発を進めるのなら，の技術を迅速に開発するというプロジェクト研究の必要組織として人員と予算を確保し，技術開発に関わる基本

性は必ずしも全面否定しませんが，あまりにもそれに偏的な問題を追究するともに各研究チームの要請に対応す

りすぎるのはいかがなものか，今少し研究現場に余裕がる研究グループを整備してもいいのではないかと思ってあってもいいのではないか常々感じています．

います．

英国の科学雑誌 “Nature” では特集記事として “Na-

独立行政法人化以降，約 20 年にわたる科学技術政策に

ture Index 2017 Japan” を掲載しています（Fuyuno, 2017;

ついていろいろな弊害が出てきており，今，tipping point

Phillips, 2017）．この中で，日本の科学論文数はこの 10
年間（2005–2015）で 14 分野中 11 分野において減少し，高品質な科学誌 68 誌に限れば 5 年間で 8.3 % 減少してを迎えているのではないかと思います．今後，予算や人員などより良い研究環境を確保できるように改善されていくことを願ってやみません．

おり，この 10 年間で日本の科学研究は失速していると報告しています．そして，その原因としては，大学等の引用文献

運営交付金が毎年減額されていること，科学技術予算の「選択と集中」による競争的資金獲得の困難さ，若手研究者を中心に任期付任用による不安定な処遇などがあげられています．
研究環境の厳しさは，研究職の新規採用にも表れています．現在，私たちの職場での新規採用者のほとんどはポスドクで原則 3 年間の任期付採用になっています．研Fuyuno, I. (2017): What price will science pay for austerity? Nature Index 2017 Japan, 543(7646): S10–S15.
中村公人 (2018)：現地土壌における土壌水分観測について思うこと．土壌の物理性, 138：53–54.

Phillips, N. (2017): Striving for a research renaissance. Nature
Index 2017 Japan, 543(7646): S7.

究チームがチーム長を中心として競争的資金獲得に翻弄

吉本行雄，渡辺優子，喜多 愛 (2017)：畑地水分消費調査の水分計される中，新規の若手研究者も初年度は課題の設定，計値のキャリブレーションを採土法で行う場合の留意点．第測機器など研究環境の整備で費やされ，3 年目は最終報

72 回農業農村工学会中国四国支部大会講演要旨：106–108．

会 務 報 告I. 2017（H29）年度 第 3 回評議員会
月日：2018 年 2 月 17 日 ∼ 26 日
開催形式：メール会議
出

2. 報告
1）編集幹事より，特集および原稿の閲読進捗状
況について報告がなされた．

席：中辻敏朗，塚本康貴，安中武幸，冠 秀昭，
宮本輝仁，北川 巌，江口定夫，小林政広，
吉川省子，吉田修一郎，渡辺晋生，小杉賢一朗，中村公人，望月秀俊，中野恵子，波多野隆介，足立泰久，牧野知之（委員 18

III. 会員消息（2018 年 6 月 9 日まで，敬称略）

名中 18 名参加：成立）

退会 正 会 員

石川

茂雄

正 会 員

鳥山

淳平

1. 2018 年土壌物理学会大会について

正 会 員

菅原

和夫

事務局案（土壌の物理性またはホームページを

正 会 員

増永

二之

岡

香菜子

山嵜

高洋

審議事項

ご参照下さい．）が全会一致で承認された．

学生会員
入会 正 会 員

報告事項：

日本大学生物資源学部

1. 地球惑星科学連合（JpGU）2018 年大会への正 会 員

山田 百合子学生会員

安達セッション共催

北九州市立大学
国際環境工学部遥

東京農工大学生物
システム応用学府

II. 2018（H30）年度 第 1 回事務局会議
月日：2018 年 4 月 17 日
場

所：北海道大学農学部

出席委員：石黒（会長），竹内（副会長），澤本（編集委員長），柏木（事務局長），倉持（会計幹事），塚本（編集幹事）
議

題：1. 2018 年度土壌物理学会大会について
2. 学会運営に関して
3. その他

III. 2018（H30）年度 第 1 回編集委員会
月日：2018 年 5 月 8 日 ∼ 5 月 15 日
開催形式：メール会議
出席委員：澤本卓治，飯山一平，岩田幸良，北川 巌，
小林幹佳，三枝俊哉，清水真理子，鈴木伸治，中川進平，中野恵子，中村和正，丹羽勝久，笛木伸彦，渡辺晋生（委員 14 名中 14 名参加：成立）
議事内容

1. 議題
1）編集委員長より「土壌の物理性 139 号」の発行計画が提案され，全会一致で承認された．

2）編集委員長より新たな投稿区分「土壌圏」の新設については白紙とする旨の提案がされ，全会一致で承認された．

現在会員数（2018 年 6 月 9 日現在）
正会員

：

240

シニア会員 ：

16

学生会員

58

：

賛助会員
：
8
—————————————
小 計
：
322
購読会員
：
37
—————————————
合 計
：
359

データロガー ＦＴＪｒ（エフティージュニア）『ワンタイム測定機能』搭載で、フィールドを移動しながら多点測定
低価格ながら、各種のセンサと直接接続でき、様々な計測に利用できるとご好評をいただいているＦＴＪｒには、設定したインターバルによる一般的なデータロガーとしての機能のほか、ボタン操作時の測定値を記録する『ワンタイム測定機能』が搭載されています。
電源・測定スイッチを【ＳＥＴ】にして、【ＥＮＴＥＲ】キーを長押しするだけの簡単記録。
ワンタイム測定機能を利用することで、従来のインターバルによる計測はもちろんのこと、フィールドを移動しながらの多点測定にご活用いただけます。
記録機能を持たない表示器を利用している場合など、表示された値を野帳に手書きするなどの作業が不要になり、手間の低減と転記時の誤記入の防止を実現します。

圃場を移動しながら、土壌水分の多点計測を行う場合の例▲
ACCESSSAMPLECANCEL

COPY

ENTER

▼

MEAS
SET
OFF

ENTERボタンを押すだけで簡単記録【ワンタイム測定機能】● 本体のボタン操作で測定実行、測定値をメモリに記録● 測定値以外に日付、時刻、測定番号も記録● 直前の測定値をキャンセルできるので、不要データの取り消しが簡単● 記録データはＵＳＢメモリで簡単回収(データファイルはＣＳＶ形式)● ワンタイム測定データと通常の測定データは回収後に識別可能

設置中の土壌水分計

例えば土壌水分計測では、標準のアルカリ乾電池４本利用時において、IMKO社TRIME-PICOの場合
で約1,100回、Delta-T社SM150Tで約1,400回、A･R･P社WD-3シリーズで約800回の測定が可能です。
※

内蔵電池での測定回数は、１回の測定及び表示時間を１分間としたときの計算上の回数です。
動作間隔、周囲温度などにより測定回数は大きく変動しますので目安としてください。

ワンタイム測定機能は、土壌水分計測に限った機能ではなく、他のセンサとの組み合わせにおいても利用できます。また、SDI-12 出力タイプの METER 社 5TE（旧 Decagon Devices 社）など、RS232C、RS485、SDI-12 のシリアル出力タイプのデジタルセンサに対応する FTJr-Digital （エフティー
ジュニア デジタル）にもワンタイム測定機能が搭載されています。
お手持ちのセンサとの接続の可否など、お気軽にお問い合わせください。
詳しくはFTJr ワンタイム で 検索
〒060-0063 札幌市中央区南３条西８丁目７番地４ 遠藤ビル５Ｆ
TEL011-596-0201 FAX 011-596-0234
E-mail info@mcs-fs.com URL http://www.mcs-fs.com

27—編 集 後 記—
土壌の物理性 139 号をお届けします．いかがだったでしょうか？次回，140 号からは寒冷地特集がスタートします．楽しみにしていてください．

入会し，最新の知見を収集するようにしています．
そのような中，昨年度，今年度と多湿黒ボク土畑で融雪時の暗渠排水の流量観測を行う機会がありました．農

さて私ですが，土壌物理学会の編集委員を担当するの耕期間が積雪で限られている北海道ですから，融雪時のは実は 2 度目で，1 度目の任期は 2011 ∼ 2012 年度でし

迅速な排水は大変重要です．しかし，昨年度は融雪時にた．編集委員の顔ぶれを見ますと，1 度目の時も 2 度目

全く暗渠から水が流出しませんでした．その後普通に営の時も大学や公設農試等，研究機関の方々が中心で，私農が開始され，そもそもこの圃場に暗渠が必要なのかとのような民間企業に所属している者が編集委員を担当すいう疑問を持ちました．しかし，今年度は融雪直後からるのはまれだと思います．構成メンバーが研究機関中心

暗渠排水からどんどん水が流出し，流出が終わったのはの学会なので，「なぜ民間企業の会社員が土壌物理学会

観測を始めてから約 4 週間後でした．確かに今年度は 3に？」と疑問を持つ方もいらっしゃると思います．そこ

月上旬に大雪が降ったのですが，そのことが，この極端で簡単に私の会社の紹介をしたいと思います．

な融雪時の暗渠排水の年次差を説明できるほどの要因に私の所属する（株）ズコーシャは，帯広に本社を構え
る建設コンサルタント会社です．北海道十勝という土地

思えません．改めて，土壌の物理性の奥深さを感じているところです．

柄，農業に関連する仕事が多く，畑地や草地の測量や設後 1 年弱，編集委員に従事する期間が残されています．

計を行う部門に加え，土壌 · 作物調査等の農業調査部門

期間中に，この奥深き土壌の物理性を誌面上で少しでもを持っています．私は調査部門に所属しており，暗渠排紐解いていければ良いなと思っています．

水や土層改良等，土壌の物理性に関係する調査に取り組
むことが多いです．このような経緯から土壌物理学会に

丹羽勝久（編集委員）

土壌物理学会
事務局構成

会
長
副 会 長
事務局長
編集幹事
会計幹事
会計監査

編集委員会

委 員 長
委
員

石黒
宗秀
竹内
晴信
柏木
淳一
塚本
康貴
倉持
寛太
志賀
弘行
横濱
充宏
澤本
卓治
飯山
一平
岩田
幸良
北川
巌
小林
幹佳
三枝
俊哉
清水 真理子
鈴木
伸治
中川
進平
中野
恵子
中村
和正
丹羽
勝久
笛木
伸彦
渡辺
晋生

（北海道大学）（北海道立総合研究機構）（北海道大学）（北海道立総合研究機構）（北海道大学）（北海道立総合研究機構）（土木研究所寒地土木研究所）（酪農学園大学）（宇都宮大学）（農業 · 食品産業技術総合研究機構）（農業 · 食品産業技術総合研究機構）（筑波大学）（酪農学園大学）（土木研究所寒地土木研究所）（東京農業大学）（秋田県農業試験場）（農業 · 食品産業技術総合研究機構）（土木研究所寒地土木研究所）（（株）ズコーシャ）（北海道立総合研究機構）（三重大学）

土壌物理学会会則

2）シニア会員への変更は，資格を有する会員からの申告を受け，事務局が承認する．原則として申告された年度か
第 1 条 本学会は土壌物理学会（Japanese Society of Soil Physics）
ら適応するが，当該年度中に有資格となった場合は，次年度から適応する．
と称する．
第 2 条 本学会は土壌物理に関する研究の進歩と普及を図り，農業 3）シニア会員の権利は正会員に準ずる．ただし，会則第 6条で定める役員および第 9 条で定める委員会委員に就く
技術および環境科学の発展に貢献することを目的とする．
ことはできない．
第 3 条 本学会はその目的を達成するために次の事業を行う.
（2011.10.28 制定 2017.10.14 改正）（1）研究発表会，討論会および見学会などの開催（2）学会誌「土壌の物理性（Journal of the Japanese Society「土壌の物理性」投稿規定
of Soil Physics）」並びにその他の印刷物の発行
1. 筆頭著者は本学会会員であること．
（3）内外の研究，技術の交流および他の学会，諸団体との協力
2. 投稿原稿の区分と規定ページは下記による．
（4）優れた論文に対する表彰
1）「論文」（Original papers）：独創性があり，土壌（5）その他，本学会の目的を達成するために必要な事業
第 4 条 本学会の会員は，正会員，シニア会員（これまで正会員の物理性に関する研究および技術の進歩に寄与するであって 60 歳以上の常勤職に就いていない者），学生会と見なされるもの．一編ごとに論文としての構成を
員および賛助会員，購読会員の 5 種とする．学会への入整えていて，他誌に未発表のものに限る．規定ペー会または退会の際は，本学会に届け出なければならない．
ジを刷り上がり 8 ページ以内とする．
第 5 条 会員は次の会費を所定の期日までに納めるものとする．
2）「研究ノート」（Notes）：土壌の物理性に関する新
しい事実や研究方法の改良などの短い報告で，独創
正会員
年額
5,500 円性があり他誌に未発表のもの．規定ページを刷り上
シニア会員
年額
3,000 円がり 6 ページ以内とする．
学生会員
年額
3,000 円3）「総説」（Reviews）：土壌の物理性に関する主題に（大学院生を含む）
ついて，それまでの研究を総括し，今後の発展方向
賛助会員
一口年額
22,500 円を展望するもの．規定ページを刷り上がり 10 ペー
購読会員
会誌年額
7,500 円ジ以内とする．
4）「解説」（Lectures）：土壌の物理性に関する諸事項
広告料賛助会員
15,000 円の理解を計るための平易な解説，ならびに研究技術賛助会員以外22,500 円の普及交換を進めるための紹介など．規定ページを
第 6 条 本学会に次の役員をおく．任期は 2 年とする．ただし，刷り上がり 8 ページ以内とする．
3 期連続の重任は認めない．選出方法は次による．
5）「研究紹介」（My research）：現在行っている研究，（1）会長 1 名，副会長 1 名，編集委員長 1 名
他の学会誌等に掲載された論文，研究上の工夫など
正会員の中から評議員会によって選出される．
についての平易な紹介．規定ページを刷り上がり 6（2）評議員
ページ以内とする．
イ 15 名 正会員の中から選挙によって選出される．
6）「資料」（Miscellaneous）：広く会員の参考となる
ロ 3 名以内を会長が正会員の中から委嘱する．
報告や情報．規定ページを刷り上がり 6 ページ以内（3）会計監査 2 名とする．
正会員の中から評議員会によって選出される．
7）「土粒子」（Readers’ column）：「土壌の物理性」（4）事務局長および幹事若干名
周辺のことについての随想や意見など．規定ページ
事務局長および若干名の庶務，会計，編集の幹事を会長を刷り上がり 2 ページ以内とする．
が委嘱する．
8）「講座」（Lecture series）：「解説」と本質的な違い
第 7 条 会長は毎年 1 回以上総会並びに評議員会を招集する．
はないが，特に編集委員会が企画して複数回にわ（1）総会は本学会の最高議決機関であって，事業計画，予算 ·たって行うものを指す．規定ページを刷り上がり 8
決算，会則の変更等を決定する．総会は正会員，シニア
ページ以内とする．
会員および学生会員で構成される．
9）「特集」（Topics）
：
「解説」
，
「総説」
，
「論文」
，
「研究
（2）評議員会は，会員の選挙によらない役員の選出，総会議
ノート」が混在したもので，特に編集委員会が企画案，その他の重要事項を検討する．評議員会の議長は会して複数回にわたって行うものを指す．規定ページ
長がこれにあたる．
は上記 1) ∼ 4) に同じ．
（3）総会および評議員会の議決は委任状を含む多数決とする．
10）「書評」（Book reviews）：規定ページを刷り上がただし，会則の変更は，総会出席者の 3 分の 2 以上の賛
り 2 ページ以内とする．
成を要する．
11）「巻頭言」（Foreword）：規定ページを刷り上がり
第 8 条 本学会の会務は，評議員会および総会の議を経て，会長，
2 ページ以内とする．
副会長，編集委員長および幹事がこれを執行する．会務
なお，上記区分において，編集委員会が会員または執行に必要な規程は別に定める．
会員以外に投稿依頼することがある．
第 9 条 本学会に次の委員会をおく．
3. 原稿の採否は編集委員会が決定する．編集委員会は，「論（1）選挙管理委員会
文」および「研究ノート」については 2 名，その他の投稿
正会員の中から評議員会によって選出され，本学会の評原稿には 1 名の閲読者を選んで閲読を依頼する．ただし，議員選挙を管理する．
「土粒子」
，
「書評」
，
「巻頭言」については，この限りでない．
（2）編集委員会
4. 投稿原稿が 2．に記す規定ページを越える場合には，1
正会員の中から評議員会によって選出される委員によっページあたり 10,000 円を著者負担とする．ただし，編て構成され，学会誌その他の印刷物の編集に当たる．
集委員会が会員または会員以外に投稿依頼した原稿の場（3）学会賞選考委員会
合はこの限りでない．また，カラー印刷の場合は実費を
正会員の中から評議員会によって選出される委員によっ著者負担とする．
て構成され，学会賞を選考する．
5. 投稿する場合は、投稿原稿（本文，図表含む）および原
第 10 条 本学会の経費は会費その他の収入を持ってあてる．本学稿送り状（別紙書式）の PDF ファイルを編集委員長宛会の会計年度は毎年 4 月 1 日に始まり翌年 3 月 31 日に提出する．ただし、編集委員会が別途指示した場合はに終わる．
この限りではない．
附則
6. 投稿原稿を受け付けた場合，編集委員会は速やかに受付番（1）本学会の連絡先は，これを当分の間次に置く．
号と受付日付を明示した電子メールまたは書状を投稿者に〒060-8589 北海道札幌市北区北 9 条西 9 丁目
送る．また，投稿原稿の掲載を決定した場合，最終原稿を北海道大学大学院農学研究院 生物環境工学分野受領したのちに受理日付を明示した書状を著者に送付する．
土壌保全学研究室内（2017.10.14 改正） 7. 著者に対して修正を求めた原稿が，特別な事情なく 6 ヶ月以内に再提出されないときは取り下げと見なす．
細則
8.
著者校正は初校のみとする．また，特に申し出がない限（1）シニア会員
り印刷ずみの原稿は返さない．
1）シニア会員は，これまで正会員であって 60 歳以上の常9.
掲載された論文等の PDF ファイルを著者に贈呈する．ま勤職に就いていない者
た，著者は別刷りを印刷業者から購入することができる．

10. 掲載された論文等は土壌物理学会ホームページ上にも掲

引用文献載する．論文，研究ノート，総説，解説，資料，講座に 24. 本文中で文献を引用する場合は，著者名（西暦），まついては，掲載誌発行後一年間は土壌物理学会員のみがたは（著者名，西暦）とする（例：Tanaka（1994），
閲覧できる．ただし，著者はオープンアクセス（50,000（Nakamura and Suzuki, 1981））．著者が 3 名以上の場
円：ページ制限なし）を選択することで，ホームページ
合は第 1 著者名のあとに—ら，または et al. をつける
掲載直後から誰もが閲覧可能な状態とすることができる．
（例：山田ら（1978），（Tayler et al., 1945））．
11. 本誌に掲載された論文等の著作権は土壌物理学会に帰属する． 25. 引用文献は，本文のあとに著者名の ABC 順に書く．未発表，私信は引用文献として記載しない．
（2017.10.14 改正）
26. 引用文献の書き方は以下のようにする．
「土壌の物理性」原稿執筆要領< 雑誌 >
森本 聡, 永田 修, 川本 健, 長谷川周一 (2009): 泥炭
全般的事項
林土壌の温室効果ガスの生成と消失．土壌の物理性，
1. 原稿は，日本語または英語で執筆する．
113: 3-12.
2. 日本語の原稿は，A4 判，25 字，26 行とし，上下左右Fujimaki, H., Ando, Y., Cui, Y. and Inoue, M. (2008):に 30 mm 以上の余白をとる（4 枚で 1 ページになる）．
Parameter estimation of root water uptake model
3. 日本語の投稿原稿（ただし，書評は除く）は，はじめに
under salinity stress. Vadose Zone Journal, 7: 31日本語で表題，著者名を書き，そのあとに英語で表題と38.
著者名（フルネーム）を書いてから本文を書き始める．
Millington, R.J. and Quirk, J.P. (1959): Permeabil英語の表題は，冒頭の単語の頭文字のみを大文字とし，
ity of porous media. Nature, 183: 387-388, doi:
それ以降は固有名詞等，特に大文字が必要な場合を除き，10.1038/183387a0.
全て小文字とする．英語の著者名は名の 1 文字目およびBitteli, M., Flury, M. and Campbell, G.S. (2003): A
姓の全てを大文字で記す．また，第 1 ページ下の脚注
thermo dielectric analyzer to measure the freezに，著者の所属機関と所在地を英語で，Corresponding
ing and moisture characteristic of porous meAuthor とする著者の氏名と所属を日本語で記す．
dia. Water Resour. Res., 39(2): W1041, doi:
4. 日本語の論文および研究ノートは，英語による表題 · 著
10.1029/2001WR000930.
者名と本文との間に英語で 300 words 程度の Abstract< 単行本 >と 5 個以内の Key words を入れる．
Tinker, P.B. and Nye, P.H. (2000): Solute movement
5. 日本語の投稿原稿（ただし，資料と土粒子，書評を除く）
in the rhizosphere, p.308. Oxford University Press,
は，引用文献の後に日本語で 350 字程度の要旨，およびNew York.
5 個以内のキーワードを書く．
土壌環境分析法編集委員会編 (1997): 土壌環境分析法,
6. 英語の原稿は，A4 判，ダブルスペースで作成する（約pp.129-138. 博友社, 東京.
600 words で 1 ページになる）．
松森堅治 (2005): 地理情報システムを用いた窒素負荷予7. 英語の投稿原稿（ただし，書評は除く）は，Manuscript
測モデル. 波夛野隆介 · 犬伏和之編 続 · 環境負荷を
title，Full name of authors を書いてから本文を書き始予測する, pp. 60-79. 博友社, 東京.
める．著者名は First name の 1 文字目および Family
Jury, W.A. (1996): Stochastic solute transport modname の全てを大文字で記す．また，第 1 ページ下
eling trends and their potential compatibility withの脚注に，Address of institutions of authors およびGIS. In Corwin, D.L. and Loague, K. (ed.) AppliCorresponding author を記す．
cation of GIS to the modeling of non-point source
8. 英語の論文および研究ノートは，Full name of authors
pollutants in the vadose zone, pp.57-67. SSSA Speと本文の間に 300 words 以内の英文 Abstract と 5 個
cial Publication No. 48, Madison.
以内の英語の Key words を入れる．
<Web ページ >9. 英語の投稿原稿（ただし，資料と土粒子，書評は除く）
FAO-AGL (2000): Global extent and location of sodic
は，引用文献の後に，日本語で表題と著者名，所属機関，
soils. Available at http://www.fao.org/ag/AGL
所在地，350 字程度の要旨，5 個以内のキーワードをつ/agll/prosoil/sodic.html.（確認年月日）
ける．なお，著者が外国人のみのために日本語部分の作
野菜茶業研究所 (2006): 野菜の硝酸イオン低減化マ
成が困難な場合，英語部分からの和訳により，編集委員
ニュアル. http://vegetea.naro.affrc.go.jp/joho/manual会の責任でこれらを作成する．
/shousan/index.html.（確認年月日）
10. 英語の適正化はすべて著者の責任において行うこと．自（雑誌のページには
p. や pp. をつけない．本の単数
然科学論文の執筆や校正に経験のある英語のネイティブ
ページを示す場合には
p.，複数ページを示す場合にはスピーカーによってチェックされている事が望ましい．
pp. をつける）
11. 閲読終了後に最終原稿を提出するときには，編集委員会
その他が指示する形式の電子ファイルを送付する．
27. 本要項のほか，科学技術情報流通技術基準（SIST）に準
本文
拠する．
12. 原稿には，閲読作業の効率化を図るため，本文に行番号（2012.12.20 改正）
を 5 行毎に書き入れる．
13. 章の見出しは 1.，2.，3.，節の見出しは 1.1，1.2，1.3，[土壌物理分野の SI 使用例]
項の見出しは 1.1.1，1.1.2，1.1.3 のようにし，字体は国際単位系（SI）については多くの出版物で解説されてゴシック体とする．
いるので，ここには，土壌物理に特に関連の深い単位の使
14. 略語，略号を使うときには，はじめにそれが出たときに正用例を示す．
式名称とともに示す（例：陽イオン交換容量（CEC））．
推奨される単位認められる単位15. 単位の使い方は，国際単位系（SI）による．
含水比
kg kg−1*
図表
体積含水率
m3 m−3
16. 日本語および英語の原稿ともに，図または写真の番号は
体積
m3
L, cm3 , mL
Fig. 1，表は Table 2 のように書く．
密度kg m−3 , Mg m−3
17. 日本語の論文，研究ノートの場合，図，表，写真の題名
濃度kg m−3 , Mg m−3
g L−1 , mg L−1には英語を併記する．
モル濃度
mol m−3
mol L−1
18. 図，表，写真は必要最小限度とし，同一事項を表と図に
力N重複して表示する事は避ける．
水分ポテンシャル J kg−1
kPa, m, cm, **19. 図，表，写真は 1 ページに 1 つずつ貼り付け，余白に図
透水係数
m s−1
cm s−1−1−1
イオン交換容量 molc kg , cmolc kg
表等の番号を記す．
s−1
min−1
20. 図と写真はそのまま製版されるので，鮮明なものを提出する． 回転数
収量kg m−2kg ha−1, Mg ha−1, t ha−1
21. 図は 1/2 倍などとおよその縮小率を指定し，2 倍くらい
* % 表示は，原則として分母と分子の内容が同じものの場に大きく書く．図の線の太さや文字の大きさは，刷り上合にだけ用いる．
がりの大きさを考慮して決める．
22. 地図や顕微鏡写真には定尺をつけ，縮尺や倍率はつけない． ** cm H2 O は使わない．
23. 図，表，写真の番号と題名は，別紙にまとめて一覧表とする．

その他，疑問のあるものについては意見を聞いて編集委
員会で検討する．

（2000. 3.10 改正）「土壌の物理性」投稿案内

1. 投稿原稿に対する基本方針

編集委員会は常に会員からの自主投稿を歓迎すると共に，投稿された原稿をなるべく早く，できるだけ掲載する方向で処理したいと考えています．ただし，「土壌の物理性」に掲載する原稿は読者にとって有益かつ完成したものであることが必要です．編集委員会は，投稿原稿が「土壌の物理性」に掲載するのにふさわしい内容を持つ原稿であるかどうかを基準として，掲載の可否を判断します．なお，他誌との二重投稿は認められません．
2. 投稿原稿に求められるもの
投稿に当たっては，「投稿案内」と共に「投稿規定」「原稿執筆要領」をよく読んでそれに従うと共に，閲読者が多忙な中で時間を割くことをよく理解して最大限に努力して完成した原稿としてから投稿して下さい．投稿原稿に求められるものは下の通りです．
（1） 内容について
1）オリジナリティーがあること（論文，研究ノートの場合）．
2）本誌の読者にとって有益であること．
3）導入，研究手法から結論に至る過程に論理性がある
こと．
（2） 表現について
1）原稿全体の構成が整っていること．
2）表題が適切であること．
3）文章表現が適切であること．
4）英文アブストラクトや和文要旨が適切であること．
5）SI 単位や記号の使い方が適切であること．
6）キーワードが適切であること．
7）英語の表現が適切であること．
3. 投稿原稿の受付
編集委員会は投稿原稿が体裁上最小限必要とする条件を満たしていることを確認した後，原稿を受け付けます．
条件を満たしていない場合には，著者への問い合わせ，もしくは原稿を受け付けずそのまま返稿する場合があります．
4. 閲読の目的と姿勢
編集委員会は受け付けた投稿原稿が「土壌の物理性」
へ掲載するのにふさわしいかどうかを判定する資料を得るために閲読を行います．編集委員会は適任の閲読者を選定し，閲読を依頼します． 編集委員会は閲読者に対して特に次の点について注意をお願いしています．
1）原稿の内容に対する責任は著者にあること．また，原稿の価値は読者が判断するものであること．
2）閲読は内容の批判や討議を行うものではないこと．
従って，見解の相違は掲載不適の理由にならないこと．
3）閲読者は投稿者に対して研究指導する立場にないこと．従って，新たな実験や計算の追加要求などは極力避けること．
4）閲読依頼期限を厳守すること．
閲読者から出された問題点について，編集委員会
は原稿が投稿原稿に求められるものを満たすために
必要に応じて著者に修正を求めます．ただし，これは著者に対する指導や助言ではありません．原稿の内容に対する責任はあくまで著者にあります．
5. 閲読の手順と投稿原稿の採否（1） 編集委員会は論文，研究ノートについては 2 人，その他の区分の投稿原稿については 1 人の適任者へ閲読を依頼します．
（2） 閲読者は「2」の投稿原稿に求められる各項目について投稿原稿を評価し，閲読結果を全体意見および個別指摘事項と共に編集委員会へ報告します．
（3） 編集委員会は閲読結果を審議し，投稿原稿の採否を決定します．この際，論文・研究ノートについては原則として下のように取り扱います．
1）閲読者 2 名とも総合判定が A（このまま掲載可）の場合には，閲読結果を受け取った日付で原稿を受理（掲載決定）します．
2）1 名の総合判定は A ，もう 1 名は B（指摘の点が修正されれば掲載可）の場合，または 2 名とも Bの場合には，指摘事項について著者に比較的小さな修正を求めます．修正原稿が提出され，編集委員会
が修正内容を確認した時点で受理します．

3）1 名の総合判定は B ，もう 1 名は C（指摘の点が修正された後，再閲読する必要あり）の場合，または 2 名とも C の場合には，指摘事項について著者にやや大きな修正を求めます．提出された修正原稿
は再び閲読を行います．
4）1 名の総合判定が D（掲載不適）の場合には，第 3の閲読者に判定を依頼することを含め，編集委員会
が対応を判断します．
5）2 名の総合判定とも D の場合には，その理由を付けて著者へ返稿します．
6）閲読は 2 回までとし，再閲読においても C または
D の判定がある場合には，編集委員会が対応を判断します．
（4）閲読結果の総合判定を D とする理由は次のようなものです．
1）既発表とみなされるもの．
2）重要な部分に根本的な誤りがあるもの．
3）事実に基づいた内容ではなく，単なる主観が述べられているに過ぎないもの（論文，研究ノートの場合）．
4）通説が述べられているだけで，新しい知見が全くな
いもの（論文，研究ノートの場合）．
5）修正を要する根本的な指摘事項を多く含んでいるもの．
6）連載形式で構成されており，1 編ごとに独立したものと認めがたいもの．
7）全体構成や文章が粗雑で，内容の判定ができないもの．
8）内容が本誌の趣旨に合わないもの．
9）その他，本誌に掲載する原稿としてふさわしくない
もの．
6. その他
1）共著者は全員が原稿に責任を負って下さい．共著者がいるにもかかわらず，その共著者が原稿を読んでないと判断されるような基本的なミスがある場合には，閲読を中止して返稿することがあります．
2）修正原稿を返送する場合には，修正結果がすぐに分かるような回答書（指摘事項に対応した修正箇所・内容一覧）を作成し，添付して下さい．
3）原稿は一度印刷公表されると末永く評価と批判の対象になります．閲読意見に対しては極力冷静かつ慎重に対処して下さい．誤解があった場合には，その誤解を解くだけでなく，誤解を受けないような書き方に改める努力が求められます．
4）「掲載不適」で返稿された原稿を再投稿する場合は新規投稿扱いとなります．「掲載不適」とした理由を十分に確認・検討の上，再投稿して下さい．
（2004. 6.29 改正）「土壌の物理性」閲読の手引き

1. 閲読依頼期限の厳守をお願い致します．万一，期限内に閲読が困難な場合には，速やかに編集委員会事務局へ連絡して下さい．
2. 「投稿案内」において，下のものを投稿原稿に求めています．閲読においても，これらをその基準とします．
（1）内容について
1）オリジナリティーがあること（論文，研究ノートの場合）．
2）本誌の読者にとって有益であること．
3）導入，研究手法から結論に至る過程に論理性がある
こと．
（2）表現について
1）原稿全体の構成が整っていること．
2）文章表現が適切であること．
3）表題が適切であること．
4）英文アブストラクトや和文要旨が適切であること．
5）SI 単位や記号の使い方が適切であること．
6）キーワードが適切であること．
7）英語の表現が適切であること．
3. 「掲載不適」とする理由としては，下のようなものが挙げられます．
1）既発表とみなされるもの．
2）重要な部分に根本的な誤りがあるもの．
3）事実に基づいた内容ではなく，単なる主観が述べられているに過ぎないもの（論文，研究ノートの場合）．
4）通説が述べられているだけで，新しい知見が全くな
いもの（論文，研究ノートの場合）．
5）修正を要する根本的な指摘事項を多く含んでいるもの．

6）連載形式で構成されており，1 編ごとに独立したものと認めがたいもの．
7）全体構成や文章が粗雑で，内容の判定ができないもの．
8）内容が本誌の趣旨に合わないもの．
9）その他，本誌に掲載する原稿としてふさわしくない
もの．
4. 原稿の内容に対する責任はあくまで著者にあります．また，原稿の価値は読者が判断するものです．
5. 閲読は内容の批判や討議を行うものではありません．
従って，見解の相違は掲載不適の理由にはなりません．
6. 閲読者は投稿者に対して研究指導する立場にありません．
よって，新たな実験や計算の追加要求などは極力避けて下さい．
7. 閲読者は著者と直接折衝をしないで下さい．著者との交渉はすべて編集委員会の責任で行います．
8. 閲読者として知り得た印刷前の原稿の内容はいかなる形でも利用・引用等出来ません．厳にご注意下さい．
9. 著者に修正を求めた後に再閲読で総合判定を D（掲載不適）とすることは，閲読者にとっても著者にとっても非生産的です．完成度の低い原稿などにおいては，修正を求め再閲読を行うよりも「掲載不適」として練り直してから新たに投稿するよう促した方がよい場合があります．
その場合には，「掲載不適」の理由（「3」参照）を明示して下さい．個別指摘事項への意見は不要です．
10. 閲読者から編集委員会へ意見等がある場合は，必ず「編集委員会への意見」として別紙への記入をお願いします．
（2003. 5.28 改正）
土壌物理学会学会賞規定

1. 種類

学会賞に論文賞と優秀ポスター賞をおく．

2. 論文賞の選考

論文賞は，年 1 回、以下の手順により論文賞選考委員会「以下，選考委員会」で審査し，評議員会で決定する．
（1）賞の対象とする論文は，「土壌の物理性」に，原則として前年度（前年 4 月 1 日から翌年の 3 月 31 日まで）に掲載された原著「論文」とする．
（2）正会員（1 名以上）は, 学会賞受賞候補として（1）に該当する論文を自薦もしくは他薦することができる.
（3）（2）に記す論文の推薦は，学会誌「土壌の物理性」に公示された募集を受け，所定の書類を必要部数揃え，

期日までに学会事務局に提出しなければならない．
（4）選考委員会では，（1）に該当する論文の中から，（2）の推薦も考慮して受賞候補論文（以下，候補論文）を審査 · 選考し，評議員会へ付議する．
（5）評議員会では，候補論文が受賞に適格かどうかを判定する．
3. 優秀ポスター賞の選考
優秀ポスター賞（以下，ポスター賞）は年 1 回，以下の手順により決定する．
（1）選考の対象は，当該年度の土壌物理学会大会において発表されたポスターであって，筆頭者は会員とする．
ただし，ポスターの筆頭者が学生の場合に限り，会員外であっても選考の対象に含める．
（2）ポスター賞は（1）に該当するポスターの中から，選考委員会および会員の投票により選考する．（5 件以下）
4. 表彰
（1）論文賞は，通常総会において表彰を行う．また会誌上に公表する．
（2）ポスター賞は，土壌物理学会大会において表彰を行う．また，会誌上に公表する．
（2002.11.23 制定 2013.10.26 改正）
土壌物理学会学会賞選考委員会規定

1. 目的

学会賞選考委員会（以下委員会という）は，論文賞および優秀ポスター賞の選考を行うことを目的とする．
2. 業務
委員会はつぎの業務を行う．
（1）学会賞規定に基づく論文賞および優秀ポスター賞の選考．
（2）土壌物理学会学会賞選考に関するその他の事項．
3. 構成
委員会は委員長 1 名および委員 5 程度を持って構成する．
（1）委員会は正会員より構成される．
（2）委員は評議員会の推薦に基づき総会で決定する．
（3）任期は 2 年とする．ただし，再任はさまたげない．
（4）委員長は委員の互選により選出する．
（5） 若干名の専門委員をおくことを認める．
4. 開催
委員会は年 3 回程度開催する．
（2002.11.23 制定 2003.11.22 改正）

原 稿 送 り 状
受付番号（本会で記入）：
発送年月日

年

区 分

論文

月

日

研究ノート

受付年月日（本会で記入）
総説

解説

資料

年

研究紹介

土粒子

月

日

書評

表題（和文）

表題（英文）

キーワード（日本語：５語以内）

キーワード（英語：５語以内）

著者名 1
ローマ字 ∗

所属

本会

会員 · 非会員

著者名 2
ローマ字 ∗

所属

本会

会員 · 非会員

著者名 3
ローマ字 ∗

所属

本会

会員 · 非会員

著者名 4
ローマ字 ∗

所属

本会

会員 · 非会員

著者名 5
ローマ字 ∗

所属

本会

会員 · 非会員

Corresponding
author の〒住 所・氏 名
電話 ·FAX 番号

TEL. ：

E-mailE-mail ：本文原稿枚数
表
別刷

FAX.：

A4：（
枚

図
（カラー図番号．

枚）
枚

希望しない · 希望する（）

写真（カラー写真番号．
部）

付記：投稿の宛先は土壌物理学会ホームページを参照のこと（http://js-soilphysics.com/）
著者には pdf ファイルを贈呈します．別刷を希望する方は印刷業者から購入することになります．
*著者名は漢字とローマ字の両方を記入してください．

枚）土壌物理学会会員登録用紙 [入会 · 退会 · 登録事項変更]
必要事項をご記入の上，このまま下記宛先に（郵便，FAX，E-mail の添付ファイルなどで）お送
りください。
届出内容

入会

申込年月日

西暦
正 ·

会員種別

· 退会 · 登録事項変更
年

シニア

月

日

· 学生 ·

購読 ·

賛助

氏名
同上ローマ字読み性別

男

生年月（西暦年 · 月）

· 女

西暦

年

勤務先

· 自宅

月

勤務先
勤務先所在地〒

勤務先電話番号 / FAX 番号
勤務先 e-mail自宅住所〒

自宅電話番号 / FAX 番号
自宅 e-mail
会誌の送付先本学会以外の主な所属学会
学生会員のみ記載（本人に連絡が取れないときの連絡先として）指導教官名指導教官 e-mail

本会記入
受付年月日：
承認年月日：
送付先 · 問い合わせ先
土壌物理学会事務局（庶務幹事）〒 060-8589 北海道札幌市北区北 9 条西 9 丁目
北海道大学大学院農学研究院 生物環境工学分野

Tel. 011-706-3641 Fax 011-706-2494
E-mail spsyomu@ml.affrc.go.jp

土壌保全学研究室内

柏木淳一

宛

複写される方に本誌に掲載された著作物を複写したい方は，（社）日本複写権センターと包括複写許諾契約を締結されている企業の従業員以外は，著作権者から複写権等の行使の委託を受けている次の団体から許諾を受けて下さい．著作物の転載 · 翻訳のような複写以外の許諾は，直接本会へご連絡下さい．

〒 107-0052

東京都港区赤坂 9-6-41

TEL：03-3475-5618

乃木坂ビル

FAX：03-3475-5619

学術著作権協会

E-mail：kammori@msh.biglobe.ne.jp

Notice about Photocopying
In order to photocopy any work from this publication, you or your organization must obtain permission from thefollowing organization which has been delegated for copyright for clearance by the copyright owner of this publication.
Japan Academic Association for Copyright Clearance (JAACC)
41-6 Akasaka 9-chome, Minato-ku, Tokyo 107-0052, Japan
TEL：81-3-3475-5618

FAX：81-3-3475-5619

E-mail：kammori@msh.biglobe.ne.jp

土壌の物理性 第 139 号（会員配布）
発行 土 壌 物 理 学 会

2018 年 7 月 20 日発行〒 060-8589

北海道札幌市北区北 9 条西 9 丁目 北海道大学大学院農学研究院
生物環境工学分野 土壌保全学研究室内 柏木淳一 宛

電011-706-3641
FAX 011-706-2494
spsyomu ＠ ml.affrc.go.jp
http://js-soilphysics.com/ent
北洋銀行北七条支店（店番 312）
普通 4093179 土壌物理学会 会計 倉持寛太
口座番号：01350-2-40943

話

E-mail
URL
銀行口座
郵便振替

加入者名：土壌物理学会
編集委員会事務局（投稿原稿送付先）
所在地はホームページ（http://js-soilphysics.com/ent）で最新の物を確認して下さい．

E-mail印

kibyosi@ml.affrc.go.jp
刷

創文印刷工業株式会社〒 116-0011 東京都荒川区西尾久 7-12-16

ISSN 0387-6012

Journal of the Japanese Society of Soil Physics
No. 139

July, 2018

Contents........................................................

Foreword

H. TAKEUCHI...

1

K. KAMEYAMA, Y. IWATA, T. MIYAMOTO and K. SASAKI

...

3

...

15

Special section for Upland Irrigation
Original Papers
Optimization of irrigation regime in a sand dune soil field using HYDRUS-1D............

Notes
Handmade Becquerel meter using commercial scintillation survey meter
S. YOSHIDA, S. SHIOZAWA, T. ZAITSU, H. YAMANO and K. NISHIDA

Readers’ Column............................................................

...

23

Announcements..............................................................................

25

Editor’s Postscript..............................................................................

27

H. INOUE

Published by

Japanese Society of Soil Physics
Research Faculty of Agriculture, Hokkaido University
Kita 9, Nishi 9, Kita-ku, Sapporo, Hokkaido 060-8589 Japan

http://js-soilphysics.com/