International Comparison CCQM-K41. Final Report. Page 1 of 55

Similar documents
Development of DMS and Acetonitrile, and Formaldehyde gas standards. Gwi Suk Heo, Yong Doo Kim, Mi-Eon Kim, and Hyunjin Jin

Final Report of APMP.QM-K46 Ammonia in Nitrogen at 30 µmol/mol Level

COOMET TECHNICAL COMMITTEE 1.8 PHYSICAL CHEMISTRY

Comparison APMP.QM-S2.1 Oxygen in nitrogen at atmospheric level

International comparison CCQM-K82: Methane in Air at Ambient level ( ) nmol/mol. (Final report)

METROLOGICAL ASSURENCE OF O3, CO2, CH4 AND CO CONTROL IN ATMOSPHERE

COOMET TECHNICAL COMMITTEE 1.8 PHYSICAL CHEMISTRY. CООМЕТ project 576/RU/12

Gases&Technology. Measurement of Impurities in Helium Using the Dielectric Barrier Discharge Helium Ionization Detector. FEATURE.

ENVIRONMENTAL AND POLLUTANTS GAS ANALYZERS

Report on Co-production of CH 4 /air within ACRM framework

Sangil LEE, Mi Eon KIM, Sang Hyub OH, Jin Seog KIM. Center for Gas Analysis Korea Research Institute of Standards and Science (KRISS) Daejeon, Korea

Final Report. International Key Comparison CCQM-K94

International Comparison CCQM-K111.1 Propane in nitrogen

Laboratory Hardware. Custom Gas Chromatography Solutions WASSON - ECE INSTRUMENTATION. Engineered Solutions, Guaranteed Results.

Improved Reliability in Natural Gas Energy Measurement with the Revised ISO Standard. Gerard Nieuwenkamp - VSL

Generating Calibration Gas Standards

Laboratory Hardware. Custom Gas Chromatography Solutions WASSON - ECE INSTRUMENTATION. Custom solutions for your analytical needs.

APPLICATION NOTE. GC Integrated Permeation Device

Gas mixtures. Individual solutions specifically for your application

Date of Shipment: Xxxxx 00, 20xx 3402c National Institute of Advanced Industrial Science and Technology National Metrology Institute of

International Comparison CCQM-K119 Liquefied Petroleum Gas

METHOD 3C - DETERMINATION OF CARBON DIOXIDE, METHANE, NITROGEN, AND OXYGEN FROM STATIONARY SOURCES

This test shall be carried out on all vehicles equipped with open type traction batteries.

International comparison Refinery gas (CCQM-K77) Final Report

Gravimetric preparation of NO 2 primary reference gas mixtures at the NMISA. Tshepiso Mphamo 03 September 2012

Title: Standard Operating Procedure for Measurement of Ethylene (C 2 H 4 ) in Ambient Air by Reduced Gas Detection (RGD)

LINEAR TRANSFORMATION APPLIED TO THE CALIBRATION OF ANALYTES IN VARIOUS MATRICES USING A TOTAL HYDROCARBON (THC) ANALYZER

Application Note AN-107

VIC offers a variety of Calibrated Gas Leaks

Carrier Gases in Capillary GC

White Paper. Chemical Sensor vs NDIR - Overview: NDIR Technology:

Research and Development. Money in the Pipeline

ASTM WK Standard Test Method for Dissolved Gases. Anne Jurek Applications Chemist

GCMSD-Headspace Analysis SOP

DETERMINATION OF TETRAHYDROTHIOPHENE IN AMBIENT AIR BY GAS CHROMATOGRAPHY WITH A PFPD DETECTOR COUPLED TO A PRECONCENTRATION TECHNOLOGY

International Comparison CCQM K53 Oxygen in Nitrogen

Automated Determination of Dissolved Gases in Water Anne Jurek. Abstract: Discussion:

Flare Gas Composition Analysis and QA/QC Lessons Learned and Lessons Lost SPECTRUM ENVIRONMENTAL SOLUTIONS, LLC 1

The analysis of complex multicomponent

Gas Mixture Two Components. Gas Mixtures Two Components

UNITY 2 TM. Air Server Series 2 Operators Manual. Version 1.0. February 2008

International Comparison CCQM-K111 Propane in nitrogen

Determination of Air Density with Buoyancy Artefacts

International comparison CCQM-K84. Carbon monoxide in Synthetic air at ambient level

Gas Metrology. SI-traceable: Primary Reference Materials Certified Reference Materials Calibrated Gas Mixtures Interlaboratory Comparisons

Detector Carrier Gas Comments Detector anode purge or reference gas. Electron Capture Nitrogen Maximum sensitivity Nitrogen Argon/Methane

You should be able to: Describe Equipment Barometer Manometer. 5.1 Pressure Read and outline 5.1 Define Barometer

Three Columns Gas Chromatograph Analysis Using Correlation between Component's Molecular Weight and Its Response Factor

Commercial Practice Test Method Internal Vapor Analysis of Hermetic Devices

APPLICATION NOTE. Fast Analysis of Coal Mine Gas Using the INFICON 3000 Micro GC ABSTRACT

OMCL Network of the Council of Europe QUALITY ASSURANCE DOCUMENT

Multiple Gas#5 GC configuration Jan 2016

Optimizing Vial Pressurization Parameters for the Analysis of <USP 467> Residual Solvents Using the 7697A Headspace Sampler

Predicted Dispense Volume vs. Gravimetric Measurement for the MICROLAB 600. November 2010

Ultimate Performance in Gas Chromatography

FPG8601 Force Balanced Piston Gauge

Speciality Gases Practice 4. Our quality, your safety. Accredited Test and Calibration Laboratory Speciality Gases Centre in Hörstel.

High Automation of Thermo Scientific FlashSmart CHNS/O Analyzer using the MultiValve Control (MVC) Module

Supporting Information for Micro-Collection of. Gases in a Capillary Tube: Preservation of Spatial

METHOD 21 - DETERMINATION OF VOLATILE ORGANIC COMPOUND LEAKS. 1.2 Scope. This method is applicable for the

CORESTA RECOMMENDED METHOD N 6

International Comparison CCQM-K76: Sulfur Dioxide in Nitrogen

CHAPTER 16 %UHDWKLQJ*DV0L[LQJ3URFHGXUHV

Demystifying the Process of Calibrating your Thermal Desorption Gas Chromatography System using Compressed Gas Standards.

Carrier Gases in Capillary GC

Advancements in Gas Chromatography Analyzers - Keeping up with New Technology. Chuck Runkle Gas Phase Product Specialist ASTS June 2013

HiQ laboratory gas generators.

ISO INTERNATIONAL STANDARD

Rapid and Reliable Detection of Dissolved Gases in Water

Gas flow calibrations performed at the National Metrology Institute of South Africa (NMISA)

Specific Accreditation Criteria Calibration ISO IEC Annex. Mass and related quantities

Appendix D: SOP of INNOVA 1412 Photoacoustic Multi-Gas Monitor. Description and Principle of Operation

International Key Comparison CCQM-K26b and Pilot Study CCQM P50b (SO 2 )

METHOD 25A - DETERMINATION OF TOTAL GASEOUS ORGANIC CONCENTRATION USING A FLAME IONIZATION ANALYZER

Thermo Scientific Model 146i Principle of Operation

MASS FLOW SYSTEMS MASS FLOW MEASURING, CONTROLLING AND BLENDING SYSTEMS

29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

Technical Data Sheet MF010-O-LC

(NIST Special Publication )

Retention Time Locking: Concepts and Applications. Application

METHOD 204F--VOLATILE ORGANIC COMPOUNDS CONTENT IN LIQUID INPUT STREAM (DISTILLATION APPROACH) 1.1 Applicability. This procedure is applicable for

Title: Standard Operating Procedure for R&R Environmental Devices Model MFC201 Gas Dilution Calibrator

EN ISO/IEC Technical Requirements

Title: Standard Operating Procedure for Dasibi Model 5008 Gas Dilution Calibrator

PERFORMANCE SPECIFICATION PROPELLANT, HYDROGEN

CORESTA RECOMMENDED METHOD Nº 67

AC : MEASUREMENT OF HYDROGEN IN HELIUM FLOW

Gases and welding distributors continue

Major Design Features

DF-310E PROCESS ANALYSERS APPLICATIONS FEATURES

A New GC/MS System Designed for Helium Carrier Gas Conservation

Operation of the Perkin Elmer TGA-GC/MS

NMISA New Gas Flow Calibration Facility

Zebron Gas Management

The Preparation of Low Concentration Hydrogen Sulfide Standards Gulf Coast Conference September 12, 2002 Galveston, TX Paper 050,12:45 pm Moody F

Gas Clean. Filters. Delivering Clean Gases for GC and GC/MS Operation

Unit 2 Kinetic Theory, Heat, and Thermodynamics: 2.A.1 Problems Temperature and Heat Sections of your book.

CHEMISTRY - CLUTCH CH.5 - GASES.

States of Matter. Q 7. Calculate the average of kinetic energy, in joules of the molecules in 8.0 g of methane at 27 o C. (IIT JEE Marks)

Applications of a Magnetic Sector Process Mass Spectrometer to the Analysis of Variable Vacuum Samples

Transcription:

International Comparison CCQM-K41 Final Report Franklin R Guenther 1, Walter R Miller 1, David L. Duewer 1, Gwi Suk Heo 2, Yong-Doo Kim 2, Adriaan M.H. van der Veen 3, Leonid Konopelko 4, Yury Kustikov 4, Nina Shor 4, Chris Brookes 5, Martin Milton 5, Florbela Dias 6, Han Qiao 7 1 National Institute of Standards and Technology (NIST), Chemical Science and Technology Laboratory, 100 Bureau Drive Stop 8393, Gaithersburg, MD USA 2 Korea Research Institute of Standards and Science (KRISS), Division of Chemical Metrology and Materials Evaluation, P.O. Box 102, Yusong, Taejon, Republic of Korea 3 NMi Van Swinden Laboratorium B.V. (NMi VSL), Schoemakerstraat 97, 2628 VK Delft, the Netherlands 4 D.I.Mendeleyev Institute for Metrology (VNIIM), Research Department fot the State Measurement Standards in the Field of Physico-Chemical Measurements, 19, Moskovsky Prospekt, 190005 St-Petersburg, Russia 5 National Physical Laboratory (NPL), Teddington. Middlesex, TW11 OLW, UK 6 Instituto Português da Qualidade, IPQ, Monte da Caparica, Portugal 7 National Research Center for Certified Reference Materials (NRCCRM), No. 7, District 11, Heping Street, 100013 Beijing, P.R. Chine Field Amount of Substance: Gas Standards Subject Comparison of gas standards containing 10 μmol/mol hydrogen sulfide in a balance of nitrogen. Participants Table 1: list of participants in this key comparison. Acronym Country Institute KRISS KR Korea Research Institute of Standards and Science (KRISS), Division of Chemical Metrology and Materials Evaluation, P.O. Box 102, Yusong, Taejon, Republic of Korea IPQ PT Instituto Português da Qualidade, IPQ, Monte da Caparica, Portugal NIST US National Institute of Standards and Technology (NIST), Chemical Science and Technology Laboratory, 100 Bureau Drive Stop 8393, Gaithersburg, MD USA NMi-VSL NL NMi Van Swinden Laboratorium B.V., Delft, the Netherlands NPL UK National Physical Laboratory (NPL), Teddington. Middlesex, TW11 OLW, UK NRCCRM CN National Research Center for Reference Materials, Beijing, P.R. China VNIIM RU D.I.Mendeleyev Institute for Metrology (VNIIM), Research Department for the State Measurement Standards in the Field of Physico-Chemical Measurements, 19, Moskovsky Prospekt, 190005 St-Petersburg, Russia Page 1 of 55

Organizing body: CCQM Gas Analysis Working Group Background This key comparison was intended to compare the capabilities for the preparation and value assignment of gas standards for hydrogen sulfide in nitrogen (subsequently referred to as H 2 S), maintained at the participating national metrology institutes. The range of the nominal amount-ofsubstance fractions of the comparison standard is 10 μmol/mol, which is close to regulatory levels in most countries. Conduct of the Comparison The National Institute of Standards and Technology (NIST) procured ten 6-liter gas mixtures of 10 μmol/mol H 2 S in nitrogen from Scott Marin Gases, Riverside CA, USA. These gas standards were monitored for 6 months prior to the start of the key comparison exercise in order to assure stability. Each standard was analyzed and referenced to a nominal 10 μmol/mol NIST stable reference cylinder (SRC), which had been observed to be stable of many years. Prior to the start of the comparison, seven of the ten standards were selected for distribution to the participants based on stability over the 6-month evaluation period. A comparison to the NIST stable reference cylinder was performed prior to shipment of the cylinders to the participants. Each participating laboratory was shipped one 6-liter cylinder. This cylinder was to be analyzed by the receiving laboratory and then returned to NIST. NIST then reanalyzed the cylinder contents, again referenced to the SRC, to assure that the compounds remained stable throughout the comparison. The comparison was officially declared closed on April 1, 2005. All participating laboratories had submitted complete reports by this date. The report from NRCCRM had originally been reported as not received; however review of email records show that the report was indeed received prior to April 1, 2005. As there is no gravimetric value for these comparison cylinders, a consensus reference value is calculated from the results received from the participants. NIST comparisons to the stable reference cylinder are used to normalize the participants data to remove cylinder to cylinder concentration differences. In effect the Key Comparison reference value is the calculated concentration of the SRC. NIST Instrumentation Non-dispersive Ultra-Violet Spectrometer (NDUV): This instrument (Western Research Series 900 Analyzer, model 922) was operated on the 0 25 mol/mol setting with a span concentration of 10.24 mol/mol. The internal pump was adjusted to deliver a constant flow of 500 ml/min. An automated stream selection device was used to randomly select the flow from one of five cylinders. The regulator on each cylinder was adjusted to give a flow of approximately 600 ml/min. This flow was swept past the inlet of the instrument, where 500 ml/min was directed into the instrument port and 100 ml/min exhausted to vent. The sample line was purged for a minimum of 3 minutes. The analog signal (10 V output) from the instrument was averaged over two minutes for each reading. Figure 1 shows the reproducibility achieved using this instrument. A series of measurements of the SRC was done over a time period of one day. Figure 1 shows the response from 90 measurements from the NDUV. From this data it was determined that the combined uncertainty Page 2 of 55

(k=1) of the ratio was 0.25 % relative. The uncertainty bars for 0.25 % relative are shown in Figure 1 with red lines. This uncertainty was applied to all ratios obtained using this instrument. Figure 1: NDUV Reproducability Signal (V) 2.055 2.050 2.045 2.040 2.035 2.030 2.025 2.020 2.015 2.010 2.005 0 10 20 30 40 50 60 70 80 90 Run Electrochemical analyzer: This instrument (Interscan, model RM-17) was operated on the 0 50 mol/mol range with the span potentiometer set to 2.55. This was done to adjust the instrument sensitivity such that the 10 mol/mol samples gave a response at 60 % of range. The internal sample pump was adjusted to deliver a constant flow rate of 100 ml/min. An automated stream selection device was used to randomly select the flow from one of five cylinders. The regulator on each cylinder was adjusted to give a flow of approximately 200 ml/min. This flow was swept past the inlet of the instrument, where 100 ml/min was directed into the instrument port and 100 ml/min exhausted to vent. The sample line was purged a minimum of 3 minutes. The analog output (10 mv) was averaged over two minutes for each reading. To measure the reproducibility of the Interscan, a series of measurements of the SRC was completed, and is show in Figure 2. This instrument does exhibit short term drift that requires correction. The sampling protocol used by NIST removes this short term drift by use of a control sample. For this project the control was the SRC. Figure 3 shows the method used to correct for this drift. In this figure five of the comparison cylinders are shown along with the SRC (labeled Control). As can be seen significant drift is being exhibited by the Interscan instrument. However, when the ratio of the sampled cylinder to the control is calculated, the drift is minimized and very good reproducibility is achieved (Table 2). A conservative estimate of 0.25 % relative was used for all ratios measured using the Interscan instrument. Figure 2: Interscan Reproducibility 0.022400 0.022350 Response (mv) 0.022300 0.022250 0.022200 0.022150 0.022100 0.022050 0 10 20 30 40 50 60 70 80 90 Run Page 3 of 55

Figure 3: Interscan Reproducibility data 0.023600 0.023400 0.023200 Response (mv) 0.023000 0.022800 Control WS-10A FA02278 FA02282 FA02277 FF21074 0.022600 0.022400 0.022200 0 5 10 15 20 25 30 35 Run Table 2: Interscan Reproducibility Data Ratio Ratio Ratio Std Sample / Sample / Sample / Average Dev % RSD Sample WS-10 WS-10 WS-10 FA02277 1.0127 1.0133 1.0134 1.01314 0.00040 0.04% FA02278 1.0177 1.0186 1.0191 1.01845 0.00068 0.07% FA02282 1.0367 1.0376 1.0379 1.03742 0.00061 0.06% FF21074 0.9923 0.9939 0.9942 0.99348 0.00105 0.11% WS-10A 1.0300 1.0296 1.0298 1.02981 0.00020 0.02% To test the linearity of the Interscan, a permeation system was adjusted to give concentration between 8 and 12 mol/mol hydrogen sulfide in nitrogen. In figure 4, the data shows that the Interscan is linear in this concentration range. Figure 4: Interscan Linearity 1.10000 1.05000 y = 0.100332x + 0.003184 R 2 = 0.999961 Response Ratio 1.00000 0.95000 0.90000 0.85000 0.80000 8 8.5 9 9.5 10 10.5 11 Concentration ( μ mol/mol) NIST Stability Measurements NIST used a non-dispersive ultra-violet spectrometer (NDUV) and an electrochemical analyzer (Interscan) to measure the hydrogen sulfide in each cylinder. The response from the sampled Page 4 of 55

cylinder (r S ) and the response from the SRC (r R ) were measured and the response ratio calculated according to the following formula: r S r i = (1) rr This was done in order to remove any short term drift the instrument may exhibit and to reference each cylinder to the SRC. A sampling sequence of SRC, sample, SRC was used. The two measurements of the SRC were averaged to calculate r R. The response ratio was measured periodically prior to shipment of the cylinders to the participant in order to monitor stability. The response ratio was again measured once the cylinder was returned from the participants. Graphs of the response ratios over the duration of this study for each cylinder are shown in Appendix A, graphs 1 through 8. The calculated average of the ratios determined by NIST is shown in the graphs along with the expanded uncertainty (k=2). The position of the average along the x-axis is the actual date the participant analyzed the cylinder, as reported by the participant. The identification of the cylinder each NMI received and analyzed along with the average response ratio ( r i ) and associated expanded uncertainty (k=2) measured by NIST is given in Table 3. For each cylinder the average ratio was used in the calculation of the Key Comparison Reference value. No stability effects could be positively determined, as all ratios were within their uncertainty bars. Some minor stability issues may be apparent in the cylinders that traveled to IPQ and NRCCRM. However in all cases the average ratio compensates for any stability issues. An extra cylinder in the study was sent to NMi and returned without analysis at NMi. This cylinder was then reanalyzed at NIST to determine if any effects from travel could be detected. No travel effects were detected. Table 3: Study Cylinders Institute Cylinder Number Average Ratio to SRC ( r i ) KRISS FA02288 0.9986 ± 0.0025 IPQ FA02267 0.9938 ± 0.0025 NIST FF21074 0.9872 ± 0.0022 NMi-VSL FA02295 1.0128 ± 0.0025 NPL FA02283 1.0091 ± 0.0023 NRCCRM FF17172 1.0005 ± 0.0029 VNIIM FA02270 0.9999 ± 0.0025 Travel Cylinder FF02282 1.0357 ± 0.0023 Results All participating laboratories sent in their complete report prior to the April 1, 2005 deadline except VNIIM, which sent only their results table on April 1. A full report from VNIIM was subsequently received on May 13, 2005. Since the results table was received before the deadline, the full report has been accepted. The report from KRISS was received on March 30, 2005, however a subsequent report was received on April 1, 2005 with corrected values. Since this report with corrected values was received before the deadline, it has been accepted as the final report from KRISS. All final reports from the participants are attached as received in Appendix B. A summary of the reported results is given in Table 4, along with the NIST determined average response ratio to the SRC. The instruments and methods used by the participants are summarized in Table 5. Since all the cylinders contained a different concentration of hydrogen sulfide, and a gravimetric value was not available, a consensus value was determined for the Key Comparison Reference value (KCRV). To calculate the KCRV, the reported value from each participant was divided by Page 5 of 55

the NIST determined average response ratio. This calculation gives the concentration of the SRC, as would be determined by the participant. The equation used was: X = i SRC X r i i, (2) X i is the concentration reported by the participant, ri is the average response ratio for where cylinder the participant analyzed, and X i, SRC is the calculated concentration of the SRC. The uncertainty of this concentration u ) was calculated using: ( X i,src ( X i, SRC ) = X i, SRC u( X i) X i 2 u( ri ) + ri u (3) where ( X i ) the average response ratio measured by NIST. u is the participants reported uncertainty divided by k, and r ) 2 u is the uncertainty of Table 4 summarizes the results obtained from the NIST stability analyses and the reported values from the participants. The final two columns in the table lists the calculated value for the stable reference cylinder (SRC) using the participants and NIST data. This data is also graphically represented in Figure 5, with the reported values on the left in red and the calculated values for the SRC on the right. Table 4: Summary of values Cylinder NMI X i Reported Concentration ( mol/mol) U ( X i ) NIST Ratio Determination r i U r ) ( i ( i Calculated SRC Value ( mol/mol) X i, SRC U (,SRC ) (k=2) (k=2) (k=2) FA02288 KRISS 10.012 0.150 0.9986 0.0025 10.03 0.15 FA02267 IPQ 9.95 0.07 0.9938 0.0025 10.01 0.08 FF21074 NIST 9.909 0.044 0.9872 0.0022 10.04 0.05 FA02295 NMi-VSL 10.14 0.06 1.0128 0.0025 10.01 0.06 FA02283 NPL 10.07 0.17 1.0091 0.0023 9.98 0.17 FF17172 NRCCRM 10.5 0.5 1.0005 0.0029 10.49 0.50 FA02270 VNIIM 9.96 0.20 0.9999 0.0025 9.96 0.20 X i Page 6 of 55

Table 5: Instruments and methods NMI Instrument Standards/Calibration Repeat Measurements 2 gravimetrically prepared primary gas Agilent 6890 GC-G2350 5 sets of measurements, each standards, single point calibration. KRISS Atomic Emission Detector, consisting of 6 subsets each with 8 Standard within 0.05 μmol/mol of CP-SIL 5CB 25m column measures of PSM and sample sample IPQ NIST LIMAS11 NDUV Model RM17 electrochemical voltammetric sensor (Interscan) 3 primary gas standards prepared by NMi-VSL (25, 50, 150 mol/mol) and one primary gas standard prepared by IPQ (9.43 μmol/mol) 20 cm permeation tube, contained within thermostatically controlled environment. Manual weightings using Mettler AT261. Flow calibration using Molblocs. 3 sets of measurements, each consisting of an unknown number averaged results. 3 measurement periods, each consisting of 5 repeat cycles between 4 permeation generated concentrations, the sample, and two controls. NMi-VSL Model RM17 electrochemical voltammetric sensor (Interscan) Gravimetric primary standards at 10, 20, 40, 60, 80, and 100 μmol/mol. Quadratic fit using ISO 6143 3 sets of measurements, each consisting of average of 90 measurements of each cylinder. NPL Varian GC with Sulfur Chemiluminescence Detection, RH-1 60m column. Two methods of sequencing samples. 4 gravimetrically prepared primary standards in BOC Spectraseal and Scott Acculife cylinders (9.7 to 10.14 mol/mol) Multiple replicates of single point calibrations against three primary standards. Each measurement run consisted of 12-20 injections at 15 s intervals NRCCRM Agilent 4890 GC with flame photometric detection, GS- GASPRO 60m column 6 gravimetrically prepared primary gas standards. All standard between 9.45 and 9.58 mol/mol 3 sets of measurements, each consisting of 6 concentration determinations. VNIIM Crystal-5000.1 GC unknown detector, DB-1 60m column 2 sets of 3 permeation tubes certified using thermo-gravimetric instrument. Unknown flow calibration with nitrogen dilution. 7 independent measurements, each consisting of 5 sub-measurements of sample to permeation system. Calculation of Key Comparison Reference Value Figure 5a displays the values and their 95% uncertainty intervals for both the reported and SRCadjusted concentrations. The empirical probability density function (PDF) for the SRC-adjusted values, calculated assuming that the X i,src and the U 95 ( X i,src ) specify N(X i,src, U 95 ( X i,src )/2) normal kernel densities, is displayed at the right edge of the Figure. The solid horizontal line represents the median of the PDF (half of the PDF area is above the line, half is below). The dashed horizontal lines enclose 95% of the area of the PDF (2.5% of the area is above the upper dashed line, 2.5% of the area is below the lower dashed line). The remarkable overlap of six of the seven results and the relatively large uncertainty on the less-concordant NCCRM result strongly suggests that the consensus location for these data is about within the region 9.96 mol/mol to 10.03 mol/mol. However, these data raise more general questions regarding the philosophical basis for assigning the consensus Key Comparison Reference Value (KCRV) and its associated 95% confidence interval, U 95 (KCRV), that could be significant for future comparisons. A number of scenarios need to be considered: 1) Ignore the U 95 ( X i,src ), treat as a Normal distribution.. The unambiguously most efficient estimates for the location and dispersion of a putatively representative, normally distributed set of data are the mean and the standard deviation. Under this assumption, the KCRV = 10.075 mol/mol and U 95 (KCRV) = (t 0.05,5 )(0.187)/ 7 = (2.57)(0.187)/2.65 = 0.173 mol/mol. 2) Ignore the U 95 ( X i,src ), identify and exclude statistical outlier X i,src. By inspection, the X NCCRM,SRC is not in the same population as the other six. Excluding this value, the mean and the standard deviation of the remaining six values give a KCRV = 10.005 mol/mol and U 95 (KCRV) = (2.57)(0.029)/2.45 = 0.031 mol/mol. Page 7 of 55

3) Ignore the U 95 (X i,src ), estimate the Normal component of a mixture distribution. There are many different robust estimates for the location and dispersion of a set of values that contain a modest percentage of discordant values. 1 Three of the simplest are: a. The median and the median absolute deviation from the median (MAD), adjusted to estimate the standard deviation when the distribution is truly normal (MADe). 2,3 These estimators give KCRV = 10.012 mol/mol and U 95 (KCRV) = (2.57)(0.038)/ 7 = 0.035 mol/mol. b. Huber s Estimate Number 2, H15 and S15. 3 These estimators give KCRV = 10.014 mol/mol and U 95 (KCRV) = (2.57)(0.042)/ 7 = 0.039 mol/mol. c. The Least Power 1.5, L1 and sl1. 4 These estimators give KCRV = 10.028 mol/mol and U 95 (KCRV) = (2.57)(0.074)/ 7 = 0.068 mol/mol. There is no agreement on which of these (or other) robust estimators is best for any given set of data. Here, it appears that Huber s method gives slightly more influence to X NCCRM,SRC than the median and the MADe but somewhat less than does the Least Power 1.5. 4) Use the U 95 (X i,src ) as weighting functions, exclude statistical outliers. There are many different weighted mean estimates for the location and dispersion, all based on the assumption that there exists a functional relationship between the estimated measurement uncertainty and the bias between the reported and the true value. 1 These methods typically estimate the variability of the location estimate, not of the population of values (and there is little agreement on how best to estimate this variability). 5 The simplest method gives influence to values as the inverse-square of the uncertainties, giving for the six concordant values a KCRV =10.020 mol/mol and a U 95 (KCRV) = (2.57)(0.012) = 0.022 mol/mol. The Mandel-Paule method limits the influence of under-estimated uncertainties; for these data, it gives KCRV =10.019 mol/mol and a U 95 (KCRV) = (2.57)(0.012) = 0.022 mol/mol. 5) Use the U 95 (X i,src ) as weighting functions, do not exclude statistical outliers. The weighted mean estimates are not robust to outliers; however, with this data, the U 95 (X NCCRM,SRC ) is large and so will not be given much influence. Using all seven values, both the inverse-square and Mandel-Paule methods give a KCRV =10.022 mol/mol and a U 95 (KCRV) = (2.57)(0.017) = 0.042 mol/mol. However, it should be realized that should the U 95 (X NCCRM,SRC ) be increasingly under-estimated, the value of the KCRV would be increasingly pulled towards the X NCCRM,SRC value = 10.49 mol/mol and the U 95 (KCRV) towards (and beyond) the (uniformly-weighted) standard deviation, 0.173 mol/mol. 6) Use the U 95 (X i,src ) as kernel densities, estimate the Normal component of a mixture distribution.. Assuming that the measurements can be validly represented as N(U 95 (X i,src ),U 95 (X i,src )/2) normal kernels, location and dispersion estimates can be 1 D.L. Duewer. A Robust Approach for the Determination of CCQM Key Comparison Reference Values and Uncertainties. Working document CCQM/04-15, BIPM, 2004. www.bipm.fr/cc/ccqm/allowed/10/ccqm04-15.pdf 2 Analytical Methods Committee of the Royal Society of Chemistry, Technical Brief 6 (April 2001), Robust statistics: a method of coping with outliers, www.rsc.org/membership/networking/interestgroups/analytical/amc/technicalbriefs.asp 3 S. Ellison, Robust Statistics Toolkit (RobStat.xla) Excel add-in, www.rsc.org/membership/networking/interestgroups/analytical/amc/software/robuststatistics.asp 4 F. Pennecchi and L. Callegaro. Between the mean and the median: the Lp estimator. Metrologia 2006;43:213-219. 5 N.F. Zhang, The uncertainty associated with the weighted mean of measurement data. Metrologia 2006;43:195-204. Page 8 of 55

defined that are simultaneously robust to discordant X i,src and under-estimated U 95 (X i,src ). 1 The PDF-median, shown in Figure 5a, is itself a location estimate, giving a KCRV of 10.024 mol/mol. While the empirical 95% confidence interval for the PDF provides a direct estimate of the dispersion of the population, 6 a more robust estimate is provided by the PDF-interquartile range (spdf-iqr), the interval spanning the central 50% of the PDF area, adjusted to estimate the standard deviation when the distribution is truly normal. 7 The spdf-iqr for these data is 0.072 mol/mol, giving U 95 (KCRV) = (2.57)(0.072)/2.65 = 0.067 mol/mol. However, the PDF-based dispersion estimates include both the variability among the X i,src and the U 95 (X i,src ). The U 95 (KCRV) estimates provided in the scenarios above only estimate the variability among the X i,src. To make the two types of estimate concordant, either the non-pdf dispersion estimates need to be expanded to include the U 95 (X i,src ) contribution or the PDF estimates need to be reduced. Given that the median U 95 (X i,src ) for these data is 0.079 mol/mol (that is, larger than the spdf-iqr), it is clear that the U 95 (X i,src ) dominate the total uncertainty and reduction of the PDF-based estimates does not provide useful estimates of just the among-x i,src. The appropriate comparison would then be to the among-x i,src estimates of dispersion augmented with 2 2 the median U 95 (X i,src ): s + u ( X ) s total among i, SRC =. Note that all U 95 (KCRV) calculated from such augmented values will be somewhat larger than that calculated from the spdf- IQR. As a general principle, it is inappropriate to exclude particular data from a consensus estimate unless there is a defensible technical reason for doubting that it is in fact a valid member of the majority distribution of the data. However, it is also generally inappropriate to bias consensus estimate towards extreme values. While weighted estimates are defensible for these particular data since the least-concordant X i,src has an appropriately large U 95 (X i,src ), it seems preferable to use methods expected to perform adequately under more general conditions. Also, given that it is desired to use the information contained in the U 95 (X i,src ) estimates, is seems appropriate to include them in the dispersion calculation. The PDF-median and spdf-iqr fully utilize the dispersion and location information provided in reported U 95 (X i,src ) and are robust to extreme X i,src and to both over- and under-estimated U 95 (X i,src ). The recommended consensus KCRV for the SRC results is therefore the PDFmedian, 10.024 mol/mol hydrogen sulfide in nitrogen. The recommended consensus U 95 (KCRV) is based upon the spdf-iqr, 0.067 mol/mol hydrogen sulfide in nitrogen. This value is used in Table 6 to determine the degrees of equivalence. In Figure 5b the KCRV is plotted (black line) along with the 95% confidence interval for the KCRV (red lines) and the 95% confidence interval for the population (blue lines). The participants reported values are plotted on the left of the figure, and the calculated values for the SRC on the right. Degree of Equivalence The degree of equivalence for each participating laboratory was calculated with the equation: D i = X, z (6) i SRC c The uncertainty of the degree of equivalence was calculated using the equation: 6 P. Ciarlini, M.G. Cox, F. Pavese. G. Regoliosi. The use of a mixture of probability distributions in temperature interlaboratory comparisons. Metrologia, 2004, 41, 116. 7 Stuart, A., Ord, J.K. Kendall s Advanced Theory of Statistics. Volume 1. Distribution Theory. Sixth Edition, Chapter 10. Edward Arnold, London, 1994. Page 9 of 55

U + 2 2 ( Di ) 2 u( X i, SRC ) u( zc ) = (7) The degree of equivalence and uncertainty is summarized in Table 6, and graphically represented in Figure 6. The KCRV lies within the uncertainty bounds of all the participants value for the SRC. For all participants, except for NRCCRM, the agreement is excellent for this reactive compound, and demonstrates excellent comparability. While the value and uncertainty from NRCCRM does show that the KCRV falls within the extremes of the uncertainty bounds, the result shows evidence of a bias of approximately 0.5%. Table 6: Degree of Equivalence Calculated SRC Value KCRV Degree of Equivalence Cylinder NMI Value uncertainty Value uncertainty D U(D) FA02270 VNIIM 9.96 0.10 10.024 0.034-0.06 0.21 FA02283 NPL 9.979 0.085 10.024 0.034-0.04 0.18 FA02295 NMi 10.012 0.032 10.024 0.034-0.01 0.09 FA02267 IPQ 10.013 0.037 10.024 0.034-0.01 0.10 FA02288 KRISS 10.026 0.076 10.024 0.034 0.00 0.17 FF21074 NIST 10.038 0.025 10.024 0.034 0.01 0.08 FF17172 NCCRM 10.49 0.25 10.024 0.034 0.47 0.51 Figure 5a: Graph of PDF median and 95% interval mol/mol (VNIIM) (NPL) (NMI/VSL) (IPQ) (KRISS:) (NIST) 9.8 10.0 10.2 10.4 10.6 10.8 11.0 Reported Concentration (NRCCRM) VNIIM NPL NMI/VSL IPQ KRISS NIST NRCCRM SRC-Adjusted Concentration Page 10 of 55

Figure 5b: Graph of Results CCQM-K41: Hydrogen Sulfide 10 mol/mol nominal (VNIIM ) (NPL ) (NMi-VSL ) (IPQ ) (KRISS ) (NIST ) (NRCCRM ) 9.8 10.0 10.2 10.4 10.6 10.8 11.0 Values as Reported VNIIM NPL NMi-VSL IPQ KRISS NIST NRCCRM KCRV Calculated Values Figure 6: Degrees of Equivalence 1.20 CCQM-K41: Degree of Equivalence (10 μmol/mol Hydrogen Sulfide) 1.00 Degree of Equivalence ( μ mol/mol) 0.80 0.60 0.40 0.20 0.00 VNIIM NPL NMi IPQ KRISS NIST NCCRM -0.20-0.40 Laboratory Page 11 of 55

Limits of Claims This Key comparison may be used to assess calibration and measurement capability (CMC claims for hydrogen sulfide in nitrogen or air, at a concentration between 1 mol/mol and 500 mol/mol. This Key comparison may also be used as evidence for CMCs of hydrogen sulfide in methane at the above concentrations, as long as the analytical method and instrumentation are similar to that used for nitrogen balance hydrogen sulfide gas mixtures. Coordinator Franklin R. Guenther National Institute of Standards and Technology (NIST) Gaithersburg, Maryland 20899-8393 USA Project Reference: CCQM-K41 Completion Date: April 2005 Participant Contact List Florbela Dias Instituto Português da Qualidade (IPQ) Monte da Caparica Portugal Han Qiao National Research Center for Certified Reference Materials (NRCCRM) No. 7, District 11 Heping Street, 100013 Beijing, P.R. Chine Dr. Jin Seok Kim Korea Research Institute of Standards and Science (KRISS) Division of Chemistry and Radiation P.O. Box 102 Yusung Taejon, 305-600 Korea Adriaan M.H. van der Veen Nederlands Meetinstituut (NMi) Schoemakerstraat 97 Postbus 654 2600 AR DELFT The Netherlands Dr. Leonid Konopelko D.I. Mendeleyev Institute for Metrology (VNIIM) 19, Moskovsky Prospekt 190005 St. Petersburg Russia Page 12 of 55

Dr. Martin J.T. Milton National Physical Laboratory (NPL) Environmental Standards Section Teddington Middlesex TW11 0LW England Dr. Franklin R. Guenther National Institute of Standards and Technology (NIST) Chemical Science and Technology Laboratory 100 Bureau Drive Gaithersburg, MD 20899-8393 USA Page 13 of 55

Appendix A FA02288 (KRISS) Ratio 1.00300 1.00200 1.00100 1.00000 0.99900 0.99800 0.99700 0.99600 0.99500 0.99400 0.99300 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average FA02267 (IPQ) 1.00000 0.99800 0.99600 Ratio 0.99400 0.99200 0.99000 0.98800 0.98600 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average Page 14 of 55

FF21074 (NIST) 0.99600 0.99400 0.99200 0.99000 Ratio 0.98800 0.98600 0.98400 0.98200 0.98000 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average FA02295 (NMi-VSL) 1.01800 1.01600 1.01400 Ratio 1.01200 1.01000 1.00800 1.00600 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average Page 15 of 55

FA02283 (NPL) 1.01400 1.01300 1.01200 1.01100 1.01000 Ratio 1.00900 1.00800 1.00700 1.00600 1.00500 1.00400 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average FF17172 (NRCCRM) 1.00600 1.00400 1.00200 Ratio 1.00000 0.99800 0.99600 0.99400 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average Page 16 of 55

FA02270 (VNIIM) 1.00600 1.00500 1.00400 1.00300 1.00200 Ratio 1.00100 1.00000 0.99900 0.99800 0.99700 0.99600 0.99500 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average FA02282 (Travel) 1.04200 1.04000 1.03800 Ratio 1.03600 1.03400 1.03200 1.03000 1.02800 April-04 August-04 November-04 February-05 May-05 September-05 Date Stability Data Average Page 17 of 55

Appendix B Participant: KRISS Confirmation form for the received cylinder Laboratory : KRISS Name of the contact person : Gwi Suk Heo e-mail address : heo@kriss.re.kr Date of reception : December, 2004 Cylinder number : FAO2288 Initial inner pressure of Cylinder as received * : 1100 PSIG Other remarks on the conditions of the cylinder (paint, packaging, and other any damages to the cylinder) : NO NOTICIABLE DAMAGE. IT COMES GOOD CONDITION. * Information on an inner pressure is important to confirm that no sample loss occurred in the course of transport. Page 18 of 55

Measurement report: Hydrogen sulfide in nitrogen Laboratory Cylinder number : KRISS : FAO2288 NOMINAL COMPOSITION : 9. 10-6 12. 10-6 mol/mol Measurement Date Result (umol/mol) stand. uncertainty (% relative) number of submeasurements No. 1 05-03-07 9.999 0.08 8 x 4 No. 2 05-03-08 10.005 0.06 8 x 2 No. 3 05-03-10 9.988 0.03 8 x 3 No. 4 05-03-10 10.033 0.15 8 x 3 No. 5 05-03-11 10.035 0.11 8 x 4 Note: Please copy this table as many times as needed for reporting additional measurements Result: with purity correction of H2S (99.6 ± 0.3%) Gas mixture Result (assigned value) Coverage factor Assigned expanded uncertainty (*) Hydrogen sulfide 10.012 umol/mol 2.0 0.150 umol/mol Reference Method: Describe your instrument(s) (principles, make, type, configuration, data collection etc.): HP6890 GC-G2350 AED (Agilent, USA) were used for analysis. Gas switching valve at 100 o C temp was used with 1 ml sample loop. Sample gas was continuously passed through sampling loop during the analysis. Sample flow-rate was controlled using MFC located at front of sampling valve. CP-SIL 5CB(Varian) GC column was used(25 m x 0.53 mm id x 5 um thickness with dimethylsilicon). GC oven temperatue was set at isothermal condition of 35 o C (4 min). Column flow was 5 ml/min and EPC used for constant column flow. High purity Helium was used as carrier gas. Sample gas was split 5:1 at VI inlet before transferring to GC column. AED transferline temp was set 250 o C and cavity temp was 250 o C. High purity Helium was used as AED discharge gas, which was supplied with EPC control. AED cavity pressure was maintained at 1.5 psi. High purity hydrogen used for ADE with EPC controlled pressure at 8.6 psi. High purity oxygen was used for AED with EPC controlled pressure at 17.5 psi. AED signal was acquired by Agilent Chemstation and GC-AED data was processed using Chemstation software. Sulfur emission line of 181 nm was monitored for sulfur analysis by AED. Calibration Standards: Describe your Calibration Standards for the measurements (preparation method, purity analyses, estimated uncertainty etc.): Page 19 of 55

Purity of hydrogen sulfide (from Aldrich) was checked by several measurement techniques. Gas MS was used for impurity analysis of moisture (4390 ppm), helium (<1 ppm)and argon(23 ppm), GC-MS was used for identification of sulfur impurities, isopropylthiol (64 ppm) and diisopropyl sulfide(2.4 ppm), which were quantified by GC-AED. GC-AED also used for analysis of other impurities, nitrogen (54 ppm), hydrogen (2 ppm), oxygen (<1 ppm), carbon dioxide(30 ppm), carbon monoxide(<1 ppm), COS(32 ppm), methane(<1 ppm), and total sulfur(except the two sulfur compounds, 1 ppm). GC-FID was used for analysis of propane (108 ppm), propylene (70 ppm), and i-butane (9.5 ppm). FTIR was used for analysis of moisture (2155 ppm) to confirm the gas MS data. The two results gave different values for the moisture content. Gas MS data of moisture was used for purity calculation of H2S with uncertainty of 0.2%. Other uncertainty contribution to purity measurement of H2S was counted as 0.!%. Therefore, total of 0.3% of uncertainty was assigned to H2S purity result of 99.6%. KRISS primary reference gases(cylinder number, MD2547, MD2611, MD5850, MD5875) with 2% concentrations of H2S had been prepared by gravimetry; MD2547 (1.99346%), MD2611 (1.99659%), MD5850 (1.99828%), MD5875 (2.00293%). The 2% H2S in cylinder MD 2547 and MD 2611 were diluted to 10 ppm at Acculife cylinder which has specially treated inner surface for sulfur compounds. Scott cylinder number D518937 and D518946 were filled with 10 ppm H2S from dilution of MD2547 of 2% H2S. Cylinder number D518935 and D518944 were filled with 10 ppm H2S from dilution of MD2611 cylinder of 2% H2S; D518937 (9.996 ppm), D518946 (10.074 ppm), D518935 (10.061 ppm), D518944 (10.004 ppm). D518935 and D518946 were selected as KRISS reference standard gas for K-41 inter-comparison work. Dilution system for making 10 ppm level H2S standard gases had been check for its cleanness (H2S contamination free) with flushing the system with N2 and collecting and analyzing the flushed N 2 which was collected in 6 L of Silco canister with pressure of 15 psig. The amounts of nitrogen used for balance gas were measured by gravimetry using high precision gas balance (Mettler, Switzerland) with readability of 1 mg and capacity of 10 kg. Luxfer cylinder (Australia) with fine-polished internal surface with stainless steel (SS) valve was used for the preparation of 2% level of H2S gas standards. For 10 ppm level of H2S gas standard, Acculife cylinder from Scott (Netherlands) cylinders were used. High purity nitrogen was used as balance gas (H2S impurities in the high purity nitrogen were checked using GC-AED) and The micro gas balance and gas balance was calibrated by E 2 grade calibration weight before starting measurement. Tare cylinder was used for buoyancy correction of gas weight measurement. Homogeneity of 2% level H2S gas standards has been check by GC-SCD (GC column used, 30 m x 0.32mm x 4 um, 5 ul of injection loop used for gas switching injection valve, 100 o C, sample split 30 : 1 at split injector, oven temp of 40 o C (4 min, isothermal), SCD burner temp 787 o C(pressure 252 psi), reactor temp 800 o C, SCD pressure 7.4 torr, column flow 2 ml/min, sample gas flow-rate of 60mL/min). Homogeneity of 10 ppm level H2S gas standards has been check by GC-AED at the same condition described above. Stability of 2% H2S standard gas was evaluated by comparing other KRISS H2S PRMs which was prepared in 2000 and 1999. New PRMs and old PRMs gave the same AED response factors with difference of 0.2%. Stability of 10 ppm H2S standard gas(prepared 2002. 2. 1) was evaluated by comparing with newly prepared 10 ppm PRM(prepared 2005. 3. 14). Old PRMs and new PRM gave the same AED response factors with difference of 0.2%. Loss of H2S due to adsorption to inner surface of cylinder was evaluated by distributing equal amount of 10 ppm H2S to other empty cylinder, then second cylinder again was distributed to another empty cylinder. The three cylinders were analyzed and compared their AED response factors to check the adsorption loss of H2S. Result showed that 0.015% loss at first distribution, 0.017% loss at second distribution. Instrument Calibration: Describe your Calibration procedure (mathematical model/calibration curve, number and concentrations of standards, measurement sequence, temperature/pressure correction etc.): Page 20 of 55

Two primary reference gas prepared as above was selected as standard for CCQM-K41 work, and one point calibration was used for the analysis since the prepared concentrations are very close to CCQM-K41 sample s concentration. At the same time, GC-AED analysis of H2S gave good linearity with passing origin always. KRISS H2S CRM conc.(umol/mol) k Uexp(umol/mol) Rel U(%) D518935 10.061 2.1 0.097 0.96 D518946 10.074 2.1 0.097 0.96 Analysis of KRISS primary reference gases and CCQM sample were repeated 8 times alternatively, and these are counted as one subset of measurement. Total of 6 subsets of measurement had been repeated which are counted as one set of measurement. Total of 5 sets of measurements had been repeated. Each analytical result of CCQM sample was obtained from the bracketed measurement of KRISS CRM to correct drift of GC-AED response. Lab temperature were kept 24 o C +_ 1 o C and atmospheric pressure were between 1000 ~ 1015 hpa. Since no appreciable change in temp and pressure were noticed, no correction for the temperature and pressure were made for the measurement. Sample Handling: How were the cylinders treated after arrival (stabilized) and how were samples transferred to the instrument (automatic, high pressure, mass-flow controller, dilution etc).: After receiving sample cylinder, box cartoon was opened, and cylinder was stored at room temperature before analysis (for 2 months). The sample cylinder was connected to GC-AED. The sample gas was transferred (with flow-rate of 100 ml/min) to gas sampling valve of GC. MFC was mounted between sample inlet-line and gas sampling valve to control constant flow of sample to the valve. Uncertainty: There are potential sources that influence the uncertainty of the final measurement result. Depending on the equipments, the applied analytical method and the target uncertainty of the final result, they have to be taken into account or can be neglected. Describe in detail how estimates of the uncertainty components were obtained and how they were combined to calculate the overall uncertainty: In support of this action, a list of potential uncertainty sources is given. This list may not be complete and is compiled from ISO-Standards, ISO-6142 and ISO-6143. a. Uncertainty related to the balance and mass pieces Readability of balance. 1 mg for gas balance, 0.01 mg for micro gas balance Accuracy of balance including linearity. 2 mg for gas balance, 0.04 mg for chemical balance Incorrect zero point. Negligible Drift (thermal and time effects). Negligible Instability due to draught. Negligible Location of cylinder on the balance pan: Negligible Always keep the cylinder positioned at the center of weighing pan. Balance has automatic centering function Errors in the mass pieces used. Not used Buoyancy effects on the weights used: Buoyancy was corrected using tare cylinder at each weighing procedure b. Uncertainties related to the gas cylinder Page 21 of 55

Mechanical handling of cylinder due to: loss of metal, paints or labels from surface of cylinder; Negligible loss of metal from threads of valve/fitting; Negligible dirt on cylinder, valves or associated fitting. Negligible Adsorption/desorption effects on the external cylinder surface. Negligible Weighing room was well air conditioned and kept temperature and humidity constant so that no absorption/desorption effects on the external cylinder surface. Buoyancy effects resulted from: the cylinder itself; differences in temperature of the cylinder from surrounding air due to e.g. filling with gas; After filling the cylinder with N2, the cylinder was stood at room temp until it cool down to room temp, and also kept at the weighing room for more than 4 hours before weighing change of cylinder volume during filling; Negligible Change of cylinder volume was checked by measuring the diameter of cylinder, but very little change was observed. change of density of air due to changes in temperature, air pressure, humidity, and carbon dioxide content. Temperature and humidity were constant during the weighing of cylinder. Buoyancy effect by air pressure change was corrected by using tare cylinder with same shape and weight. CO2 content of weighing room was also analyzed, and CO2 effect was negligible. Uncertainty in determination of external cylinder volume. Negligible c. Uncertainties related to the component gases Residual gas in cylinder. Cylinder used for KRISS CRMs for H2S were evacuated under oil-less high vacuum system (consisting turbo-molecular pump with diaphram pump backup). Leakage resulting from; leakage of air into the cylinder after evacuation; Negligible Leakage of the cylinder was check by reweighing the cylinder after evacuation (one day later). No leak was noticed leakage of gas from the cylinder valve during filling; Negligible Leakage of cylinder also was check by pressurizing the filling system and monitoring pressure change for a day. No pressure change was observed. escape of gas from cylinder into transport lines. Negligible Gas remaining in transfer system when weight loss method is used. Weight loss method does not used. Absorption/reaction of components on internal cylinder surface. Loss of 10 ppm level of H2S by adsorption or reaction on the surface of cylinder has been checked by filling to new evacuated cylinder with same specification. Absorption less than 0.1% could not be checked by analysis. Therefore, 0.1% uncertainty was counted as a loss due to absorption. Reaction between components. : No No other component except H2S and N2 was not found from H2S CRM. Impurities in the balance gase, N2 used. Less than 1 ppb of H2S High purity nitrogen with 99.9999% purity was obtained from gas company, and the purity was checked by KRISS. Moisture, CO2, O2, Ar, CH4 were analyzed. Using the same analysis system, GC-AED, H2S was check down to ppb level. Impurities in H2S, 0.4% Insufficient homogenization. Negligible KRISS CRMs were sufficiently mixed after preparation using rolling mixer, then stood for a day at room temp before use. Uncertainty of molecular mass. Negligible Page 22 of 55

d. Uncertainties related to the analysis Repeatability and selectivity of the analyzer. Repeatability of GC-AED for 10 ppm level of H2S analysis was 0.2 ~ 0.4 %. Appropriateness of the calibration curve: response model and its residuals; Negligible mismatch of the sample gas and the calibration gas. Negligible Uncertainty related the corrections for the sampling: change in atmospheric pressure, temperature, sample flow rate, and sampling time. 0.1% Uncertainty sources with small contribution are not listed on following tables although the sources had been counted for uncertainty budget calculaltions. Uncertainty table: Hydrogen sulfide in N2 Uncertainty source Estimate X I Uncertainty of analysis, analytical response of sample, A sample Uncertainty of analysis, analytical response of standard, A x Factor related to purity of pure H2S chemical, fimp Factor related to adsorption loss in preparation of 10 ppm standard gas, fads10 Factor related to preparation homogeneity of 10 ppm standard gas, fhomo10 Factor related to stability of 10 ppm standard gas, fstab10 Factor related to preparation homogeneity of 2% standard gas, fhomo2 Factor related to stability of 2% standard gas, fstab2 x I 2592 2606 Weighing uncertainty of 2% 0.3170 standard gas in preparation of 10ppm std gas, W g balh2s210 Uncertainty related to micro gas balance, Wug Uncertainty related to gas balance, Wg Uncertainty related to handling of sampling cylinder in micro balance measurement, Whadug2 Uncertainty related to handling of cylinder in weighing of std gas cylinder, Whad2 Assumed distributio n t- distrubuti on t- distrubuti on Standard uncertainty u(x i) Sensitivity coefficient c I Contribution to standard uncertainty u I(y) umol/mol 2.78 3.87 x 10-3 0.0108 1.22 3.85 x 10-3 0.0047 0.996 normal 0.003 9.80 0.0294 1 normal 0.001 10.0 0.0100 1 t- distrubuti on 0.0023 10.0 0.0231 1 normal 0.001 10.0 0.0100 1 t- distrubuti on 0.002 9.80 0.0196 1 normal 0.001 9.80 0.0098 t- distrubuti on 43 x 10-6 g 31.7 0.0014 0 mg normal 0.05mg 32 x 10-6 0.0016 0 mg normal 10 mg 522 x 10-6 0.0052 0 mg normal 0.05 mg 32 x 10-6 0.0016 0 mg normal 5 mg 522 x 10-6 0.0026 Page 23 of 55

Coverage factor or degree of freedom: k=2.0 or dof=193 Expanded uncertainty: 0.15 umol/mol Optional You may provide additional data like the raw measurement data, information on your measurement procedure etc: Preparation of 10 umol/mol level of H2S CRM and Analysis of K-41 sample. Model equation for CRM preparation and analyis C sample = (A sample /A x * x H2S10 ) * f press ; x H2S10 = (n H2S210 / n tot10 * 1000000) * (f ads10 * f homo10 * f stab10 ); n tot10 = n H2S210 + n N212 + n N2dil ; n H2S210 = (m H2S210 / M toth2s10 ) * (x H2S2 /100) ; n N212 = (m H2S210 / M toth2s10 ) * (x N21 /100) ; n N2dil = m N2dil /M N2 ; M toth2s10 = M H2S * (x H2S2 /100) + M N2 * (x N21 /100) ; m H2S210 = (W balh2s210 + W ug / 1000000 + W hadug2 / 1000000) ; m N2dil = (W baln2dil + W g / 1000 + W had2 / 1000) ; x H2S2 = (n H2S2 / n tot2c * 100 ) * f purity * f homo2 * f stab2 ; x N21 = (n N21 / n tot2c * 100 ) ; n tot2c = n H2S2 + n N21 ; n H2S2 = m H2S2 / M H2S ; n N21 = m N21 / M N2 ; m H2S2 = (W balh2s2 + W g / 1000 + W had2 / 1000) ; m N21 = (W baln21 + W g / 1000 + W had2 / 1000) ; Page 24 of 55

List of quantities: Quantity Unit Definition C sample ppm conc of CCQM-41 from analysis A sample A x peak area from sample analysis peak area from standard analysis x H2S10 ppm conc of KRISS H2S std gas f press factor for pressure difference in sample introduction n H2S210 mole mole of H2S in 10 ppm std gas (kriss) n tot10 mole total mole of gases in 10 ppm H2S std gas f ads10 f homo10 f stab10 factor related to loss of H2S by adsorption in 10 ppm std gas uncertainty related to prep (homogeneity) of 10 ppm std gas uncertainty related to stability of 10 ppm H2S std gas n N212 mole mole of N2 contained 2% H2S used in preparation of 10 ppm H2S n N2dil mole mole of N2 used for preparation of 10 ppm H2S(as dilution gas) m H2S210 g wt(g) of 2% H2S gas taken for prep of 10 ppm gas M toth2s10 MW of 2% H2S x H2S2 % mole conc of H2S at 2% H2S gas x N21 % mole conc of N2 at 2% H2S gas m N2dil g wt of N2 dil gas in prep of 10 ppm gas M N2 M H2S MW of N2 MW of H2S W balh2s210 g wt of 2% H2S std gas taken for prep of 10 ppm gas W ug ug uncertainty related to chemical balance measurement W hadug2 ug uncertainty related to handling of cylinder in balance measurement W baln2dil g wt of N2 taken for prep of 10 ppm gas W g mg uncertainty related to chemical balance measurement W had2 mg uncertainty related to handling of cylinder in balance measurement n H2S2 mole mole of H2S used in 2% std gas n tot2c mole total mole of gases in 2% std gas f purity f homo2 f stab2 uncertainty of H2S purity uncertainty related to prep (homogeneity) uncertainty related to stability of 2% H2S std gas n N21 mole mole of N2 used in prep of 2% std gas m H2S2 g wt of H2S in prep of 2% std gas m N21 g wt of N2 used in prep of 2% std gas W balh2s2 g wt of H2S in prep of 2% std gas W baln21 g wt of N2 used in prep of 2% std gas Page 25 of 55

Uncertainty budget: Quantity Value Standard uncertainty Degrees of freedom Sensitivity coefficient Uncertainty contribution Corr.- coeff. Index A sample 2592.00 2.78 7 3.87 10-3 0.0108 ppm 0.22 0.049 A x 2606.00 1.22 13-3.85 10-3 -4.69 10-3 ppm -0.10 0.009 x H2S10 10.0977 ppm 0.0463 ppm f press 1.00000 1.00 10-3 50 10.0 0.0100 ppm 0.21 0.043 n H2S210 225.903 10-6 mole n tot10 838 10-9 mole 22.371636 mole 399 10-6 mole f ads10 1.00000 1.00 10-3 50 10.0 0.0100 ppm 0.21 0.043 f homo10 1.00000 2.30 10-3 3 10.0 0.0231 ppm 0.48 0.226 f stab10 1.00000 1.00 10-3 50 10.0 0.0100 ppm 0.21 0.043 n N212 0.01104116 mole n N2dil 3.07 10-6 mole 22.360369 mole 399 10-6 mole m H2S210 0.3170000 g 83.0 10-6 g M toth2s10 28.13511 2.56 10-3 x H2S2 2.00498 % 7.60 10-3 % x N21 97.99502 % 1.21 10-3 % m N2dil 626.3918 g 0.0112 g M N2 28.01348000 5.00 10-6 50 0.359 1.79 10-6 ppm 0.00 0.000 M H2S 34.08000000 5.00 10-6 50-0.295-1.47 10-6 ppm 0.00 0.000 W balh2s210 0.3170000 g 43.5 10-6 g 1.36 10 6 31.7 1.38 10-3 ppm 0.03 0.001 W ug 0.0 ug 50.0 ug 50 31.7 10-6 1.58 10-3 ppm 0.03 0.001 W hadug2 0.0 ug 50.0 ug 50 31.7 10-6 1.58 10-3 ppm 0.03 0.001 W baln2dil 626.3917500 g 43.5 10-6 g 21.8 10 6-0.0160-697 10-9 ppm 0.00 0.000 W g 0.0 mg 10.0 mg 50 522 10-6 5.22 10-3 ppm 0.11 0.012 W had2 0.0 mg 5.00 mg 50 522 10-6 2.61 10-3 ppm 0.05 0.003 n H2S2 0.521339 mole 328 10-6 mole n tot2c 26.002191 mole 727 10-6 mole f purity 0.99600 3.00 10-3 50 9.80 0.0294 ppm 0.61 0.366 f homo2 1.00000 2.00 10-3 3 9.80 0.0196 ppm 0.40 0.163 f stab2 1.00000 1.00 10-3 50 9.80 9.80 10-3 ppm 0.20 0.041 n N21 25.480852 mole 399 10-6 mole m H2S2 17.7673 g 0.0112 g m N21 713.8073 g 0.0112 g W balh2s2 17.767250000 g 310 10-9 g 2040 0.552 171 10-9 ppm 0.00 0.000 W baln21 713.807333333 g C sample 10.0435 ppm 0.0486 ppm 34 310 10-9 g 1940 0.0 0.0 ppm 0.0 0.0 Page 26 of 55