Dit artikel presenteert een protocol en software tool voor de kwantificering van onzekerheden in de kalibratie en data-analyse van een semi-continue thermische optische organische/Elemental Carbon Analyzer.
Onderzoekers uit talloze gebieden zijn vaak gericht op het kwantificeren en classificeren van concentraties van koolzuurhoudende aerosols als organische koolstof (OC) of elementair koolstof (EC). Dit wordt vaak bereikt met behulp van thermische optische OC/EC-analyzers (TOAs), die metingen mogelijk maken via gecontroleerde thermische pyrolyse en oxidatie onder specifieke temperatuur protocollen en binnen beperkte atmosferen. Er bestaan verschillende commerciële Toa’s, waaronder een semi-continu instrument dat on-line analyses in het veld mogelijk maakt. Dit instrument maakt gebruik van een in-test kalibratieprocedure die een relatief frequente kalibratie vereist. Dit artikel beschrijft een kalibratieprotocol voor deze semi-continue TOA en presenteert een open-source software tool voor data-analyse en strenge Monte Carlo kwantificering van onzekerheden. Met name de software tool bevat nieuwe middelen om te corrigeren voor instrument drift en identificeren en kwantificeren de onzekerheid in de OC/EC split point. Dit is een aanzienlijke verbetering ten opzichte van de onzekerheids schatting in de software van de fabrikant, die een gesplitste punt onzekerheid negeert en anderszins vaste vergelijkingen gebruikt voor relatieve en absolute fouten (in het algemeen leidend tot ondergeschatte onzekerheden en vaak niet-fysieke resultaten opleveren, zoals wordt gedemonstreerd in verschillende voorbeeldgegevenssets). Het Gedemonstreerde kalibratieprotocol en het nieuwe softwareprogramma zorgen voor een nauwkeurige kwantificering van gecombineerde onzekerheden van kalibratie, herhaalbaarheid en OC/EC-splitspunt worden gedeeld met de intentie om andere onderzoekers te assisteren bij het bereiken van betere metingen van OC, EC en totale koolstof massa in aërosolmonsters.
De mogelijkheid om de atmosferische concentraties van koolzuurhoudende soorten nauwkeurig te meten is voor veel onderzoekers uitermate belangrijk. Koolzuurhoudende soorten in zwevende deeltjes (PM, de grootste milieurisico factor voor vroegtijdige dood1) zijn voorgesteld als het belangrijkste onderdeel van PM dat verantwoordelijk is voor nadelige gezondheidseffecten en resultaten2,3 ,4. De deeltjes koolstof in de atmosfeer is een kritische klimaat vervuilend, waar bekend is dat verschillende koolzuurhoudende soorten variabel, zelfs tegenovergesteld, effect hebben. Zwarte koolstof is mogelijk de tweede sterkste directe radiatieve forcer in de atmosfeer van de aarde5,6,7,8. Bij het storten op sneeuw en ijs vermindert Black Carbon ook de reflectie van het Arctische landschap, verbetert de absorptie van zonlicht en verhoogt het de smelt snelheid9,10,11,12 . Contrastingly fungeren hygroscopische organische koolstofdeeltjes als condenscondensatie kernen, waardoor de gemiddelde reflectiviteit van de aarde toeneemt en een verkoelend effect13ontstaat. Nauwkeurige classificatie van gesamplede koolstofhoudende materialen en gelijktijdige kwantificering van meetonzekerheden zijn dus essentiële aspecten van de metingen van deeltjes.
Differentiëren tussen organisch en elementair koolstof in een monster met een deeltjes lading kan worden bereikt met behulp van een thermische-optische analyse14. Commerciële, laboratorium-gebaseerde systemen voor thermische-optische koolstof analyses zijn gemaakt op15,16,17 , met inbegrip van een on-line, semi-continue Analyzer18 die het mogelijk maakt de uitvoering van thermische optische analyses in het veld. Het huidige werk beschrijft een gedetailleerde procedure voor het kalibreren van dit laatste OCEC-instrument (Zie tabel met materialen) en deelt een open-source software tool voor de strikte kwantificering van de onzekerheden van kalibratie en analyse. Hoewel de eerste release van de open-source software is ontworpen voor het uitvoerbestandsformaat van het semi-continue instrument, kan de software tool in de toekomst gemakkelijk worden uitgebreid door anderen om te werken met outputs die door andere instrumenten worden gegenereerd.
De semi-continue thermische-optische organische/Elemental Carbon Analyzer (OCEC) kwantificeert organische koolstof (OC) en elementair koolstof (EC) in een monstervolume. De analyseprocedure bevat vier fasen die in Figuur 1worden beschreven. Ten eerste wordt een monstervolume getrokken door het instrument, waar deeltjes worden afgezet, en gasfase Organics worden geadsorreven door een kwarts filter. Bij het beëindigen van de bemonstering wordt het kwarts filter verwarmd door middel van een voorgeschreven temperatuur protocol in een inerte, helium (he) atmosfeer. Tijdens deze procedure wordt een deel van het koolzuurhoudende materiaal thermisch gepyrolyseerd uit het kwarts filter. Gasvormig uitlaatgas wordt doorgestuurd naar een vaste temperatuur, mangaan oxide (MnO2) oven die gepyrolyseerd koolzuurhoudende soorten omzet in koolstofdioxide (co2). De Tijdopgeloste concentratie van de vrijgekomen CO2 wordt vervolgens gemeten door een niet-dispersieve infrarooddetector (NDIR). Na de initiële Verwarming in de he omgeving wordt het monster verwarmd via een soortgelijk protocol in een oxiderende (OX) omgeving. In aanwezigheid van zuurstof worden refractaire koolzuurhoudende soorten die op het kwarts filter blijven, geoxideerd en vervolgens via de MnO2 -oven en NDIR-detector op dezelfde manier doorgestuurd. Eenmaal gesamplede koolzuurhoudende soorten zijn volledig geëvolueerd uit het kwarts filter, een laatste in-test kalibratieprocedure wordt uitgevoerd. Een vaste hoeveelheid (nominaal 0,8 mL) van een 5% methaan (h4)-helium mengsel wordt ingebracht in het instrument, geoxideerd in de MnO2 oven, omgezet in co2, en vervolgens gemeten door de NDIR. Het geïntegreerde NDIR-signaal tijdens deze in-test kalibratie fase (aangeduid als de CH4-lus) komt overeen met de bekende koolstof massa (geïntroduceerd als CH4) en kwantificeert daarom de gevoeligheid van de NDIR, die in de tijd kan drijven. Deze meting van de NDIR-gevoeligheid wordt vervolgens gebruikt om koolstof massa’s uit het NDIR-signaal te afleiden tijdens de eerdere he-en Ox-fasen van de analyse.
Figuur 1: thermische-optische analyseprocedure. Analyseprocedure van het thermisch optische OCEC-instrument. Na de overname van een monster op het kwarts filter (stap 0) worden drie belangrijke analysestappen uitgevoerd. Twee thermische protocollen eerst in een helium atmosfeer (hij-fase, stap 1) en vervolgens in een oxiderende atmosfeer (Ox-fase, stap 2) worden uitgevoerd, waar koolstofhoudende componenten worden pyrolyzed/geoxideerd uit het kwarts filter, verder geoxideerd in CO2 binnen een secundaire katalyserende oven, en vervolgens gemeten door een NDIR detector. Een laatste in-test kalibratieprocedure (CH4-loop, stap 3) wordt uitgevoerd, waarbij de oxidatie van een bekende massa methaan een mate van gevoeligheid van de NDIR-detector biedt. Detector gevoeligheid is gekoppeld met NDIR-gemeten CO2 tijdens de he-en Ox-fasen om koolstof massa belasting op het kwarts filter te kwantificeren. Klik hier om een grotere versie van dit cijfer te bekijken.
De massa van de geïnjecteerde koolstof tijdens de CH4-lus kan gevoelig zijn voor operationele omstandigheden, zodat intermitterende kalibratie vereist is. Deze kalibratie maakt gebruik van een waterige sucrose-oplossing met een bekende concentratie (ongeveer 0,99%m/m) als een externe standaard. Herhaalde tests worden uitgevoerd waarbij verschillende bekende volumes van de sacharoseoplossing in willekeurige volgorde in het instrument worden geïntroduceerd en de thermische optische analyse wordt uitgevoerd. Elke herhaalde test (d.w.z. resultaten van elke injectie en daaropvolgende analyse) levert een geïntegreerd NDIR-signaal op tijdens de CH4-lus (“kalibratiegebied”) en een geïntegreerd NDIR-signaal voor totale koolstof (d.w.z. signaal tijdens de he-en Ox-fasen; Doorverwezen als “totale oppervlakte”), wat overeenkomt met de bekende massa koolstof in sacharose. Lineaire regressie van de bekende koolstof massa met instrument-gerapporteerde “totale oppervlakte” biedt een meting van de gemiddelde NDIR-gevoeligheid. Deze gevoeligheid wordt vervolgens gekoppeld aan het gemiddelde “kalibratiegebied” om gekalibreerde kennis van de koolstof massa te leveren die tijdens de CH4-lus wordt geïnjecteerd.
Naast kalibratie is een belangrijke uitdaging bij de interpretatie van de resultaten van het OCEC-instrument de bepaling van de relatieve fracties OC en EC in het gemeten monster. Als OC pyrolyzes tijdens de he-fase van het temperatuur protocol, een breuk tekens op het filter in plaats van wordt vrijgegeven, geoxideerd in de MnO2 oven, en gedetecteerd door de NDIR. In theorie blijft deze verkroerde refractaire OC (genaamd gepyrolyseerd Carbon, PC) op het filter tot de Ox-fase, wanneer het wordt geoxideerd naast de EC. Een naïef etikettering van alle geëvolueerde koolstof tijdens de he-fase als OC en geëvolueerde koolstof in de Ox-fase, aangezien de EG leidt tot partijdige schattingen van de ware Fractie van OC en EG. Een gemeenschappelijk middel om het splitspunt te definiëren (d.w.z. het moment waarop alle voorgaande koolstof evolutie wordt beschouwd als OC en alle daaropvolgende koolstof evolutie als EC) is de thermische/optische-doorlatendheid (tot) methode19. Hier wordt een laser door het kwarts filter geleid tijdens de thermische analyse en wordt het vermogen (optisch stroomafwaarts van het kwarts filter) gedetecteerd door een foto detector. Door te veronderstellen dat OC niet optisch actief is op de laser golflengte (d.w.z. dat OC een verwaarloosbaar licht absorbeert) en PC deelt optische eigenschappen met EC, kan het splitspunt worden geschat. Het uitgangspunt is om eerst het verzwakte Laser vermogen te meten aan het begin van de analyse. Naarmate OC evolueert (deels in lichtabsorberende PC), wordt de demping van de laser zodanig verbeterd dat het fotodetector signaal daalt. Als de Ox-fase wordt ingevoerd en EC/PC zijn co-geëvolueerd, verzwakking wordt verminderd, en de foto detectorsignaal begint te verhogen. Het splitspunt wordt gedefinieerd als het exemplaar in de tijd waarin het gemeten Laser vermogen terugkeert naar de beginwaarde. Hoewel de logica van deze aanpak robuust is, is het resultaat gebaseerd op de bovengenoemde veronderstellingen. Als zodanig is het gebruikelijk om te verklaren dat de gerapporteerde OC-en EG-resultaten “operationeel gedefinieerd” zijn, d.w.z. dat ze specifiek zijn voor de toegepaste techniek om het splitspunt14,20,21te evalueren.
Hoewel het waar is in theorie dat OC is geëvolueerd in de he-fase en PC/EC zijn geëvolueerd in de Ox-fase, is geconstateerd dat de evolutie van PC/EC inderdaad kan optreden tijdens de he-fase als gevolg van verschillende mechanismen22,23,24 ,25, zodat het echte gespleten punt kan voorkomen vóór de invoering van zuurstof. Deze onduidelijkheid in de voorspelling van waar het splitspunt moet liggen, in combinatie met onzekerheid in de veronderstellingen van de optische eigenschappen van OC, PC en EG, suggereert dat onzekerheid in het splitspunt een dominante bron van onzekerheid kan zijn in de gemeten koolstof Massa. Gelukkig maakt de methodische schatting van het splitspunt via de tot-methode een objectieve raming van de onzekerheid van de gesplitste punten mogelijk. Voor de kennis van de auteurs is er echter geen directe raming (en daaropvolgende vermeerdering) van de onzekerheid van de gesplitste punten in de software van de fabrikant; gerapporteerde totale onzekerheden worden in plaats daarvan berekend met vaste relatieve en absolute componenten26. Als onderdeel van dit werk wordt een nieuwe techniek voor het inschatten van de onzekerheid van het gesplitste punt gepresenteerd − de “verzwakking daling”-techniek. In deze techniek wordt de onzekerheid in het splitspunt gekwantificeerd als het verschil tussen de geëvolueerde koolstof massa bij de nominale splitsing (via de tot-methode) en de geëvolueerde koolstof massa op een volgend punt waar Laser demping verder is gedaald dan een kritische hoeveelheid, een voorgeschreven fractie van de initiële verzwakking. Een kritische afname daling wordt geschat op basis van de onzekerheid in laser demping ten opzichte van de initiële waarde; in theorie vangt deze aanpak onzekerheid op in het sleutel dempings principe van de tot-methode. Bovendien, om te overwegen (ten minste gedeeltelijk) splitpunt onzekerheid als gevolg van veronderstelde optische eigenschappen van PC en EG, de voorgestelde kritische verzwakking daling wordt uitgebreid met een factor van twee.
Dit artikel presenteert een gedetailleerd protocol voor het kalibreren van het OCEC-instrument naast een software tool om de onzekerheden van kalibratie en analyse grondig te kwantificeren. Ten eerste, paragrafen 1 tot en met 3 van het protocol schetsen instructies voor het maken van de waterige sucrose-oplossing, het voorbereiden van het instrument voor kalibratie en het verkrijgen van kalibratiegegevens. Sectie 4 maakt gebruik van de roman, open-source, software tool (Zie tabel met materialen) om kalibratiegegevens te analyseren via de grafische gebruikersinterface van de software. Paragraaf 5 bevat overwegingen voor het verkrijgen van een monster met het OCEC-instrument en sectie 6 beschrijft het gebruik van de bovengenoemde software voor de berekening van koolstof massa’s en bijbehorende onzekerheden, met inbegrip van bijdragen uit de raming van de gesplitste Punt. Nieuwe technieken om de verwerking van OCEC-gegevens te verbeteren ― met inbegrip van de hierboven geïntroduceerde “dempings daling”-techniek worden beschreven in de online documentatie van de software.
Binnen de gepresenteerde software tool worden kalibratie constanten, gemeten koolstof massa’s en bijbehorende onzekerheden berekend met behulp van een Monte Carlo (MC)-methode. Met deze procedure worden fouten doorgegeven die niet, naar de kennis van de auteurs, momenteel worden overwogen in de eigen software van de fabrikant. Voor de kalibratie omvatten deze foutbronnen onzekerheid in de concentratie van sucrose in de waterige oplossing, nauwkeurigheid in het volume van de toegepaste sacharoseoplossing (instrument nauwkeurigheid, reproduceerbaarheid tussen gebruikers en herhaalbaarheid binnen de gebruiker), en onzekerheid in de lineaire regressie. Met betrekking tot gegevensanalyse, beschouwd als bronnen van fouten omvatten de kalibratie onzekerheid en herhaalbaarheid en, belangrijk, de geschatte onzekerheid van het splitspunt. Uiteindelijk, de software stelt een gebruiker in staat om nauwkeurig kwantificeren onzekerheid in de kalibratie van het instrument en propageren deze onzekerheid naast die van splitpunt schatting in de berekening van de koolstof massa’s. Dit is een opmerkelijke verbetering ten opzichte van het Protocol van de fabrikant, door rechtstreeks rekening te houdend met belangrijke foutbronnen in de meting, in plaats van gedeeltelijke onzekerheden te schatten met behulp van een vaste vergelijking.
Tabel 3 toont de bijdrage van de specifieke bronnen van onzekerheid aan de massa-kalibratie constante voor de in tabel 1 pt Figuur 2beschreven voorbeeld case. Cumulatieve kalibratie onzekerheid als gevolg van bias-fout in de NDIR-detector, bias-fout in de concentratie van sucrose, en precisie-en bias-fout in pipettted volumes worden vermeld. Bias fout in de NDIR detector (dat wil zeggen, variantie in het “kalibratiegebied”) neigt te domineren, met vooroordelen in de Pipetteer procedure wordt de tweede belangrijkste (hoewel vrij klein in het representatieve voorbeeld). Een juiste schatting van de Pipetteer fout is dus van cruciaal belang voor een nauwkeurige kwantificering van de totale kalibratie onzekerheid; verwijzend naar protocol stap 4.3.2, wordt daarom voorgesteld om de herhaalbaarheid van de intra-gebruiker en de reproduceerbaarheid van de Inter-gebruiker te beoordelen voor elke groep gebruikers en pipetten. Daarentegen is de onzekerheid als gevolg van de concentratie van sucrose in de externe standaard verwaarloosbaar klein. Bovendien lijkt er een verwaarloosbare bijdrage van de regressie onzekerheid, een waarschijnlijke gevolg van de goede lineariteit van het instrument − de determinatiecoëfficiënt (R2) voor een lineaire pasvorm van de kalibratiegegevens over het algemeen meer dan 99,95% . Als kalibratiegegevens niet voldoende lineair zijn, waarschuwt de software automatisch de gebruiker, die vervolgens in staat is om de gegevensset op te lossen via het Boolean-besturingselement dat is genoteerd in Protocol stap 4,2; de gebruiker kan vervolgens de kalibratie-gegevensset wijzigen door indien nodig vervangende gegevens te verkrijgen.
Beschouwd als onzekerheden | Volledige (6-punts) kalibratie |
NDIR bias | ± 2,61% |
+ Sacharoseoplossing | ± 2,61% |
+ Pipet | ± 2,78% |
Nominaal resultaat [μgC] | 18,49 |
Tabel 3: onzekerheden in massa kalibratie constante. Representatieve bijdrage van onzekerheden bij de kalibratie van het OCEC-instrument voor de voorbeeld-kalibratie met zes punten (Zie tabel 1). De totale kalibratie onzekerheid wordt gedomineerd door bias in de NDIR-detector met fout als gevolg van de nauwkeurigheid van het pipetteren van de sacharoseoplossing, met inbegrip van menselijke fouten (reproduceerbaarheid tussen gebruikers en herhaalnauwkeurigheid binnen de gebruiker), waarbij de tweede belangrijkste is, gevolgd door onzekerheid in de lineaire regressie en de sacharose concentratie (die beide verwaarloosbaar zijn).
Kalibratie van het OCEC-instrument is een tijdrovende procedure, waarbij doorgaans 2 − 3 uur nodig is om te voltooien, afhankelijk van de lengte van het gebruikte thermische protocol. Een snellere kalibratieprocedure is wenselijk. Hiertoe werd de werkzaamheid van een gewijzigd, afgeknotte kalibratieprotocol geanalyseerd met de gepresenteerde software tool. De ontwikkelde MC-procedure werd uitgevoerd met behulp van alle mogelijke subsets van de voorbeeld kalibratiegegevens vermeld in tabel 1 — beperkt tot gevallen met drie of meer gegevens en een minimum van één blanco meting. Alle resulterende massa kalibratie constanten uit deze analyse worden uitgezet in Figuur 4 als functie van het aantal gebruikte kalibratiegegevens, waarbij kalibratie constanten zijn genormaliseerd door het volledige kalibratie resultaat (6 punten). Het is niet verrassend dat de onzekerheid in de kalibratie constante toeneemt naarmate de beschikbare kalibratiegegevens afnemen. Belangrijk is echter dat de middelen van alle afgeknotte kalibraties binnen de 2σ CI van het volledige kalibratie resultaat vallen, wat een gevolg is van de hierboven vermelde lineariteit van het instrument. Deze consistentie in het MC-gemiddelde suggereert dat een gewijzigde, snellere kalibratie bestaande uit weinig kalibratiegegevens kan worden gebruikt als een “bumptest”-controle van de kalibratie van het OCEC-instrument. Dat wil doen, als het MC-gemiddelde van een 3-punts kalibratie gegevensset binnen de 2σ CI van de bestaande kalibratie ligt, is het waarschijnlijk dat het OCEC-instrument geen herkalibratie vereist. In Figuur 4 wordt ook duidelijk dat de kalibratie onzekerheid met meer kalibratiegegevens wordt gereduceerd, maar de afname van de onzekerheid lijdt aan afnemende rendementen. Verwijzend naar tabel 3 en de bovenstaande bespreking, aangezien kalibratie onzekerheden worden gedomineerd door NDIR-bias (gekwantificeerd met de standaardfout van de “kalibratie gebieden”), de marginale afname van de kalibratie onzekerheid door de nth gegevenspunt kan worden geschat met de factor √ (1 − 1/n). Bijgevolg is in het representatieve voorbeeld de marginale afname van de onzekerheid groter bij het verplaatsen van een drie-naar-vier-punt kalibratie dan bij het verplaatsen van een vijf-naar een zespunts kalibratie. De ontwikkelde software tool, die kan worden uitgevoerd na verwerving van elk kalibratiegegevens punt (d.w.z. na elke herhaling van protocol stap 3,3), stelt de gebruiker in staat om de kalibratie onzekerheid te kwantificeren gedurende de gehele overname van gegevens. Kritisch kan deze mogelijkheid de gebruiker niet alleen beslissen over een adequate convergentie van de kalibratie in het kader van zijn onzekerheid, maar ook de aanwezigheid van valse gegevens detecteren — dat wil doen, een incrementele afname van de kalibratie onzekerheid die met name anders dan verwachte hoogtepunten voor de gebruiker dat de meest recent verworven kalibratiegegevens punt kan worden gebrekkig.
Figuur 4: beoordeling van de steekproefgrootte bij kalibratie onzekerheid. Berekende massa kalibratie constante voor alle mogelijke combinaties van kalibratiegegevens die worden vermeld in tabel 1 (waarbij een minimum van drie gegevens, inclusief ten minste één blanco, is genormaliseerd door het resultaat met behulp van alle zes gegevens. De relatieve onzekerheid in de kalibratie constante neemt af met een toename van het aantal kalibratiegegevens. Het blauwe gearceerde gebied in de figuur komt overeen met de 2σ CI van de kalibratie constante, berekend op basis van alle kalibratiegegevens. Het is duidelijk dat alle nominale resultaten binnen dit CI zijn, wat suggereert dat – hoewel onzeker – een afgeknotte kalibratieprocedure van slechts drie kalibratiegegevens punten kan worden gebruikt als een “bumptest”-controle van de kalibratie van het instrument. Klik hier om een grotere versie van dit cijfer te bekijken.
Berekende koolstof massa’s en onzekerheden voor vier voorbeeld gegevenssets worden gedetailleerd beschreven in tabel 4; deze gegevens zijn afkomstig van metingen van koolstofhoudende emissies uit koolwaterstof vlammen27,29, gasturbines29, en Fine-mode (< 2 μm) koolzuurhoudende deeltjes verkregen uit sub-Ocean sedimentmonsters30. Berekende (MC-gemiddelde) OC-, EC-en TC-massa’s worden weergegeven in tabel 4 naast de berekende kritieke afname daling en de EC/TC-ratio voor elke gegevensset, waarbij de breedte van de voorbeeldgegevens wordt weergegeven in de context van de bemonsterde koolstof samenstelling. Samengevat in tabel 4 zijn de bronnen en algemene onzekerheden in koolstof massa met behulp van de gepresenteerde software tool in vergelijking met die gerapporteerd door het instrument. Binnen de software tool levert de combinatie (productdistributie) van kalibratie onzekerheid en herhaalbaarheid de totale onzekerheid op in MC-berekende TC-massa (gekwantificeerd als − 8.32/+ 8.40% in het huidige werk), die onafhankelijk is van onzekerheid in de splitspunt en fungeert daarom als een ondergrens in de onzekerheid van OC en EG-massa’s. Deze representatieve kalibratie onzekerheden worden toegepast op elke voorbeeldgegevensset met behoud van de nominale massa kalibratie constante die in de oorspronkelijke analyses wordt gebruikt.
Tabel 4: onzekerheden in gegevensanalyse. Bijdrage tot onzekerheid in de door de ocec gemeten koolstof massa’s voor vier voorbeeld gegevenssets uit een groot aantal bronnen en uitgevoerd door verschillende laboratoria27,29,30. (a) numerieke resultaten van de voorbeeldgegevenssets: OC-, EC-en TC-massa’s, kritische afname daling voor de kwantificering van een gesplitste-punt onzekerheid en een elementair-tot-totaal koolstof verhouding. b) een samenvatting van de onzekerheden in de berekende koolstof massa’s. Bijdragende bronnen van onzekerheid omvatten die van de constante van de massa kalibratie, de reproduceerbaarheid van de kalibratieprocedure en de onzekerheid in het splitspunt (ten opzichte van de TC-massa) die overeenkomt met de in (a) vermelde kritische verzwakking van de verzwakking. Onzekerheden in de koolstof massa’s met behulp van de vaste vergelijkingen (EQ. (1)) die door het instrument worden gebruikt, worden ook weergegeven in (b). De rode en gele markering van de gegevens komt overeen met onder-en overschatting van de onzekerheid, respectievelijk bij het gebruik van de vaste vergelijkingen ten opzichte van de huidige methode. In de meeste gevallen kan het instrument de onzekerheids onzekerheid van de koolstof massa echter, indien gemeten OC of EG-massa klein is, mag het instrument de onzekerheid in vergelijking met de huidige software te zeer inschatten. Klik hier om dit bestand te downloaden.
Voor deze voorbeelden werd de techniek van de verzwakking van het gesplitste punt gebruikt om de onzekerheid te kwantificeren. Ocecgo-berekende kritische waarden van afname daling varieerden van 1,342% tot 2,059% resulterend in splitspunt onzekerheden van 0,10% tot 4,50% van TC Mass. Hoewel de toegepaste waarden van de verzwakking van de afname inderdaad enigszins subjectief zijn-met name de factor in loondienst om de onzekerheid in het splitspunt te schatten als gevolg van die van optische eigenschappen-deze voorbeelden wijzen op de afhankelijkheid van gespleten punt onzekerheid op de specifieke analysegegevens. Bijvoorbeeld, splitpunt onzekerheid is gevoelig voor de helling van de AVEC plot in de nabijheid van het nominale splitspunt. Overweeg de gegevensset ‘ A ‘ die overeenkomt met de voorbeeldgegevens in afbeelding 3 en gegevensset ‘ D ‘; Ondanks het feit dat vergelijkbare kritische demping daalt, geven de relatief ondiepe en steile hellingen van hun respectieve AVEC plots (bijv. Zie Figuur 3b voor dataset “A”) de grootste en kleinste split point onzekerheden van 4,50% en 0,10% van TC Mass. Bovendien blijkt uit de voorbeeldgegevens dat de invloed van onzekerheid in het splitspunt grotendeels afhangt van de schaal van de nominale OC-en EG-massa’s. Overweeg voorbeelden van gegevenssets “B” en “C”, die bijna identieke splitspunt onzekerheden hebben (≈ 1,22% van TC Mass); gegevensset “C” bevat ≈ 43% OC terwijl de gegevensset “B” ≈ 8% bevat; de lagere relatieve hoeveelheid OC in de laatste resulteert in een bijna verdubbeling van overklokmassa onzekerheid. Kritisch wijzen deze resultaten op de eis om een gesplitste punt onzekerheid direct in de context van de analyse ‘ AVEC data en algemene koolstof massa’s ‘ te beschouwen.
De door de fabrikant gerapporteerde onzekerheden in koolstof massa’s worden ook weergegeven in tabel 4. Deze schattingen houden niet direct rekening met onzekerheden in het kalibratie-en splitspunt, maar worden berekend aan de hand van de vaste relaties die worden weergegeven in EQ. (1)26, waarbij mi de nominale massa van de specifieke koolstof component vertegenwoordigt.
1
Deze vaste relaties maken het mogelijk dat geschatte onzekerheden in OC en/of de EG-massa kunstmatig minder zijn dan die van TC Mass-deze aandoening treedt op wanneer OC of EC-massa minder dan een derde van TC Mass is, zoals het geval is voor gegevenssets “A”, “B” en “D”. Dit resultaat is niet-fysiek, aangezien de relatieve onzekerheid in OC en de EG-massa’s hieronder moet worden begrensd door de die van TC Mass, een gevolg van de vermeerdering van een gespleten punt onzekerheid in de berekende OC-en EG-massa’s. Cellen die rood en geel in de tabel worden gemarkeerd, komen overeen met de onder-en overschattingen van de koolstof massa-onzekerheid bij gebruik van EQ. (1) van de fabrikant. De vaste vergelijkingen onderschatting TC-massa onzekerheid voor alle vier de voorbeelden, een gevolg van berekende TC-massa’s die voldoende groot zijn. In de meeste gevallen, de vaste vergelijkingen ook onderschatting EC en OC massa onzekerheid, behalve waar OC (dataset “B”) en EC (dataset “D”) waren voldoende klein om te overschatting via EQ. (1). Deze asymptotische toename van onzekerheid via EQ. (1) is het eens met de huidige software, aangezien onzekerheid in OC en EG-massa’s als gevolg van die van het splitspunt afhankelijk zijn van hun absolute grootte; onzekerheden via de vaste vergelijkingen kunnen echter niet worden bijgehouden met die van de huidige software, die in de context van specifieke analysegegevens rechtstreeks rekening houdt met de onzekerheid van de Splits punten en deze verspreidt.
Het gebruik van een MC-Framework in de gepresenteerde software tool is noodzakelijk om de component onzekerheden nauwkeurig te propageren via de niet-lineaire algoritmen van thermische optische OC/EC-analyses. Het is belangrijk op te merken echter, door hun inherent willekeurige aard, dat MC methoden niet deterministisch zijn en de neiging om inconsistente resultaten opleveren als het werk aantal van MC tekent/herhalingen (zie protocol stappen 4.3.3 en 6.5.2) is onvoldoende ― vergelijkbaar met een statistisch ondermaatse monster. Daarom is er een inherente consistentie versus deze berekening tijd trade-off worden overwogen bij het verwerken van gegevens met behulp van ocecgo. Het wordt dus aanbevolen voor gebruikers om voorlopige verwerking en probleemoplossing van de gegevens uit te voeren met behulp van een klein aantal MC-trekkingen (bijv. 104). Zodra de berekeningen bevredigende resultaten opleveren, moet de gebruiker vervolgens het aantal MC-trekkingen (naar 106-108) verhogen om een resultaat te opleveren dat minder wordt beïnvloed door de discrete en willekeurige aard van de MC-methode. Naast de noodzakelijke “operationele definitie” van thermische optische OC/EC-analyses, zijn er andere beperkingen in de verwerking van deze gegevens die ook moeten worden herkend bij het gebruik van Ocecgo voor het berekenen en rapporteren van OC/EC-gegevens. Ten eerste, NDIR-gebaseerde instrumenten (zoals die vermeld in de tafel van de materialen) lijden aan drift in het NDIR-signaal dat moet worden gecorrigeerd. In dit Protocol (zie stap 6.3.2 en de Ocecgo -documentatie) kan een nieuwe benadering voor het corrigeren van drift in de NDIR-detector optioneel door de gebruiker worden gebruikt. Hoewel dit in de ervaring van de auteurs verbeterde resultaten oplevert ten opzichte van de standaard lineaire NDIR-correctie van de fabrikant, moet worden opgemerkt dat onzekerheid in deze NDIR-correctie uitdagend is als het niet onmogelijk is om te kwantificeren en blijft dus een onadministratief onderdeel van de onzekerheid bij de berekening van koolstof massa’s. In dezelfde geest is het ook uitdagend om onzekerheid te kwantificeren in de noodzakelijke veronderstelling dat PC en EG optische eigenschappen delen. Als deze optie is geselecteerd (zie protocol stap 6.4.1), probeert de kritieke dempings methode het effect van deze veronderstelling door middel van een subjectieve expansie factor conservatief te binden. Belangrijk echter, dit is noodzakelijkerwijs slechts een schatting, en de gebruiker wordt voorgesteld om te beoordelen van het effect van deze expansie factor (dat wil zeggen, de kritische verzwakking daling) op hun specifieke gegevens. Ocecgo wordt geleverd als een open-source tool, zodat het gemakkelijk kan worden uitgebreid door de auteurs en andere geïnteresseerde medewerkers om niet alleen met andere instrumenten samen te vallen, maar ook andere nuttige, Veldspecifieke functionaliteiten te omvatten. In het algemeen is de ontwikkelde open source software tool in combinatie met de gedetailleerde kalibratieprocedure bedoeld om te helpen bij het bereiken van nauwkeurigere metingen van OC, EC en TC Mass in aërosolmonsters, terwijl de robuuste berekening van de meting Onzekerheden.
The authors have nothing to disclose.
Dit werk werd gesteund door de natuurwetenschappen en ingenieurs Onderzoekraad van Canada (NSERC) FlareNet strategisch netwerk (Grant # 479641), NSERC Discovery Research subsidies (Grant #06632 en 522658), en natuurlijke hulpbronnen Canada (project manager, Michael Layer). De auteurs zijn dankbaar voor degenen die RAW-gegevensbestanden hebben gedeeld voor gebruik als representatieve voorbeelden in dit werk.
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 |