BASE


====== Distúrbio e Coexistência ======
{{noiteestrelada.jpg?200  |}}
Vamos investigar nesse exercício modelos de sucessão e incorporar neles distúrbio. No exercício de coexistência de duas populações simulamos distúrbio através da taxa de extinção. Vamos partir desse modelo e avaliar o que acontece com a coexistência conforme aumento da intensidade do distúrbio.  
Vamos relembrar o modelo de coexistência e o exercício que fizemos [[ecovirt:roteiro:metap_duas:metap_coexrcmdr|Coexistência em Metapopulações]]. Tínhamos duas espécie, sendo a primeira  melhor competidora pois podia colonizar manchas já ocupadas pela outra, enquanto a segunda só colonizava manchas vazias. A variação na ocupação de manchas era dada por:

$$  (df_1)/dt  =  i_1f_1(1-f_1)-p_e f_1 $$

$$  (df_2)/dt = i_2f_2(1-f_1-f_2)-i_1f_1f_2-p_e f_2   $$

onde:
  * f = fração de manchas ocupadas
  * pe = probabilidade de extinção por mancha
  * i = taxa de incremento da probabilidade de colonização com o aumento de f
  * A taxa de colonização, portanto, é o produto i.f, e varia com a fração de manchas ocupadas (quanto mais ocupação, mais propágulos).

Para que a espécie dois persistisse no sistema era necessário que:

$$ e/i_{(1)} > i_{(1)}/i_{(2)} $$
 
 
====== EcoVirtual ======



Segue abaixo a descrição dos parâmetros do modelo para relembrarmos:


^ opção ^  parâmetro ^definição ^
^ data set |objeto no R | guarda os resultados|
^ Maximum time |tmax |Número de iterações da simulação |
^ coluns | cl |número de colunas de habitat da paisagem |
^ rows | rw |número de linhas de habitat da paisagem |
^ Best Competitor ^ Parâmetros para a melhor competidora ^^
^ initial occupance | f01 |proporção de manchas ocupadas no inicio pela sp1|
^ colonization coef. | i1 |coeficiente de colonização i das sp1 |
^ Inferior Competitor ^ Parâmetros para a pior competidora ^^
^ initial occupance |f02 |proporção de manchas ocupadas no inicio pela sp2|
^ colonization coef. | i2 |coeficiente de colonização i da sp2| 
^ Both Species ^ Parâmetros comuns para as duas espécies ^^
^ prob. extinction | pe |probabilidade de extinção |
^ Habitat Destruction | d |proporção de manchas não disponibilizadas|
^ Show simulation frames | //anima=TRUE// |mostra cada tempo simulado | 

===== Simulando Distúrbio =====

Vamos simular um aumento crescente da probabilidade de extinção fixando os seguintes parâmetros:
<code>
tmax=100;
cl=100;
rw=100;
f01=0.1;
f02=0.1;
i1=0.4;
i2=0.8;
d=0;
anima=FALSE
</code>

Agora vamos variar a probabilidade de extinção para simular uma aumento de distúrbio:

  - //pe// = 0.1
  - //pe// = 0.2
  - //pe// = 0.25 
  - //pe// = 0.3
  - //pe// = 0.4
  - //pe// = 0.5 


  - O que está acontecendo com o sistema conforme aumentamos a intensidade do distúrbio (pe)?
  - Em que cenário o sistema é mais diverso? 
  - Interprete o comportamento da espécie pior competidora no início das quatro primeiras simulações. O que está acontecendo?