Top Banner
49

N t +1 /N t = r t = exp[ R m (1 - N t /K )] (1.6) =λ exp(- cN t ) (1.6a)

Dec 31, 2015

Download

Documents

daphne-daniel

N t +1 /N t = r t = exp[ R m (1 - N t /K )] (1.6) =λ exp(- cN t ) (1.6a) c = R m /K, R m =lnλ. . ロジスチック式の拡張( k (競争による増殖率の低下率)の密度による変化) ロジスチック式では,ある個体の周りの面積 s にいる競争者が i 匹いたときの増殖率は r i =r 0 k i である。 k i を第 i 番目の競争者が加わることによって生じる r 0 の減少率とすると, r i =r o k 1 k 2 ・・・ k i 。 - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 2: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 3: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Nt+1/Nt = rt

= exp[Rm(1 - Nt   /K)]    (1.6)

=λexp(-cNt) (1.6a)

c = Rm/K, Rm =lnλ

Page 4: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

.ロジスチック式の拡張( k(競争による増殖率の低下率)の密度による変化) ロジスチック式では,ある個体の周りの面積 s にいる競争者が i 匹いたとき

の増殖率はri=r0ki である。

ki を第 i 番目の競争者が加わることによって生じるr 0 の減少率とすると, ri=rok1k2 ・・・ ki 。1) 1 > ki-1 > ki > 0 の場合 これは,競争者が増えるにしたがって競争が激しくなることを示す。2) 0 < ki-1 < ki < 1 の場合 これは,競争者が増えるにしたがって競争が緩和することを示す。1)と2)の効果を組み込む一つの方法は,1 .6a(ただし λ=r0 , Nt = xt

)で xt を xta ( a は正の定数)で置き換えることである.

a > 1のときは1)の場合を表わし, a < 1の場合は2)の場合を表わす。   ( ln r0 ≡ Rm )

)15.1()exp(0att cxrr

Page 5: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

N-N増殖曲線ある世代の密度と次世代の密度との関係を表わしたもの(同世代内のある発育段階とそのあとのある発育段階の関係を表わしたものを競争曲線と呼ぶ)。

R-N増殖曲線密度と増殖率の関係を表わした曲線

増殖曲線(再生産曲線, reproduction curve)

Page 6: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

コンテスト型(勝ち残り型)

スクランブル型(共倒れ型)

N-N増殖曲線(再生産曲線, reproduction curve)

Page 7: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

X t

R t  =Xt+1/X t

R-N増殖曲線

Page 8: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

  rt = r0 /(1 + cxta) (1.16)

      (Maynard Smith & Slatkin,

1973)

   rt = r0 /(1 + cxt)a (1.17) 

       (Hassell, 1975)

今まで述べた式では,コンテスト型の増殖曲線を描くことが出来ない

Page 9: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

競争者の増加に伴い競争の激しさが加速度的に増える

競争者の増加に伴う競争の激しさの増加率は減少する

Page 10: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

rt = r0 /(1 + cxt)a (1.17) を利用 

昆虫個体群

Page 11: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

rt = r0 /(1 + cxta) (1.16)を利用

昆虫個体群

Page 12: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

平衡密度一定( 昆虫個体群の場合餌量で決まる)

制御 regulation

個体群の維持 population persistence  1.密度依存要因による個体群の維持  

Page 13: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

平衡密度一定

平衡密度変化(トレンド)

制御 regulation

制御

  平衡密度変化(トレンド)の制御とトレンドに対する個体群のの制御 ( 動的平衡)

( 昆虫個体群の場合餌量で決まる)

Page 14: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

密度依存性による変動要因の分類

密度依存要因 密度逆依存要因

密度独立要因

Page 15: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

危険の分散 (spreading the risk)

X将来、今いるパッチ、生息地が維持される保証がなければ、密度の高い低いに拘わらず移動分散する必要がある。

Page 16: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

オレンジあたりの密度変化

ハダニ

カブリダニ

カブリダニの移動分散をワセリンを使って制限する。 ハダニは,糸を吐きそれをパラシュート代わりにオレンジ間を移動。この分散を助けるため棒を立てる。

危険の分散

Page 17: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

密度依存要因と危険の分散の相対的重要性は、種によって、時には同種内の個体群の間でも異なっていると考えられる。

Page 18: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

生命表を基に個体群変動の機構を解明することを生命表分析と言う。

特に,相関と回帰分析を使った,変動を起こしている要因の解明と個体群制御要因の解明を,変動主要因分析 (key factor analysis) と言う。

III .個体群変動の解析 (Analysis of Population Fluctuations)

生命表分析Life table analysis

Page 19: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Winter moth ナミスジフユナミシャク

Page 20: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Winter moth ナミスジフユナミシャク

Page 21: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

生命表

Page 22: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

年 75 76 77成虫 X(1,75) X(1,76) X(1,77) 平均 V卵 X(2,75) X(2,76) X(2,77) 平均 (V 分散)

幼虫 X(3,75) X(3,76) X(3,77) 平均 V蛹 X(4,75) X(4,76) X(4,77) 平均 V

成虫 X(5,75) X(5,76) X(5,77) 平均 V

成虫から卵X(2,75)-X(1,75)

X(2, )-76X(1, )76

X(2, )-77X(1, )77 平均 V

卵から幼虫X(3,75)-X(2,75)

X(3, )-76X(2, )76

X(3, )-77X(2, )77 平均 V

幼虫から蛹X(4,75)-X(3,75)

X(4, )-76X(3, )76

X(4, )-77X(3, )77 平均 V

蛹から成虫X(5,75)-X(4,75)

X(5, )-76X(4, )76

X(5, )-77X(4, )77 平均 V

変動主要因分析  Key factor analysis

a.対数変換及び密度と変化率の平均値と分散の算出

Page 23: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Egg Larva 1 Larva 5 Pupa Adult

Developmental stage

0

1

2

3

4

5V

ari

an

ce

DenistyDenisty

r

密度撹乱要因 密度安定要因

変化率

Page 24: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Varley-Gradwell法

K= -Rk i = -ri

b.変動主要因の検出

Page 25: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Podoler-Rogers法

R

rx

Page 26: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

回帰分析による密度依存性の検出:普通の回帰分析では密度依存性がないのに密度依存性が検出されることがある。

Page 27: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

b,通常の計算法から得る回帰係数n,データー数5%水準 : μγ  = 1.9601%水準: μγ  = 2.576

μγ をいれて計算してbcの範囲が負であればそれぞれの水準で有意。

Kuno(1991)の方法

Page 28: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Bulmerの第1検定法

B = V / U を計算。

B < Bα(B0.05 = 0.25 + (n - 2)0.0366; B0.01= 0.25 +

(n - 2)0.0248;B0.001 = 0.25 + (n - 2)0.0170) のとき,危険率 αで密度依存性があると言える。

97 98 99 2000 01 02 03

X97 X98 X99 X00 X01 X02 X03

2

1

1

1

21

)(

)(

XXV

XXU

n

tt

n

ttt

Page 29: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

生命表分析を行なう上での注意点・ R は,同じ発育段階の個体数の世代変化率を用いるべきである。・年に2世代以上を経過する昆虫では, R として1年間の総変化率を用いるべきである。・草食昆虫の場合,変動要因のほとんどが, 1m2 当たりの密度といった絶対密度に反応するのでなく,植物当たりの密度に反応することが多い。そこで,植物当たりの密度の変化を変動主要因分析することのほうが,変動機構をより理解できることがある。このとき,植物側の変化を独立した変動要因として扱い,その r 値を log(xi/Pi) - log(xi/Pi-1)(Pi,  世代 i の植物量 )とすれば,植物の変化が及ぼす影響を評価出来る(Yamada 1996) 。

Page 30: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

Pi

Pi+1 Pi+2Plant

xi

Xi+1

X i +2Insect

1

11

1

logloglog

i

i

i

i

i

i

i

i

P

P

x

x

PxP

x

R

Page 31: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

• 植物量の年次変化が変動主要因となることがある。

• 昆虫 / 植物の変動は昆虫だけの変動をみるより,密度依存性が検出されやすい。

• 昆虫 / 植物の変動は,昆虫だけの変動より,過剰密度依存が起こり,連続振動,カオスになることが多い。

Page 32: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

生命表分析を行なう上での注意点(続き)

・通常は, R に密度依存性が発見された場合,第 1変動主要因に密度依存性が発見されるが,密度依存性の様式によっては,第 1 変動主要因に密度依存性が見られず,第 2 以下の変動主要因に密度依存性が見られる場合が理論的に考えられる (Royama 1996)。

これが起こるのは, Royamaでは R と x の関係を非線形モデルで表しているからである。すでに示したように, R と x の関係は線形(直線関係)ではなく,非線形と考えられる。密度依存に対する貢献度を評価する新しい解析法の誕生が待たれる。少なくとも, R―N 増殖曲線の形状を調べることが重要である。

Page 33: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

ヒトリガ Tyria jacobaeae

Page 34: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

ヤコブボロギク Senecio jacobaea

Page 35: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 36: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 37: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

ヒトリガのR-N増殖曲線

絶対密度

植物あたりの相対密度

密度は振動が続くことが予測される

Page 38: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

ヤマトアザミテントウ(コブオオニジュウヤホシテントウ) Epilachna niponica

Page 39: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 40: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

産卵の抑制 ー> 植物の回復待ち、新鮮な餌植物を求めての分散、第2回目の越冬

Page 41: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

ヤマトアザミテントウ のR N曲線ー越冬成虫

成虫と幼虫の間に餌をめぐる競争 成虫密度が多いと、次年度の成虫密度を過度に減らす

絶対密度

植物あたりの相対密度

*洪水

Page 42: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

シジュウカラ

                                                

    

Page 43: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 44: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

冬期から翌春の繁殖開始までの死亡

変動主要因越冬期の死亡:ヒナが巣立ってから,翌年の繁殖期迄の変動=冬期までの死亡 +冬期から翌春の繁殖開始までの死亡

他の要因 一腹卵数 孵化率 ひなの死亡

Page 45: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

3 日後

縄張り (territory):

取り除く

Page 46: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

縄張り (territory): 1個体あるいは複数個体によって独占的に使用される空間。通常,異種個体あるいは同種の他個体に対して防御行動をとる。縄張りのメリット:1)資源の独占使用。2)間おき効果による病気,捕食の軽減。3)縄張り内の熟知による捕食からの逃避。4)個体群の安定。

4)は結果として起こるのであって,このために縄張り制が進化したのではないと考えられている。よく似た用語に,行動圏 (home range )があるが,他個体を排除しない点で,縄張りと異なる。

Page 47: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)
Page 48: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

被食者の増加↓捕食者にとっての餌の増加↓捕食者の増加↓捕食圧の増加↓被食者の減少↓捕食者の減少

Page 49: N t +1 /N t  =  r t = exp[ R m (1 -  N t /K )]  (1.6) =λ exp(- cN t )         (1.6a)

植物(潅木若芽)の質悪化

カンジキウサギ

ヤマネコ

ライチョウ