Denne artikkelen presenterer en protokoll og programvareverktøy for kvantifisering av usikkerhet i kalibrering og dataanalyse av en semi-kontinuerlig termisk-optiske organiske/elementær karbon analysator.
Forskere fra utallige felt ofte søker å kvantifisere og klassifisere konsentrasjoner av karbonrike aerosoler som organisk karbon (OC) eller elementær karbon (EC). Dette oppnås vanligvis ved hjelp av termisk-optiske OC/EC-analyserer (Passerernær), som muliggjør måling via kontrollert termisk pyrolyse og oksidasjon under spesifikke temperatur protokoller og innenfor begrensede atmosfærer. Flere kommersielle Passerernær finnes, inkludert et halvt kontinuerlig instrument som muliggjør on-line analyser i feltet. Dette instrumentet bruker en test kalibrerings prosedyre som krever relativt hyppig kalibrering. Denne artikkelen beskriver en kalibreringsprotokoll for denne halv kontinuerlige TOA og presenterer et programvareverktøy med åpen kildekode for dataanalyse og streng Monte Carlo-kvantifisering av usikkerhet. Spesielt, inneholder programvaren verktøyet romanen betyr å korrigere for instrumentet drift og identifisere og kvantifisere usikkerheten i OC/EC delt punkt. Dette er en betydelig forbedring av usikkerhet estimering i produsentens programvare, som ignorerer delt punkt usikkerhet og ellers bruker faste ligninger for relative og absolutte feil (vanligvis fører til under-anslått usikkerhet og ofte gir ikke-fysiske resultater som vist i flere eksempeldatasett). Den viste kalibrerings protokollen og nytt programvareverktøy som muliggjør nøyaktig kvantifisering av kombinert usikkerhet fra kalibrering, repeterbarhet og OC/EC delingspunkt deles med den hensikt å bistå andre forskere med å oppnå bedre målinger av OC, EC og total karbon masse i aerosol prøver.
Evnen til å nøyaktig måle atmosfæriske konsentrasjoner av karbonrike arter er overmåte viktig for mange forskere. Karbonrike arter i Omgivelses partikler (PM, den største miljørisiko faktoren for tidlig død1) har blitt foreslått å være den viktigste komponenten i PM ansvarlig for uønskede helseeffekter og utfall2,3 ,4. Partikler karbon i atmosfæren er et kritisk klima forurensende, hvor ulike karbonrike arter er kjent for å ha variabel, selv motsatt, virkninger. Black Carbon er potensielt den nest sterkeste direkte strålingspådrivet Forcer i jordens atmosfære5,6,7,8. Når avsatt på snø og is, Black Carbon også reduserer Reflektivitet av det arktiske landskapet, forsterke absorpsjonen av sollys, og øke frekvensen av smelte9,10,11,12 . Contrastingly, hygroskopisk organiske karbon partikler fungere som Sky kondens kjerner, øke gjennomsnittlig Reflektivitet av jorden, og forårsaker en avkjølende effekt13. Nøyaktig klassifisering av samplet karbonrike materiale og samtidige kvantifisering av målings usikkerhet er dermed vesentlige aspekter ved målinger av partikler.
Skille mellom organisk og elementær karbon i en partikkel-Laden prøven kan oppnås ved hjelp av en termisk-optisk analyse14. Kommersielle, laboratorium-baserte systemer for termisk-optiske karbon analyser har blitt opprettet15,16,17 inkludert en on-line, semi-kontinuerlig analysator18 som gjør det mulig å utføre optiske analyser i felten. Det nåværende arbeidet beskriver en detaljert prosedyre for kalibrering av dette siste OCEC-instrumentet (se tabell over materialer) og deler et programvareverktøy med åpen kildekode for den strenge kvantifisering av kalibrering og analyse usikkerhet. Selv om den første utgivelsen av åpen kildekode programvare er utviklet for output filformatet på semi-kontinuerlig instrument, kan programvaren verktøyet lett utvides av andre i fremtiden for å arbeide med utganger generert av andre instrumenter.
Semi-kontinuerlig termisk-optisk organisk/elementær karbon analysator (OCEC) kvantifiserer organisk karbon (OC) og elementær karbon (EC) i et utvalg volum. Analyse prosedyren inneholder fire faser som er skissert i figur 1. For det første, et prøvevolum trekkes gjennom instrumentet, der partikler er avsatt til, og gass-fase organiske er adsorberes av, et kvarts filter. Ved opphør av prøvetaking, er kvarts filteret varmes opp gjennom en foreskrevet temperatur protokoll i en inert, helium (HE) atmosfære. Under denne prosedyren er en del av det karbonrike materialet termisk pyrolysert fra kvarts filteret. Gass eksos er videresendt til en fast temperatur, mangan oksid (MnO2) ovn som konverterer pyrolysert karbonrike arter i KARBONDIOKSID (co2). Den tids løse konsentrasjonen av frigitt CO2 måles deretter av en ikke-dispersive infrarød (NDIR) detektor. Etter første oppvarming i he miljøet, er prøven varmes opp gjennom en lignende protokoll i en oksiderende (Ox) miljø. I nærvær av oksygen, ildfaste karbonrike arter igjen på kvarts filter er oksidert og deretter videresendt gjennom MnO2 ovn og NDIR detektor på samme måte. Når samplet karbonrike arter har blitt fullt utviklet fra kvarts filter, en endelig in-test kalibrering prosedyren utføres. Et fast antall (nominelt 0,8 mL) av en 5% metan (CH4)-helium blandingen innføres i instrumentet, oksidert i MnO2 ovn, konvertert til co2, og deretter målt ved NDIR. Den integrerte NDIR signalet under denne in-test kalibrerings fasen (kalt CH4-loop) tilsvarer den kjente karbon massen (INTRODUSERT som ch4) og derfor kvantifiserer følsomheten av NDIR, som kan drive i tid. Dette tiltaket av NDIR følsomhet blir deretter brukt til å antyde karbon masser fra NDIR signalet i løpet av tidligere he-og Ox-faser av analysen.
Figur 1: termisk-optisk analyse prosedyre. Analyse prosedyre av termisk-optiske OCEC instrument. Etter oppkjøpet av et utvalg på kvarts filter (trinn 0), tre viktigste analyse trinn utføres. To termiske protokoller først i en helium atmosfære (HE-fase, trinn 1) og deretter i en oksiderende atmosfære (okse-fase, trinn 2) er henrettet, der karbonrike komponenter er pyrolysert/oksidert fra kvarts filter, ytterligere oksidert i CO2 innen en sekundær katalyserende ovn, og deretter målt ved en NDIR detektor. En siste in-test kalibrering prosedyre (CH4-loop, trinn 3) er utført, der oksidasjon av en kjent masse av metan gir et mål på følsomhet for NDIR detektoren. Detektor følsomheten er kombinert med NDIR-målt CO2 under he-og Ox-faser for å kvantifisere karbon masse lasting på kvarts filteret. Vennligst klikk her for å se en større versjon av dette tallet.
Massen av injisert karbon under CH4-loop kan være følsomme for operasjonelle forhold slik at intermitterende kalibrering er nødvendig. Denne kalibreringen bruker en vandig sukrose løsning av kjente konsentrasjon (ca 0,99%m/m) som en ekstern standard. Gjentatte tester utføres der ulike kjente volumer av sukrose løsningen introduseres i tilfeldig rekkefølge i instrumentet, og den termiske-optiske analysen utføres. Hver gjentatt test (dvs. resultater fra hver injeksjon og påfølgende analyse) gir et integrert NDIR-signal under CH4-loop (“kalibrerings området”) og et integrert NDIR-signal for total karbon (dvs. signal under he-og Ox-fasene; referert til som “totalt areal”), som tilsvarer den kjente massen av karbon i sukrose. Lineær regresjon av den kjente karbon massen med instrument rapporterte “totalareal” gir et mål på gjennomsnittlig NDIR-følsomhet. Denne følsomheten er så kombinert med gjennomsnittet “kalibrerings området” for å gi kalibrert kunnskap om karbon massen injisert under CH4-loop.
Utover kalibrering er en viktig utfordring i tolkningen av resultatene fra OCEC-instrumentet fastsettelse av de relative brøkdeler av OC og EC i den målte prøven. Som OC pyrolyzes under he-fasen av temperatur protokollen, en brøk tegn på filteret i stedet for å bli utgitt, oksidert i MnO2 ovn, og oppdaget av NDIR. I teorien er dette forkullede ildfast OC (betegnet pyrolysert karbon, PC) forblir på filteret til Ox-fase, når det er oksidert sammen med EC. Følgelig naivt merking alle utviklet karbon i løpet av he-fase som OC og utviklet seg karbon i Ox-fase som EF fører til partisk anslag over den sanne brøkdel av OC og EC. En vanlig måte å definere delt punkt (dvs. det øyeblikket i tid når alle tidligere karbon evolusjon regnes som OC og alle påfølgende karbon evolusjon som EC) er termisk/optisk-transmisjon (TOT) metode19. Her er en laser rettet gjennom kvarts filter under termisk analyse og dens makt (optisk nedstrøms av kvarts filter) oppdages av en photodetector. Ved å anta at OC ikke er optisk aktiv på laser bølgelengde (dvs. OC negligibly absorberer lys) og PC aksjer optiske egenskaper med EC, kan dele punktet anslås. Forutsetningen er å først måle den svekkede laser kraften ved starten av analysen. Etter hvert som OC utvikler seg (delvis i lett absorberende PC), forbedres dempingen av laseren slik at photodetector signalet synker. Som Ox-fase er angitt og EC/PC er co-utviklet, demping er redusert, og photodetector signalet begynner å øke. Dele punktet er definert som forekomsten i tid når den målte laser effekten går tilbake til sin opprinnelige verdi. Mens logikken i denne tilnærmingen er robust, avhenger resultatet på ovennevnte forutsetninger. Som sådan er det vanlig å erklære at rapporterte OC-og EF-resultater er “operativt definert” – det vil si at de er spesifikke for teknikken som brukes til å evaluere Delings punktet14,20,21.
Selv om det er sant i teorien om at OC er utviklet i he-fase og PC/EC er utviklet i okse-fase, det har blitt observert at utviklingen av PC/EC kan faktisk skje i løpet av han-fase på grunn av ulike mekanismer22,23,24 ,25, slik at den sanne delt punkt kan oppstå før innføringen av oksygen. Dette tvetydighet i prediksjon av hvor dele punktet bør ligge, kombinert med usikkerhet i forutsetningene for de optiske egenskapene til OC, PC og EC, tyder på at usikkerhet i delt punkt kan være en dominerende kilde til usikkerhet i målte karbon Massene. Heldigvis muliggjør metodisk estimering av dele punktet via TOT-metoden et objektiv anslag over brøk punkt usikkerhet. Til forfatternes kunnskap derimot, er det ingen direkte estimering (og påfølgende forplantning) av delt punkt usikkerhet i produsentens programvare; rapportert total usikkerhet er i stedet beregnet med faste relative og absolutte komponenter26. Som en del av dette arbeidet presenteres en ny teknikk for å anslå at det er usikkerhet knyttet til punkt usikkerheten − “dempings nedgang”-teknikken. I denne teknikken, usikkerhet i Splitt punktet er kvantifisert som forskjellen mellom utviklet karbon massen ved nominell splitt (via TOT metoden) og utviklet karbon massen på et etterfølgende punkt der laser demping har sunket utover noen kritiske kvantitet, en foreskrevet brøkdel av den første dempingen. En kritisk reduksjon av demping er anslått basert på usikkerheten i laser demping i forhold til den opprinnelige verdien; i teorien, denne tilnærmingen fanger usikkerhet i nøkkelen demping-samsvarende prinsippet om TOT metoden. Videre, for å vurdere (i hvert fall delvis) delt punkt usikkerhet på grunn av antatt optiske egenskaper av PC og EC, den foreslåtte kritiske demping nedgangen er utvidet med en faktor på to.
Denne artikkelen presenterer en detaljert protokoll for kalibrering av OCEC instrumentet sammen med et programvareverktøy for strengt kvantifisere kalibrering og analyse usikkerhet. For det første deler 1 til 3 av protokollen skissere instruksjoner for å lage den vandige sukrose løsning, forberede instrumentet for kalibrering, og anskaffe kalibreringsdata. Del 4 bruker romanen, åpen kildekode, Programvareverktøyet (se tabell over materialer) til å analysere kalibreringsdata via programvaren grafiske brukergrensesnitt. Del 5 spesifiserer betraktninger for å anskaffe en prøve ved hjelp av OCEC-instrumentet og Seksjon 6 beskriver bruken av den ovennevnte programvaren for beregning av karbon masser og tilhørende usikkerhet, inkludert bidrag fra vurderingen av delingen Punkt. Nye teknikker for å forbedre behandlingen av OCEC data, inkludert “dempings reduksjon”-teknikken som ble innført ovenfor, er beskrevet i programvaredokumentasjonen på Internett.
Innenfor den presenterte programvaren verktøyet, er kalibrering konstanter, målt karbon masser, og tilhørende usikkerhet beregnet ved hjelp av en Monte Carlo (MC) metode. Denne prosedyren sprer feil som ikke er, til forfatternes kunnskap, for tiden vurderes i produsentens proprietær programvare. For kalibrering inkluderer disse feilkildene usikkerhet i konsentrasjonen av sukrose i den vandige oppløsningen, nøyaktigheten i volumet av anvendt sukrose-løsning (instrument nøyaktighet, reproduserbarhet for Inter-bruker og repeterbarhet av brukeren), og usikkerhet i lineær regresjon. Med hensyn til dataanalyse, anses feilkilder inkluderer kalibrering usikkerhet og repeterbarhet, og viktigst, anslått usikkerhet i delt punkt. Til syvende og sist gjør programvaren en bruker til å nøyaktig kvantifisere usikkerhet i kalibrering av instrumentet og spre denne usikkerheten sammen med delt punkt estimering i beregningen av karbon masser. Dette representerer en merkbar forbedring av produsentens protokoll, ved å direkte vurdere viktige kilder til feil i målingen, i stedet for å estimere delvis usikkerhet ved hjelp av en fast ligning.
Tabell 3 viser bidraget fra de spesifikke kildene til usikkerhet til masse kalibrerings konstanten for eksempel saken som er beskrevet i tabell 1 og figur 2. Kumulativ kalibrering usikkerhet som følge av bias feil i NDIR detektor, bias feil i konsentrasjonen av sukrose, og presisjon og bias feil i Pipet tert volumer er oppført. Bias feil i NDIR detektor (dvs. varians i “kalibrerings området”) har en tendens til å dominere, med bias i pipettering prosedyren blir nest viktigste (men ganske liten i representative eksempel). Riktig estimering av pipettering feil er dermed avgjørende for å sikre nøyaktig kvantifisering av total kalibrering usikkerhet; med henvisning til protokoll trinn 4.3.2, er det derfor antydet at repeterbarhet og reproduserbarhet mellom brukere vurderes for hver gruppe av brukere og pipette. I kontrast, usikkerhet på grunn av konsentrasjonen av sukrose i den eksterne standarden er negligibly liten. Videre synes det å være ubetydelig bidrag fra regresjon usikkerhet, en sannsynlig konsekvens av den gode linearitet av instrumentet − koeffisient av besluttsomhet (R2) for en lineær tilpasning av kalibreringsdata generelt overstiger 99,95% . Hvis kalibreringsdataene ikke er tilstrekkelig lineære, advarer programvaren automatisk brukeren, som deretter kan feilsøke datasettet via den boolske kontrollen som er angitt i protokoll trinn 4,2. brukeren kan deretter endre kalibreringsdata settet ved å anskaffe erstatnings data etter behov.
Ansett usikkerhet | Full (6-punkts) kalibrering |
NDIR-bias | ± 2,61% |
+ Sukrose løsning | ± 2,61% |
+ Pipette | ± 2,78% |
Nominell resultat [μgC] | 18,49 |
Tabell 3: usikkerhet i masse kalibrerings konstant. Representativ bidrag av usikkerhet i kalibreringen av OCEC-instrumentet for eksempel seks-punkts kalibrering (se tabell 1). Total kalibrering usikkerhet er dominert av bias i NDIR detektor med feil på grunn av nøyaktighet i pipettering av sukrose løsningen, inkludert menneske-feil (Inter-User reproduserbarhet og intra-User repeterbarhet), som nest viktigste, etterfulgt av usikkerhet i lineær regresjon og sukrose konsentrasjon (som begge er ubetydelig).
Kalibrering av OCEC-instrumentet er en tidkrevende prosedyre, som vanligvis krever 2 − 3 timer for å fullføre avhengig av lengden på den termiske protokollen som brukes. En raskere kalibrerings prosedyre er ønskelig. For dette formål ble effekten av en modifisert, avkortet kalibreringsprotokoll analysert med det presenterte Programvareverktøyet. Den utviklede MC prosedyren ble utført ved hjelp av alle mulige delsett av eksempel kalibreringsdata oppført i tabell 1 -begrenset til tilfeller med tre eller flere data og minst én blank måling. Alle resulterende massekalibrering konstanter fra denne analysen er plottet i Figur 4 som en funksjon av antall kalibrering data som brukes, der kalibrering konstanter har blitt normalisert av full (6-punkts) kalibrering resultat. Det er ikke overraskende at usikkerheten i kalibrerings konstanten øker når tilgjengelige kalibreringsdata reduseres. Viktigere imidlertid, betyr alle avkortet kalibreringer faller innenfor 2σ CI av full kalibrering resultatet, som er en konsekvens av ovennevnte linearitet av instrumentet. Denne konsistensen i MC-gjennomsnittet tyder på at en modifisert, raskere kalibrering bestående av få kalibreringsdata kan benyttes som en “bump test” kontroll av OCEC instrumentets kalibrering. Det vil si at hvis MC-gjennomsnittet av et 3-punkts kalibreringsdata sett er innenfor 2σ CI av den eksisterende kalibreringen, er det sannsynlig at OCEC-instrumentet ikke krever ny kalibrering. Det er også tydelig i Figur 4 at kalibrerings usikkerheten reduseres med flere kalibreringsdata, men reduksjonen i usikkerhet lider av avtagende avkastning. Med henvisning til tabell 3 og dens diskusjon ovenfor, siden kalibrering usikkerhet er DOMINERT av NDIR bias (kvantifisert med standard feil av “kalibrerings områder”), den marginale reduksjonen i kalibrering usikkerhet ved å inkludere nth datapunkt kan anslås med faktor √ (1 − 1/n). I representant eksempelet er derfor den marginale reduksjonen i usikkerhet større når du flytter fra en tre-til en fire-punkts kalibrering enn ved flytting fra en fem-til en seks-punkts kalibrering. Den utviklede Programvareverktøyet, som kan utføres etter oppkjøpet av hvert kalibreringsdata punkt (dvs. etter hver gjentakelse av protokollen trinn 3,3), tillater brukeren å kvantifisere kalibrering usikkerhet gjennom oppkjøpet av data. Kritisk, gjør denne evnen brukeren til å ikke bare bestemme på tilstrekkelig konvergens av kalibreringen i sammenheng med sin usikkerhet, men også oppdage tilstedeværelsen av falske data-det vil si en trinnvis reduksjon i kalibrering usikkerhet som er spesielt annerledes enn forventet høydepunkter for brukeren at det sist mottatte kalibreringsdata punktet kan være feil.
Figur 4: vurdering av utvalgsstørrelsen på kalibrerings usikkerheten. Beregnet massekalibrering konstant for alle mulige kombinasjoner av kalibreringsdata oppført i tabell 1 (som krever minimum tre data, inkludert minst én blank) normalisert av resultatet ved hjelp av alle seks data. Relativ usikkerhet i kalibrerings konstanten reduseres med en økning i antall kalibreringsdata. Det blå skyggelagte området i figuren tilsvarer 2σ CI av kalibrerings konstanten som er beregnet ved hjelp av alle kalibreringsdata. Det er tydelig at alle nominelle resultater er innenfor dette CI-tyder på at-selv om usikkert-en avkortet kalibrering prosedyre av bare tre kalibrering datapunkter kan brukes som en “bump test” kontroll av instrumentets kalibrering. Vennligst klikk her for å se en større versjon av dette tallet.
Beregnede karbon masser og usikkerheter for fire eksempeldatasett er beskrevet i Tabell 4; disse dataene er Hentet fra målinger av karbonrike utslipp fra hydrokarboner, som er27,29, gassturbiner29og fin modus (< 2 μm) karbonrike partikler innhentet fra sub-Ocean sediment prøver30. Beregnede (MC-gjennomsnitt) OC-, EC-og TC-masser vises i Tabell 4 ved siden av den beregnede kritiske dempings NEDGANGEN og EC/TC-forholdet for hvert datasett, som viser bredden på eksempeldataene i en samplet karbon komposisjon. Oppsummert i Tabell 4 er kildene og generell usikkerhet i karbon massen ved hjelp av presentert Programvareverktøyet i forhold til de som rapporteres av instrumentet. I Programvareverktøyet gir kombinasjonen (produkt distribusjon) av kalibrerings usikkerhet og repeterbarhet den totale usikkerheten i MC-beregnet TC-masse (kvantifisert som − 8.32/+ 8.40% i det foreliggende arbeidet), som er uavhengig av usikkerhet i og derfor fungerer som en nedre grense i usikkerheten av OC og EC massene. Disse representative kalibrerings usikkerheten brukes på hvert eksempeldatasett og samtidig opprettholde den nominelle masse kalibrerings konstanten som brukes i de opprinnelige analysene.
Tabell 4: usikkerhet i dataanalyse. Bidrag til usikkerhet i OCEC karbon masser for fire eksempeldatasett fra et bredt spekter av kilder og utført av ulike laboratorier27,29,30. (a) numeriske tall resultater fra eksempel datasettene: OC, EC og TC-masser, kritisk dempings nedgang for kvantifisering av usikkerhet i splittet punkt, og elementær-til-total-karbon-forhold. (b) et sammendrag av usikkerhet i de beregnede karbon massene. Bidra kilder til usikkerhet omfatter at av massen kalibrering konstant, repeterbarhet av kalibreringen prosedyren, og usikkerhet i delt punkt (i forhold til TC masse) som tilsvarer den kritiske demping nedgangen oppført i (a). Usikkerhet i karbon massene ved hjelp av de faste formlene (EQ. (1)) som brukes av instrumentet, vises også i (b). Rød og gul utheving av dataene tilsvarer under-og over-estimering av usikkerhet, henholdsvis når du bruker den faste ligninger i forhold til dagens metode. I de fleste tilfeller, instrumentet under-estimater karbon masse usikkerhet men hvis målt OC eller EC massen er liten, kan instrumentet over-anslå usikkerhet i forhold til dagens programvare. Vennligst klikk her for å laste ned denne filen.
For disse eksemplene ble metoden for reduksjon av demping av kvantifisere brukt. OCECgo-beregnede kritiske verdier av demping nedgang varierte fra 1,342% til 2,059% resulterer i delt punkt usikkerhet fra 0,10% til 4,50% av TC masse. Mens de næringsdrivende verdiene for reduksjon av demping faktisk er noe subjektiv − særlig den sysselsatte faktoren for å anslå usikkerhet i dele punktet på grunn av at optiske egenskaper − disse eksemplene fremhever avhengigheten av del punkt usikkerhet på de spesifikke analyse dataene. For eksempel, delt punkt usikkerhet er følsom for skråningen av AVEC tomten i nærheten av den nominelle delt punkt. Vurder datasettet “A” som tilsvarer eksempeldataene i Figur 3 og datasettet “D”. til tross for å ha lignende kritisk demping avtar, de relativt grunne og bratte bakkene i sine respektive AVEC tomter (f. eks, se Figur 3b for datasett “A”) gi den største og minste delt punkt usikkerhet på 4,50% og 0,10% av TC Mass. I tillegg viser eksempeldataene at påvirkningen av usikkerhet i dele punktet avhenger i stor grad av skalaen i forhold til de nominelle OC-og EC-massene. Betrakt eksempeldatasett “B” og “C”, som har nesten identiske delt punkt usikkerhet (≈ 1,22% av TC masse); datasettet “C” inneholder ≈ 43% OC, mens datasettet “B” inneholder ≈ 8%. den lavere relative mengden OC i sistnevnte resulterer i en nesten dobling av OC masse usikkerhet. Kritisk, disse resultatene markere kravet til direkte vurdere delt punkt usikkerhet i sammenheng med analysen ‘ AVEC data og samlede karbon masser.
Produsent-rapportert usikkerhet i karbon masser er også vist i Tabell 4. Disse anslagene tar ikke direkte hensyn til usikkerhet i kalibrerings-og dele punktet, men beregnes i stedet ved hjelp av de faste relasjonene som vises i EQ. (1)26, der mjeg representerer nominell masse av den spesifikke karbon komponent.
1
Disse faste forbindelsene tillater estimert usikkerhet i OC og/eller EF-massen til å være kunstig mindre enn TC-massen, denne tilstanden oppstår når OC eller EF-massen er mindre enn en tredjedel av TC-massen, slik tilfellet er for datasett “A”, “B” og “D”. Dette resultatet er ikke-fysisk, siden relativ usikkerhet i OC og EC massene må være avgrenset under av det av TC masse, en konsekvens av utbredelsen av delt punkt usikkerhet i den beregnede OC og EC massene. Celler markert rødt og gult i tabellen tilsvarer under-og over-estimater av karbon masse usikkerhet ved bruk av EQ. (1) fra produsenten. Den faste ligninger under-anslag TC masse usikkerhet for alle fire eksempler, en konsekvens av beregnede TC massene blir tilstrekkelig store. I de fleste tilfeller, de faste ligninger også under-estimat EC og OC masse usikkerhet, unntatt der OC (datasett “B”) og EC (datasett “D”) var tilstrekkelig liten til å forårsake over-estimering via EQ. (1). Denne asymptotiske øker usikkerheten via EQ. (1) er enig med den nåværende programvaren siden usikkerheten i OC og EC massene på grunn av at delt punkt er avhengig av deres absolutte størrelse; usikkerhet via de faste formlene klarer imidlertid ikke å spore med de aktuelle programvare, som direkte vurderer og overfører del punkt usikkerhet i sammenheng med spesifikke analysedata.
Bruk av en MC rammeverk i presentert Programvareverktøyet er nødvendig for å nøyaktig spre komponenten usikkerhet gjennom ikke-lineære algoritmer for termisk-optiske OC/EC analyser. Det er viktig å merke seg imidlertid gjennom sin iboende tilfeldig natur, at MC-metoder ikke er deterministisk og har en tendens til å gi inkonsekvente resultater dersom sysselsatt antall MC-uavgjorte/repetisjoner (se protokoll trinn 4.3.3 og 6.5.2) er utilstrekkelig – beslektet med en Statistisk undersized prøven. Derfor er det en iboende konsistens kontra beregning tid trade-off som skal vurderes ved behandling av data ved hjelp av OCECgo. Det er derfor anbefalt for brukere å utføre foreløpige behandling og feilsøking av data ved hjelp av et lite antall MC-trekninger (f. eks 104). Når beregninger gir tilfredsstillende resultater, bør brukeren da øke antall MC trekker (til 106-108) for å gi et resultat som er mindre påvirket av diskret og tilfeldig karakter av MC-metoden. Utover den nødvendige “operative definisjon” av termisk-optiske OC/EC analyser, er det andre begrensninger i behandlingen av disse dataene som også må anerkjennes når du bruker OCECgo å beregne og rapportere OC/EC data. For det første, NDIR-baserte instrumenter (slik som det som er oppført i tabellen av materialer) lider av drift i NDIR signalet som må rettes opp. I den nåværende protokollen (se trinn 6.3.2 og OCECgo dokumentasjon), kan en ny tilnærming til riktig drift i NDIR detektoren være frivillig ansatt av brukeren. Mens, i forfatterne erfaring dette gir bedre resultater over produsentens standard lineær NDIR korreksjon, må det bemerkes at usikkerhet i denne NDIR korreksjon er utfordrende om ikke umulig å kvantifisere og derfor forblir en savnet komponent av usikkerhet i beregningen av karbon masser. I en lignende blodåre, er det også utfordrende å kvantifisere usikkerhet i den nødvendige antakelsen om at PC og EC dele optiske egenskaper. Hvis dette er valgt (se protokoll trinn 6.4.1), forsøker den kritiske dempingen for tilbakefall å konservativt bundet effekten av denne antakelsen gjennom en subjektiv Utvidelses faktor. Dette er imidlertid nødvendigvis bare et anslag, og brukeren foreslås å vurdere effekten av denne Utvidelses faktoren (dvs. den kritiske dempings nedgangen) på deres spesifikke data. OCECgo er gitt som en åpen kildekode-verktøyet slik at det kan lett bli utvidet av forfatterne og andre interesserte samarbeidspartnere til ikke bare grensesnitt med andre instrumenter, men også å inkludere andre nyttige, felt-spesifikke, funksjonalitet. Samlet er utviklet åpen kildekode Programvareverktøyet kombinert med detaljert kalibrering prosedyren ment å bidra til å oppnå mer nøyaktige målinger av OC, EC, og TC masse i aerosol prøver mens effektivisere den robuste beregning av måling Usikkerhet.
The authors have nothing to disclose.
Dette arbeidet ble støttet av Natural Sciences and Engineering Research Council of Canada (NSERC) FlareNet strategiske nettverk (Grant # 479641), NSERC Discovery Research Grants (Grant #06632 og 522658), og Natural Resources Canada (prosjektleder, Michael Lag). Forfatterne er takknemlige for de som delte rå datafiler for bruk som representative eksempler i dette arbeidet.
10% oxygen gas in helium | Local gas supplier | – – – | Primary or certified standard preferred |
5% methane gas in helium | Local gas supplier | – – – | Primary or certified standard preferred |
Distilled, de-ionized water | Harleco | 6442-85 | ASTM D1193-91 Type II or Type I (preferred) |
Filter punch tool | Sunset Laboratories Inc. | – – – | Included with carbon analyzer |
Filter removal tool | Sunset Laboratories Inc. | – – – | Included with carbon analyzer |
Glass jar (4 oz.) | ULINE | S-17982P-BL | Or suitable equivalent; borosillicate glass preferred |
Helium gas | Local gas supplier | – – – | Ultra-high purity (> 99.999%) or better preferred |
High-accuracy thermal gas mass flow meter | Bronkhorst | EL-FLOW Prestige | For accurate measurement of sample volume (see Protocol step 5) |
High-purity sucrose | Sigma Aldrich | S9378 | Purity ≥ 99%m/m or higher preferred |
Lint-free tissues | Kimtech | 34155 | Or suitable equivalent |
MatLab Runtime (R2016a or newer) | MathWorks Inc. | mathworks.com | Search "runtime compiler" and install appropriate version for the operating system |
Non-powdered, disposable, plastic gloves | VWR | 89428-752 | Or suitable, properly-sized equivalent |
OCECgo software | Carleton University, Energy and Emissions Research Lab. | GitHub Repository | Source and build distributions of the software are available on GitHub |
Oxygen trap | Supelco | 22449 | Or suitable GC-quality equivalent |
Pipette | Eppendorf | 3120000020 | Model: Research® Plus 0.5 – 10 μL – Or any single-channel, adjustable volume, manual pipette |
Pipette tips | Eppendorf | 022492012 | Model: epT.I.P.S.® Standard, 0.1 – 20 μL |
Precision scale / balance | AND | FX-3000IWP | Precision balance with capacity > 1 kg |
Quartz filters | Pall | 7202 | Model: Tissuquartz 2500 QAT-UP – 47 mm |
Semi-continuous thermal-optical organic/elemental carbon analyzer | Sunset Laboratories Inc. | – – – | Model 4 semi-continuous analyzer |