Neozephyrus japonicus

Ecosystem Conservation Society - Japan

Ecosystem Conservation Society - Japan

CONTENTS

1． ハビタット利用に関する既存文献情報
1 概要

12 食物

13 水

14 繁殖

15 カバー

16 行動圏

1

2． ハビタット適性指数モデルの構築
2

1 方法(a) ハビタット変数候補の定義

2(b) ハビタット変数の選択

3(c) HSI結合式の選択

32 結果

4(a) ハビタット変数の選択

4(b) ハンノキ林縁部5m帯面積（ｍ ）のSI関数

5(c) 周辺150m広葉樹林林縁部面積割合（％）のSI関数

5(d) HSIの統合式

62

3 考察および今後の課題

7

3． ハビタット適性指数モデル
1 モデルの適用範囲

8(a) 地理的範囲

8(b) 季節8(c) 最小ハビタット面積

82 モデルの構造

83 適性指数グラフ

94 HSI算出のための統合式

94． 謝辞105． 引用文献

10Neozephyrus japonicus
ミドリシジミのH S I モデル ver.1.0

ハビタット利用に関する既存文献情報

3水
成虫による吸水は確認されている（筆者

1 概要
ミドリシジミ（ Neozephyrus japonicus ）は中国東北部、西シベリア、アムール、ウスリー、サハリン、日本に分布する。国内では北海道・本州・四国・九未発表）

が、生息の上での水の必要性について述べられた文献はない。なお、幼虫の餌であるハンノキ類は湿潤な環境に生育する。

4 繁殖

州および周辺離島（利尻島・佐渡・淡路島・小豆島）

ハンノキ、ヤマハンノキ、ミヤマハンノキの樹幹、に分布し、平地や低山地の湿地のハンノキ群落、あ枝、特に休眠芽付近に多く産卵する（福田ほか

るいは山地のヤマハンノキやミヤマハンノキの自生

1 9 8 4 ）。 一 般 に 卵 は 林 の 周 辺 部 の ハ ン ノ キ に 多 く 、する林に生息する（福田ほか 1984）。年1回の発生特に少し離れた独立樹に集中することが知られていで6月中旬から羽化を始め、7月に最盛期に達する。

る（竹内 1992）。蛹化場所については、食樹の根元寒冷地では7〜8月に出現する。産卵活動のピークは付近の落ち葉の裏で蛹化することが知られている(市8月で、卵で越冬し、関東地方平地では4月上旬に孵川 1978；福田ほか 1984)。繁殖行動に関してオス化し、5月末から6月上旬ごろ蛹化する（福田ほか

成虫は強い占有性を示す。オスは夕方、午後4時〜61984）。本種は、東京都、埼玉県、千葉県、愛知県、時ごろにかけて活発に樹上を飛び、他個体を追飛す

大阪府、兵庫県の都道府県版レッドデータブックにる（福田ほか 1984）。朝も弱い活動性を示すが、昼掲載されている（東京都 1998；大阪府 2000；千間 は 不 活 発 で 、 食 樹 ※1 や 付 近 の 下 草 に 静 止 し て い る

葉県 2000；愛知県 2002；埼玉県 2002；兵庫県

ことが多い。メスも夕方を中心に活動するが、オス2003）。また、埼玉県では県の蝶に指定されているほど活発ではなく、昼間は食樹付近の葉上に静止し（江村 1992）。ていることが多い（福田ほか 1984）。

2 食物
幼虫は4月から5月にかけてハンノキ、ヤマハンノ

5 カバー
本種のカバーの必要性について述べられた文献はない。

キ、ミヤマハンノキの葉を寄主植物として利用する（福田ほか 1984；猪又 1990）。成虫は活動期に開花する花の花蜜（クリ・ノリウツギ・チダケサシ・

6 行動圏
行動圏について述べられた文献はない。

シシウド・ノダケ・ヤマゼリ・ソバ・オカトラノ
オ・ヒトツバハギ・シラキなど）、樹液（クヌギ・ナラガシワ・ヤナギ類など）、鳥糞を餌として利用している（福田ほか 1984）。

【本モデルの引用例】（財）日本生態系協会ハビタット評価グループ (2004) ミドリシジミのHSIモデル ver.1.0．（財）日本生態系協会編．ハビタット評価モデルシリー
ズ３．（財）日本生態系協会, 東京．
※１食樹：幼虫の食物となる植物。樹木である場合、食樹と呼ぶ。ここでは、ハンノキ、ヤマハンノキ、ミヤマハンノキ。

1Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

2

ハビタット適性指数モデルの構築(a) ハビタット変数候補の定義
前節の文献調査の結果より、ミドリシジミのハビタット適性は繁殖条件と食物条件によって決定され、

1 方法

水は制限要因となっていないと考えられた。カバー

本モデルは、2000年から2004年に埼玉県さいたま市、川島町、鴻巣市の荒川河川敷で実施された、条件については、その重要性が不明であったため、本モデルでは検討対象としなかった。

樹林とミドリシジミの分布調査データの一部を用い繁殖条件については、ハンノキ林分の面積、あるて 構 築 し た （ 埼 玉 県 生 態 系 保 護 協 会

いはハンノキ林分の林縁部が増えるほど、本種にと2001,2002,2003,2004）。調査方法の概要を以下にって適性が高い条件となることが予想された。
食物条件については、成虫の食物資源（花蜜・樹示す。
樹林の分布調査として、まず、既存データや空中液）を供給する広葉樹林の面積、あるいは吸蜜植物

写真、現地踏査などによって、調査地域内の全てのとなる草本類が繁茂しやすい林縁部が増えるほど、ハンノキ林分（87林分）を地図化した（最終的に本種にとって適性が高い条件となることが予想されはGISデータ化した）。次に、ミドリシジミの生息確た（食物資源を供給する広葉樹林は、産卵場所であ認調査を、これら87箇所のハンノキ林分において実るハンノキ林分そのものと、そのハンノキ林分から施した。2000年の6月下旬には、飛翔個体の観察お一定範囲内に存在する広葉樹林とした）。

※２

よび叩きだし法によって成虫の確認を行った。2001

なお、繁殖や食物条件の適性に影響する可能性の年の5月中旬には踏査により幼虫確認調査を、同年6ある「林縁部の幅」の具体的な値については、文献

月中下旬には前年同様に成虫確認調査を行った。

調査では特定できなかったため、ここでは、2mとハンノキ林以外の樹林も含めた全樹林の分布状況5mの幅を変数候補として設定することとした。同様

は、既存データや空中写真、IKONOS画像、現地踏査に、「ハンノキ林分から一定範囲内に存在する広葉樹などにより、2000年から2004年にかけて調査し、林」といった場合の「範囲」についても、距離が明GISデータ化した。なお、河川敷以外の区域（提内地）

確に示された文献はなかったため、20、40、50、についても、必要に応じて樹林の分布調査を行った。

60、80、100、150、200mという範囲を変数候補

パフォーマンスメジャー ※３ は、ハンノキ林分におけるミドリシジミの生息の有無とし、成虫あるいは幼虫の生息が確認できたハンノキ林分に対しては1として設定した。
以上の検討を踏まえ、最終的に、38変数をハビタット変数候補として設定した（表1）。

を、確認できなかったハンノキ林分に対しては0の値を割り当てた。
HSIモデルに用いるハビタット変数は、以下の手順で選定した※2基本的に、道路等で区切られている場合は別の林分としたが、隣接し樹冠部が接している場合は同一林分として扱った。
※３パフォーマンスメジャー：HSIと対応する、個体群の具体的な指標値。

2

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

表１.ハビタット変数候補一覧(b) ハビタット変数の選択
分類木（classification tree）を用いて、38のハビ不純度の最小変化量=0.0001
カテゴリ目的変数の不純度=Gine

タット変数候補の中から、ミドリシジミの生息の有無に影響する重要な変数を選択した。解析ソフトは、
Answer Tree 2.0J（SPSS社）を使用し、以下の条件で樹木を成長させた後、剪定を行った。

(c) HSI統合式の選択
HSIは、単一のSI関数あるいは複数のSI関数の統合式によって求められる。SI関数の統合に際しては、いくつかのタイプの統合式（例えば、相加平均、相3

成長方法：CART乗平均など）でHSIを計算し、実際の各ハンノキ林分

樹木の最大深さ：10におけるミドリシジミの確認の有無と比較すること

ケースの最小数：親ノード=10、子ノード=5により最も適切な統合式を選定した。

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

辺1 5 0 m広 葉樹林林縁部面積割合（％）が1 0 .3％以

2 結果

上の林分ではミドリシジミが生息せず、10.3％未満の林分でミドリシジミが生息すると予想された。

(a) ハビタット変数の選択
樹形モデルを成長させた結果、第1段階の分岐では、ハ ン ノ キ 林 縁 部 5 m 帯面積（ｍ ２ ）が1100.8ｍ ２ 以 上以上の解析結果より、ハンノキ林縁部5m帯面積の林分でミドリシジミが生息すると予想された（図（ｍ）（繁殖・食物条件）を第1段階の説明変数とし、1）。1100.8ｍ 未満の林分については、周辺150m広その一部に対し周辺150m広葉樹林林縁部面積割合

葉樹林林縁部面積割合（％）によって、さらに分岐（ ％ ）（ 食 物 条 件 ） を 第 2 段 階 の 説 明 変 数 、 ミ ド リ シ

した。この第2段階の分岐では、周辺150m広葉樹林ジミの確認の有無を目的変数として、モデルの構築林縁部面積割合（％）が6.8％以下の林分では、ミドを行うこととした。

２

リシジミが生息しないと予想された。6.8％より大きい林分については、同じく周辺150m広葉樹林林縁部面積割合（％）によって第3段階の分岐が発生し、周図１

樹形モデルの結果

4

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

9.08×10 -4 ×V1とした。これらをハンノキ林縁部5m(b) ハンノキ林縁部5ｍ帯面積（ｍ２）のSI関数樹形モデルの第1段階目では、ミドリシジミの確認帯面積のSI関数とした（図3）。の有無は、ハンノキ林縁部5m帯面積（ｍ ）が、２1100.8ｍ ２ の値で分岐した。ハンノキ林縁部5m帯面積が1100.8ｍ ２ 以上の15林分では、全ての林分でミドリシジミが確認できた。一方、1100.8ｍ ２ 未満の72林分ではミドリシジミが確認できなかった林分数は47であり、残りの25林分は確認された林分となった（図2）。

図3．ハンノキ林縁部5m帯面積（V1）と適性度（SIV1）の関係(c) 周辺150ｍ広葉樹林林縁部面積割合（％）のSI関数
樹形モデルの第2段階目では、ミドリシジミの確認の有無は、周辺150m広葉樹林林縁部面積割合（％）
が6.8％の値で分岐した。6.8％より大きい36林分の
図2．ハンノキ林縁部5m帯面積（ｍ ）とミドリシジミの確内、ミドリシジミの生息が確認できた林分数は20、認の有無との関係確認できなかった林分数は16であった。一方、

２6.8％以下の36林分の内、ミドリシジミの生息が確
図2を元に、以下のようにSI関数を定義した。

認できた林分数は5、確認できなかった林分数は31

ハンノキ林縁部5m帯面積（以下、V1とする）がとなった（図4）。さらに、6.8％より大きい林分のう

1100.8ｍ ２ 以上の林分については、V1に対応する適ち、10.3％以上に該当する21林分では、ミドリシジ性度（以下、SIV1とする）を1とした。V1が0の場

ミの生息が確認できた林分数は6、確認できなかった合は、産卵可能なハンノキが全くない状態となるた

林分数は15であった（図4）。

め、SIV1も0とした。0〜1100.8の区間については、ミドリシジミが確認されている林分と、されていない林分が混在することから、SIV1は0から1へと徐々に増加するものとした。増加方法は不明であったた
め、ここでは直線的な増加、すなわち、SIV1＝5

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

図4．周辺150m広葉樹林林縁部面積割合とミドリシジミの確認の有無との関係

図4を元に、以下のようにSI関数を定義した。

図5．周辺150m広葉樹林林縁部面積割合（V2）と適性度（SIV2）の関係(d) HSIの統合式

周 辺 1 5 0 m 広 葉 樹 林 林 縁 部 面 積 割 合 （ ％ ）（ 以 下 、繁殖条件と食物条件はどちらも本種の生息には欠V2 とする）が6.8％より大きい林分についてはかせないものであるが、本モデルで選択されたハビ

SIV2＝1とした。V2が0％の場合は食物資源が全くタット変数（V1,V2）を見ると、前者（V1）は両条ない状態となるため、SIV2＝0とした。0〜6.8の区件を兼ねていることから、SIV1とSIV2は、必ずしも間は、SIV2は0から1へと直線的に増加するものとし制限的に作用しなくても良いと考えられる。従って、た（SIV2＝1.471×10 -1 ×V2）。また、樹形モデル

ここでは乗法：SIV1×SIV2、幾何平均；(SIV1×の第3段階の分岐では、10.3％以上の林分はミドリ

SIV2) 1/2 、三乗根：(SIV1×SIV2) １/3、最小：Min(SIV1，シジミの生息の可能性が低いと予想された。しかし、SIV2)、算術平均：（SIV1＋SIV2）/2の5式をHSIのこの分岐を、ミドリシジミの生態から合理的に説明統合式候補として、検討を行った。

することは困難であり、今回、データが得られなかった他のハビタット変数（林齢など）が影響してい各統合式候補により算出された値が0.5未満の場合る可能性も考えられた。このため、本モデルでは、を0、0.5以上の場合を1で示したカテゴリ（「予測」）10.3％以上の林分についてもSIV2＝1とした。これと、実際の生息の有無（生息有り＝1、生息無し＝0）

らを周辺150m広葉樹林林縁部面積割合のSI関数としを 示 し た カ テ ゴ リ （「 実 際 」） と の 間 で 、 誤 分 類 率た（図5）。

（％）を比較したところ、算術平均が最も低くなった（表2）。
また、HSIは、パフォーマンスメジャーがハビタッ
トにより規定される上限を示すため、「予測」＝1で「 実 際 」 ＝ 0 は 許 容 さ れ る 。 し か し 、「 予 測 」 ＝ 0 で

6

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0「実際」＝1となる割合が多い場合は、不適切なモデルといえる。統合式候補を比較すると、「予測」＝0で「実際」＝1となる割合についても、全体の誤分類率と同様に、算術平均が最も低くなった（表2）3 考察および今後の課題
ハンノキ林分におけるミドリシジミの確認の有無とハンノキ林縁部5m帯面積、周辺150m広葉樹林林縁部面積割合を分析することにより、HSIモデルを構

表2.各統合式候補における誤分類率（％）

築することができた。第１に良好な産卵場所、食物資源が存在すると考えられるハンノキ林の林縁部が広ければ、本種が生息する可能性が高くなると考えられた。そして、ハンノキ林縁部が小さくても、成虫の食物資源を供給する場所と考えられる周辺樹林地、特に林縁部が十分に存在することで、本種の生息の可能性が高くなることが考えられた。
ただし、周辺150m広葉樹林林縁部面積割合に関して、第3段階では10.3%未満の林分の方が、10.3％以上の林分よりも、本種の生息が確認された林分の割合が多かった。今回、周辺150m広葉樹林林縁部面積割合が10.3％以上あり、かつ、本種が確認されなかったハンノキ林は、ハンノキ以外の木本種の構成割合や林齢が高いという特徴があった。従って、10.3％以上の林分については、今回、検討できなか以上の検討から、結合式候補の中から算術平均をHSIった他のハビタット変数（すなわち、ハンノキ以外の結合式とした。の木本種の構成割合や林齢など）が影響していた可すなわち、能性もあるだろう。
また、周辺樹林地については、今回の調査地域でHSI＝（SIV1＋SIV2）/2

は針葉樹林が存在しなかったため、広葉樹林のみで解析を行った。食物資源の供給場所として考える場

ただし、SIV1＝1.0のときは、SIV2の値に関わらずHSI＝1とする。

合、針葉樹林は広葉樹林に比べ、特に木本からの供給は少ない可能性がある。従って、今後は、食物供給における針葉樹林の役割など、周辺樹林地についても、より詳細な評価を行う必要がある。

7

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

3

ハビタット適性指数モデル

2 モデルの構造
ミドリシジミのHSIは、繁殖条件と食物条件によっ1 モデルの適用範囲て規定されるものとした。なお、水による制限はな(a) 地理的範囲あるため考慮しないものとした。本モデルにおいていと仮定した。また、カバー条件に関しては不明で本モデルは、埼玉県平野部の荒川本川（埼玉県生は、本種の繁殖条件を一般に産卵場所に適している

態系保護協会 2001；2002；2003；2004）におけといわれているハンノキ林の林縁部面積から評価し、る調査データに基づいて作成した。そのため、本モ

食物条件を蜜源と考えられるハンノキ林の林縁部面デルの適用範囲としては、主に関東地方の平野部が積あるいはハンノキ林の周辺部に存在する広葉樹林推奨される。の林縁部面積の割合により評価する。

(b) 季節周辺部とはハンノキ林から150ｍ以内とした。そのハンノキ林の林縁部は内側に5ｍ幅とした。また、本モデルは、ミドリシジミの通年のハビタット適性を評価するものとする。

範囲内に存在する広葉樹林（以下、周辺広葉樹林とする）の林縁部の内側に5ｍのバッファーを設定し、変数はその合計値のハンノキ林の外側に設定した150mバッファーの面積に対する割合とした。

(c) 最小ハビタット面積
ミドリシジミの最小ハビタット面積に関する情報が記載された文献はない。このため、最小ハビタッ

ハビタット変数、生存必須条件、およびミドリシジミのHSIの関係については図7に示した。

ト面積は設定しないものとする。

＜ハビタット変数 ＞
ハンノキ林縁部5ｍ面積（V1）
周辺150ｍ広葉樹林林縁部面積割合（V2）＜ 生存必須条件＞
繁殖・食物
食物HSI

図7．ミドリシジミのモデルにおける、ハビタット変数、生存必須条件、HSIの関係

8

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

3 適性指数グラフ
V1：ハンノキ林縁部5m帯面積（ｍ ２）

1.0

Ｖ1≧1100.8の場合

9.08×10 -4×V1

V1＜1100.8の場合

SIV1＝

V2：周辺150m広葉樹林林縁部面積割合（％）

1.0

V2＞6.8の場合

1.47×10 -1 ×V2

V2≦6.8の場合

SIV2＝

４ HSI算出のための統合式
1.0

SIV1＝1.0の場合

HSI＝（SIV1＋SIV2）/29

SIV1＜1.0の場合

Neozephyrus japonicus
ミドリシジミのH S I モデル v e r . 1.0

4物編．

謝辞

埼玉県生態系保護協会（2001）平成12年度 荒川生
態系モニタリング調査報告書．国土交通省荒川上

本モデルの作成に当たり、国土交通省荒川上流河流工事事務所委託．

川事務所には調査データの使用を快く許可していた

埼玉県生態系保護協会（2002）平成13年度 荒川生

だいた。山梨県環境科学研究所動物生態学研究室の態系モニタリング調査報告書．国土交通省荒川上北原正彦主任研究員(当協会専門研究委員)には、ご流工事事務所委託．

多忙のところ、モデルの草稿をご確認いただき有益

埼玉県生態系保護協会（2003）平成14年度 荒川生

なご助言を賜った。ここに記して厚く御礼申し上げ

態系モニタリング調査報告書．国土交通省荒川上る。

流河川事務所委託．
埼玉県生態系保護協会（2004）平成15年度 荒川生

5態系モニタリング調査報告書．国土交通省荒川上引用文献

流河川事務所委託．
白水隆（1961）ミドリシジミ類における幼虫食性の進化，蝶と蛾，12．

愛知県（2002）愛知県の絶滅のおそれのある野生生物 レッドデータブックあいち−動物編−．
千葉県（2000）千葉県の保護上重要な野生生物−千竹内崇夫（1992）ミドリシジミの産卵位置について．
寄せ蛾記，63：1317-1322．
東京都（1998）東京都の保護上重要な野生生物種．

葉県レッドデータブック−動物編
江村薫（1992）県蝶（埼玉県）決定のいきさつ．昆虫と自然，27（2）：2-7．
福田晴夫・浜栄一・葛谷健・高橋昭・高橋真弓・田中蕃・田中洋・若林守男・渡辺康之（1984）原色
日本蝶類生態図鑑III．保育社，大阪．
兵庫県（2003）改訂・兵庫の貴重な自然 兵庫県版
レッドデータブック2003．
市川和夫（1978）ミドリシジミの蛹化場所はハンノキ林の樹下落葉中．昆虫と自然，13（13）：2930．
猪又敏男（1990）原色蝶類検索図鑑．北隆館，東京．
大阪府（2000）大阪府における保護上重要な野生生物−大阪府レッドデータブック−．
埼玉県（2002）埼玉県レッドデータブック2002 動

10

ハビタット評価モデル等利用規約
１．本モデルの著作権は(公財)日本生態系協会に帰属し、著作権法によって保護されています。当協会の許可なく本モデルをウェブサイトや印刷媒体に転載することはできません。
２．非営利の学術研究または教育を目的として利用する場合は、出典を明記した上でご利用ください。
営利目的などその他の目的で利用する場合は、事前に当協会の許可が必要となります
ので、利用申請書に必要事項を記載の上、当協会まで郵送してください。
３．利用者が本モデルの利用や利用不能により被った直接的または間接的損害に対し、(公財)日本生態系協会は一切の責任を負いません。

お問合せ・送付先(公財)日本生態系協会 生態系研究センター
ハビタット評価グループ〒330-0802 埼玉県さいたま市大宮区宮町 1-103-1
TEL 048-649-3860
FAX 048-649-3859

YK ビル 6F

ハビタット評価モデル等利用申請書
平成年月日(公財)日本生態系協会会長 殿
団体名
申請者

代表者

印

担当者
住 所
TEL
E-mail
利用規約および利用条件に同意の上、下記のとおり利用を申請します。

記

利用を希望するモデル

利用目的
および事業名・発注者名

利用期間

利用条件
１．上記の目的以外に利用しないこと。
２．利用結果を公表した場合は、速やかに当協会へそのコピーを提出するか、公表資料の入手方法を報告すること。

(財)日本生態系協会は、国内や海外の情報を広く集め、自然と共存する豊かな国づくり、まちづくりを進める専門集団です。自然科学、社会科学の両側面から調査研究を行い、持続可能な発展を目指す国土計画や地域再生プロジェクトに、市民とともに取り組んでいます。また、各種法制度に、自然との共存という視点を加えるための提言を行うとともに、環境教育活動、国際シンポジウムやセミナーの開催、ビオトープ管理士制度の設置、書籍の企画、編集など、さまざまな普及啓発活動や技術指導も行っています。こうした活動は、会員をはじめとする多くのボランティアによって支えられています。
主な著書としては、『日本を救う「最後の選択」』（情報センター出版局）、『ビオトープネットワーク』（ぎょうせい）、『ビオトープネットワークII』（ぎょうせい）、『環境を守る最新知識』（信山社サイテック）、『環境の時代を迎える世界の農業』（(財)日本生態系協会）、『学校ビオトープ』（講談社）、『環境アセスメントはヘップ（HEP）でいきる』（ぎょうせい）などがあります。
連絡先〒171-0021

東京都豊島区西池袋2-30-20 音羽ビル

TEL:03-5951-0244

FAX:03-5951-2974

http://www.ecosys.or.jp/eco-japan/
コメントの送付について本モデルに関するご意見等がございましたら、所属、氏名、連絡先をご記入の上、下記送付先までコメントをお寄せ下さい。貴重な情報として、モデル改訂時の参考とさせていただきますので、どうぞよろしくお願い申し上げます。
【コメント送付先】(財)日本生態系協会 ハビタット評価グループ
住所 〒171-0021

東京都豊島区西池袋2-30-20 音羽ビル

FAX 03-5951-2974
E-mail habitat@ecosys.or.jp※電子メールはテキスト形式にてご送付下さい（ファイルの添付、html形式のメールはご遠慮下さい）。
※個別のご質問についてはお答え致しかねますので、予めご了承下さい。

モデル作成（財）日本生態系協会ハビタット評価グループ

ミドリシジミのHSIモデル ver.1.0
2004年6月 発行
編集

財団法人日本生態系協会

発行

財団法人日本生態系協会〒171-0021
東京都豊島区西池袋2-30-20 音羽ビル

＊禁無断転載・複製c（財）日本生態系協会 2004