Volendo allora inserire in STRATA i dati di sorgente del terremoto dell'AQ, apriamo il file -source, andiamo a 'motions' e clicchiamo su outcrop 2A sotto 'type'. In questa maniera possiamo modificare i parametri sismologici.
Decidiamo di prendere coem riferimento iniziale il mainshock del 06/04/2009.
Mw è noto = 6.3 R=11 km è noto, a partire dalla più recente determinazione dell'ipocentro, ossia a poca distanza dalla stazione ferroviarai sulel pendici del Colle di Roio
Depth= 9 km è noto
Parameter region: custom
Stress drop = 92 bars (da articolo di Bindi et al., segnalato da Mucciarelli)
Geometric attenuation = 1/R = 0.09
"Data speak for themselves" -Reverend Thomas Bayes 1702-1761 P(Ai|E)=(P(E|Ai)P(Ai))/P(E)
path duration = 0.16 (valore indicato nel manuale, eq. 2.35
path attehuation: valori indicati nell'articolo di bindi et al.: a= 48=~50, b=0.60
shear velocity alla rottura = 3 km/s (stimata) density : 2.7 g/cc (stimata per basamento a 9 km) site attenuation: lascio di default=0.04, non avendo dati a disposizione
Duration: la durata non posso variarla ed è = 7.04 secondi, molto simile alla durata calcolata dalla I arias del segnale mainshock
"Data speak for themselves" -Reverend Thomas Bayes 1702-1761 P(Ai|E)=(P(E|Ai)P(Ai))/P(E)
a questo punto, dò l'OK e vedo che in PGA mi apapre il valore 0.18 g.
Il valore<valore di normativa=0.26 g, per cui devo manipolare qualche parametro.
Sceglo di variare il parametro ceh in realtà può essere davvero soggetto a variazione, la Mw e la pongo= 6.7, che in realtà è perfettamente legittimo (la struttura di cui si tratta può arrivare a generare Mw 7 e poco oltre).
a questo punto noto che PGA dello spettro generato dal programam = 0.26 = spettro di normativa.
Coem ultimo step devo solo controllare la compatibilità alle altre frequenze con lo spettro di normativa.
"Data speak for themselves" -Reverend Thomas Bayes 1702-1761 P(Ai|E)=(P(E|Ai)P(Ai))/P(E)
Copio i dati presenti nelal tabella 'RSdata', li incollo sul foglio di calcolo, lasciando una colonna per lo spettro di normativa cat 'A'. Il paragoen diretto, come evidente dal grafico allegato,è uno di ottima compatibilità. In questo caso abbiamo un modello di sorgente perfettamente realistico e uno spettro perfettamente norma-compatibile.
"Data speak for themselves" -Reverend Thomas Bayes 1702-1761 P(Ai|E)=(P(E|Ai)P(Ai))/P(E)
Seguo il discorso ma a questo punto mi viene un altro dubbio. Ma se hai scelto di usare l'eq 2.35 non sarebbe corretto usare come crustal amplification quello Eastern NA?
path duration = 0.16 (valore indicato nel manuale, eq. 2.35
mccoy ma l'equazione 2.35 non dice 0.16R?
saluti, agi
In effetti la posizioen della virgola suggerirebbe quello che dici, ma se provi a variare il valore si ferma al massimo di 0.2. Prova anceh tu, se confermi credo proprio che quel valore sia un 0.16
"Data speak for themselves" -Reverend Thomas Bayes 1702-1761 P(Ai|E)=(P(E|Ai)P(Ai))/P(E)
Seguo il discorso ma a questo punto mi viene un altro dubbio. Ma se hai scelto di usare l'eq 2.35 non sarebbe corretto usare come crustal amplification quello Eastern NA?
saluti, agi
Il tuo ragionamento è logico, ma se ci provi esce una PGA troppo bassa e non si riesce a portarla a 0.26 neanche con un Mw=7.06 A questo punto dovremmo forse pensare che il modello di amplificazione eastern NA non è adatto all'Appennino, mentre i corrispondenti Tp possono essere usati rimanendo entro ouputs realistici
"Data speak for themselves" -Reverend Thomas Bayes 1702-1761 P(Ai|E)=(P(E|Ai)P(Ai))/P(E)
In effetti la posizioen della virgola suggerirebbe quello che dici, ma se provi a variare il valore si ferma al massimo di 0.2. Prova anceh tu, se confermi credo proprio che quel valore sia un 0.16
E' così mccoy non va oltre lo 0.2 anche provando a inserirlo manualmente. "Errore" nel manuale mi verrebbe da dire..
A questo punto dovremmo forse pensare che il modello di amplificazione eastern NA non è adatto all'Appennino, mentre i corrispondenti Tp possono essere usati rimanendo entro ouputs realistici
la cosa strana è che utilizzando il modello western della formula 2.35 il tp=0.05R e ad es. nel tuo caso verrebbe 0.55 che è maggiore a 0.2.
doppio "errore" nel manuale? mi sa di si perchè guardando la path duration di default è proprio 0.05 per western e 0.16 per eastern.
