Baseado na apostila de H.T. Odum, Univ. of Florida,
e E. Odum, Santa Fe Community College.
Minimodelos para simular o comportamento de ecossistemas em computador. Gainesville, Fl., EUA, 1994. http://www.unicamp.br/fea/ortega/eco/simul/indice-p.htm |
Apresentamos agora
os seguintes modelos básicos:
(a) Crescimento a partir de uma fonte renovável de reposição lenta. (b) Fatores limitantes externos (c) Fatores limitantes internos (reciclagem de nutrientes) (d) Sistema com produção e consumo (respiração) diários |
Principais tipos de Fontes e
formas de Captura e Consumo de Energia
7. Crescimento a partir
de uma fonte renovável de reposição lenta
Programa slowren Este sistema no início aproveita um estoque grande (E) resultante de um acúmulo prévio de J sem uso até o momento. O consumidor Q cresce exponencialmente aumentando seu estoque de energia. Quando o estoque E mingua o crescimento de Q para e depois decresce, ele aproveita durante o declínio seu próprio estoque acumulado (Q). Depois surto ou pico de crescimento o sistema se ajusta a um "modelo sustentável" de menor energia que corresponderia a taxa de acumulo em E do fluxo externo (J).
Não ha interações que afetem
o fluxo de energia que entra, por tanto, o estoque E muda
diretamente com o fluxo J. Assim sendo a entrada de energia
no intervalo DT fica como:
O estoque E tem um dreno (perda simples, por exemplo:
mortalidade) que varia em proporção direta com esse
estoque e também tem a saída que alimenta Q cujo fluxo varia
em proporção dos dois estoques (Q*E) envolvidos nessa
interação..
A variação do acumulo de E seria:
No caso da mudança do estoque Q teríamos,
no intervalo DT:
Se o intervalo de tempo fosse unitário
teríamos DT = 1.
Inicio do programa
Imprimir ou plotar os valores dos estoques e o tempo respectivo Calcular o incremento dos estoques no intervalo DT: DE = J*DT - K4*Q*DT - KO*E*Q*DT DQ = K1*Q*E*DT - K3*Q*DT Calcular os novos valores dos estoques, somando ao valor anterior o incremento: E = E + DE Q = Q + DQ Incrementar o valor de T: T = T + DT Testar o valor de T para saber se estamos no do intervalo de tempo permitido T<TMAX? Em caso positivo fique dentro do laço senão saia dele Fim do programa Um exemplo deste modelo seria o crescimento de
peixes em uma represa construída em um local onde havia uma floresta
No início a pesca é abundante devido ao uso do estoque de
matéria orgânica residual depois diminui e se ajusta em um
patamar que corresponde ao regime sustentável em recursos naturais
renováveis.
|
10
REM PC: SLOWREN (Slowly renewable)
20 SCREEN 1, 0: COLOR 0, 0 35 LINE (0, 0)-(319, 180), 3, B 40 Q = 3 50 J = 2 60 E = 159 70 K0 = .001 80 K1 = .001 90 K3 = .03 100 K4 = .01 110 Q0 = 1 |
120 E0 = 1.2
130 DT = 1 140 T0 = 1 200 PSET (T * T0, 180 - Q * Q0), 1 210 PSET (T * T0, 180 - E * E0), 2 220 DQ = K1 * E * Q - K3 * Q 230 DE = J - K0 * E * Q - K4 * E 240 E = E + DE * DT 250 Q = Q + DQ * DT 260 T = T + DT 270 IF T * T0 < 319 GOTO 200 |
8. Fatores limitantes externos
|
J = 5
I = 1 |
Z = 2
DI = 2 |
K1 = .08
K2 = .05 |
K3 = .01 |
10 REM PC: FACTORS (External
limiting factors)
20 SCREEN 1, 0: COLOR 7, 0 40 LINE (0, 0)-(240, 180), 3, B 50 J = 5 60 E = 1 80 K1 = .08 90 K2 = .01 110 DE = 2 115 E0 = 2: REM fator de escala 118 P0 = 1: REM fator de escala |
120 M = J / (1 + K2 * E)
130 P = K1 * E * M 140 E = E + DE 150 PSET (E / E0, 180 - P / P0), 1 160 IF (E / E0) < 240 GOTO 120 170 J = J + 5 180 E = 0 190 IF J < 25 GOTO 120 |
9. Fatores limitantes internos
|
Outra equações:
Produção P = K1*I*N Saída = K2*Q Acúmulo DQ = P - K2*Q |
No estado estacionário quando
a taxa de acumulo se reduz totalmente temos DQ=0, nesse caso:
K1*I*N = K2*Q K1*I(TN-F*Q) = K2*Q K1*I*TN-K1*I*F*Q = K2*Q K1(I*TN) = K2*Q+K1*I*F*Q K1(I*TN) = Q(K2+K1*F) Q = (K1*I*TN)/(K2+K1*F) |
TN = .5
I = 0 |
Q = 5
Z = .05 P0 = .1 |
K1 = 2
K2 = .5 |
F = .05
DI = .05 |
3 REM RECYCLMT (Recycle
Limiting Factor)
4 SCREEN 1, 0: COLOR 0, 0 6 LINE (0, 0)-(319, 180), 1, B 10 TN = .5 15 I = 0 17 Q = 5 20 Z = .05 22 P0 = .1 30 K1 = 2 40 K2 = .5 45 F = .05 47 DI = .05 |
50 Q = (K1 * TN * I) / (K2 + K1 * F * I)
55 N = TN - F * Q 57 P = K1 * I * N 60 I = I + DI 70 PSET (I / Z, 180 - P / P0), 2 80 IF (I / Z) < 319 GOTO 50 81 BEEP: SLEEP: END 82 REM type CONT for more material to recycle 85 Q = 1 90 TN = TN + .5 95 I = 0 100 Q = 1 110 IF TN < 40 GOTO 50 |
10. Produção e Consumo. Programa daypr![]() Nos processos ecossistêmicos onde ocorre fotossíntese temos uma oscilação diária de energia externa, o que origina um ciclo de produção e consumo (ou respiração). A energia solar é uma fonte renovável limitada cujo potencial é estabelecido pela energia não utilizada JR e varia como uma senoide truncada. Um pico que cresce e decresce durante o dia e de noite é nula. As equações seriam as seguintes: J = 40*SIN(T/15.9) IF (J<)) THEN J=0 J = JR + K1*JR*N A biomassa e o oxigênio são produzidos de dia e consumidos de noite. As plantas geram um excedente que é pela cadeia trófica. No caso de um ecossistema florestal teríamos: TN = quantidade total de nutrientes N = nutrientes no solo disponíveis para uso pela plantas K3 = proporção da biomassa Q composta por nutrientes Equação do balanço de nutrientes: TN = N + K3*N Acumulo em Q: DQ = K1*JR*N - K2*Q
|
10 REM DAYPC (Daily production
and consumption)
20 SCREEN 1, 0: COLOR 0, 0 35 LINE (0, 0)-(319, 60), 3, B 40 LINE (0, 70)-(319, 130), 3, B 50 LINE (0, 140)-(319, 190), 3, B 60 Q = 5 80 NT = 1 90 K = .9 100 K1 = .2 110 K2 = .04 120 K3 = .01 130 Q0 = .5 140 N0 = 50: REM .02 |
150 T0 = 1
160 DT = 1 170 J = 40 * SIN(T / 15.9) 180 IF J < 0 THEN J = 0 190 R = J / (1 + K * N) 200 N = NT - K3 * Q 210 DQ = K1 * R * N - K2 * Q 220 Q = Q + DQ * DT 230 PSET (T * T0, 130 - Q * Q0), 2 240 PSET (T * T0, 190 - N * N0), 1 250 LINE (T, 60)-(T, (60 - J)), 2 260 T = T + DT 300 IF T * T0 < 319 GOTO 170 |