Modello per
Introduzione
Gli insetti sono animali "eterotermi" o a "sangue freddo". In termini concreti ciò significa che essi risentono fortemente delle variazioni termiche ambientali che condizionano, rallentandolo o accelerandolo, il loro sviluppo. Tale caratteristica è comune anche alle piante e ciò determina uno stretto legame tra sviluppo dell'ospite e del fitofago. Lo sviluppo avviene in un intervallo di temperature favorevoli, limitate da una soglia termica inferiore e da una superiore, al di sotto ed al di sopra delle quali si ha l'arresto dello sviluppo. Allontanandosi ulteriormente dall'intervallo termico favorevole allo sviluppo si incontrano le soglie termiche letali, al di sopra ed al di sotto delle quali è impossibile la sopravvivenza degli individui. Queste soglie possono essere determinate attraverso opportune indagini di laboratorio o con osservazioni in campo.
Modelli
di dinamica di popolazione
Si definisce popolazione un insieme di individui appartenenti ad una singola specie che esplicano la loro attività in un determinato luogo e che sono tra loro legati sia in termini riproduttivi che in termini funzionali (interazioni trofiche, competitive o cooperative). Le popolazioni di insetti sono caratterizzate da una struttura demografica in stadi di sviluppo, pertanto la loro dinamica deve essere descritta in termini di variazione della densità relativa ad ogni stadio di sviluppo. Le popolazioni sono soggette a variazioni dell'abbondanza del numero di individui dovute a fenomeni di natura demografica (natalità e mortalità) e a fenomeni migratori (immigrazione ed emigrazione). La variazione del numero di individui di una popolazione è descritta in termini matematici dalla seguente relazione per intervalli discreti di tempo:
(1.1)
dove: n(t) è il
numero di individui al tempo t; b(t) è il tasso di
natalità; m(t) è il tasso di mortalità; i(t) è il
tasso di immigrazione; e(t) è il tasso di emigrazione e Dt rappresenta
l'intervallo temporale.
La variazione nel tempo della struttura demografica di una popolazione dipende
dalla variazione del numero di individui presenti in
ciascuno stadio del suo ciclo di vita. La variazione del numero di individui in un generico stadio s al tempo t
è data dagli individui che si sviluppano provenienti dallo stadio precedente
meno quelli che, a causa dello sviluppo, escono dallo stadio in questione,
sommati agli individui che muoiono. In termini matematici, escludendo la
variazione dovuta a fenomeni migratori, il fenomeno è
descritto dalla relazione:
(1.2)
dove: ns(t) è il numero degli individui dello stadio s al tempo t; t è l'intervallo temporale; rs-1(t) è il tasso di reclutamento allo stadio s degli individui dallo stadio s-1 al tempo t; rs(t) è il tasso di sviluppo degli individui dello stadio s al tempo t e ms(t) è il tasso di mortalità dello stadio s.
Per una
corretta analisi della dinamica di popolazione di
specie singole è inoltre necessario considerare la diversità genetica esistente
tra gli individui, che è una delle componenti all'origine della variabilità
biologica. Tale diversità risulta evidente quando
individui che sperimentano uguali condizioni ambientali manifestano una diversa
risposta fisiologica e comportamentale.
La variabilità genetica interviene influenzando tutti i processi biologici che
sono implicati nella dinamica della popolazione,
andando in definitiva a modificare:
- la natalità: il numero di uova che ciascuna femmina depone nel corso del periodo riproduttivo;
- la mortalità: la differente probabilità di morte di ciascun individuo dovuta sia ad avversità biotiche che abiotiche;
- il tasso di sviluppo: la diversa velocità (o tasso) di accrescimento di ciascun individuo;
- i tassi di emigrazione e di immigrazione.
Abbiamo visto che uno dei parametri che influenza lo sviluppo degli insetti è la temperatura. I parametri della relazione che descrive il tasso di sviluppo in funzione della temperatura sono stimati da dati ottenuti, in genere, con esperimenti di laboratorio allevando gli insetti a temperatura costante. In questo modo si ottengono dapprima i tempi di sviluppo D alle differenti temperature di laboratorio T e, dai loro reciproci, le velocità di sviluppo in rapporto alla temperatura di allevamento, per cui si avrà che: v=1/D. Nell'intervallo termico in cui si può ritenere che il tasso di sviluppo sia linearmente dipendente dalla temperatura si può eseguire un fitting dei dati sperimentali e costruire l'equazione di regressione lineare v=a+bT in grado di rappresentare il fenomeno. Per il calcolo della soglia termica inferiore di sviluppo Tinf si risolve l'equazione di regressione per v=0, quindi Tinf = -a/b. L'accumulo termico espresso in gradi-giorno GDD, necessario al completamento dello sviluppo, viene calcolato come differenza tra il valore di T per v=1 e quello di T per v=0; quindi GDD =1/b. Con questo metodo di calcolo, i gradi-giorno non sono dipendenti dalla temperatura e quindi il modello fenologico, anche se realizzato con esperimenti di laboratorio a temperature costanti, può essere adottato per simulare i tempi di sviluppo a temperature di campo. La sommatoria termica viene calcolata di norma fra l'inizio ed il termine di una generazione, tranne che per la generazione diapausante (quando è presente), per la quale viene presa come data di inizio del calcolo il 1 gennaio o una data successiva stabilita empiricamente. Nel metodo di calcolo usato tradizionalmente vengono sommate le differenze tra la temperatura media giornaliera e la soglia termica inferiore di sviluppo. E' evidente che il valore della somma termica dipenderà dalle soglie di temperatura prescelte per il calcolo ed inoltre, a parità di soglia, le somme termiche potranno risultare differenti per osservazioni fenologiche riguardanti anni diversi. Pertanto, per un uso operativo dei modelli gradi-giorno per la simulazione fenologica dello sviluppo degli insetti, è necessaria un'ampia sperimentazione in campo, in particolare nelle aree climatiche diverse da quelle in cui sono state compiute le osservazioni biologiche e meteorologiche.
·
Gli adulti volano al crepuscolo (attività
crepuscolare) e sono favoriti da temperature superiori ai
·
La durata degli stadi pre-immaginali
è massima a
·
La mortalità della Lobesia è crescente a valori
superiori ai
·
La fertilità, risultata
da esempi di laboratorio, variava fra 9 uova, ad una temperatura di ovideposizione di
In estrema
sintesi le condizioni abiotiche avverse sono: uova,
larve e crisalidi resistono a
Il modello
Il modello descritto in questa sezione costituisce la parte fondamentale di uno strumento operativo in grado di elaborare previsioni dello sviluppo fenologico delle popolazioni di Lobesia Botrana. Tali previsioni verranno utilizzate per l'emissione di bollettini fitopatologici diffusi via Internet e destinati ai tecnici che operano nell'assistenza alle aziende. Il modello permette di stabilire le epoche di emergenza dei diversi stadi del fitofago utilizzando due sub-modelli. Il primo simula la fenologia degli stadi preimmaginali sulla base di esperimenti di dinamica di sviluppo condotti in laboratorio. Il secondo simula la fenologia dello stadio adulto attraverso l'uso delle sommatorie termiche calcolate a partire da dati di campo. Entrambi i sub-modelli richiedono in ingresso le temperature dell'aria massime e minime giornaliere per il periodo considerato.
Fenologia
degli stadi preimmaginali
Nella figura 1 viene riportata la velocità di sviluppo delle varie fasi del ciclo vitale della Lobesia Botrana: uovo, larva e crisalide, in funzione della temperatura.

fig. 1
Il calcolo della velocità di sviluppo degli stadi di uovo, larva e crisalide viene effettuato attraverso il modello proposto da Logan (1976):
![]()
(1.3)
dove: v
rappresenta la velocità di sviluppo in funzione della temperatura dell'aria T. Nel modello si è introdotta l'approssimazione di
considerare il tasso medio giornaliero di sviluppo come pari alla velocità di sviluppo calcolata con l'equazione (1.3) ad una temperatura
T pari alla temperatura media giornaliera. I coefficienti a, b e c (Tabella 1),
sono stati calcolati a partire da esperimenti sullo
sviluppo dell'insetto condotti in laboratorio a temperature costanti (Rapagnani et al., 1988). Le
soglie termiche di sviluppo Tsup e Tinf sono state determinate sperimentalmente e sono pari,
rispettivamente, a
Tabella 1: Coefficienti dell'equazione della velocità di sviluppo (1.3) per gli stadi preimmaginali di Lobesia Botrana.
|
Stadio |
a |
b |
c |
|
Uovo |
0,297370 |
0,183374 |
0,187975 |
|
Larva |
0,225958 |
0,195681 |
0,197009 |
|
Crisalide |
0,429051 |
0,211920 |
0,212915 |
Fenologia
degli adulti
Lo sfarfallamento degli adulti viene simulato impiegando un modello proposto da Arca (1993) già sottoposto a validazione in alcune aree della Sardegna (Nurra e Campidano). In questo modello l'accumulo dei gradi-giorno (GDD), utili ai fini dello sviluppo degli adulti, viene calcolato tenendo conto della diversa posizione del ciclo delle temperature giornaliere rispetto alle soglie termiche di sviluppo di Lobesia Botrana. In relazione a 6 diverse situazioni individuate, si definiscono le seguenti 6 equazioni per il calcolo dei GDD, proposte da Allen (1976):
caso 1: Tmin
> Tsup
(A)
caso 2: Tmax < Tinf
(B)
caso 3: Tmin
> Tinf;
Tmax < Tsup
(C)
caso 4: Tmin
< Tinf
; Tinf<
Tmax < Tsup
(D)
caso 5: Tmax
>Tsup;
Tinf<
Tmin < Tsup
(E)
caso 6:Tmax > Tsup ; Tmin < Tinf
(F)
dove:
GDD
= gradi-giorno;
Tmax = temperatura
massima giornaliera;
Tmin = temperatura
minima giornaliera;
;
;
;
;
= soglia termica inferiore di sviluppo (
= soglia termica superiore di sviluppo (
I
valori di sommatoria termica individuati per le principali fasi del volo degli
adulti delle diverse generazioni, stimati per le diverse aree in cui è stata
fatta la sperimentazione, vengono riportati nella
Tabella 2.
Tabella 2: Valore medio e
deviazione standard in GDD delle sommatorie termiche, relative ai voli di L. botrana nella Nurra di Sassari e nel Campidano di Oristano
e Cagliari, per le diverse generazioni (Arca et al.,
1993).
|
Generazione |
Inizio volo |
Picco volo |
Fine volo |
|
1 |
351 ± 28,1 |
442 ± 23,8 |
518 ± 28,9 |
|
2 |
871 ± 44,2 |
1048 ± 61,28 |
1150 ± 67,7 |
|
3 |
1570 ± 61,4 |
1865 ± 35,4 |
1948 ± 42,53 |
|
4 |
2096 ± 61,8 |
2444 ± 8,061 |
2528 ± 2,404 |
Simulazioni
sulla Calabria
Di seguito vengono riportati i grafici ottenuti dalle simulazioni dello sviluppo della Lobesia Botrana nella regione Calabria. Questi grafici sono stati prodotti, per validare il modello, utilizzando i dati di temperatura delle stazioni di Palmi, Spezzano Albanese e Stalettì, e le temperature per cui si hanno i passaggi di fase nel ciclo biologico della Lobesia Botrana, ricavate dai dati trovati in bibliografia. In questo modo è stato possibile ricavare il giorno giuliano dello sviluppo delle diverse generazioni dell’insetto.
L’uscita del modello prevede, per ogni stazione e per ogni anno, un grafico ed una tabella. Sul grafico vengono riportate le gdd in funzione del giorno giuliano (linea nera continua). Vengono riportate, inoltre, le gdd relative alle diverse fasi di volo delle diverse generazioni di Lobesia Botrana secondo quanto riportato nella tabella 2. Quando si hanno delle intersezioni tra le gdd ottenute dai dati di stazione e quelle di tabella 1 viene riportata la data per la quale si ha l’intersezione corrispondente a quella particolare generazione (fase del volo). Nella tabella vengono riportate, se applicabili, le date per le quali si ottengono le interazioni tra le gdd misurate dalle stazioni e quelle di tabella 2.






BIBLIOGRAFIA
meridionale. Informatore Fitopatologico
n. 5: 8-11