Uma população em que as taxas de nascimento e mortalidade são constantes tem um crescimento ou decréscimo independente da densidade dela própria. Essa situação é geralmente relacionada à ausência de restrição ao crescimento, quando os recursos são ilimitados, mas pode também estar associada à extinção de populações.
Vamos imaginar agora uma população com taxas constante de crescimento e mortalidade e sem migrações. Vamos imaginar também que essa população cresce a intervalos regulares. O tamanho da população no próximo passo de tempo ($N_{t+1}$) é o número de indivíduos da geração anterior ($N_t$) mais o número de nascimentos (B), e menos o número mortes (D) no intervalo de tempo:
$$N_{t+1} = N_t + B - D $$
O número de mortes e nascimentos são resultado de taxas per capita multiplicadas pelo tamanho da população:
onde: $b$ = taxa de nascimento per capita a cada geração ; $d$ = taxa de mortalidade per capita a cada geração. Note que a taxa não muda com o tamanho da população, e que o número de nascimentos e mortes é proporcional ao tamanho populacional. Vamos apenas deixar claro mais uma premissa, para fins didáticos: os nascimentos e mortes ocorrem simultaneamente na população, no intervalo de tempo $t$. Podemos então dizer que :
se definimos um fator de crescimento discreto: $r_d = b-d$
Como $ 1+r_d $ é uma constante, vamos designá-la como $\lambda$, um número positivo que expressa o aumento proporcional da população de uma geração a outra. Portanto:
$$ \lambda=\frac{N_{t+1}}{N_t} \implies N_{t+1} = \lambda N_t$$
Podemos então projetar a nossa população a cada passo de tempo $t$. Por exemplo:
Se uma população com 100 indivíduos tem uma taxa per capita de natalidade de 0,8/ano e de mortalidade de 0,75/ano, qual o tamanho esperado da população no próximo ano?
$$N_{t+1}=100 \times (1+0,8-0,75) = 100 \times 1,05 = 105$$
Podemos também projetar a população para outras gerações, usando iterações:
prosseguindo e tomando o tamanho da população no tempo zero ($N_0$):
Generalizando:
$$N_{t}=N_0 \lambda^t $$
Assim, para nosso exemplo a projeção para 10 intervalos de tempo é
$$100 \times {1,05}^{10} = 162,8895$$
Com um pouco de manipulação algébrica a equação para tempo discreto
$$N_{t+1} \ = \ (1+r_d)N_t$$
Pode ser reescrita como
$$ N_{t+1} - N_t \ = \ \Delta N \ = \ r_dN_t $$
o que explicita que a velocidade de crescimento $\Delta N = N_{t+1} - N_t$ é proporcional ao tamanho poulacional $N_t$. Essa é característica essencial do crescimento populacional sem limites a uma taxa constante:
O fato básico da reprodução faz com que a variação do número de indivíduos de uma população no tempo seja proporcional ao número de indivíduos, em um ambiente constante.
A velocidade de crescimento é a variação no número de indivíduos $\Delta N$ dividida pelo tempo em que se deu esta variação, $\Delta t$. No modelo de tempo discreto $\Delta t = 1$, por definição. É um intervalo que pode ser tão grande como uma geração. Isso faz sentido se as mudanças no tamanho populacional se dão em intervalos discretos, como por exemplo em espécies semélparas com estação reprodutiva sincronizada.
No entanto, o tamanho populacionais da maioria das espécies varia continuamente. O cálculo diferencial foi criado para descrever esse tipo de dinâmica. Como o problema é que podem ocorrer mudanças a qualquer instante, o conceito-chave aqui é o de taxa instantânea, ou derivada.
Se nascimentos e mortes podem acontecer a todo momento, faz mais sentido pensarmos em uma velocidade instantânea do tamanho populacional. Isto equivale a reduzir $\Delta t$ tanto que pode ser considerado um instante. Esta velocidade instantânea é a derivada do tamanho populacional, que representamos com ${dN}/{dt}$, para diferenciar da velocidade a intervalos grandes e arbitrários $\Delta N/\Delta t$ 1).
Agora podemos expressar que velocidade instantânea de uma população é proporcional ao tamanho populacional com a equação:
$$\frac{dN}{dt} = rN$$
Que é o modelo de crescimento populacional a taxa constante em tempo contínuo. O parâmetro $r$ é chamado taxa instântanea de crescimento per capita. Essa taxa $r$ expressa o número médio de filhotes que cada indivíduo da população produz num intervalo de tempo muito curto. Por isso, muito livros de ecologia indicam que a unidade de $r$ é indivíduos/indivíduo.tempo. Físicos e matemáticos são mais rigorosos e lembram que a expressão correta da unidade é 1/tempo.
Para prever o tamanho de uma população que cresce a uma taxa constante em tempo contínuo usamos a equação
$$N(t) = N_0e^{rt}$$
E a relação entre a taxa de crescimento instantânea e a taxa de crescimento do modelo discreto é
$$ r \ = \ ln(\lambda)$$
As duas equações acima são deduzidas da equação $dN/dt=rN$ com técnicas de cálculo numérico. Se quiser entender um pouco mais sobre isso veja o roteiro sobre taxas, derivadas e função exponencial.
Para prosseguir você deve ter o ambiente R com os pacotes Rcmdr e Ecovirtual instalados e carregados. Se você não tem e não sabe como ter, consulte a página de Instalação.
Caso já tenha o R e pacotes instalados
Carregue o pacote principal RcmdrPlugin.EcoVirtual pelo menu do R Pacotes > Carregar Pacotes, ou pela linha de comando com o código:
library("RcmdrPlugin.EcoVirtual")
Vamos usar o EcoVirtual para projetar o tamanho de populações que crescem a taxas constantes sem limite, em tempo contínuo e discreto.
Na janela do Rcmdr clique na opção de menu EcoVirtual, e depois em One population em seguida em Exponential growth. Uma janela como esta irá se abrir:
Clique em “aplicar” para gerar um gráfico como este:
Arranje as janelas de modo que o gráfico esteja sempre visível, mesmo quando você alterar valores na janela de opções.
Os parâmetros da nossa função de crescimento denso-independente são:
Opção | parâmetro | definição |
---|---|---|
Enter name for last simulation data set | objeto no R | nome para salvar os resultados da simulação em um objeto no R |
Maximum time | tmax | número de interações |
Interval time size | intt | divisões do intervalo de tempo para o modelo discreto |
Initial population size | N0 | tamanho inicial da população |
Population growth rate (lambda) | lambda | taxa de crescimento do modelo discreto |
O resultado da função será um gráfico com o tamanho da população em função do tempo previstos pelos modelos. Os pontos azuis indicam os tamanhos populacionais previstos pelo modelo de crescimento em tempo discreto:
$$N_t = N_0\lambda^t$$
usando os valores de $N_0$ e $\lambda$ da caixa de opções. A linha preta indica os tamanhos populacionais previstos pelo modelo de crescimento em tempo contínuo:
$$N(t) = N_0e^{rt}$$
usando os mesmos valores de parâmetros. Para isso, o $\lambda$ é usado para calcular a taxa de crescimento instantânea per capita correspondente do modelo contínuo pela relação:
$$r=ln(\lambda)$$
Os valores para o modelo discreto são pontos, porque este modelo prevê o tamanho da população a intervalos discretos. Já o modelo em tempo contínuo prevê o tamanho da população a qualquer momento, e por isso é representado por uma linha contínua.
Os pontos se sobrepõem à linha porque o EcoVirtual usa taxas de crescimento equivalentes para tornar as projeções comparáveis. Abaixo do eixo X do gráfico está o valor de $\lambda$ e de $r$ usados na simulação.
Em suas memoráveis aulas de dinâmica de populações, o físico Roberto Kraenkel costuma dizer que não sabe se o paraíso existe, mas em caso afirmativo sabe que só lá as populações crescem sem limites. Apesar disso, prossegue, os modelos para este tipo de crescimento são essenciais para entender a dinâmica de populações reais, assim como a irrealista primeira lei de Newton é essencial para entender o movimento dos corpos no espaço.
Então vamos usar o EcoVirtual para explorar o comportamento dos modelos de crescimento a taxas constantes!
Experimente diferentes valores da taxa de crescimento discreto $\lambda$ com a opção Population growth rate (lambda)
. Isso vai alterar também a taxa de crescimento contínuo $r$. Veja o efeito no gráfico e use uma calculadora para conferir a correspondência entre as duas taxas, indicadas abaixo do eixo X do gráfico.
Qual valor ou intervalos de valores de $\lambda$ e de $r$ sob o quais a população:
O parâmetro Interval time size (intt)
define um novo intervalo de tempo para o modelo de tempo discreto. O novo valor de intervalo é uma fração do original de valor um. Assim, se você escolher $0.5$ para esta opção as projeções do modelo discreto são recalculadas para intervalos que correspondem à metade do intervalo original. Por isso, você verá no gráfico pontos azuis a cada meia unidade de tempo. Alterando o intervalo para $0.25$ você terá projeções para intervalos que são um quarto da unidade original, e o gráfico terá pontos a cada 0,25 unidade de tempo.
Note que à medida que você diminui o intervalo de tempo do modelo discreto os pontos se aproximam, até que parecem formar uma linha contínua. Ou seja, a projeção a intervalos discretos tende à projeção em tempo contínuo à medida que os intervalos são reduzidos!
De fato, o modelo exponencial pode ser visto como um limite do modelo discreto. Os detalhes estão no já indicado roteiro taxas, derivadas e função exponencial. Mas para que isso funcione o EcoVirtual recalcula o valor de $\lambda$ e o correspondente $r$ para cada novo intervalo de tempo discreto. E aqui vai nossa pergunta: como são feitos estes cálculos? Mais precisamente:
Dado uma razão de crescimento discreto para um intervalo de tempo de valor de uma unidade, $$\lambda_1=\frac{N_{t+1}}{N_t}$$
como calcular a razão de crescimento para um intervalo fracionário
$$\lambda_{1/n}=\frac{N_{t+{1/n}}}{N_t}$$
de modo que ao final de uma unidade de tempo a razão de crescimento permaneça $\lambda_1$? Verifique sua solução contra os valores que o EcoVirtual retorna quando você reduz o intervalo de tempo discreto pela metade (intt=0.5).
Esse exercício utiliza os dados de censos do IBGE para modelar e fazer predições sobre o crescimento da população brasileira. Antes de continuar, baixe os arquivos de dados necessários:
Os modelos que iremos utilizar nesse exercício são:
$$N_t=N_0\lambda^t$$
$$N_t=N_0e^{rt}$$
E a transformação entre eles:
$$r=\ln(\lambda)$$ $$\lambda=e^r$$
Abra o arquivo censo90.csv em uma planilha eletrônica.
Os dados estão estratificados por gênero (Homens, Mulheres) e por local de residência (Urbano ou Rural). Para esse exercício só utilizaremos os dados de Homens e Mulheres (colunas 2 e 3) somados e representando o tamanho total da população. O primeiro passo é, portanto, calcular esse valor em uma nova coluna da tabela.
Para a população total, calcule:
Compare sua estimativa com o tamanho populacional observado no censo IBGE 2010 abaixo:
População no censo 2010: 169872856
Discuta com os colegas da bancada (trios) as possíveis fontes da diferença entre a estimativa do modelo e os dados do censo 2010, baseado nos pressupostos que estruturam o modelo. Anotem as explicações que o grupo julga plausíveis.
Vamos agora usar os dados de uma série temporal mais longa e ver como as taxas de crescimento anual estão se comportando ao longo do tempo.
Abra o arquivo censoDecadas.csv em uma planilha eletrônica. Os dados estão estruturados da mesma maneira, de forma que também será necessário calcular a população total somando homens e mulheres.
Discuta com os colegas de bancada (trios) como poderíamos modelar a demografia da população brasileira para incorporar o(s) problema(as) diagnosticado(os) e fazer previsões mais plausíveis. Anote uma proposta para ser discutida com a turma.