Source: http://docplayer.it/30873-Protocollo-per-la-definizione-dei-valori-di-fondo-per-le-sostanze-inorganiche-nelle-acque-sotterranee.html
Timestamp: 2017-12-11 13:36:23+00:00

Document:
Protocollo per la Definizione dei Valori di Fondo per le Sostanze Inorganiche nelle Acque Sotterranee - PDF
Download "Protocollo per la Definizione dei Valori di Fondo per le Sostanze Inorganiche nelle Acque Sotterranee"
1 Servizio Interdipartimentale per le Emergenze Ambientali Settore Siti Contaminati * * * Protocollo per la Definizione dei Valori di Fondo per le Sostanze Inorganiche nelle Acque Sotterranee * * * Aprile 009
2 Elaborato da: E. Bartolucci (1), M. Bussettini (), N. Calace (1), L. D Aprile (1), M. Fratini (3), M. Guerra (1), L. Marangio (1), G. Pirani (1), A. Vecchio (1) (1) ISPRA - Servizio Interdipartimentale per le Emergenze Ambientali Settore Siti Contaminati () ISPRA - Dipartimento Acque Servizio Monitoraggio ed idrologia acque interne- Settore Idrologia (3) ISPRA - Dipartimento difesa del suolo Servizio Istruttorie, Piani di Bacino, Raccolta Dati Condiviso da: ARPA Toscana: L. Balocchi, E. Baldini, S. Cavalieri, L. Gori, S. Menichetti ISS: L. Musmeci _
3 INDICE Acronimi PREMESSA E OBIETTIVI DEL PROTOCOLLO... 5 Inquadramento Normativo Normativa Europea Normativa Nazionale Acque Bonifiche e Danno Ambientale Documenti di riferimento Progetto Bridge Risultati di altri studi condotti in Europa PROCEDURA PER LA DETERMINAZIONE DEI VALORI Di FONDO Definizione del Modello concettuale Raccolta dati Ricostruzione dell assetto geologico e idrogeologico Valutazione delle pressioni antropiche Organizzazione della banca dati Completezza del set di dati Revisione e selezione dei dati raccolti Pianificazione indagini ex novo Analisi dei dati Numerosità campionaria Trattamento dei non detect Individuazione e trattamento degli outliers Definizione della distribuzione dei dati Determinazione del valore del fondo Confronto tra le popolazioni dei dati del fondo e di un area di interesse BIBLIOGRAFIA APPENDICI: A - ANALISI STATISTICA DEI DATI
4 ACRONIMI DQA DAS MATTM VF VFN VS Direttiva Quadro Acque Direttiva Acque Sotterranee Ministero dell ambiente e della tutela del territorio e del mare Valore di fondo (Background level) Valore di fondo naturale (Natural background level) Valore soglia (Threshold value) Pagina 4 di 17
5 1 PREMESSA E OBIETTIVI DEL PROTOCOLLO Il presente protocollo è stato elaborato da ISPRA per rispondere alla richiesta pervenuta dalla Direzione Qualità della Vita del MATTM con nota Prot /QdV/DI del 7 Maggio 008, acquisita al Prot. APAT del 13 Maggio 008. In tale nota il MATTM ha richiesto al Servizio Interdipartimentale per le Emergenze Ambientali dell ex APAT (ora ISPRA), l elaborazione di una metodologia per la determinazione del fondo naturale nelle acque di falda nei Siti di Interesse Nazionale di Massa Carrara e Livorno con particolare riferimento ai seguenti parametri: - per il SIN di Massa: Mn - per il SIN di Livorno: Fe, Mn, Solfati e Cloruri. Il documento, ancorché predisposto sulla base di una richiesta specifica, trova una sua generale applicazione alla determinazione del valore di fondo di composti inorganici con particolare riferimento ai metalli e metalloidi Occorre premettere che, in considerazione delle problematiche ambientali relative ai Siti di Interesse Nazionale, caratterizzati da un elevato impatto antropico sulle matrici ambientali, non è possibile determinare veri e propri valori di fondo naturale, ovvero dovuti esclusivamente a processi naturali che caratterizzano il corpo idrico sotterraneo, non essendo distinguibili gli apporti di inquinanti connessi alle attività antropiche da quelli connessi ai processi naturali stessi. Pertanto nel presente protocollo si intende per Valore di Fondo (background level), VF, la concentrazione di una sostanza o il valore di un indicatore in un corpo idrico sotterraneo corrispondente all'assenza di alterazioni antropogeniche, o alla presenza di alterazioni estremamente limitate, rispetto a condizioni inalterate, così come indicato all Art.5 della Direttiva Acque Sotterranee 006/118/CE (DAS) e nel testo del decreto nazionale di recepimento della medesima direttiva, attualmente in fase di emanazione. I VF possono essere il risultato di vari fenomeni di origine geochimica, chimica e biologica che hanno luogo nella zona insatura e/o satura; anche la piovosità, eventuali interconnessioni tra acquiferi nonché lievi alterazioni antropogeniche possono influenzare tali valori. Pertanto ogni acquifero è caratterizzato da un chimismo unico che può presentare forti variazioni spaziali. E possibile definire range di valori tipici di ciascun acquifero per ogni analita ricercato, dal quale ricavare un valore rappresentativo da utilizzare convenzionalmente nelle procedure di bonifica dei siti contaminati. INQUADRAMENTO NORMATIVO.1 Normativa Europea La Direttiva Quadro Acque 000/60/CE (DQA) obbliga gli Stati Membri a pianificare, in un ottica di integrazione multisettoriale, l utilizzo sostenibile delle risorse idriche e del territorio alla scala di bacino, impedendone il deterioramento, al fine di raggiungere il buono stato ambientale per tutti i corpi idrici entro il 015. In particolare, l art. 17 della DQA sottolinea la necessità di misure specifiche per il raggiungimento del buono stato chimico dei corpi idrici sotterranei. Tali misure comprendono (art.17, comma ) criteri per valutare il buono stato chimico delle acque sotterranee, secondo l'allegato II, punto. e dell'allegato V, punti.3. e.4.5 e criteri per individuare tendenze significative e durature all'aumento e per la determinazione di punti di partenza da utilizzare per le inversioni di tendenza secondo l'allegato V, punto.4.4. La Direttiva 006/118/CE, avente come oggetto la tutela delle acque sotterranee dall inquinamento, attua quanto richiesto dall art.17 della DQA e stabilisce i criteri per la valutazione dello stato chimico e dei trend. In particolare, l articolo 3 definisce lo stato chimico delle acque sotterranee sulla base del confronto delle concentrazioni misurate con dei valori di conformità, o valori soglia (threshold Pagina 5 di 17
6 values) VS. L allegato II, parte A, fornisce indicazioni per la fissazione dei valori soglia da parte degli Stati Membri, sottolineando la necessità di considerare le caratteristiche idrogeologiche comprendenti informazioni sui livelli di fondo e sul bilancio idrico. Ulteriori indicazioni tecniche sono fornite dalle linee guida della Commissione Europea Guidance on Groundwater Status and Trend Assessment che, ancorché non cogenti, sono state approvate ufficialmente dai Direttori alle Acque di tutti gli Stati Membri e recepite nelle norme tecniche nazionali. La rilevanza dell analisi idrogeochimica era già presente all art..3 dell allegato II alla DQA, ove è consentito agli Stati di Membri di utilizzare tipologie di caratterizzazione delle acque sotterranee all'atto di stabilire i livelli di fondo naturale per questi corpi idrici sotterranei.. Normativa Nazionale..1 Acque La normativa nazionale vigente (D.Lgs. 15/06) e il decreto attuativo in fase di emanazione, in continuità con la precedente (D.Lgs. 15/99) ed in conformità con le direttive comunitarie 000/60 e 006/118, impone il raggiungimento, entro il 015, del buono stato ambientale dei corpi idrici sotterranei, definito come il più scadente tra lo stato quantitativo e lo stato chimico. La classificazione dello stato chimico si basa sul non superamento di valori soglia VS. Tali valori, fissati a livello nazionale su base ecotossicologica, possono essere rivisti a scala sito-specifica laddove il fondo naturale delle acque sotterranee assuma delle concentrazioni superiori ai VS. In tal caso i VF vengono assunti quali VS. I VS costituiscono anche il riferimento per l attivazione delle misure necessarie ad invertire i trend (Decreto di recepimento direttiva 006/118/CE, allegato 6, parte B, all art.5, comma1): Il punto di partenza per attuare misure atte a provocare l'inversione delle tendenze significative e durature all'aumento è stabilito quando la concentrazione di inquinanti raggiunge il 75 % dei valori parametrici degli standard di qualità o dei valori soglia delle acque sotterranee La determinazione dei VF assume pertanto una rilevanza prioritaria al fine di non classificare le acque di scarsa qualità come in cattivo stato (o di valutare improbabili punti di inversione dei trend) con conseguente attivazione di impossibili misure di ripristino. Infine, occorre ricordare che già la precedente normativa (D. Lgs. 15/99) prevedeva un apposita classe di stato, 0 o classe naturale per le acque caratterizzate da presenza naturale di sostanze con concentrazioni superiori ai VS fissati per quelle sostanze a livello nazionale... Bonifiche e Danno Ambientale La necessità di determinare valori di fondo per il suolo e le acque sotterranee ai quali riferire gli obiettivi degli interventi di bonifica e ripristino ambientale era stata stabilita già nel Decreto ministeriale 471 del 1999 (DM 471/99), laddove all Articolo 4 (Obbligo di bonifica e ripristino ambientale), comma, si riportava: Per ogni sostanza i valori di concentrazione da raggiungere con gli interventi di bonifica e ripristino ambientale sono tuttavia riferiti ai valori del fondo naturale nei casi in cui, applicando le procedure di cui all'allegato, sia dimostrato che nell'intorno non influenzato dalla contaminazione del sito i valori di concentrazione del fondo naturale per la stessa sostanza risultano superiori a quelli indicati nell'allegato 3 In particolare, per le acque sotterranee, il DM 471/99 proponeva l adozione di obiettivi di bonifica e ripristino ambientale più restrittivi in caso di aree sensibili o situazioni di particolare vulnerabilità degli acquiferi, privilegiando la necessità di tutela della risorsa per l uso potabile (Articolo 4, comma 3): I valori di concentrazione da raggiungere con la bonifica. ed il ripristino ambientale possono essere più restrittivi di quelli previsti per la specifica destinazione d'uso qualora il corpo idrico Pagina 6 di 17
7 ricettore compreso, anche parzialmente, nel sito da bonificare sia classificato come area sensibile ai sensi della normativa sulla tutela delle acque dagli inquinamenti, ovvero ricorrano situazioni accertate di particolare vulnerabilità delle acque all'inquinamento ovvero sia necessario tutelare la qualità delle acque destinate ad uso potabile. Anche il Decreto Legislativo n.15/06 Norme in materia ambientale", pubblicato nella Gazzetta Ufficiale n. 88 del 14 aprile Supplemento Ordinario n. 96, all art. 40, che con la Parte IV, Titolo V, sostituisce il DM 471/99, prevede l utilizzo dei valori di fondo, laddove riporta, nella definizione di Concentrazione Soglia di Contaminazione (CSC), (Art. 40, comma 1, lettera b): CSC:i livelli di contaminazione delle matrici ambientali che costituiscono valori al di sopra dei quali è necessaria la caratterizzazione del sito e l'analisi di rischio sito specifica, come individuati nell'allegato 5 alla parte quarta del presente decreto. Nel caso in cui il sito potenzialmente contaminato sia ubicato in un'area interessata da fenomeni antropici o naturali che abbiano determinato il superamento di una o più concentrazioni soglia di contaminazione, queste ultime si assumono pari al valore di fondo esistente per tutti i parametri superati; I valori di fondo, quindi, sono sostitutivi dei valori di riferimento per terreni e acque sotterranee, al di sopra dei quali è necessaria l elaborazione dell analisi di rischio sito-specifica. In riferimento a quest ultima occorre inoltre ricordare le modifiche apportate alla Parte IV, Titolo V del DLgs 15/06 dal Decreto Legislativo 16 gennaio 008, n. 4 "Ulteriori disposizioni correttive ed integrative del decreto legislativo 3 aprile 006, n. 15, recante norme in materia ambientale", pubblicato nella Gazzetta Ufficiale n. 4 del 9 gennaio Suppl. Ordinario n. 4. L Art. 43 del suddetto decreto riporta infatti:.il punto di conformità per le acque sotterranee rappresenta il punto a valle idrogeologico della sorgente al quale deve essere garantito il ripristino dello stato originale (ecologico, chimico e/o quantitativo) del corpo idrico sotterraneo, onde consentire tutti i suoi usi potenziali, secondo quanto previsto nella parte terza (in particolare articolo 76) e nella parte sesta del presente decreto (in particolare articolo 300). Pertanto in attuazione del principio generale di precauzione, il punto di conformità deve essere di norma fissato non oltre i confini del sito contaminato oggetto di bonifica e la relativa CSR per ciascun contaminante deve essere fissata equivalente alle CSC di cui all'allegato 5 della parte quarta del presente decreto. Valori superiori possono essere ammissibili solo in caso di fondo naturale più elevato o di modifiche allo stato originario dovute all'inquinamento diffuso, ove accertati o validati dalla Autorità pubblica competente, o in caso di specifici minori obiettivi di qualità per il corpo idrico sotterraneo o per altri corpi idrici recettori, ove stabiliti e indicati dall'autorità pubblica competente, comunque compatibilmente con l'assenza di rischio igienico-sanitario per eventuali altri recettori a valle. Sulla base di quanto disposto dal DLgs 04/08, quindi, i valori di fondo stabiliti per le acque sotterranee validati dall Autorità pubblica competente costituiscono gli obiettivi di bonifica sitospecifici da rispettare al punto di conformità. Occorre inoltre ricordare quanto disposto dal Dlgs 15/06 Parte VI, in merito al risarcimento del danno ambientale. L art. 300 riporta infatti: comma 1: E' danno ambientale qualsiasi deterioramento significativo e misurabile, diretto o indiretto, di una risorsa naturale o dell'utilità assicurata da quest'ultima comma. Ai sensi della direttiva 004/35/CE costituisce danno ambientale il deterioramento, in confronto alle condizioni originarie, provocato: a). b) alle acque interne, mediante azioni che incidano in modo significativamente negativo sullo stato ecologico, chimico e/o quantitativo oppure sul potenziale ecologico delle acque interessate, quali definiti nella direttiva 000/60/CE, ad eccezione degli effetti negativi cui si applica l'articolo 4, paragrafo 7, di tale direttiva; Pagina 7 di 17
8 c) ; d).. e ancora, all art. 30, comma 9: Per "ripristino", anche "naturale", s'intende: nel caso delle acque, delle specie e degli habitat protetti, il ritorno delle risorse naturali o dei servizi danneggiati alle condizioni originarie; La definizione di un danno alle acque sotterranee è quindi basata sull accertamento delle condizioni originarie delle stesse che non può prescindere dalla definizione del valore di fondo. Tabella I: Riferimenti Normativi e Modalità di Utilizzo dei Valori di Fondo Riferimento Normativo DM 471/99 (sostituito dal DLgs 15/06) D.Lgs 15/99 DLgs 15/06 (modificato dal DLgs 04/08) DLgs 04/08 (correttivo del DLgs 15/06) Modalità di Utilizzo dei Valori di Fondo Valori di intervento al di sopra dei quali è necessaria la bonifica (sostituitivi valori tabellari Allegato 3) Definizione dei valori soglia VS superati i quali è necessario porre in atto misure di ripristino; definizione dello stato chimico naturale per il quale non sono previsti interventi di risanamento. Parte III (Acque): decreto correttivo in fase di emanazione valori di riferimento per la fissazione dei valori soglia, con i quali coincidono quando VF>VS. Al di sopra dei VS occorre attuare misure di ripristino. Parte IV, Titolo V (Bonifiche): valori di riferimento (screening, CSC) al di sopra dei quali è necessaria l applicazione dell analisi di rischio sito specifica Parte VI (Danno Ambientale): valori di riferimento per la determinazione delle condizioni originarie della risorsa Obiettivi di bonifica sito-specifici da rispettare al Punto di Conformità 3 DOCUMENTI DI RIFERIMENTO 3.1 Progetto Bridge Il Progetto Europeo Backgoround criteria for the IDentification of Groundwater thresholds BRIDGE ha proposto una metodologia di derivazione dei VF che può essere schematizzata come segue: determinazione della tipologia dell acquifero sulla base di una tipizzazione idrogeochimica degli acquiferi europei operata tenendo in considerazione la litologia, il grado di intrusione salina, parametri idrodinamici, condizioni redox, contenuto di sostanza organica, ossidi e solfuri minerali, nonché l età geologica. identificazione dell approccio da utilizzare sulla base della tipologia, dalla quantità e dalla affidabilità dei dati disponibili. Nei casi in cui si ha un livello medio di conoscenza dei processi geochimici che interessano l acquifero e sono disponibili dati relativi al monitoraggio delle acque sotterranee i VF possono essere: Pagina 8 di 17
9 stimati sulla base dei risultati di elaborazioni statistiche (uso del 90 percentile dei valori elaborati statisticamente in altri Paesi per acquiferi della stessa tipologia); derivati sulla base di un approccio semplificato (pre-selezione). La pre-selezione è basata sull esclusione dei campioni che contengono indicatori di contaminazione antropica (ad es: nitrati, potassio, ammoniaca)al di sopra di un certo valore. Ulteriori limitazioni sull uso dei dati da utilizzare per la determinazione dei VF sono costituite dall incompletezza di informazioni circa i punti di campionamento (es. profondità rappresentativa del campione) oppure quando gli stessi afferiscono a situazioni geologiche/geochimiche specifiche (es. presenza di sorgenti idrotermali, intrusioni saline, ecc.). I VF possono essere derivati dopo la esclusione dei campioni secondo i criteri sopra accennati, attraverso semplici metodi statistici. In particolare possono essere utilizzati i valori relativi al 90 percentile della distribuzione dei dati. 3. Risultati di altri studi condotti in Europa Il Servizio Geologico Francese (BRGM) ha proposto un metodo finalizzato alla determinazione del fondo geochimico di un acquifero, in funzione del grado delle conoscenze acquisite in merito al sistema fisico in studio (Guide technique qualité naturelle des eaux souterraines: Méthode de caractérisation des états de référence des aquifères français, Laurence Chéry, BRGM édition, 006). Lo schema proposto prevede quattro livelli di conoscenza: 1. informazioni di base non geochimiche : - carta dei complessi idrogeologici; - carte idrogeologiche; - log stratigrafici;. informazioni geochimiche di base: - analisi chimiche da punti di captazione; 3. informazioni geochimiche specifiche: - dati acquisiti a seguito di studi di dettaglio; - dati a scala locale; 4. informazioni geochimiche specifiche complementari: - elementi in tracce, analisi isotopiche da acquisire mediante indagini specifiche. Ai quattro livelli di conoscenza sopra riportati corrispondono tre possibili approcci: Descrizione della facies geochimica teorica a partire dalle conoscenze geologiche di base. Si basa sullo studio della correlazione geologia-fondo geochimico in base al quale è possibile definire l ordine di grandezza della concentrazione di un determinato elemento in funzione della tipologia litologica dell acquifero. Analisi dei dati geochimici esistenti ed esclusione di quelli in cui sono presenti indicatori di contaminazione antropica. Questo approccio utilizza alcuni elementi chimici caratteristici di apporti antropici come i cloruri, i nitrati, ecc. Analisi dei dati geochimici esistenti e definizione di quelli rappresentativi del fondo geochimico naturale mediante distinzione tra componente naturale e componente antropica sulla base dei rapporti isotopici. Pagina 9 di 17
10 Livello di conoscenza dell acquifero SCARSA ELEVATA Tabella II: Schema BRGM APPROCCIO APPLICABILE LIVELLO DI CONFIDENZA ATTESO SEMPLICE = COMPLESSO SCARSO = ELEVATO Descrizione della facies geochimica teorica Analisi dei dati esistenti ed esclusione di quelli antropizzati Ricostruzione del fondo geochimico mediante correzione dell apporto antropico Base non geochimica qualitativo semiquantitativo quantitativo Base geochmica Geochimica specifica Geochimica specifica da acquisire 4 PROCEDURA PER LA DETERMINAZIONE DEI VALORI DI FONDO Nella Figura seguente è schematizzata la procedura per la determinazione dei valori di fondo per le acque sotterranee. Nelle sezioni successive sono descritti, nel dettaglio, i singoli passaggi. DEFINIZIONE DEL MODELLO CONCETTUALE Raccolta dati Acquisizione nuovi dati ORGANIZZAZIONE DEL DATA BASE REVISIONE DEI DATI Definizione dei criteri di esclusione NO set di dati sufficiente? SI ANALISI DEI DATI - non detect - outliers - distribuzione DETERMINAZIONE DEL VALORE DI FONDO CONFRONTO POPOLAZIONI FONDO / AREA DI INTERESSE Figura 1. Schema della procedura per la determinazione dei valori di fondo Pagina 10 di 17
11 4.1 Definizione del Modello concettuale In accordo a quanto riportato nell allegato 1 dello Schema di Decreto legislativo di recepimento della Direttiva 006/118/CE del Parlamento Europeo e del Consiglio, del 1 dicembre 006, sulla protezione delle acque sotterranee dall inquinamento e dal deterioramento, il modello concettuale rappresenta il sistema delle acque sotterranee sulla base delle conoscenze delle caratteristiche naturali (tipo di acquifero, struttura tridimensionale, condizioni idrauliche ed al contorno) e delle pressioni e degli impatti. La definizione del modello concettuale, ottenuta combinando le informazioni relative all assetto geologico/idrogeologico e alla valutazione delle pressioni antropiche, fornisce gli elementi necessari all identificazione dei punti di indagine (pozzi esistenti, pozzi di monitoraggio, ecc.) idonei per lo studio delle concentrazioni di fondo Raccolta dati La prima fase della ricostruzione del modello concettuale, consiste nella raccolta delle informazioni disponibili riguardanti sia gli aspetti geologici/idrogeologici sia la presenza di potenziali fonti di contaminazione. Per questo scopo saranno analizzati anche gli studi e le ricerche condotte da Istituti pubblici, Università, soggetti privati che contengano informazioni sullo stato ambientale dell area come, ad esempio: studi riguardanti l assetto geologico-strutturale e le conoscenze idrogeologiche ed idrochimiche; indagini geologiche - geognostiche (carotaggi, prove penetrometriche, ecc.); censimento pozzi; dati idrogeologici (misure piezometriche, prove idrogeologiche in pozzo, analisi idrochimiche, ecc.); dati meteoclimatici; analisi chimiche dei suoli e delle falde; carte dell uso del suolo; dati riguardanti le attività industriali. Tutti i dati raccolti verranno inseriti ed organizzati in un data-base georeferenziato per le successive elaborazioni Ricostruzione dell assetto geologico e idrogeologico Ai fini della ricostruzione del modello geologico e idrogeologico del sito dovranno essere presi in considerazione i seguenti aspetti: ricostruzione dell assetto stratigrafico, anche mediante sezioni geologiche interpretative; identificazione degli acquiferi; delimitazione spaziale degli acquiferi - spessore, profondità della superficie piezometrica, eventuali variazioni connesse alle maree; direzioni di deflusso delle acque sotterranee e gradiente idraulico; definizione dei parametri idrogeologici trasmissività, permeabilità, ; valutazione di eventuali flussi tra i differenti acquiferi; valutazione dello stato qualitativo delle acque sotterranee; presenza e caratteristiche di acque superficiali (fiumi, laghi ) e interazioni con l acquifero. Pagina 11 di 17
12 La ricostruzione potrà essere presentata mediante relazioni ed elaborazioni grafiche, quali ad esempio: ubicazione punti d acqua censiti; carta delle curve isopiezometriche; sezioni idrogeologiche esplicative; carta delle facies idrochimiche e della qualità delle acque. I criteri definiti nel presente documento non si applicano agli acquiferi non significativi, come, ad esempio, falde superficiali di estensione e spessore limitati, soggette a importanti variazioni stagionali, potenzialmente prossime a sorgenti dirette di contaminazione. Per definire la significatività di un acquifero può essere utilizzato il criterio proposto nell allegato 1 dello schema di decreto legislativo di recepimento della Direttiva 006/118/CE, sulla protezione delle acque sotterranee dall inquinamento e dal deterioramento. Secondo quanto riportato nella norma citata, l identificazione degli acquiferi è effettuata sulla base di due criteri: flusso significativo e quantità significativa. Se uno o entrambi i criteri sono soddisfatti, gli strati geologici sono da considerarsi acquifero, come illustrato nello schema seguente. Figura. Schema per l identificazione degli acquiferi Valutazione delle pressioni antropiche Per la valutazione delle pressioni che possono causare una variazione dello stato qualitativo dell acqua si procederà alla verifica della presenza e della ubicazione delle aree maggiormente antropizzate. In questa fase saranno valutate: ubicazione di tutte le sorgenti attive o potenziali, presenti all interno dell area in studio ovvero all esterno ma che impattino sull area stessa; tipologia ed estensione delle sorgenti di contaminazione; profondità delle sorgenti di contaminazione rispetto alla superficie topografica e alla superficie piezometrica; Pagina 1 di 17
13 i parametri chimici e fisici (es. ph, Eh, T) che influenzano la diffusione ed il trasporto delle specie chimiche di interesse. 4. Organizzazione della banca dati Nei successivi paragrafi sono descritti i criteri da seguire per la composizione del set di dati da utilizzare nella determinazione dei valori di fondo. I requisiti di rappresentatività, omogeneità e qualità proposti, potranno essere utilizzati sia per la validazione di dati di bibliografia, acquisiti in seguito ad indagini già effettuate, sia per la programmazione di campagne di indagine appositamente predisposte. Al fine di verificarne la rappresentatività i dati dovranno essere soggetti ad una attenta revisione secondo quanto prospettato nelle sezioni successive. E necessario sottolineare tuttavia che in funzione dell importanza e/o complessità del corpo idrico sotterraneo di studio, dello stato dell arte delle conoscenze ad esso relative e degli obiettivi dello studio, la completezza del set dei dati disponibili dovrà essere valutata con un certo grado di tolleranza Completezza del set di dati La complessità di un sistema idrogeologico è tale che maggiore è il numero di informazioni disponibili, maggiore è la possibilità di identificare i processi idraulici in senso stretto (es. ricarica e deflusso), i processi idrochimici (es. interazione acqua/roccia, processi di mixing) nonché eventuali fenomeni di contaminazione diffusa e/o puntuale che caratterizzano il sistema in esame. Per tali motivi la banca dati dovrà essere organizzata in modo tale da fornire le seguenti informazioni: coordinate, tipologia punto monitoraggio (pozzo, piezometro, sorgente), ricostruzione stratigrafica, profondità del tratto filtrato; profondità di campionamento, data di campionamento, metodo di campionamento, pretrattamento del campione (ovvero se le analisi sono state eseguite sul campione tal quale o sul filtrato a 0,45 µm, aggiunta di reagenti e/o di preservanti), conservazione e metodi di analisi per ciascun parametro; dati chimico-fisici: temperatura, conducibilità elettrica, ph, Eh, ossigeno disciolto, solidi sospesi totali. Tali parametri permettono di valutare gli equilibri chimici a cui possono partecipare gli analiti in esame e le loro eventuali variazioni caratterizzanti il corpo idrico (es. mobilizzazione dei metalli pesanti): il parametro solidi sospesi totali, per esempio, può evidenziare, laddove non fosse esplicitamente riportato, le condizioni di campionamento (es spurgo del pozzo); chimismo principale: contenuto in K +, Na +, Mg +, Ca +, Cl -, SO4 -, HCO 3-, NO 3-. Tali dati consentono di attribuire al campione una facies idrochimica risultante dai processi che avvengono nel sistema in studio (rappresentazione tramite diagramma di Piper o simili); composti organici: (es. idrocarburi totali) nella maggior parte dei casi essi possono essere considerati indicatori di contaminazione dovuta a pressione antropica; parametri di interesse: sono i composti inorganici per i quali si vogliono determinare i valori di fondo nelle acque sotterranee. 4.. Revisione e selezione dei dati raccolti Il primo livello di revisione e selezione dei dati raccolti consiste nell individuazione di quelli pertinenti al sistema di interesse, con particolare riferimento anche alle profondità e alla litologia in corrispondenza delle quali è stata prelevata l acqua. Tali dati dovranno pertanto essere relativi a punti di campionamento (pozzi, piezometri, sorgenti) che presentino le seguenti caratteristiche: interessino lo stesso sistema di quello oggetto di indagine; Pagina 13 di 17
14 siano ubicati in aree distanti da sorgenti di contaminazione puntuali e/o specifiche attive nel presente o nel passato. In linea teorica i punti di campionamento più idonei per la acquisizione di dati analitici, saranno quelli ubicati in aree sopra gradiente rispetto a potenziali sorgenti di contaminazione. Occorre tener presente che potrà essere valutata l opportunità di utilizzare anche i dati provenienti da punti di campionamento ubicati sottogradiente rispetto a potenziali sorgenti, nel caso in cui questi fossero ubicati all esterno della zona interessata dalla dispersione laterale dei contaminanti di interesse. A tal fine si suggerisce di rappresentare graficamente (es. post map) la distribuzione spaziale dei parametri indagati correlandole con lo scenario idrogeologico ed antropico. Per la verifica della estensione della zona contaminata si potrà far ricorso a modelli numerici o programmi di calcolo specifici. La valutazione dell opportunità di escludere i dati, sia perché non ritenuti rappresentativi del sistema in esame, sia perché potenzialmente affetti da significativi input antropogenici, sarà basata sui seguenti criteri: i campioni che mostrano concentrazioni di contaminanti organici maggiori di 3 volte il limite di rilevabilità (assimilabile sperimentalmente al limite di quantificazione); il limite di rilevabilità dovrà comunque essere accettabile rispetto ai metodi analitici attualmente diffusi; i campioni prelevati in pozzi/piezometri di cui non è nota la profondità del tratto filtrato (campioni per i quali non è possibile, con i dati a disposizione, verificare l appartenenza all acquifero in esame); i campioni con concentrazioni di NO 3 > 10 mg/l; i campioni con concentrazioni di NaCl > 1000 mg/l. Tale vincolo deve tuttavia essere valutato caso per caso, in funzione delle interazioni, anche naturali, che il corpo idrico di interesse ha con sistemi esterni (es. prossimità alla linea di costa, risalita di fluidi di origine profonda). In generale, le difformità dai requisiti sopra enunciati dovranno essere opportunamente documentate in modo da poter effettuare, in fase di elaborazione dei dati, una valutazione ad hoc, anche mediante l applicazione di opportuni strumenti statistici. Qualora i dati ritenuti idonei alla valutazione del fondo sulla base dei criteri descritti in precedenza risultino non coerenti da un punto di vista temporale, si potrà tener conto dei seguenti criteri: se i dati non presentano significative correlazioni rispetto alla variabile tempo, ovvero non si registrano dei trend, si possono trattare come se fossero temporalmente coerenti quindi includere in un dataset unico. se i dati presentano un trend temporale non è possibile trattarli, in linea di massima, come dataset unico. In tal caso la valutazione del fondo potrà tener conto di tali trend o, qualora necessario, si valuterà la possibilità di reperire nuovi dati (es. nuova campagna di monitoraggio). Nel caso in cui i composti inorganici di interesse mostrino evidenti correlazioni con le caratteristiche chimico fisiche e che queste ultime mostrino una certa variabilità all interno del set di dati (es. il parametro Eh che individua porzioni della falda a diverso stato di ossidoriduzione) si consiglia di individuare dei subset in cui dette caratteristiche siano sostanzialmente omogenee. Nelle situazioni più complesse, e se il set di dati ha una certa consistenza, questo orientamento può indirizzarsi verso una cluster analysis più articolata. Pagina 14 di 17
15 4..3 Pianificazione indagini ex novo Nel caso in cui i dati raccolti non fossero idonei per la determinazione del fondo (numerosità non sufficiente, eccessiva variabilità temporale, ) lo studio potrà prevedere la necessità di effettuare ulteriori campagne di monitoraggio utilizzando punti di misura esistenti e riconosciuti idonei per la determinazione dei valori di fondo. In funzione dell importanza e/o complessità del sistema in esame nonché degli obiettivi dello studio, si potrà valutare l opportunità di effettuare più campagne di monitoraggio per individuare eventuali variazioni stagionali. Nel caso in cui non fossero disponibili pozzi di monitoraggio ubicati in zone idonee e con caratteristiche adeguate alla caratterizzazione del sistema in studio, lo studio potrà prevedere la necessità di realizzarne appositi. Le caratteristiche costruttive dei pozzi, le modalità di prelievo dei campioni e le metodiche analitiche dovranno essere tali da garantire la rappresentatività del dato. Per rispettare i criteri di rappresentatività e riproducibilità dei dati prodotti le analisi chimiche devono essere effettuate sul campione filtrato (possibilmente in campo). Qualora le indagini ex novo vadano ad integrare un set di dati sui quali non è stata eseguita la filtrazione del campione, si potrà valutare di eseguire le analisi sia sul filtrato sia sul tal quale al fine di meglio interpretare i dati pregressi e di verificare l opportunità di impiegarli per la determinazione dei VF. Le fasi di campionamento devono essere accompagnate dalla raccolta di campioni di bianco (bianco di campo etc.) al fine di valutare l accuratezza dei dati chimici prodotti (APAT, 006). Ciascun campione prelevato dovrà essere etichettato, e dovrà essere compilata una scheda del campione in cui siano riportate le seguenti informazioni: coordinate geografiche, profondità e spessore di acquifero campionato, data di campionamento e tecniche di pretrattamento del campione se operate in campo (come ad esempio campione filtrato, aggiunta di reattivi, ecc), conservazione del campione (contenitori utilizzati, temperatura di conservazione, ecc.). Dal punto di vista analitico i parametri da determinare sono (Hem, 1985): temperatura, potenziale redox, conducibilità elettrica, ph, ossigeno disciolto da effettuarsi in campo; solidi sospesi totali; ioni maggiori: sodio, potassio, calcio, magnesio, cloruri, solfati, nitrati, bicarbonati; contaminanti organici, scelti in funzione delle problematiche sito-specifiche analiti di interesse (ad esempio metalli pesanti). In funzione della complessità sito specifica potranno essere acquisite ulteriori informazioni geochimiche (es:elementi in tracce, analisi isotopiche). Per tutti gli analiti da determinare dovranno essere utilizzati metodi di analisi ufficiali riconosciuti a livello nazionale e/o internazionale e che tali metodi facciano riferimento alle più avanzate tecniche di impiego generale. E auspicabile inoltre che i dati prodotti siano validati da un Ente di Controllo nell eventualità in cui lo studio sui valori di fondo sia effettuato da Aziende private. 4.3 Analisi dei dati Numerosità campionaria L indicazione del numero minimo di record su cui basare la procedura per la determinazione del fondo delle acque sotterranee dipende dalle condizioni al contorno quali ad es. l estensione del corpo idrico, dal tipo di distribuzione dei valori e dal livello di accettabilità dell errore definito a priori dal decisore. In linea di massima il numero di stazioni di campionamento necessario a garantire la significatività statistica del valore determinato è compreso fra 10 e 30 (APAT-ISS, 006), in funzione dello scenario sito specifico (es. rilevanza del sistema idrico in esame, difficoltà Pagina 15 di 17
16 di individuare punti di campionamento idonei, ecc) Come ovvio un numero di campioni maggiore aumenta la significatività del valore determinato, a patto che siano comunque rispettati i criteri descritti nel paragrafi precedenti 4.3. Trattamento dei non detect Le concentrazioni di alcuni parametri possono risultare inferiori al limite di rilevabilità (detection limit, d.l.) del metodo analitico con il quale sono stati analizzati. I metodi con cui si trattano questi non detect (n.d.) sono funzione del tipo di distribuzione che caratterizza la popolazione, della percentuale con cui essi si registrano nel data set, e della variabilità dei limiti di rilevabilità associati ai dati. In prima istanza si ritiene opportuno associare ai non detect un valore di concentrazione pari al corrispondente limite di rilevabilità (n.d. = d.l.); laddove siano presenti differenti limiti di rilevabilità per un singolo parametro si dovrà associare ai non detect il limite di rilevabilità più basso Individuazione e trattamento degli outliers Il set di dati ottenuto dovrà essere analizzato al fine di individuare la presenza di eventuali outliers, distinguendo in particolare i veri outlier dai falsi outlier. I veri outlier possono derivare da errori di trascrizione, di codifica dei dati o da una qualsiasi inefficienza degli strumenti del sistema di rilevazione dei dati. I falsi outlier sono valori estremi reali anche definibili come hot spot. E dunque necessario identificare e differenziare i tipi di outlier, in modo da rimuovere i primi e gestire i secondi. Per ulteriori dettagli si rimanda alla sezione 4 della Appendice A Definizione della distribuzione dei dati Lo scopo è quello di individuare la distribuzione di probabilità che approssimi meglio l insieme dei dati disponibili. L individuazione del tipo di distribuzione che meglio approssima il campione di dati serve a definire i descrittori statistici più appropriati per stimare il valore del fondo. Dal tipo di distribuzione dipendono inoltre i test statistici da applicare per il confronto tra due set di dati (ad esempio il set relativo ai valori di fondo e quello relativo alle concentrazioni riscontrate in uno specifico sito). Le caratteristiche delle distribuzioni suddette e i test da applicare per la selezione delle stesse sono descritti nel dettaglio in Appendice A. 4.4 Determinazione del valore del fondo Una rassegna dei descrittori statistici, proposti a livello nazionale e internazionale, per rappresentare il VF è contenuta nel documento Protocollo operativo per la determinazione dei valori di fondo di metalli /metalloidi nei suoli dei siti di interesse nazionale (APAT-ISS, 006). Tale rassegna ha evidenziato l assenza di un criterio condiviso sulle modalità di determinazione della concentrazione del fondo in un area. Con lo scopo di fornire una procedura basata su presupposti scientifici sufficientemente rigorosi, che permetta al momento stesso un applicazione ai siti di interesse nazionale, si suggerisce, a valle delle procedure descritte nelle sezioni precedenti, il seguente approccio: verifica dell adeguatezza del numero di dati utili disponibili (minimo 30); applicazione di un test statistico per la verifica del tipo di distribuzione; calcolo e presentazione dei descrittori statistici; costruzione della curva cumulativa di frequenza e individuazione di eventuali punti di discontinuità; Pagina 16 di 17
17 il valore di fondo potrà essere individuato, fra gli indicatori statistici, anche in funzione del tipo di distribuzione (es. media +σ, mediana, percentili) ovvero tramite l individuazione di punti di discontinuità della curva cumulativa di frequenza. In analogia con il Protocollo operativo per la determinazione dei valori di fondo di metalli /metalloidi nei suoli dei siti di interesse nazionale (APAT-ISS, 006) si propone, nel caso di valutazione basata sui percentili, di riferirsi al Confronto tra le popolazioni dei dati del fondo e di un area di interesse In modo complementare rispetto all identificazione del VF, il processo decisionale potrà basarsi anche sul metodo comparativo. Tale metodo prevede il confronto tra la distribuzione dei dati rappresentativi del fondo con quella dei dati sito specifici, con lo scopo di verificare se le due popolazioni da cui sono originati i campioni siano statisticamente uguali o una mostri dei descrittori significativamente diversi (più alti o più bassi) dell altra. Il criterio comparativo può basarsi su metodi statistici più o meno complessi e raffinati di cui, in Appendice A, si riporta una descrizione dettagliata. L applicazione di test statistici si concretizza con l accettazione o il rifiuto di ipotesi statistiche circa l uguaglianza o meno dei set di campioni in esame. Con questi test si ha anche il controllo sulla probabilità di commettere errori decisionali. Il metodo comparativo, nelle diverse forme, richiede un certo grado di conoscenza del sito e di un consistente set di dati rappresentativi del fondo ovvero dell area di riferimento (non interessata da contaminazione). 5 BIBLIOGRAFIA APAT: Manuale per le indagini ambientali nei siti contaminati, Manuali e Linee guida 43/006, APAT-ISS: Protocollo operativo per la determinazione dei valori di fondo di metalli/metalloidi nei suoli dei siti di interesse nazionale, 006, Muller D., Blum A., Hart A., Hookey J., Kunkel R., Scheidleder A., Tomlin C., Wendland F., 006 Final Proposal for methodology to setup groundwater treshold values in Europe, Deliverable D18, BRIDGE project, 63 p., Chéry L.: Qualité naturelle des eaux souterraines. Méthode de caractérisation des états de référence des aquifères français, Guide technique. BRGM édition, 006, p. 38 European Commission, WGC, Guidance on Groundwater Status and Trend Assessment, Office for Official Publications of the European Communities, 008 in pubblicazione, Hem J.D., 1985, Study and interpretation of the chemical characteristics of natural water (3d ed.): U.S. Geological Survey Water-Supply Paper 54 Pagina 17 di 17
18 APPENDICE A ANALISI STATISTICA
19 SOMMARIO 1 INTRODUZIONE... 1 TEST STATISTICI NUMEROSITÀ CAMPIONARIA... 4 GLI OUTLIER I NON-DETECT ANALISI MULTIVARIATA DISTRIBUZIONE DEI DATI DESCRITTORI STATISTICI DEFINIZIONE DEL TIPO DI DISTRIBUZIONE CONFRONTO FRA I VALORI DI FONDO E I VALORI SITO SPECIFICI Indice delle Tabelle Tabella 1: Confidenza e Potenza del Test... Tabella : Criteri di selezione del test per la identificazione degli outlier... 4 Tabella 4: Test per la selezione del tipo di distribuzione Tabella 6: Tipologia dei Test Statistici Tabella 8: Applicabilità dei Test Statistici Tabella 7 Esempio di calcolo del Rango Tabella 1 Valori critici per la distribuzione t di Student;... 0 Tabella 10 Valori critici per il Wilcoxon Rank Sum Test... 1 Indice delle Figure Figura 1 Esempio di distribuzione normale... 8 Figura 3 Esempio di distribuzione lognormale Figura 5: Parametri statistici rappresentati nel box plot... 1 Figura 7 Curva cumulativa di frequenza
20 1 INTRODUZIONE La presente Appendice contiene indicazioni e riferimenti bibliografici per l applicazione dei criteri finalizzati alla determinazione dei valori di fondo di metalli e metalloidi nei terreni dei siti d interesse nazionale. Nel seguito sono descritti i principali elementi che è necessario prendere in considerazione per stabilire l applicabilità di criteri statistici atti ad individuare la distribuzione dei valori di fondo dall insieme di dati a disposizione. Si presuppone che i dati analitici a disposizione siano stati già validati, ossia sia stata verificata la loro attendibilità. Nella redazione della presente appendice si è fatto riferimento ai contenuti dell APPENDICE H al documento Criteri metodologici per l'applicazione dell'analisi assoluta di rischio ai siti contaminati rev, disponibile sul sito dell APAT, (ISPRA, marzo 008). TEST STATISTICI Il test statistico può essere visto come un mezzo per verificare in maniera quantitativa la validità di un ipotesi. In statistica, l ipotesi da verificare si chiama ipotesi nulla e si indica con H 0, mentre con H a si indica l ipotesi alternativa. Come esempio, nel caso in cui fosse necessario verificare se un set di dati ha una distribuzione normale, è possibile prevedere le seguenti ipotesi: H 0 : il set di dati presenta una distribuzione normale H a : il set di dati non presenta una distribuzione normale Se, applicando il test, risultasse scartata l ipotesi H 0 in favore di H a, si potrebbe concludere che il set di dati, non supportando l ipotesi nulla, deve derivare da un tipo di distribuzione non normale. Occorre però tener presente che l applicazione di un test statistico comporta sempre un rischio di errore. Nella pratica statistica si individuano due tipi di errori errore di primo tipo; è quello che porta a rifiutare H 0 quando è vera; errore di secondo tipo; è quello che porta ad accettare H 0 quando è falsa. Con la lettera α si indica la probabilità di commettere un errore di primo tipo. Con 100(1 α)% si indica il livello di confidenza del test. Se il test non scarta H 0, ovvero conferma l ipotesi nulla, può significare che le informazioni del set di dati non sono sufficienti per scartare H 0 con quel livello di confidenza. Ad esempio, fissato α = 0,05(5%), con i dati campionari si esegue il test e si valuta se il suo valore cade nella regione di rifiuto e nella regione di accettazione. Se, ad esempio, cade nella regione di rifiuto si dice che il test è significativo al 5%. Con la lettera β i indica la probabilità di commettere un errore del secondo tipo. Con 100(1 β)% si indica la potenza del test, cioè la probabilità di scartare correttamente l ipotesi nulla. 1
21 DECISIONE BASATA SU CONDIZIONE DEL SITO ATTUALE UN CAMPIONE DI DATI H 0 è vero H 0 non è vero H 0 non è rigettata H 0 è rigettata Decisione corretta: (1-α) Errore di tipo I: Falso positivo (α) Errore di tipo II: Falso negativo (β) Decisione corretta: (1-β) Tabella 1: Confidenza e Potenza del Test Nell ambito di indagini ambientali che comportino l applicazione di test statistici, i limiti di tolleranza della probabilità di commettere errori (del primo o del secondo tipo) dovrebbero essere specificati in fase di progettazione. Sono stati elaborati differenti metodi statistici finalizzati alla verifica di ipotesi, la scelta del più appropriato è funzione di una serie di fattori tra i quali il tipo di distribuzione dei dati (come definito al capitolo 7), ha un peso determinante. In funzione del tipo di distribuzione dei dati è possibile distinguere i Metodi Parametrici e i Metodi Non-parametrici: Parametrici: si tratta di metodi statistici che si basano su distribuzioni probabilistiche quale, ad esempio, la distribuzione normale. I test statistici parametrici sono utilizzati per la valutazione di ipotesi che riguardano i parametri della distribuzione; Non-parametrici: metodi la cui applicazione prescinde dalla conoscenza del tipo di distribuzione della popolazione. In generale i test non parametrici dovrebbero essere preferiti quando i dati non si distribuiscono secondo una distribuzione normale, o comunque non si è in grado di dimostrarlo, ad esempio per numerosità ridotta. Per riepilogare i concetti sopra riportati, nel seguito sono descritti i passaggi da seguire nella applicazione di test statistici: definizione dell ipotesi nulla e dell ipotesi alternativa; scelta del test da adottare; decisione del livello di significatività; esecuzione dei calcoli previsti nel test; decisione se accettare o meno la validità dell ipotesi nulla, in genere confrontando il valore ottenuto nel test con un valore tabulato. 3 NUMEROSITÀ CAMPIONARIA Per ogni data-set, il numero di dati a disposizione non può essere inferiore ad un valore minimo. L ampiezza del data set è di particolare importanza soprattutto nei casi in cui si abbia una grande variabilità della distribuzione dei dati. Affinché l analisi statistica sia significativa, ovvero il campione di dati sia rappresentativo della popolazione, si fa generalmente riferimento ad un numero minimo di dati, che i diversi testi consultati, riportano variabile tra 10 e 30. La scelta del numero di campioni rappresentativi è funzione, in primo luogo, dello scopo dell indagine: posto l obiettivo dell indagine di caratterizzazione, viene formulata una ipotesi da verificare mediante l applicazione di un determinato test statistico; lo stesso test serve ad effettuare la stima del numero di campioni necessari.
22 Ad esempio, si vuole verificare (ipotesi nulla) se la media delle concentrazioni di un composto in un area è superiore ad un valore soglia (ad esempio la concentrazione massima ammissibile CMA). Una delle possibili equazioni utilizzate per determinare il numero minimo di campioni necessari per la verifica della media di una distribuzione nei confronti di un valore soglia di intervento (US EPA, 006), mediante l applicazione del t-test, è: dove: n è il numero minimo di campioni ( Z1 α + Z1 β ) s n = + Z 0,5 1α s è la stima della varianza totale vera (σ ); nel caso in cui (molto frequente) non fosse nota la σ della distribuzione, dovrà essere inserita nella formula una stima di detto parametro (ad esempio il valore della deviazione standard derivata da un campione noto della popolazione di riferimento desunto da dati bibliografici o da precedenti studi). α è la probabilità accettabile che il test, applicato sul numero n di dati, indichi in maniera errata che la media delle concentrazioni non supera la CMA (in poche parole che un sito contaminato venga definito pulito ) β è la probabilità accettabile che il test, applicato sul numero n di dati, indichi in maniera errata che la media delle concentrazioni supera la CMA (in poche parole che un sito pulito venga definito contaminato) definito come la minima differenza rilevabile, ovvero, se l obiettivo dello studio è quello di confrontare la media di concentrazioni di un area con le CMA, rappresenta la massima differenza tra la media delle concentrazioni e le CMA, che è importante rilevare con una probabilità paria a 1- β. Z è il valore, per una distribuzione di dati normale, per il quale la proporzione della distribuzione a sinistra di Z 1- α è paria a 1- α. I valori di Z 1- α sono riportati in numerosi testi di statistica (es. Gilbert 1987, Table A1, pag. 54). Ad esempio se l ipotesi nulla è quella che le concentrazioni misurate superano le CMA, Z 1- α rappresenta una quantificazione del livello di accettabilità di evitare di considerare contaminato un sito che in realtà è pulito. Le assunzioni alla base della equazione sopra riportata sono che la distribuzione dei dati sia di tipo normale, i dati siano rappresentativi della popolazione, che i dati non siano correlati nel tempo e nello spazio. Il livello accettabile di errore viene definito dal decisore e viene espresso tramite il livello di confidenza (1-α) e potenza (1-β). Per gli scopi del presente documento, il numero minimo di campioni necessari per la determinazione della distribuzione di concentrazioni di fondo è posto pari a GLI OUTLIER Gli outlier sono quei valori di un data set che non sono rappresentativi dell insieme di dati nel suo complesso. Non sono rappresentativi perché, in genere, sono quantitativamente in numero estremamente ridotto e qualitativamente assumono dei valori molto grandi o molto piccoli 3
23 rispetto al resto del data set. In campo ambientale di inquinamento dei suoli, valori di concentrazione molto alti in genere corrispondono ai picchi (hot spot) locali di concentrazione. Comunque, in generale, tali valori estremi possono costituire dei veri outlier o dei falsi outlier. I primi possono derivare da errori di trascrizione, di codifica dei dati o da una qualsiasi inefficienza degli strumenti del sistema di rilevazione dei dati. I secondi sono quei valori estremi reali, spesso presenti in questo tipo di indagini soprattutto, come gia detto, in campo ambientale. La rimozione dei secondi e/o la mancata rimozione dei primi può condurre ad una visione errata del data set (EPA 000b, QA/G-9). Infatti è di fondamentale importanza tener conto e quindi non rimuovere i falsi outlier dal data set (OSWER , EPA 00). Se il data-set a disposizione è stato già validato si esclude automaticamente la presenza di veri outlier. L identificazione degli outlier può essere condotta attraverso le seguenti fasi (EPA 000b, QA/G-9). 1. Identificazione dei valore estremi che potranno essere potenziali outlier. Questo può essere fatto mediante rappresentazione grafica dell insieme dei valori rilevati: è possibile così individuare velocemente quei punti che corrispondono a valori più elevati o più ridotti rispetto agli altri. Una volta identificati i potenziali outlier, è necessario procedere a ulteriori indagini, applicando uno dei test statistici disponibili.. Applicazione di un opportuno test statistico. Esistono molti test statistici atti a verificare se un outlier statistico, cioè un potenziale vero outlier, sia tale o meno. I principali test statistici utili a tale scopo sono quattro: - Extreme value test (Dixon s Test) - Discordance Test - Rosner s test - Walsh s test Questi, descritti dettagliatamente nel seguito, si differenziano per le dimensioni del data set da considerare, il numero di potenziali outlier da analizzare e la necessità o meno di una distribuzione di tipo normale dei dati raccolti. In particolare la guida raccomanda l uso del Rosner s test quando il data set contiene un numero di elementi maggiore di 5; in caso contrario suggerisce quello dell Extreme Value test. Se si ha un solo valore sospetto outlier il Discordance test può essere sostituito a uno di questi test. Se però i dati non seguono una distribuzione normale si deve considerare un test non parametrico, come il Walsh s test (Tabella ). Per la descrizione di dettaglio dei test si rimanda al documento (EPA 000b, QA/G-9) DIMENSIONE DEL DATA SET TEST DISTRIBUZIONE NORMALE 5 n Estreme Value Test si 50 n Discordance Test si 5 n Rosner s Test si 50 n Walsh s Test no Tabella : Criteri di selezione del test per la identificazione degli outlier 4
24 3. Studio scientifico degli outlier identificati per la scelta di trattamento del dato. I test statistici (fase ) da soli non permettono di stabilire se comprendere o escludere il dato dall insieme considerato. Le scelte possibili sono: - Effettuare ulteriori approfondimenti e indagini al fine di correggere il valore di outlier. - Utilizzare il data set comprensivo dei valori di outlier. - Escludere l inserimento di tali valori dal data set. Tale scelta può avvenire solo se è possibile accompagnare i risultati dei test statistici (fase ) con valide giustificazioni scientifiche. 4. Nel caso di esclusione degli outlier dal data set, conduzione della successiva analisi statistica dei dati sia sull insieme dei dati comprensivo di outlier, sia su quello rivisto con l eventuale soppressione degli outlier. 5. Documentazione dell intero procedimento, con la descrizione di tutti i passaggi e le scelte effettuate. 5 I NON-DETECT Tutte le tecniche analitiche di laboratorio hanno un Detection Limit (DL) (limite di rilevazione): i valori cosiddetti non-detect (ND) sono quelle concentrazioni realmente o virtualmente pari a zero, o comunque maggiori di zero, ma al di sotto delle possibilità di misurazione del laboratorio. Il DL dipende dalla sensibilità della metodica di estrazione ed analisi. Un data set contenente non-detect viene definito in letteratura censored a indicare la sua incompletezza, che può essere più o meno influente a seconda del DL del laboratorio che ha condotto il campionamento: per questo motivo è opportuno che il laboratorio alleghi, alla documentazione dello studio, le informazioni sul Quantitation Limit (limite di misura) che dipenderà dalla strumentazione di cui si è servito. Il Quantitation Limit può essere definito come il livello più basso al quale una sostanza chimica può essere misurata con precisione, generalmente pari al DL dello strumento moltiplicato per un fattore compreso fra tre e cinque, ma comunque variabile a seconda della sostanza considerata e del tipo di campione (RAGS Part A, EPA 1989). La presenza di ND in un insieme di dati può influire pesantemente sul calcolo della media, della varianza, sull orientamento dei dati e su vari altri parametri, pregiudicando quindi il procedimento statistico nel caso in cui questo risulti applicabile nonostante la loro presenza. I laboratori di analisi riportano questi valori come non-detect (ND), oppure li pongono pari a zero o come dati less-than (LT) cioè minori di una certa quantità, in genere pari proprio al DL, o ancora capita di trovarli indicati come valori pari ad una frazione del DL (in genere a ½ DL). E comunque preferibile, qualora le tecniche di misurazione lo consentano, che siano riportate le loro misure esatte, benché minime, per non perdere informazioni utili all analisi dei dati. Nel seguito è riportato quanto proposto dai testi bibliografici presi quali riferimento. Il documento (OSWER , EPA 00) descrive quattro possibili approcci per la trattazione dei non-detect, finalizzati all applicazione di analisi statistiche dell insieme dei dati e alla conseguente individuazione di un valore rappresentativo. 5
25 1. Riesame del modello concettuale del sito: da questo riesame potrebbe risultare una distribuzione dei valori di concentrazione tali da permettere l individuazione di aree a maggior grado di contaminazione e aree a minor grado di contaminazione. In tal caso, il sito oggetto di indagine potrebbe essere suddiviso in sotto-aree, alcune delle quali presenteranno una maggiore e altre una minore concentrazione di non-detect. In tale caso potrebbe risultare necessario raccogliere un maggior numero di campioni per permettere una migliore caratterizzazione del sito.. Metodo della sostituzione semplice ( Simple Substitution Methods ): questo metodo prevede l assegnazione di un valore costante ai dati non-detect. Tale valore potrà essere: - pari a zero; - pari al Detection Limit; - pari alla metà del DL. L incertezza associata a questo metodo aumenta all aumentare del valore del DL e all aumentare del numero di non-detect. Quindi si consiglia di scegliere, senza un preciso criterio, il valore costante da attribuire tra i tre proposti solo nel caso in cui il numero dei non-detect costituisce al massimo il 15% di tutto il data set (EPA 000b, QA/G-9). 3. Metodo degli estremi ( Bounding Methods ): Tale metodo propone di calcolare il valore di concentrazione rappresentativo della distribuzione attribuendo, di volta in volta, uno dei valori costanti elencati sopra. Questi metodi forniscono una stima del limite superiore e di quello inferiore, calcolati sulla base dell intero range di valori dei non-detects possibili (da 0 fino al DL). 4. Metodi della distribuzione ( Distributional Methods ): Si basano sull ipotesi che la forma della distribuzione dei dati non-detects sia simile a quella delle concentrazioni misurate che superano il DL. Tra questi metodi il più utilizzato è il Metodo di Cohen ( Cohen s Method ). Metodo di Cohen ( Cohen s Method) (EPA 000b, QA/G-9): è applicabile se i nondetect costituiscono il 15-50% del data set disponibile, se la forma della distribuzione dei dati senza i valori non-detect sia di tipo normale e che il DL sia sempre lo stesso. Questo metodo adatta la media e la deviazione standard per valori al di sotto del DL, basandosi sulla tecnica statistica della stima più probabile della media e della varianza, in modo che sarà possibile applicare i vari test statistici al data set. Nella applicazione di questo metodo i non-detect non si assumono mai pari a zero. Le stime derivanti dai campionamenti sono X 1, X X 3,, X n di cui i primi m valori rappresentano i dati sopra il DL. Quelli sotto il DL saranno dunque n-m. La scelta del metodo più appropriato dipende dal grado di incompletezza del data set, dalle sue dimensioni e dalla distribuzione più idonea a rappresentare i campioni. Inoltre, sempre il documento (OSWER , EPA 00) fornisce cinque raccomandazioni su come trattare un insieme di dati in cui siano presenti dei non-detects: - I Detection Limits devono sempre essere specificati e i non-detects riportati con il valore osservato se possibile. - I non-detects non devono mai essere riportati come valori zero senza specifiche giustificazioni. 6
26 - Se un analisi condotta con un Bounding Method rivela che gli effetti quantitativi della presenza di non-detects nel data set è trascurabile non sono necessari ulteriori esami. - Se si vuole procedere ad ulteriori analisi è consigliabile usare un metodo per una specifica distribuzione. - Se la quantità dei non-detects nel data set è alta (>75%) oppure se il numero di campioni è basso (n<5) nessun metodo funzionerà bene. In tal caso si può riportare la percentuale di valori al di sotto del DL, ricorrere ancora ad un Bounding Method nel quale i non-detects saranno sostituiti dal DL nel calcolo del fondo, che sarà riportato come un numero probabile considerevolmente maggiore della media reale. Per l applicazione delle presenti linee guida, si ritiene opportuno porre, in ogni caso e per qualsiasi distribuzione dell insieme dei dati, i non-detect pari al detection limit (n.d.=d.l.); laddove siano presenti differenti detection limit per un singolo parametro si dovrà associare ai non detect il limite di rilevabilità più basso. In ogni caso andrà sempre riportata la percentuale di non detect che caratterizza il data set. 6 ANALISI MULTIVARIATA In alcuni casi il sistema in esame potrebbe risultare estremamente complesso tale da rendere molto laboriosa l applicazione dell analisi statistica univariata; in questi casi si può ricorrere a metodi esplorativi del data set basati sull analisi multivariata. Gli obiettivi principali delle metodologie di analisi multivariata sono riassumibili nella semplificazione della struttura dei dati (riduzione del numero delle variabili), nell'ordinamento e nel raggruppamento (classificazione) dei campioni in esame, nello studio delle interdipendenze tra le variabili, nella formulazione e verifica di ipotesi operative. Tra le più utili metodologie di analisi multivariata si può annoverare l analisi dei cluster (cluster analysis). Tale tecnica, partendo dall analisi di un insieme di campioni multidimensionale, cerca di individuare dei sottogruppi di campioni omogenei al loro interno ed eterogenei tra loro (ad esempio sottogruppi caratterizzati dalle stesse condizioni di ph, potenziale redox, ecc.). L obiettivo ultimo è dunque quello di ottenere una vera e propria classificazione dei campioni in esame. 7 DISTRIBUZIONE DEI DATI Quando si ha a che fare con dati ambientali (in particolare, concentrazioni di specie chimiche nei comparti ambientali: suolo, acqua, aria), le distribuzioni di probabilità più comunemente utilizzate per la loro rappresentazione sono: - distribuzione gaussiana o normale - distribuzione lognormale - distribuzione gamma - distribuzione non parametrica. Nel seguito sono descritte sinteticamente le caratteristiche delle distribuzioni e i test utili per identificare quale di queste distribuzioni rappresenti al meglio l insieme di dati in esame. Distribuzione Gaussiana o normale La distribuzione Gaussiana, o normale, è una distribuzione di tipo simmetrico la cui tendenza centrale è data dal calcolo della media aritmetica dei valori X 1, X X 3,, X n delle grandezze considerate. 7
27 La forma della distribuzione normale è descritta dalla funzione Densità di Probabilità, definita da due parametri: la media aritmetica e la varianza del campione, che è indice della dispersione dei dati rispetto al valor medio. Funzione f ( x) _ 1 1 = exp x x σ π σ Media _ n 1 x i i= 1 x = n Varianza σ = 1 1 n n i= 1 _ x i x dove n è il numero di valori considerati. In Figura 1 è riportato un esempio di distribuzione normale. Figura 1 Esempio di distribuzione normale Distribuzione lognormale La distribuzione lognormale è un tipo di distribuzione asimmetrica, derivante dal calcolo della media geometrica dei valori. La sua forma è più pendente di quella di una distribuzione normale ed è delimitata a sinistra dallo zero, mentre la parte finale all altra estremità risulta avere una specie di coda più lunga di quella normale. Quindi, la distribuzione lognormale è caratterizzata da una asimmetria positiva (coda a destra) dovuta al fatto che ad un elevata frequenza di valori bassi si associa una coda di valori molto meno frequenti ma, allo stesso tempo, molto elevati. La distribuzione lognormale è generalmente definita da due parametri y e (media e varianza della variabile trasformata y = lnx). Funzione f ( x) σ y > 0 _ 1 1 = exp ln x y xσ y π σ y _ σ y x > 0, < y <, _ 8
28 Media _ y = 1 n n i= 1 y i _ = 1 n n i= 1 ln x i = Varianza σ n 1 y = n i= 1 y i _ y dove n è il numero di valori considerati. In Figura è riportato un esempio di distribuzione lognormale. Figura Esempio di distribuzione lognormale. Distribuzione Gamma - Molti data set che presentano asimmetrie possono essere rappresentati sia mediante una distribuzione lognormale che da una distribuzione di tipo gamma, specialmente nei casi in cui il numero di campioni n è inferiore a La distribuzione gamma è generalmente definita da due parametri: k (parametro di forma) e θ (parametro di scala); il loro prodotto è pari alla media aritmetica Funzione f ( x k, θ ) x 1 ( ) ( k 1) θ, x e k = x>0, k>0, θ>0 θ Γ( k) Distribuzione non parametrica Nel caso in cui non sia possibile dimostrare che i valori di un data set seguano una tra le suddette distribuzioni (ad esempio a causa dello scarso numero di campioni) o qualora risulti, dalla applicazione dei test statistici, che nessuna distribuzione approssimi bene l insieme dei dati, allora si parla di data set non parametrici. In tal caso esistono delle procedure specifiche, per l individuazione del valore rappresentativo dell insieme dei dati, indipendenti dai parametri statistici e dal tipo di distribuzione dei dati. 8 DESCRITTORI STATISTICI Nel seguito sono descritte altre grandezze statistiche utili per lo studio del tipo di distribuzione dei dati. 9

References: Art.5
 art. 17
 art.17
 articolo 3
 art.5
 Articolo 4
 art. 40
 Art. 43
 articolo 76
 articolo 300
 art. 300
 art. 30