Silviase/jgov_v1.4_full_ann huggingface dataset:
### 1. はじめに

土壤には微生物とともに多数の土壤動物が生息している $ ^{1)2)} $ 。土壤動物は陸域生態系の動物群集の5から9割程度のバイオマスを占めており、陸域の物質循環の大きなフローを制御している。たとえば、土壤動物のバイオマスが増えることで、植物の生長が増大することが明らかにされている $ ^{2)3)4)} $ 。土壤微生物と土壤動物は有機物の分解により栄養塩を植物に再利用できるようにしているが、土壤動物は土壤構造を物理的に変えることにより、植物の成長にも影響している（図(5)-1）。また、土壤動物は地上動物の餌として分解系と生食連鎖系を結びつけている。たとえば、森林性の鳥類は繁殖期の餌としてミミズに大きく依存しており、両生類やは虫類、そして昆虫類もミミズを餌資源として利用している。地下部と地上部の生物群集の応答の研究は、陸域生態系の動態を把握するために緊急に必要とされている $ ^{5)} $ 。

日本列島における環境変動が生態系に与える影響を明らかにするには、土壌動物の地理分布を調べ、環境変動に対する応答を予測する必要がある。そこで、大型ミミズ類の地理分布に関する整理を行い、地上部食物網との関係を評価するための基礎データを得た。また、詳細で精度の高い既存の地理分布データがあるカマアシム類（原尾目、Protura）の地理分布を環境条件との関係で解析した。

陸域生態系の気候変動への応答

図(5)-1 本研究の概念図

## 2 ．研究開発目的

## （1） ミミズ群集と地理的分布

モニ1000（森林）では地表徘徊性甲虫が森林の状態を表現する生物指標であると考え、コアサイトでピットホール調査が実施されている。これまでに、「日本では北へいくほど地表徘徊性甲虫が多く捕獲され、種数も増大する」という地理分布パタンが得られている $ ^{6} $ 。これは、多くの生物群で北（高緯度）になるほど多様性が低下することが知られている（図(5)-2）のに対して、対照的なパタンとなっている。そこで、この地理分布に多くの甲虫類の潜在的な餌であ
るミミズの分布が大きく影響するという仮説を立て、モニ1000の調査地でミミズ群集調査を行った。

## （２）日本列島におけるカマアシムシ類の分布と環境条件との関係

カマアシムシ類は、体長1～2mmの六脚類(Hexapoda)に属する節足動物で、土壌では中型土壌動物に区分されている。おなじ中型土壌動物のトビムシ、ササラダニ類と比べると個体数は少ないが、世界各地の土壌に出現する。食性は主に土壌中の糸状菌を食べると考えられており、小型でクチクラの発達が十分ではない体を持つため、地表面を活発に移動することはない。日本では、北海道から与那国島まで、日本全国の広い範囲から多数の記録が得られ、すでに公表されており、シベリア、中国、台湾との種の共通性から、日本列島への侵入経路が検討されてきた。種の同定はカマアシムシの分類研究者である故今立源太郎博士と中村修美氏の2名のみが行ったため、分類学上の見解の相違や、同定ミスは少ないことが期待できる。土壌では、糸状菌を摂食し、菌根菌を選択的に摂食するという見解もある。これらのことからカマアシムシ類の地理分布と植生や気候条件との関係を解析することで、気候変動が土壌生物の分布に与える影響について予測することができる $ ^{7} $ 。

### 3. 研究開発方法

（１）ピットホールトラップによる地表徘徊性甲虫調査

モニ1000の森林調査区のうち、苫小牧、カヌマ沢、佐渡（大佐渡、豊岡）、綾、田野、与那の各調査地でピットホールトラップ調査と土壌動物の掘り取り調査を行った。

ピットホールトラップ法は、一辺約2mの方形区を設定し、四隅にピットホールトラップ（口径84mm、深さ116mm）を、開口部が地表面と同じ高さになるように埋設した。トラップの底面には、直径約4mmの排水孔を2つ開けた。雨水の浸入を抑えるために、トラップの直上5cm程度の高さに15cm×15cmの板を設置した。設置から3日後に、トラップ内に落下したすべての動物を、酢酸エチルを充満させた容器に回収した。回収した動物はエタノール中に保存し、甲虫類は種まで、その他の動物は目または綱まで同定した。

掘り取りは、まず25cm×25cmの範囲のリターを採取し、乾燥重を測定した。リターと土壌（深さ10cm）に生息する大型土壌動物を現地でのハンドソーティングで採集した。ハンドソーティングは目視で、ピンセットと吸虫管を用い、すべての大型土壌動物を捕獲し、ミミズ以外はアルコールに保存した。ミミズは持ち帰り、ホルマリンで固定し、体重を測定後、すべて種名を同定した。

捕食者とその餌となる被食者との種間関係の強さは、連結度(Connectance)と呼ばれている。本調査では列島スケールでミミズのデータが得られているが、各サイトで正確にミミズとその捕食者との関係を確認することは困難である。また、ミミズ群集は限られた掘り取り調査では十分に把握できていない。そこで、各サイトで2009年から2011年の3年間にモニタリング1000によって捕獲された地表徘徊性甲虫のリストから総種数を得た。これらの甲虫の食性は、すべて肉食性もしくは雑食性であり、基本的にミミズの生体または、死体を摂食すると考えた。また、各サイトの半径50kmの地点からこれまでに報告された大型ミミズ類のリストを作成し、機能群データベースを作成した。ミミズを被食者、地表徘徊性甲虫を捕食者と考え、各地点での最大
となる種間関係数（ミミズ種数 X 捕食者種数）を、日本全体の種間関係数で割ったものを連結度とみなした $ ^{8)} $ 。

## （２）カマアシムシ類地理分布解析

## 1 ）データセット

全国のカマアシムシ類の群集解析を行うため、2007年までに出版された日本におけるカマアシムシの分類学的論文や採集記録の報告をもとに、採集場所や植生、標高が明らかである地点について、種の在/不在データを収集した。一方で、同定者自身が採集したサンプルだけではなく、他の小型節足土壌動物研究者が抽出・寄贈した標本も用いられているため、採集者・採集努力は一定ではない。ウェノカマアシムシNipponentomon uenoi Imadate et Yosii, 1959は今立(1988) $ ^{9)} $ により6つの型が認められ、それぞれの型が分類体系内のいかなる階級(taxon)に属するかを断定するだけの根拠はないものの、相補的な地理分布を示すため、独立したtaxonとして扱った。いくつかの採集記録は、同定に疑問があるために除外した。高山帯植生、浜辺、ササ原、コケ、道ばたという植生条件の記録は、地点数が少ないために除外した。また、藪地は森林や草原との定義上の境界が不明瞭であることから、記録を除外した。複数の植生タイプにまたがる地点も存在したが、このうち筆頭植生条件を植生条件とした。緯度経度は世界測地系(WGS84)を用いて記述し、緯度経度が明記されていない記録については、採集地点の住所からアドレスジオコーディングを用いて座標値を得た。

## 2 ）地域ごとの種組成の類型と類似度

採集記録の地理的な偏りを明らかにするため、行政区分と生物地理学的重要性を加味して、全国を10地域に分割した。この際、対馬・壱岐は朝鮮半島と共通の生物群が、琉球列島は固有種や中国・台湾と共通の生物群が多く独特な生物相を持つことが知られているため、九州とするのではなく、対馬、琉球として独立した地域とした。それぞれの面積あたり採集地点数を算出し、採集努力量とした。また、地域ごとの採集努力量を評価するために、EstimateS ver. 8.2.0 のChao1指数を用いて総出現種数を予測した $ ^{10} $ 。種多様性の高い地域を明らかにするため、それぞれの地域の面積あたり予測総出現種数を算出した。

1地点あたりの出現種数が多い地域を明らかにするため、採集地点数の頻度を算出した。1地点あたりの出現種数が多い地域の地理的偏在を明らかにするため、ArcGIS 9.3（ESRI Japan, Tokyo, Japan）のSpatial analystツールにより、Inverse Distance Weighted（IDW）手法を用いて各地点の出現種数からサーフェスを内挿してラスタを作成した。この時、出力セルサイズを0.005、乗数を3、検索範囲を可変、検索半径のポイント数を50とした。

全国10地域の種組成を比較するため、各地域ごとの種の在／不在データを集計した。このデータをもとに、PC-ORD $ ^{11} $ のTwinspan $ ^{12} $ を用いて、地域ごとの種組成を類型化した。類型化された地域間の種組成の類似度を算出するため、統計ソフトR ver. 2.13.0 $ ^{13} $ を用いてJaccard指数を算出した。

## 3 ）共出現種の解析

共出現種を明らかにするために、出現種数の半数である6種以上が出現した地点の種組成を元
にWard法を用いてクラスター分析を行った。6種以上出現する地点は比較的攪乱の少ない、安定した環境であると考えられ、人為的影響を排除できると考えられる。このとき、解析には統計ソフトRのvegan、ade4、gclus、cluster、RcolorBrewer、labdsvの各packageを用いた。なお、これらの地点の地理的偏りを明らかにするために、表(5)-2と同様に各地域の面積あたりの採集地点数を算出した。

類型化された各グループの指標種Indval $ ^{14)} $ を明らかにするため、Borcard et al. (2011) $ ^{15)} $ に従って統計ソフトRにより算出した。また、植生、標高、気候（最高気温、最低気温、年平均降水量）を説明変数、各グループの有無を応答変数としたロジスティック回帰を行った。植生は広葉樹林を1、針葉樹林を0とした二値化を行った。各採集地点の気象条件（最高気温、最低気温、年平均降水量）は、メッシュ気候値2000 $ ^{16)} $ より取得した1kmメッシュ気候値を利用した。応答変数の選択については、ベストモデルだけでなく、ΔAICが2未満のモデルを全て考慮した。

Total Earthworm biomass (wet)

Epigeic Earthworm biomass (wet)

図(5)-2 モニ1000サイトで得られた樹木の種数と平均気温との関係（上左）、オサムシ類の種数と平均気温との関係（上右）、および7地点で得られたミミズ全体の現存量（下左）、および表層性種の現存量（下右）と年平均気温との関係。

## 4 ．結果及び考察

## （１）森林のミミズ現存量の地理分布と捕食者との関係

平均気温と現存量との関係を求めた(図(5)-2下)。全現存量は苫小牧で最も多く、豊岡で最も少なかったが、明瞭な関係はみられなかった。ミミズ群集は地表性(Epigeic)、地中性（Endogeic）、
そして表層採食地中性（Anecic）の3つの生活型機能群に分けられている $ ^{17} $ 。地表性がもっとも甲虫による捕食を受けるが、地表性種は苦小牧と綾で採集できたが、それ以外のサイトではほとんど採集できなかった。田野、綾には大型で、おそらく地表性のシーボルトミミズ(Metaphire sieboldii)が、そして与那には表層採食地中性のヤンバルオオフトミミズ（A. yambaruensis）がそれぞれ生息しているが、個体数が少なく、現存量の推定が難しいため、今回の掘り取り調査では十分に把握されていない。与那のヤンバルオオフトミミズでは暫定値で生息密度が0.5/m²、現存量が8 g/m²、という値を得ている（未発表データ）。この値を加えても、与那の値が低く、苦小牧の値が高いことは変わらない。苫小牧だけでなく、本州の温帯林におけるデータを増やす必要があるが、少なくとも北方でミミズ現存量が低くなる傾向は認められないと言える。
表(5)-1 森林モニタリングサイト1000調査地におけるミミズ相と生活型の構成。佐渡の2カ所のサイトは場所が近いため、ひとつにまとめた

和名	学名	調查地
苦小牧	カヌマ沢	佐渡	綾	田野	与那	生活型
フスジミミズ	Amynthas vittatus	○	○	○	○	○	表層
ヒトツモンミミズ	Metaphire hilgendorfi	○	○	○	○	○	表層
ハタケミミズ	Metaphire agrestis	○	○	○	○	○	表層
セグロミミズ	Amynthas divergens	○	○	○	地中
ユノシマミミズ	Amynthas yunoshimensis	○	○	表層
メガネミミズ	Metaphire acincta	○	○	地中
フキソクミミズ	Amynthas irregularis	○	○	○	○	○	表層
イロジロミミズ	Amynthas phaselus	○	地中
パライロツリミミズ	Allolobophora rosea	○	地中
ムラサキツリミミズ	Dendrobaena octaedra	○	表層
シマミミズ	Eisenia fetida	○	表層
サクラミミズ	Eisenia japonica	○	○	○	○	地中
フクロナシツリミミズ	Bimastos parvus	○	○	○	表層
クロイロツリミミズ	Aporrectodea caliginosa	○	○	地中
ヘンイセイミミズ	Amynthas heteropodus	○	○	○	○	○	地中
クソミミズ	Amynthas hupeiensis	○	○	○	地中
ヒナナミミズ	Amynthas micronarius	○	○	○	○	地中
アオキミミズ	Metaphire aokii	○	表層
ホソスジミミズ	Metaphire striata	○	表層
フツウミミズ	Metaphire communissima	○	○	○	表層
シーボルトミミズ	Metaphire sieboldi	○	○	表層
ミタマミミズ	Metaphire soulensis	○	○	表層
チョウミミズ	Amynthas glabrus	○	地中
ハワイミミズ	Amynthas gracilis	○	地中
チビミミズ	Amynthas minimus	○	地中
モリスミミズ	Amynthas morrisi	○	地中
ヨロイミミズ	Amynthas robustus	○	地中
ヤンバルオオフトミミズ	Amynthas yambaruensis	○	地中
アカシマフトミミズ	Amynthas kunigamiensis	○	表層採食地中
メキシコミミズ	Metaphire californica	○	表層採食地中
キクチミミズ	Metaphire schmardae	○	表層
Polypheretima elongata	○	地中
Ocnerodrilus occidentalis	○	情報無し
Dichogaster bolauni	○	情報無し

モニ1000のピットホール調査では、7箇所の調査地から合計75分類群のオサムシ科とホソクビゴミムシ科の甲虫が捕獲された。一方、モニ1000サイト周辺の採集記録から、それぞれの調査地付近の潜在的な分布種のリストを作成した（表(5)-1）。調査地周辺からはミミズが34種報告されていた。北海道と本州から九州にかけて広く分布する種がいる一方、沖縄のミミズ種組成は大きく異なっていることが明らかとなった。すなわち、ほとんどが地中性であり、地表性の種は少ない。また、表層採食地中性種は、本州にも別の種が分布するが、モニ1000サイトと
なった森林では、沖縄にだけ2種分布していた。

各調査地の甲虫とミミズの連結度は0.01から0.16の値をとり、温度との間に弱い負の相関関係が見られた $ \left(R_{\mathrm{adj}}=0.3944,\mathrm{p}=0.075\right) $ （図(5)-3）。各地に生息するミミズの種数は、苫小牧と与那でやや高い値が見られたが、地点間の大きな差がなく、連結度には各地の甲虫の種数の違いが影響していた。この連結度は実際にミミズが補食されたという観察に基づくものではなく、各地においても潜在的な数値を表していると考える方がよいが、それでも温度との負の関係が見られた点で注目に値する。地中性や、大型のヤンバルオオフトミミズのような表層採食地中性種よりも地表性種が容易に捕獲されるため、地表性甲虫の潜在的な餌としては、地表性ミミズの現存量を考える必要がある。シーボルトミミズの生活型は不明であるが、生重が20gを越す大きなものになると、甲虫による捕獲は困難であろう。日本の温暖な地域では、地表性のミミズが少なく、地表徘徊性甲虫の餌資源が少ないことが示唆される。この原因としては、温暖な地域で落葉の分解速度が速く、ミミズに十分な資源量がなく、地表性が少ない可能性があるが、掘り取り調査の結果はこの仮説を支持しない。連結度の低さを、甲虫群集がより南の地点ではミミズへの依存度が低いと考えると、地表徘徊性甲虫以外の捕食者がミミズ群集を制御している可能性が指摘できる。

図(5)-3 ミミズ-地表徘徊性甲虫群集の連結度と平均気温との関係（右）

これまで、日本では大型土壌動物全体の現存量は南方ほど多くなることが知られていた $ ^{3)} $ が、ミミズの現存量の地理分布について詳しく比較したデータはない。モニ1000で採用されているピットホールによる地表徘徊性の無脊椎動物のモニタリングデータは、森林の物質循環に大きな影響をもつミミズとその捕食者との食物連鎖関係を簡便に検討するためにすぐれた方法であるといえる。

## （２）カマアシムシ分布データ

植生や標高も明記されていた3103地点、71分類群（64種3亜種4型）の情報から合計7077種・地点の採集記録を得ることでき、ファイルメーカーのデータベースとして整理できた。琉球や対馬では面積当たりの調査努力が大きかったが、出現種数は少なく、おおむね面積が広いほど多くの種が分布していた（表(5)-2）。特に、関東と中部では他の地域より多くの種が分布していた。
表(5)-2 日本の10地域におけるカマアシムシ類分布調査状況

Hokkaido	Tohoku	Kanto	Chubu	Kinki	Chugoku	Shikoku	Kyushu	Tsushima	Ryukyu	Total
Numbers of sampling localities	480	339	743	568	299	171	116	296	34	64	3110
Area ( $ 100km^{2} $ )	846	669	324	726	273	319	183	422	7	23	3792
Ratio1 (N/Area)	0.6	0.5	2.3	0.8	1.1	0.5	0.6	0.7	4.9	2.8
Total number of species	27	30	43	45	22	17	20	20	12	9	71
N of Exp Spe* $ ^{1} $	27.0	32.0	46.0	47.3	32.0	17.0	27.5	20.5	12.5	9.0
Ratio2 (N of exp species/Area)	0.03	0.05	0.14	0.07	0.12	0.05	0.15	0.05	1.79	0.39

$ ^{*} $ Expected number of species calculated using by Chao1

これを元に、森林で一箇所から6種以上出現した143地点の情報を用いてTwinspanでグループ分けを行い、大きく地理分布の異なる3ループを得た（図(5)-4）。グループⅠ北海道から中部日本まで、グループⅡは東北から西日本に、グループⅢは関東から西日本まで分布していた（図(5)-4）。これら3グループの分布が重なる関東から東北南部が、日本でのカマアシムシ類の多様性が高い地点であった。これらのグループの分布を環境要因と一般化線形モデルを用いて関係を見たところ、グループⅠは最高、最低気温が低い場所、グループⅡは降水量が少ないところ、そしてグループⅢは最高、最低気温が高く、降水量が多いところに見られた。また、針葉樹と広葉樹という樹種の違いは、どのグループにも影響していなかった（表(5)-3）。

図(5)-4 日本列島のカマアシムシ群集の分布。左上は、Twinspanによる群集のクラスター分け。右下は、6種以上採取された地点の3グループの分布。

表(5)-3 カマアシムシ群集3グループの分布と環境要因との関係のGLM（一般化線型）モデルによる解析。vegeN：広葉樹、elev：標高、te_max：月最高気温、te_min：月最低気温、preci_an：年間降水量。AICの変化が2までのモデルと、フルモデルを表示。
Group	No. of localities	Model	Constant	Standardized regression coefficient	AIC	$ \Delta $ AIC
vegeN	elev	te_max	te_min	preci_an
Gr1	47	1	-0.87	-0.81	-0.59	-0.49	153.05	0.00
2	-0.86	-1.08	-0.77	153.85	0.80
3	-0.83	-0.55	-0.93	154.49	1.44
4	-0.86	-0.15	-0.85	-0.53	-0.53	154.53	1.48
5	-0.86	-0.23	-1.09	-0.78	154.62	1.56
6	-0.87	-0.13	-0.83	-0.67	-0.48	154.82	1.76
full	-0.86	-0.14	-0.10	-0.86	-0.60	-0.52	156.39	3.34
null	-0.72	183.9	30.85
Gr2	35	1	-1.60	-1.49	137.78	0.00
2	-1.60	0.19	-1.45	138.97	1.19
3	-1.61	0.50	-0.56	-1.12	139.10	1.32
4	-1.60	0.14	-1.46	139.32	1.55
5	-1.60	0.12	-1.48	139.44	1.66
6	-1.58	0.27	0.30	-1.36	139.64	1.87
7	-1.60	-0.07	-1.46	139.71	1.93
full	-1.60	0.16	0.17	0.55	-0.54	-1.07	141.94	4.16
null	-1.14	161.72	23.94
Gr3	62	1	-0.50	0.75	1.83	1.13	120.56	0.00
2	-0.50	-0.12	0.72	1.76	1.13	122.45	1.90
3	-0.50	0.07	0.77	1.82	1.14	122.47	1.91
full	-0.51	0.08	-0.14	0.75	1.73	1.15	124.34	3.79
null	-0.28	198.84	78.28

カマアシムシ類は、主に土壌中の糸状菌を摂食していると考えられており、糸状菌を餌に飼育可能である。一方で、菌根菌を摂食しているとの主張もあるが、この解析のように広葉樹と針葉樹で分けた樹種との関係が強くないところから、特定の菌根菌との強い結びつきはない可能性が高い。したがって、この動物の地理分布は温度や降水量が重要であることがわかる。これまで土壌動物の地理分布については、極域から赤道のような大きな傾度での分布 $ ^{3} $ や、自然誌をもとにした地域ごとの分布解析から緯度の違いに対して種数が異なることが知られている $ ^{17} $ 。一方、日本列島のスケールで環境要因との関係を調べた研究はない。グループⅠは大陸から樺太を経て北海道、本州へと移動してきた種を含んでおり、低温に耐える特徴を持つ。一方、グループⅡ、Ⅲは朝鮮半島経由で九州、本州へと移動してきた種を含んでおり、これらの結果は、今立、中村によるこれまでの研究と一致する。グループⅡは降水量と負の、グループⅢが最高、最低気温および降水量と正の関係を持ち、本州のなかで異なる環境を利用していることがわかる。関東、本州中部でカマアシムシ類の種数が多くなっており、関東、本州中部でカマアシムシ類の種数が多くなっているのは3グループの分布が重なるためであることがわかった。

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

## （1） 科学的意義

モニタリングサイト1000では、コアサイトでピットホールトラップ法により地表徘徊性の無脊椎動物のモニタリングを行っている⁶。各サイトで実際に詳細な食物網関係を記述するには多大な努力が必要である。本研究で、ミミズと地表徘徊性甲虫の連結度を求めたところ、温度との間に弱いながらも負の相関がみられたことにより、地表徘徊性甲虫のモニタリングが、森林の物質循環に重要な役割を果たしているミミズの状態とモニタリングでカバーできていないトリや哺乳類といった他の捕食者の状態を指標するのに適していることを明らかに出来た。地表性甲虫の変化は、餌資源としてのミミズか、競争関係にある他の捕食者の変化を意味する。したがって、ピットホールトラップによるモニタリングは妥当な方法である。
カマアシムシは、長い歴史を経て日本列島に定着した。土壌中に生息し、地表面を積極的に移動しないので、気候変動が生じた場合に、カマアシムシは敏速に分布を変更できない動物群の代表である。カマアシムシ群集のタイプごとの分布が、気温や降水量と強い関係を持つことから、これらの生物とより移動性に富む生物、たとえば植物やカマアシムシの餌となっている糸状菌との間に従来の分布とのずれが乗じる可能性が示唆できる。また、関東が種多様性の中心であることは、過去の日本列島への侵入過程を反映していることから、生物多様性の保全に歴史的な経緯の理解が必要であることを示す新たな事例であることを明らかにできた。

## （２）環境政策への貢献

今後、環境省、農林水産省林野庁等の生物多様性モニタリング関連の委員会等を通して、本研究で得られた知見の政策への反映を提案していくとともに、成果の広報・普及に努める。

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

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

### 7. 研究成果の発表状況

## （1） 誌上発表

<論文（査読あり）>

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

＜查読付論文に準ずる成果発表＞（「持続可能な社会・政策研究分野」の課題のみ記載可。）特に記載すべき事項はない

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

2) 金子信博：時空間情報プラットフォーム-環境情報の可視化と協働、佐土原聡編．東京大学出版会. 92-108 (2010).

「生態系サービス維持のための土壌生態系保全-土壌生態（第8章）」

3） 金子信博：めぐる（生命誌年刊号 vol.61-64）、中村桂子編、東京、新曜社. 162-172 (2010). 「土はやっぱり生きている一土壌動物が育む土壌環境」

4) 金子信博：タクサ28:28-31(2010) 「土壌のbioturbation-土壌生物の多様性と土壌構造の関係」

5) 金子信博：ビオフィリア7:56-58 (2011) 「『ミミズ』と『地球温暖化』—土壌生物と地球環境問題—」.

6) 金子信博：「森のバランス－植物と土壌の相互作用－」丹下健，高橋正通，金子真司，荒木誠，金子信博，佐藤保，松浦陽次郎，森下智陽，三浦覚，戸田浩人，徳地直子（編著）東海大学出版会，pp.21-27(2012）「第3章 生物多様性」
## （2） 口頭発表（学会等）

1) Toyota, A., I. Tayasu, et al. (2009). Comparison of soil carbon accumulation in conifer and broad leaf forests using radiocarbon measurements. International Radiocarbon Conference, Hawaii.

Arai, M., I. Tayasu, et al. (2010). Changes in aggregate carbon dynamics during conversion from tillage to no-tillage system revealed by radio carbon analysis. Isotope. Kyoto.

2) Kawaguchi, T. and N. Kaneko (2010). N:P stoichiometry and life form pattern of Japanese earthworm species. 生態学研究センターセミナースペシャル. 京都大学生態学研究センター.

3) 甘楽法, 金子信博, et al. (2010). 農地におけるミミズの土壌改変効果がダニ・トビムシの多様性に与える影響. 第33回日本土壌動物学会. 同志社大学、京田辺市.

4) 甘楽法，金子信博, et al. (2010). 林地における土壌構造が小型節足動物と物質動態に与える影響-土壌填圧実験から-. 生物地球化学研究会島根大会. 松江.

5) 橋本みのり，中村修美, et al. (2010). 土壤動物の地理分布情報および機能解析データベースの構築にむけて〜カマアシムシ類のデータを用いた分布図作成〜. 第33回日本土壤動物学会. 同志社大学、京田辺市.

6) 中森泰三, 川口達也, et al. (2011). 餌質によりシマミミズの摂食量と微量元素の生物蓄積量が変化する. 第34回日本土壌動物学会, 札幌.

7) 丹羽慈, 南波興之, et al. (2011). 地表徘徊性甲虫類のトップダウン・ボトムアップ制御: 北海道の落葉広葉樹林における操作実験. 第34回日本土壌動物学会.

8) 廿楽法, 金子信博, et al. (2011). ミミズによるニッチ構築がササラダニ群集構造に与える影響. 日本生態学会第58回大会. 札幌.

9) 橋本みのり，高橋正通, et al. (2011). スギ人工林におけるミミズ土壌団粒形成量の推定. 日本生態学会第58回大会. 札幌.

10) 荒井見和 and 金子信博 (2011). 不耕起栽培圃場におけるミミズの土壌改変が土壌の炭素貯留に及ぼす影響. 日本生態学会第58回大会. 札幌.

11) Arai, M., N. Kaneko, et al. (2012). Changes in soil carbon accumulation and soil structure in the no-tillage management after conversion from conventional managements. The 5th East Asian Federation of Ecological Societies (EAFES), Ryukoku University, Otsu, Shiga, Japan.

12) Miura, T., N. Kaneko, et al. (2012). Agricultural management alters soil microbial community, hence litter decomposition: Case study of a sugarcane plantation in Lampung province, south Sumatra, Indonesia The 5th East Asian Federation of Ecological Societies (EAFES), Ryukoku University, Otsu, Shiga, Japan.

## （3） 出願特許

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

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

特に記載すべき事項はない
## （５）マスコミ等への公表・報道等

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

## (6) その他

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

### 8 \. 引用文献

1) 金子信博：地球環境学事典、総合地球環境学研究所編、東京、弘文堂、146-147 (2010). 「土壌動物の多様性と機能」

2) 金子信博：めぐる（生命誌年刊号 vol.61-64）、中村桂子編、東京、新曜社、62-172（2010）.「土はやっぱり生きている一土壌動物が育む土壌環境」

3) 金子信博：「土壤生態学入門」 199pp (2007).

4) Sackett, T.E., Classen, A.T. and Sanders, N.J.: Oikos, 119(12), 1984-1992 (2010).

“Linking soil food web structure to above- and belowground ecosystem processes: a meta-analysis.”

5) Tylianakis, J.M., Raphael K. Didham, Jordi Bascompte and David A. Wardle: Ecology Letters, 11(12), 1351-1363 (2008).

“Global change and species interactions in terrestrial ecosystems.”

6) 環境省自然環境局生物多様性センター : (2009).

「重要生態系監視地域モニタリング推進事業（モニタリングサイト1000）森林・草原調査第1期取りまとめ報告書」

7) Kaneko N, Minamiya Y, Saito M, Hashimoto M, Nishina K, Nakamura O. : Diversity (投稿中). “Biogeography and species assembly of Protura in forest soils in Japan.”

8) Sánchez-Moreno S, Ferris H, Young-Mathews A, Culman SW, Jackson LE.: Soil Biology and Biochemistry 43, 2374-2383 (2011).

“Abundance, diversity and connectivity of soil food web channels along environmental gradients in an agricultural landscape.”

9) 今立源太良：Edaphologia, 38: 17-26 (1988). 「ウェノカマアシムシの諸型」

10) Colwell, R.K., : University of Connecticut, USA (2009).

“EstimateS: Statistical estimation of species richness and shared species from samples, ver.8.2.0.”

11) McCune, B., Mefford, M.J.: MjM Software Design, Gleneden Beach, Oregon, USA. 1999. "PC-ORD. Multicariate Analysis of Ecological Data, Version 4."

12) Hill, M.O.: Cornell University, Ithaca, New York, USA. 1979.

“TWINSPAN- a FORTRAN program for arranging multivariate data in an ordered two-way table by classification of the individuals and attributes. Ecology and Systematics.”

13) R Development Core Team: Vienna, Austria. URL http://www.R-project.org/ (2011) "R: A language and environment for statistical computing. R Foundation for Statistical Computing."

14) Dufrêne, M., Legendre, P.: Ecological Monograph, 67: 345-366 (1997).

“Species assemblages and indicator species: the need for a flexible asymmetrical approach.”
15) Borcard, D., Gillet, F., Legendre, P.: Springer Science+Business Media, New York, USA 2011. “Numerical Ecology with R.”

16) 気象庁：http://data.sokki.jmbsc.or.jp/cdrom/mesh_climatic_data/index.htm 「メッシュ気候値2000」

17) Uchida T, Kaneko N, Ito MT, Futagami K, Sasaki K, Sugimoto A: Applied Soil Ecology, 27, 153-163 (2004).

“Analysis of the feeding ecology of earthworms (Megascolecidae) in Japanese forests using gut content fractionation and d15N and d13C stable isotope natural abundances.”

18) Lilleskov, E.A., Mattson, W.J., Storer, A.J.: Diversity and Distribution 14: 893-904 (2008). "Divergent biogeography of native and introduced soil macroinvertebrates in North America north of Mexico."
# Large-scale Estimations of Ecosystem Functions Using Biological Indexes and Development of Eco-informatics

Principal Investigator: HIURA Tsutom

Institution:

Field Science Center for Northern Biosphere

Hokkaido University

Takaoka, Tomakomai 053-0035 Japan

Tel:+81-144-33-2171 Fax:+81-144-33-2173

E-mail: hiura@fsc.hokudai.ac.jp

Cooperated by: National Institute for Environmental Studies, Gifu University, Tokyo University, Yokohama National University

## [Abstract]

Key Words: Global warming, Litter fall, Canopy photosynthesis, Remote sensing, Data base

We revealed responses of biological traits, community structure, and ecosystem function in the forest ecosystem to environmental change focusing on 1) Evaluation of biological indexes and its response to environmental change, 2) Development of monitoring technique for canopy photosynthetic activity, 3) Development of evaluation protocols of ecosystem functions using stable isotope ratios, 4) Development of convenient remote sensing sensors, 5) Construction of the data base of soil invertebrates distribution and leaf traits. Changes in photosynthetic resource use efficiency, secondary metabolites, and litter fall of oak trees would suggest considerable ecological consequences such as herbivory, decomposition, soil fauna, and nutrient cycling in a warmer climate. We also developed some new techniques to monitor the change of ecological features and functions of forest ecosystems.
## 生產量

## 葉量

機能量

生態系機能の広域評価

葉量推定法の開発

植生指数と

機能量の解析

簡易カメラの開発

環境省モニ1000などのネットワーク・データベースの有効活用

## 指標群

同位体による窒素起源の推定

大規模実験による

温暖化応答の検出

気候変動の影響を強く受ける

土壌動物では初のパターン発見
## 平成29年5月30日

政府統計

【照会先】

雇用均等・児童家庭局職業家庭両立課

課長源河真規子

報道関係者 各位

課長補佐 土岐 祥蔵

(代表電話) 03(5253)1111(内線 7855)

(直通電話) 03(3595)3274

## 「平成28年度雇用均等基本調査（速報版）」を公表します

～育児休業取得者の割合は、女性は81.8%、男性は3.16%でともに微増～

厚生労働省は、このたび、「平成28年度雇用均等基本調査（速報版）」の結果を取りまとめましたので、公表します。

「雇用均等基本調査」は、男女の雇用均等問題に関わる雇用管理の実態把握を目的に、毎年実施しています。

今回の速報版では、育児休業取得者割合に関する調査項目について取りまとめています。なお、その他の項目を加えた確報版は、7月末ごろに発表する予定です。

## 【育児休業取得者の割合】

女性 : 81.8% （対前年度比 0.3 ポイント上昇）

男性 : 3.16% （対前年度比 0.51 ポイント上昇）

## [調査の概要]

## ■調査の時期

平成28年10月1日現在の状況について、平成28年10月1日から10月31日までの間に実施。

## ☑調査対象数

6,092 事業所（有効回答数：4,213 事業所、有効回答率：69.2%）

## 調查方法

厚生労働省雇用均等・児童家庭局から調査対象事業所に対して、

郵送により調査票を配布・回収する方法

[别添]

· 事業所調查 結果概要

· 事業所調查 付属統計表

・調査の概要
## 第7節 障害者雇用の動向

## ●障害者の雇用者数は13年連続で過去最高

第1－（2）-36図で、障害者の雇用状況についてみていこう。2016年6月1日現在の障害者の雇用者数は47万人と、前年比で4.7%の増加となり、13年連続で過去最高となったほか、実雇用率についても、1.92%と前年比0.04%ポイントの上昇となり、5年連続で過去最高となった。前年比で増加率をみると身体障害者は2.1%、知的障害者は7.2%、精神障害者は21.3%となっており、精神障害者の伸び率が特に大きくなった。

法定雇用率達成企業の割合は、法定雇用率達成企業の割合は、50～100人未満が45.7%（前年は44.7%）、100～300人未満が52.2%（同50.2%）、300～500人未満が44.8%（同44.0%）、500～1,000人未満が48.1%（同44.6%）、1,000人以上が58.9%（同55.0%）となり、全ての規模の区分で前年より増加した。

## 第1-(2)-36図 雇用されている障害者の数と実雇用率の推移

○ 2016 年の雇用障害者数は、13 年連続で過去最高を更新した。

資料出所 厚生労働省「障害者雇用状況報告」をもとに厚生労働省労働政策担当参事官室にて作成

注）1）雇用義務のある企業（平成24年までは56人以上規模、平成25年は50人以上規模の企業）についての集計である。

2）「障害者の数」とは、次に掲げる者の合計数である。

～2005年 身体障害者（重度身体障害者はダブルカウント）

知的障害者（重度知的障害者はダブルカウント）

重度身体障害者である短時間労働者

重度知的障害者である短時間労働者

2006年～ 身体障害者（重度身体障害者はダブルカウント）

知的障害者（重度知的障害者はダブルカウント）

重度身体障害者である短時間労働者

重度知的障害者である短時間労働者

精神障害者

精神障害者である短時間労働者（精神障害者である短時間労働者は0.5人でカウント）

2011年～ 身体障害者（重度身体障害者はダブルカウント）

知的障害者（重度知的障害者はダブルカウント）

重度身体障害者である短時間労働者

重度知的障害者である短時間労働者

精神障害者

身体障害者である短時間労働者（身体障害者である短時間労働者は0.5人でカウント）

知的障害者である短時間労働者（知的障害者である短時間労働者は0.5人でカウント）

精神障害者である短時間労働者（精神障害者である短時間労働者は0.5人でカウント）

3）法定雇用率は、2012年までは1.8%、2013年4月以降は2.0%となっている。

4）2010年5月に制定の働法（短時間労働者の箇子、除くの引き下げ箇）が、本法は、2011年以降上昇を想定する。

4）2010年7月に制度改正（短時間労働者の算入、除外率の引き下げ等）があったため、2011年以降と2010年までの数値を単純に比較することは適当ではない状況である。
## 物価・消費の動向

経済の好循環を継続させていくためには、企業収益の拡大を賃金上昇や雇用拡大につなげていき、消費の拡大に結び付けることが重要である中で、消費の動向をみると、消費総合指数は緩やかながら上昇傾向で推移しており、雇用・所得環境が改善する中で持ち直しがみられる。そこで、本章では、このような2016年度の物価と消費の動向についてみていく。

## 第1節 消費者物価の動向

## ●生鮮食品及びエネルギーを除く総合は緩やかに上昇

物価の状況について、まず最終需要財の価格の動きを表す消費者物価指数の推移をみていく。第1－（4）－1図により、「総合」「生鮮食品を除く総合（いわゆるコア）」「生鮮食品及びエネルギーを除く総合」に分けて消費者物価指数の推移をみると、「総合」については、2014年4月の消費税率引上げに伴う上昇分が剥落した2015年4月以降は、前年同月比0％台で推移していた。その後、円高傾向に伴う輸入財の価格の押し下げやエネルギー価格の低下等による物価の低下が進み、2016年4月には前年同月比がマイナスに転じたが、2016年10月以降は天候不順の影響による生鮮食品の値上がりや原油等の値上がりを受け、プラスとなっている。

一方、「生鮮食品及びエネルギーを除く総合」は、2016年に入ってからも引き続き緩やかな上昇を続けていたが、2016年後半は前年同月比0％近傍と横ばいの動きとなっている。

そこで、2016年平均の消費者物価指数を10大費目別の前年比でみると、生鮮食品や外食の価格上昇を受けて「食料」「被服及び履物」などで上昇した一方、「交通・通信」「光熱・水道」などは低下した（付1-(4)-1表）。
## ●消費者マインドは持ち直しの動きがみられる

消費の動向については、所得以外の要因である消費者のマインドの影響も受ける。そこで、第1－（4）－4図により、消費者態度指数の推移をみていこう。消費者態度指数は、2014年4月の消費税率引上げ時を底としていったん回復の動きをみせた後、2016年に入ってからは、消費者マインドには足踏みがみられたが、その後は持ち直しの動きがみられた。

また、項目別に2016年の推移をみると、「雇用環境」は比較的高水準で推移している中、「暮らし向き」も緩やかな上昇傾向で推移した。

## 第1-(4)-4図 消費者態度指数の推移

○ 雇用環境は高水準を続けている中で、暮らし向きや収入の増え方は緩やかな上昇傾向で推移。

資料出所 内閣府「消費動向調査」をもとに厚生労働省労働政策担当参事官室にて作成

（注）1）二人以上の世帯、季節調整値。

2）消費者態度指数の作成方法は以下のとおり。

①「暮らし向き」「収入の増え方」「雇用環境」「耐久消費財の買い時判断」の4項目について、今後半年間の見通しについて5段階評価で回答を求める。

②5段階評価のそれぞれ「良くなる・大きくなる」に（+1）、「やや良くなる・やや大きくなる」に（+0.75）、「変わらない」に（+0.5）、「やや悪くなる・やや小さくなる」に（+0.25）、「悪くなる・小さくなる」に（+0）の点数を与え、これに各回答区分のそれぞれの構成比（%）を乗じ、乗じた結果を合計して、項目ごとに消費者意識指標（原数値）を算出する。

③4項目の消費者意識指標（原数値）それぞれについて季節調整値を算出し、それらを単純平均して消費者態度指数（季節調整値）を算出する。

3）2013年4月調査から、訪問留置調査法から郵送調査法に変更したことにより、不連続が生じている。

## ●消費総合指数は緩やかに上昇

次に、消費者以外の面にも着目して消費の動向についてみていくために、第1-(4)-5図により、需要側と供給側の要因を包括的にまとめた消費総合指数と、販売側の状況を示す実質小売業販売額指数についてみていこう。まず、実質小売業販売額指数をみると、2014年3月までは消費総合指数と同様の動きを示していたが、2014年4月の落ち込みが大きく、同年5月に回復に転じた後は横ばいで推移しており、2016年には前年比0.6%減少と2013年の水準まで回復していない状況にある。

そこで、消費総合指数をみると、2014年4月には消費税率引上げによる駆け込み需要の反動でいったん低下したが、その後は緩やかながら上昇傾向で推移していることが分かる。
# 平成29年3月27日栃木県那須町の雪崩について

※これは速報であり、数値等は今後も変わることがある。

平成29年3月29日

8時30分現在

内閣府

1 事故発生場所（消防庁情報）

栃木県那須郡那須町大字湯本字那須岳国有林内 那須温泉ファミリースキー場付近

2 事故発生日時（消防庁情報）

平成29年3月27日 8時30分（推定）

3 事故の発生状況（栃木県情報）

3月27日9時20分、那須温泉ファミリースキー場で雪崩が発生したとの110番通報（登山研修中の高校生等が雪崩に巻き込まれたもの）

## 4 気象状況（気象庁情報）

・気象概況

27日（月）は、低気圧が関東の南海上を発達しながら北東へ進み、さらに上空に寒気が流れ込んだため、栃木県の山地では大雪となった。

栃木県那須高原では、27日（月）1時から10時までの降雪量が34センチに達し、短い時間にまとまった降雪となった。

・3月27日8時30分現在での那須町への注意報発表状況

大雪、なだれ、着雪注意報（3月26日10時32分）

## 5 被害状况

(1) 人的被害（消防厅情報：3月29日08:30現在）

死者 8名（男性教員1名、男子生徒7名）

重傷 2名

軽傷 38名

※高校関係の入山者は48名（生徒40名、教員8名）、他に入山しなかった高校関係者14名（うち生徒11名）

※行方不明者なし

## 6 政府の主な対応

(1) 官邸の対応

3 月 27 日 11:19 官房長官会見（～11:36）

11:50 情報連絡室設置

(2) 関係省庁災害対策会議等の実施

3月29日10:00 平成29年3月雪崩等に係る関係省庁災害対策会議（予定）

## 7 各省庁等の対応

(1) 内閣 府の対応

3 月 27 日 10:15 内閣府情報連絡室設置

(2) 警察 庁の対応

・警察署員・県警機動隊が現場臨場

・県警へりが上空から状況確認

(3) 消防 庁の対応

3月27日10時00分 応急対策室長を長とする消防庁災害対策室設置(第1次応急体制)

12時25分 国民保護・防災部長を長とする消防庁災害対策本部に改組（第2次応
## 急体制）

１２時２５分 消防庁職員２名を現地に派遣

## 【消防機関等の活動】

・雪中の要救助者の救助活動を実施。自力下山が出来ない入山者の下山支援及び負傷者の救護活動等を実施（3月27日）

・無人へリによる上空からの救助活動現場の確認を実施（3月28日）

ア 地元消防機関の活動

3月27日 9時22分 那須地区消防組合消防本部覚知

→18隊52名が出動

17時10分 栃木県消防防災航空隊が出動（負傷者の医療機関への搬送）

3月28日　9時00分　無人へリによる上空からの救助活動現場の確認を実施（2隊6名）

県内消防応援等の活動

3月27日10時44分 県内消防応援を要請

→宇都宮市消防局、塩谷広域行政組合消防本部、南那須地区広域行政事務組合消防本部が出動（合計13隊54名）

１１時２０分 県内消防応援を追加要請

→日光市消防本部が出動（1隊3名）

１１時３６分 消防相互応援協定に基づく出動を要請

→福島県 白河地方広域市町村圏消防本部が出動（3隊8名）

## ウ 緊急消防援助隊の活動

3月27日12時15分 消防庁から埼玉県に対し、情報収集活動を任務として、さいたま市消防局に配備している消防活動用偵察システム（無人ヘリ）の出動準備を依頼

12時25分 栃木県知事から消防庁長官に対し、緊急消防援助隊（上空からの情報収集活動）の出動要請

12時25分 消防庁長官から埼玉県知事に対し緊急消防援助隊（いたま市消防局）の出動の求め

12時37分 緊急消防援助隊（さいたま市消防局）出動（3隊10名）

16時15分 緊急消防援助隊（さいたま市消防局）現場到着

救助及び救護活動を実施

3月28日　9時00分　無人へリによる上空からの救助活動現場の確認を実施

→12時00分　緊急消防援助隊（さいたま市消防局）引揚決定

## (4) 防衛省の対応

ア 災害派遣の概要

（1）要請日時平成29年3月27日（月）10時37分

（2）要請元栃木県知事

（3）要請先陸上自衛隊第12特科隊長（宇都宮）

（４）要請の概要　雪崩遭難者の人命救助

（５）発 生 場 所 栃木県那須郡那須町那須温泉ファミリースキー場

（6）撤收要請日時 平成29年3月27日（月）17時23分

イ 活動部隊 陸 自 第12特科隊（宇都宮）、第12ヘリコプター隊（北宇都宮）、第12偵察隊（相馬原）、第12化学防護隊（相馬原）、第48普通科連隊（相馬原）、第12高射特科中隊（相馬原）、東部方面後方支援隊（相馬原）、第2普通科連隊（高田）、第30普通科連隊（新発田）、第12後方支援隊（宇都宮）、東部方面ヘリコプター隊（立川）、自衛隊栃木地方協力本部

ウ 活動規模 人 員 約230名