--------------------------------------------------------------------- Stable Isotopic Composition of Atmospheric Carbon Dioxide (13C and 18O) from the NOAA ESRL Carbon Cycle Cooperative Global Air Sampling Network, 1990-2014 Version: 2015-10-26 -------------------------------------------------------------------- CONTENTS 1. Data source and contacts 2. Use of data 2.1 Citation 3. Reciprocity 4. Warnings 5. Update notes 6. Introduction 7. DATA - General Comments 7.1 DATA - Sampling Locations 7.2 DATA - File Name Description 7.3 DATA - Event with single parameter 7.4 DATA - Event with multiple parameters 7.5 DATA - QC Flags 7.6 DATA - Monthly Averages 8. Data retrieval 9. References -------------------------------------------------------------------- 1. DATA SOURCE AND CONTACTS University of Colorado Institute of Arctic and Alpine Research (INSTAAR) James White, Bruce Vaughn, and Sylvia Michel Institute of Arctic and Alpine Research Campus Box 450 University of Colorado, Boulder, CO 80309-0450 USA telephone: (303) 735-5850 facsimile: (303) 492-6388 e-mail: sylvia.englund@colorado.edu SIL web site: http://instaar.colorado.edu/sil Jim White, CU-INSTAAR Sylvia Michel, CU-INSTAAR Bruce Vaughn, CU-INSTAAR -------------------------------------------------------------------- 2. USE OF DATA These data are made freely available to the public and the scientific community in the belief that their wide dissemination will lead to greater understanding and new scientific insights. The availability of these data does not constitute publication of the data. NOAA relies on the ethics and integrity of the user to insure that ESRL receives fair credit for their work. If the data are obtained for potential use in a publication or presentation, ESRL should be informed at the outset of the nature of this work. If the ESRL data are essential to the work, or if an important result or conclusion depends on the ESRL data, co-authorship may be appropriate. This should be discussed at an early stage in the work. Manuscripts using the ESRL data should be sent to ESRL for review before they are submitted for publication so we can insure that the quality and limitations of the data are accurately represented. 2.1 CITATION Please reference these data as White, J.W.C., B.H. Vaughn, and S.E. Michel (2015), University of Colorado, Institute of Arctic and Alpine Research (INSTAAR), Stable Isotopic Composition of Atmospheric Carbon Dioxide (13C and 18O) from the NOAA ESRL Carbon Cycle Cooperative Global Air Sampling Network, 1990-2014, Version: 2015-10-26, Path: ftp://aftp.cmdl.noaa.gov/data/trace_gases/co2c13/flask/. -------------------------------------------------------------------- 3. RECIPROCITY Use of these data implies an agreement to reciprocate. Laboratories making similar measurements agree to make their own data available to the general public and to the scientific community in an equally complete and easily accessible form. Modelers are encouraged to make available to the community, upon request, their own tools used in the interpretation of the ESRL data, namely well documented model code, transport fields, and additional information necessary for other scientists to repeat the work and to run modified versions. Model availability includes collaborative support for new users of the models. -------------------------------------------------------------------- 4. WARNINGS Every effort is made to produce the most accurate and precise measurements possible. However, we reserve the right to make corrections to the data based on recalibration of standard gases or for other reasons deemed scientifically justified. We are not responsible for results and conclusions based on use of these data without regard to this warning. The following important information is CRITICAL for the successful use of this data! Every effort is made to produce the most accurate and precise measurements possible. However, we reserve the right to make corrections to the data based on recalibration of standard gases or for other reasons deemed scientifically justified. We are not responsible for results and conclusions based on use of these data without regard to this warning. ********************** Attention must be paid to the flagging scheme, because we put ALL of our data (warts and all) into this data base. It is only through effective use of the flagging field that meaningful data can be identified. Please see the section that explains flagging. *********************** Special consideration for using CO2O18 data: Users of O18 data should be aware that many flasks in the network have been sampled without drying the air. This can result in isotopic exchange between CO2 and H2O, which can significantly deplete d18O values (see Gemery et al, 1996). NOAA is working to dry ALL air at wet sites prior to sampling, but a few 'wet' sites remain. As of May 2009, we have applied a new flagging scheme that makes more of the data from potentially wet sites available for use. SEE MORE DETAILS IN "COMMENTS" SECTION BELOW. Users of the O18 data are encouraged to please contact James White (James.White@colorado.edu) prior to using the O18 data for a current update. -------------------------------------------------------------------- 5. UPDATE NOTES +++++++++++++++++++++++++++++++ Lab-wide notes: 2015-08-03 The 3-letter site identification code for Ushuaia, Argentina (TDF) was changed to USH to be consistent with the WMO GAWSIS. 2011-10-07 We have introduced the term "measurement group", which identifies the group within NOAA and INSTAAR that made the actual measurement. We can now have multiple groups measuring some of the same trace gas species in our discrete samples. Measurement groups within NOAA and INSTAAR are ccgg: NOAA Carbon Cycle Greenhouse Gases (CCGG) hats: NOAA Halocarbons and other Atmospheric Trace Species (HATS) arl: INSTAAR Atmospheric Research Laboratory (ARL) sil: INSTAAR Stable Isotope Laboratory (SIL) curl: INSTAAR Laboratory for Radiocarbon Preparation and Research (CURL) We have also changed the file naming convention (see section "DATA - FILE NAME DESCRIPTION"). +++++++++++++++++++++++++++++++ Project-specific notes: 2013-08-27 Coordinates of some of the sample locations have changed. These changes improve the specified location based on new information. Changes tend to be minor and do not necessarily reflect a change in the actual sampling location. 2011-10-01 The data file format has been modified to include the measurement group and, additionally, the sample collection and analysis times now include second information (e.g., 2011 03 15 23 06 12). See section 7.3 for details. 2010-10-01 The format of the NOAA ESRL data records has been changed to include an estimate of the uncertainty associated with each measurement. The determination of the estimate is trace gas specific and described in section 6 (INTRODUCTION). +++++++++++++++++++++++++++++++ Parameter-specific notes: ************************************* January 2013 - Second position flags moved to first and third positions: Due to confusion resulting from our second position flags, they were moved to either first or third position. All flags are applied automatically by code except for '!' hand flags. Please see flagging section for more information. November 2010 - New second position flags: In order to improve the quality of the data, we have implemented additional flags in the second position. See more in the flagging section. October 2009 Update - Correction for 17O: All data are now processed with a revised correction for the influence of 17O (since 12C17O16O is an isobaric isotopologue of 13C16O2). Prior to this reprocessing (August 2009), the data were calculated based on the Craig (1957) measurements of the ratios of 17O and 13C in PDB carbonate. Craig also assumed a fractionation factor (lambda) between 17O and 18O of 0.5. However, work by Assonov and Brenninkmeijer (2003) suggests a more appropriate value of lamda is 0.528, due to non-mass-dependent fractionation in the global hydrological cycle. Furthermore, their calculations of the isotopic ratios of PDB are slightly different than those ascribed by Craig (1957). The resultant algorithms for the 17O correction were difficult to apply; however, the equations defined in Brand, Assonov, and Coplen (2009) simplify this alogrithm, and are now applied to our dataset. This is in keeping with recommendations provided by the WMO/IAEA CO2 Experts meeting held in Jena, Germany in 2008. June 2009 Update - New and improved flagging of compromised data: Historically, flasks were sampled at sites with high humidity without being dried, during some or all seasons. This is problematic because there can be isotopic exchange in the flasks between CO2 and H2O (See Gemery et al, 1996). As of May 2009, we have implemented a data filtering method which distinguishes between samples that are compromised by water vapor from those that can be trusted as valid measurements. See more details in "Comments" section below. November 2006 Update - d18O scale change: The CO2O18 measurement scale has been adjusted by +0.82 per mil, relative to V-PDB. This was done as a result of calibrations and comparisons with multiple reference materials from other labs. See the Comments section for more details. ************************************* -------------------------------------------------------------------- 6. INTRODUCTION The files in this archive contain measurements of the stable isotopic composition (13C and 18O) of atmospheric carbon dioxide. The isotopic analyses were carried out at the Stable Isotope Laboratory, CU-INSTAAR, using samples of air provided by the NOAA CMDL Carbon Cycle Cooperative Global Air Sampling Network. For a complete listing of flask sampling site names and locations, see Section 7.1. Since 1989, the Stable Isotope Laboratory at INSTAAR has been measuring the stable isotopic composition of CO2 from weekly flask samples of air obtained from the CMDL network. Starting with a selection of samples from just six sites and two ships in 1990, the measurement effort has grown to include all sites in the NOAA/ESRL program. The degree to which isotopic measurements made on atmospheric samples are useful is seriously constrained by the precision of the analysis method used. A small change of just 0.01 per mil in delta 13C (henceforth written as d13C) could have a large effect on models of partitioning of carbon between the oceans and the biosphere. Therefore, a high precision instrument and high accuracy calibrations are desirable. In 1990 we began making measurements of flask samples with a VG SIRA Series II mass spectrometer. This instrument and extraction system produced an overall reproducibility of 0.03 per mil for d13C and 0.05 per mille for d18O (1 sigma of 9 standard cylinders within a run). (See more about the uncertainty of measurements in the comments section). In 1996 a Micromass Optima mass spectrometer was purchased and fitted with a custom manifold and extraction system, based on the methods that were proven on the SIRA instrument. This system yielded better precision: +/- 0.012 per mil d13C and +/- 0.03 1 per mil for d18O. In 2005 a GV Instruments (formerly Micromass and now Elementar) Isoprime dual Inlet mass spectrometer was brought online for analysis of Progammable Flask Packages (PFPs) which obtain samples from aircraft, towers, and also remote surface sites. The extraction system is conceptually identical to the previous two instruments (described further the Comments section). The reproducibility for the Isoprime dual inlet system is generally +/- 0.01 per mil for d13C and +/- 0.03 for d18O. This instrument can also analyze network flasks and reference cylinders. In 2006, we acquired another GV Instruments Isoprime dual Inlet mass spectrometer which has been be fitted with extraction systems to handle CO2 extracted from whole air, CO2 evolved from carbonates in phosphoric acid, and CO2 equilibrated with water. With these three sample preparation systems, we hope to better calibrate CO2 in air to both d13C and d18O measurement scales. We continue to cooperate with the many different efforts to cross calibrate between stable isotope labs, including flask inter- comparisons, round-robin high pressure cylinder exchanges, and measurement of the new Jena Air Reference Set (JRAS), produced from CO2 evolved from carbonates mixed with CO2-free air by the Max Planck Institute for Biogeochemistry. If, as a result of these calibration efforts it becomes necessary to adjust our scale, we will document it here. For more information about the Stable Isotope Lab or the instruments used for theses analyses, including a downloadable document on methods, see the SIL web site at: http://instaar.colorado.edu/sil. ------------------------------------------------------------------- 7. DATA - GENERAL COMMENTS Carbon and oxygen isotopes of carbon dioxide (d13C and d18O of CO2) The methods to measure isotope ratios of carbon and oxygen in carbon dioxide described below have been developed at the University of Colorado/INSTAAR Stable Isotope Laboratory. For more information, please refer to our web site: http://instaar.colorado.edu/sil or contact our assistant laboratory manager, Sylvia Michel: sylvia.englund@colorado.edu. Methods: To measure the stable isotopes of carbon dioxide, the CO2 must first be separated from bulk air. A vacuum pulls the air from the flask or PFP through a mass flow controller, a water trap held at -85 degrees C, and a liquid nitrogen trap at -197 degrees C. This extraction lasts for 10 minutes, sampling approximately 400 cc of air. After the CO2 sample has been isolated in the liquid nitrogen trap, a heating wire warms the sample to -20 degrees C, releasing it to the sample bellows of the dual inlet mass spectrometer. While the isotope ratio is measured by the mass spectrometer (which takes approximately ten minutes), extraction of the next sample is underway, increasing the throughput of the instrument. Mass spectrometry works on the principle that an ion in a magnetic field bends at a radius proportional to its mass-to-charge ratio. When the CO2 extracted from our sample is ionized by electrons, the ion beam is accelerated through a flight tube, separated by an electromagnetic field into masses 44, 45, and 46, and detected by Faraday cups. The currents obtained by these cups determine the ratio of heavy to light isotopes of the sample. The sample is iteratively measured against a reference gas of isotopically similar CO2 in the reference bellows. Ten ratios are obtained, with a standard error usually less than 0.005. The isotopic value of the gas in the reference bellows is not known; to calibrate to the PDB scale, the international standard for carbon and oxygen isotope ratios in CO2, we extract air from a cylinder with known isotope values of CO2 (the "working reference"). This follows the "identical treatment" principle, which states that samples and references should be handled using identical procedures to avoid introuducing bias. The isotope ratios are reported as delta values (written here as "d") where d13C = [(13C/12Csample/(13C/12Cstandard/)-1] x 1000. Delta values are reported in units of per mil. Corrections: Several corrections are applied to the raw data. To correct for the fractionation of gas in the reference bellows, we apply a drift correction using the working reference, which is run at the beginning, middle, and end of the run. Because N2O also freezes at liquid nitrogen and is isobaric with CO2, we mathematically correct for it using the concentration of N2O obtained by NOAA measurements (Mook and Jongsma 1987). We also correct for the presence of 12C17O16O, an isotopologue of 13CO2, using the Brand, Coplen, and Assonov correction (2009; see more in the updates section). Finally, we anchor the measurements to the PDB scale using the calibrated isotopic values of our working reference. Calibrations: PDB (Pee Dee Bellemnite) is a carbonate, and its link to CO2-in-air standards differs among laboratories. Ongoing intercomparison experiments are essential for assessing the comparability of isotope measurements from one lab with another (Masarie et al., 2001; Allison et al., 2003; Ghosh, P., Patecki, M., Rothe, M. and Brand, W.A. 2005). Measurement accuracy based on results from intercomparison experiments is 0.03 per mil. Uncertainty: Samples are run daily from a cylinder whose isotope value is known: this "trap" tank alerts us to any problems with the mass spectrometer or extraction system. We also use this cylinder to calculate uncertainty, which we define as the standard deviation of ten runs' worth of trap measurements. Usually four trap samples are run a day, of which all but the first are used (due to known irregularities caused by the tank regulator). For each run, the uncertainty is calculated from the trap data from that run, regardless of its flagging, and the previous nine runs of unflagged data. (For example, if the current run is flagged, that trap data will be used in the uncertainty calculation; however, the next day's run will not use the flagged data in its uncertainty calculation. Presumably, whatever caused one run to be flagged will have been repaired.) As of November 2009, uncertanty is calculated with each run, and has been back-calculated to 2004. On the Optima (flask measurement system) uncertainty averages 0.014 permil for d13C and 0.035 permil for d18O; for the Isoprime (PFP system), uncertainty averages 0.017 permil for d13C and 0.04 permil for d18O. Special considerations for d18O of CO2: The INSTAAR-NOAA scale was known to be +0.82 per mil different from other labs measuring the d18O of CO2 in air, due to the temperature of the carbonate reaction that established the original CO2-in-air reference (Masarie et al., 2001; Allison et al., 2003). The application of a +0.82 per mil scale adjustment has been made as of November 2006 based on the results of the following comparisons: the CLASSIC cylinders; standards prepared by Hitoshi Mukai's NARCIS I and II; CO2 in air samples from W. Brand; flask inter comparisons with CSIRO (Cape Grim); and IAEA standard materials. If air is sampled without drying, oxygen can exchange between carbon dioxide and water vapor in the flasks, in some cases entirely masking the desired signal (Gemery et al, 1996). While pair agreement flagging identifies most of the problem flasks, some flasks still have good pair aggreement and are within the normal range of acceptable values, even though they are biased by the influence of water vapor (Evans, 2008). Comparisons between d18O of CO2 and specific humidity, calculated from the National Climatic Data Center (NCDC), show that at tropical and semi-tropical NOAA/ESRL sites, d18O of CO2 of non-dried samples has a seasonal cycle that is strongly anti-correlated to the specific humidity, while the d18O of dried samples display a seasonal cycle that occurs 1-2 months earlier than the specific humidity seasonal cycle. The latter phasing is expected, given the seasonal phasing between climate over the ocean and land; the former is consistent with a small, but measurable isotope exchange in the flasks. Based on investigations by Evans et al (2008), a threshold of specific humidity, below which non-dried air flask samples can be considered unaltered, has been set at 12 g/kg, and this filter has been applied to the d18O of CO2 database. For example, for sites where specific humidity exceeds 12 g/kg, each non-dried flask sample has been matched to a NCDC measurement of specific humidity within a 4-hour window. If specific humidity for that sample exceeded 12 g/kg, the SIL measurement for d18O is flagged with a 'W' (for wet) in the first position. If specific humidity data was unavailable, the sample is also flagged. Unfortunately there are several sites where the specific humidity is always above 12 g/kg (including CHR, GMI, POCN00, POCN05, POCN10, POCN15, POCS05, POCS10, POCS15, RPB, SEY, and SMO); therefore all data sampled without drying are flagged. This approach segregates suspect and reliable data, providing a much-improved record of d18O of CO2 over the past two decades. Sample collection: In the case of samples acquired from ships at sea, the Pacific Ocean Cruise (POC, travelling between the US west coast and New Zealand or Australia) data have been merged and grouped into 5 degree latitude bins. For the South China Sea cruises (SCS) the = data are grouped in 3 degree latitude bins. Sampling frequencies are approximately weekly for the fixed sites and average one sample every 3 weeks per latitude zone for POC and about one sample every week per latitude for SCS. The air samples are collected by two general methods: flushing and then pressurizing glass flasks with a pump, or opening a stopcock on an evacuated glass flask. During each sampling event, a pair of flasks is filled. For more information on sampling units, see http://www.esrl.noaa.gov/gmd/ccgg/psu/index.html The Pacific Ocean Cruise (POC, travelling between the US west coast and New Zealand or Australia) data have been merged and grouped into 5 degree latitude bins. For the South China Sea cruises (SCS) the data are grouped in 3 degree latitude bins. Sampling frequencies are approximately weekly for the fixed sites and average one sample every 3 weeks per latitude zone for POC and about one sample every week per latitude for SCS. Historically, samples have been collected using two general methods: flushing and then pressurizing glass flasks with a pump, or opening a stopcock on an evacuated glass flask; since 28 April 2003, only the former method is used. During each sampling event, a pair of flasks is filled. ------------------------------------------------------------------- 7.1 DATA - SAMPLING LOCATIONS For a summary of sampling locations, please visit http://www.esrl.noaa.gov/gmd/dv/site/site_table.html. IMPORTANT NOTES: 1. Data for all species may not be available for all sites listed in the table. 2. The exact location of a sampling site recorded in our database may change or become better defined over time. The latitude, longitude, and altitude of a sample event is based on the best information available at the time of sample collection. Differences in sample position associated with a particular site may be due to the site moving or changes in technology that permit a more accurate location determination. To view near real-time data, manipulate and compare data, and create custom graphs, please visit http://www.esrl.noaa.gov/gmd/dv/iadv/. ------------------------------------------------------------------- 7.2 DATA - FILE NAME DESCRIPTION Encoded into each file name are the parameter (trace gas identifier); sampling site; sampling project; laboratory ID number; measurement group; and optional qualifiers that further define the file contents. All file names use the following naming scheme: 1 2 3 4 5 6 7 [parameter]_[site]_[project]_[lab ID number]_[measurement group]_[optional qualifiers].[file type] 1. [parameter] Identifies the measured parameter or trace gas species. (ex) co2 Carbon dioxide ch4 Methane co2c13 d13C (co2) merge more than one parameter 2. [site] Identifies the sampling site code. (ex) brw pocn30 car amt 3. [project] Identifies sampling platform and strategy. (ex) surface-flask surface-pfp surface-insitu aircraft-pfp aircraft-insitu tower-insitu 4. [lab ID number] A numeric field that identifies the measurement laboratory (1,2,3, ...). NOAA ESRL is lab number 1 (see http://www.esrl.noaa.gov/gmd/ccgg/obspack/labinfo.html). 5. [measurement group] Identifies the group within NOAA or INSTAAR that makes the actual measurement. See Section 5 (UPDATE NOTES) for details. (ex) ccgg hats arl sil 6. [optional qualifiers] Optional qualifier(s) may indicate data subsetting or averaging. Multiple qualifiers are delimited by an underscore (_). A more detailed description of the file contents is included within each data file. (ex) event All measurement results for all collected samples (discrete (flask) data only). hour_#### Computed hourly averages for the specified 4-digit year (quasi-continuous data only) HourlyData Computed hourly averages for entire record (quasi-continuous data only) DailyData Computed daily averages for entire record (quasi-continuous data only) day MonthlyData Computed monthly averages for entire record (quasi-continuous data only) month 7. [file type] We now provide some NOAA Global Monitoring Division measurements in two unique file formats (netCDF and ASCII text). The Network Common Data Form (NetCDF) is a self-describing, machine-independent data format that supports creation, access, and sharing of array-oriented scientific data. To learn more about netCDF and how to read netCDF files, please visit http://www.unidata.ucar.edu. The ASCII text file is derived directly from the netCDF file. The text file is also self-describing and can be viewed using any text editor. "Self-describing" means the file includes enough information about the included data (called metadata) that no additional file is required to understand the structure of the data and how to read and use the data. (ex) txt ASCII text file nc netCDF4 file ------------------------------------------------------------------- 7.3 DATA - EVENT WITH SINGLE PARAMETER The event data files in ftp://aftp.cmdl.noaa.gov/data/trace_gases/co2c13/flask/surface/ use the following naming scheme (see Section 7.2): [parameter]_[site]_[project]_[lab ID number]_[measurement group]_[optional qualifiers].txt (ex) CH4_pocn30_surface-flask_1_ccgg.txt contains CH4 ccgg measurement results for all surface flask samples collected on the Pacific Ocean Cruise sampling platform and grouped at 30N +/- 2.5 degrees. (ex) CO2_brw_surface-flask_1_ccgg.txt contains CO2 ccgg measurement results for all surface flask samples collected at Barrow, Alaska. The data files contain multiple lines of header information followed by one record for each atmospheric measurement of a single parameter or trace gas species. Fields are defined as follows: Field 1: [SITE CODE] The three-character sampling location code (see above). Field 2: [YEAR] The sample collection date and time in UTC. Field 3: [MONTH] Field 4: [DAY] Field 5: [HOUR] Field 6: [MINUTE] Field 7: [SECOND] Field 8: [FLASK ID] The sample container ID. Field 9: [METHOD] A single-character code that identifies the sample collection method. The codes are: P - Sample collected using a portable, battery powered pumping unit. Two flasks are connected in series, flushed with air, and then pressurized to 1.2 - 1.5 times ambient pressure. D - Similar to P but the air passes through a condenser cooled to about 5 deg C to partially dry the sample. G - Similar to D but with a gold-plated condenser. T - Evacuated flask filled by opening an O-ring sealed stopcock. S - Flasks filled at NOAA ESRL observatories by sampling air from the in situ CO2 measurement air intake system. N - Before 1981, flasks filled using a hand-held aspirator bulb. After 1981, flasks filled using a pump different from those used in method P, D, or G. F - Five liter evacuated flasks filled by opening a ground glass, greased stopcock. Field 10: [TRACE GAS NAME] Gas identifier (e.g., co2, co2c13). Field 11: [MEASUREMENT GROUP] Identifies the group within NOAA and INSTAAR making the actual measurement (e.g., ccgg, hats, arl). See Section 5 (UPDATE NOTES) for details. Field 12: [MEASURED VALUE] Dry air mole fraction or isotopic composition. Missing values are denoted by -999.99[9]. Field 13: [ESTIMATED UNCERTAINTY] Estimated uncertainty of the reported measurement value. Missing values are denoted by -999.99[9]. Field 14: [QC FLAG] A three-character field indicating the results of our data rejection and selection process, described in section 7.5. Field 15: [INSTRUMENT] A 2-character code that identifies the instrument used for the measurement. Field 16: [YEAR] The measurement date and time in LT. Field 17: [MONTH] Field 18: [DAY] Field 19: [HOUR] Field 20: [MINUTE] Field 21: [SECOND] Field 22: [LATITUDE] The latitude where the sample was collected, (negative (-) numbers indicate samples collected in the Southern Hemipshere). Field 23: [LONGITUDE] The longitude where the sample was collected, (negative (-) numbers indicate samples collected in the Western Hemisphere). Field 24: [ALTITUDE] The altitude of the sample inlet (masl). Field 25: [EVENT NUMBER] A long integer that uniquely identifies the sampling event. Fields in each line are delimited by whitespace. (ex) SUM 2009 08 21 13 53 00 2145-99 P ch4 ccgg 1857.010 1.100 ... H11 2010 05 12 16 57 00 72.5800 -38.4800 3238.00 296091 ------------------------------------------------------------------- 7.4 DATA - EVENT WITH MULTIPLE PARAMETERS On special request we can distribute a "merged" file, which includes for each sampling event, measurement results for muliple parameters or trace gas species. A merged file does not include all information found in a single parameter data file. For example, merged files exclude measurement uncertainty, analysis instrument ID and date and time for each parameter. Thus, the single parameter data file is our most comprehensive data archive. The format of a merged file is slightly different from single parameter event file. A "merged" file will have the word "merge" in the parameter field of the file name. The file name does not inform on the number of parameters included in the file. Merged data files use the following naming scheme (see Section 7.2): merge_[site]_[project]_[lab ID number]_[measurement group]_[optional qualifiers].txt (ex) merge_pocn30_surface-flask_1_ccgg.txt contains ccgg measurement results for two or more parameters for all surface flask samples collected on the Pacific Ocean Cruise sampling platform and grouped at 30N +/- 2.5 degrees. (ex) merge_brw_surface-flask_1_ccgg.txt contains ccgg measurement results for two or more parameters for all surface flask samples collected at Barrow, Alaska. The data files contain multiple lines of header information followed by one record for each atmospheric measurement of a single parameter or trace gas species. Fields are defined as follows: Field 1: [SITE CODE] The three-character sampling location code (see above). Field 2: [YEAR] The sample collection date and time in UTC. Field 3: [MONTH] Field 4: [DAY] Field 5: [HOUR] Field 6: [MINUTE] Field 7: [SECOND] Field 8: [FLASK ID] The sample container ID. Field 9: [METHOD] A single-character code that identifies the sample collection method. The codes are: P - Sample collected using a portable, battery powered pumping unit. Two flasks are connected in series, flushed with air, and then pressurized to 1.2 - 1.5 times ambient pressure. D - Similar to P but the air passes through a condenser cooled to about 5 deg C to partially dry the sample. G - Similar to D but with a gold-plated condenser. T - Evacuated flask filled by opening an O-ring sealed stopcock. S - Flasks filled at NOAA ESRL observatories by sampling air from the in situ CO2 measurement air intake system. N - Before 1981, flasks filled using a hand-held aspirator bulb. After 1981, flasks filled using a pump different from those used in method P, D, or G. F - Five liter evacuated flasks filled by opening a ground glass, greased stopcock. Field 10: [LATITUDE] The latitude where the sample was collected, (negative (-) numbers indicate samples collected in the Southern Hemipshere). Field 11: [LONGITUDE] The longitude where the sample was collected, (negative (-) numbers indicate samples collected in the Western Hemisphere). Field 12: [ALTITUDE] The altitude of the sample inlet (masl). Field 13: [EVENT NUMBER] A long integer that uniquely identifies the sampling event. There is a group of 4 fields for each parameter and measurement group included in the merge file. Field ##+1: [TRACE GAS NAME] Gas identifier (e.g., co2, co2c13). Field ##+2: [MEASUREMENT GROUP] Identifies the group within NOAA and INSTAAR making the actual measurement (e.g., ccgg, hats, arl). See Section 5 (UPDATE NOTES) for details. Field ##+3: [MEASURED VALUE] Dry air mole fraction or isotopic composition. Missing values are denoted by -999.99[9]. Field ##+4: [QC FLAG] A three-character field indicating the results of our data rejection and selection process, described in section 7.5. Fields in each line are delimited by whitespace. (ex) SGP 2002 04 02 21 15 00 551-91 D 36.6200 -97.4800 374.00 115605 co2 CCGG 378.870 -.. ch4 CCGG 1874.995 ... co CCGG 163.800 ... h2 CCGG 496.350 ... n2o CCGG 319.380 ..P sf6 CCGG 5.130 ..P co2c13 SIL -8.477 +.. co2o18 SIL -0.271 ... ------------------------------------------------------------------- 7.5 QC FLAGS NOAA ESRL uses a 3-column quality control flag where each column is defined as follows: column 1 REJECTION flag. An alphanumeric other than a period (.) in the FIRST column indicates a sample with obvious problems during collection or analysis. This measurement should not be interpreted. column 2 SELECTION flag. An alphanumeric other than a period (.) in the SECOND column indicates a sample that is likely valid but does not meet selection criteria determined by the goals of a particular investigation. column 3 INFORMATION flag. An alphanumeric other than a period (.) in the THIRD column provides additional information about the collection or analysis of the sample. WARNING: A "P" in the 3rd column of the QC flag indicates the measurement result is preliminary and has not yet been carefully examined by the PI. The "P" flag is removed once the quality of the measurement has been determined. Samples are collected in pairs, and in order for a value to be retained, the flask pair difference must meet the threshold tolerance for the specific instrument being used. For samples analyzed prior to 1996 pairs must agree to within 0.09 per mil for C13 and 0.15 per mil for O18. For post-1996 analyses, done on both Optima and Isoprime mass spectrometers, pairs must agree to within 0.06 per mil for C13 and 0.12 per mil for O18. Typically, pairs are flagged for pair rejection more frequently for d18O of CO2 than for d13C of CO2. The reason for the lower success rate in O18 pair agreement is likely related to the isotopic exchange of water vapor and CO2 in the flask during storage and analysis. This is dramatically demonstrated in the particularly low pair agreements at high humidity, low latitude sites. See Section 7, Data-General Comments for more details on the O18 data set. FLAG DEFINITIONS *********************Reject flags (1st flag character)******************** Automatic flags - applied daily during data reduction by a processing program A problems in analysis or data reduction. Typically, the A is applied when the working references run at the beginning, middle, or end of the run have a standard deviation higher than 0.04 or 0.08 permil for d13C and d18O, respectively. For measurements prior to 1996, values differ by more than 0.075 and 0.14 permil, respectively. + bad pair (high member). Values differ by more than 0.06 or 0.12 permil for d13C and d18O, respectively. For measurements prior to 1996, values differ by more than 0.09 and 0.15 permil, respectively) - bad pair (low menber). See above. H The trap tank run on the same day has an average d13C or d18O value greater than 0.15 or 0.3 permil above its long term averages, respectively (since 2005.) L The trap tank run on the same day has an average d13C or d18O value less than 0.15 or 0.3 permil above its long term averages, respectively (since 2005.) . good flask Flags applied periodically by code C flagged for CO2 mole fraction (by NOAA CO2 measurement lab) W flask sampled 'wet' (applied to o18 data only). See comments. N/n problem due to sample collection (inherited from CO2 measurements) Hand flags - applied by hand or in software (REFLAG) to selected flasks ! hand flag *******************Non-background flags (2nd flag character)**************** ground flags - applied by outlier-identification software (CCG_FILTER) X Outlier by more than 3-sigma from a CCGVU curve x Outlier by more than 3-sigma in CO2 concentration . Good flask Hand flags - applied by hand to selected flasks ! hand flag *********************Retain flags (3rd flag character)*********************** Automatic flags - applied during data reduction (by the processing program) S - single flask (flask without a pair mate) o - no trap data available for comparisons H The trap tank run on the same day has an average d13C or d18O value greater than 0.045 or 0.09 permil above its long term averages, respectively (since 2005.) L The trap tank run on the same day has an average d13C or d18O value less than 0.045 or 0.09 permil above its long term averages, respectively (since 2005.) P Data has poor precision in dual inlet analysis: greater than 0.02 permil for d13C and 0.03 permil for d18O (since 2005). T The trap tank run on the same day has high standard deviation of its (typically) three measurements: above 0.08 permil and 0.16 permil for d13C and d18O respectively (since 2005). Additional flags L linked flask (0.5-liter flask analyzed together with its mate) I flask also analyzed by another lab (aliquot taken) i same as "I" above, but displaced a previous flag in this field. . no comments ------------------------------------------------------------------- 7.6 DATA - MONTHLY AVERAGES The monthly data files in ftp://aftp.cmdl.noaa.gov/data/trace_gases/co2c13/flask/surface/ use the following naming scheme (see Section 7.2): [parameter]_[site]_[project]_[lab ID number]_[measurement group]_month.txt (ex) CH4_pocn30_surface-flask_1_ccgg_month.txt contains CH4 ccgg monthly mean values for all surface flask samples collected on the Pacific Ocean Cruise sampling platform and grouped at 30N +/- 2.5 degrees. (ex) CO2_brw_surface-flask_1_ccgg_month.txt contains CO2 ccgg monthly mean values for all surface flask samples collected at Barrow, Alaska. Monthly means are produced for each site by first averaging all valid measurement results in the event file with a unique sample date and time. Values are then extracted at weekly intervals from a smooth curve (Thoning et al., 1989) fitted to the averaged data and these weekly values are averaged for each month to give the monthly means recorded in the files. Flagged data are excluded from the curve fitting process. Some sites are excluded from the monthly mean directory because sparse data or a short record does not allow a reasonable curve fit. Also, if there are 3 or more consecutive months without data, monthly means are not calculated for these months. The data files contain multiple lines of header information followed by one line for each available month. Fields are defined as follows: Field 1: [SITE CODE] The three-character sampling location code (see above). Field 2: [YEAR] The sample collection year and month. Field 3: [MONTH] Field 4: [MEAN VALUE] Computed monthly mean value In these files a monthly mean value of -999.99 denotes months where there are insufficient data to calculate a monthly mean. ------------------------------------------------------------------- 8. DATA RETRIEVAL To transfer all files in a directory, it is more efficient to download the tar or zipped files. To transfer a tar file, use the following steps from the ftp prompt: 1. ftp> binary ! set transfer mode to binary 2. ftp> get filename.tar.gz ! transfer the file 3. ftp> bye ! leave ftp 4. $ gunzip filename.tar.gz ! unzip your local copy 5. $ tar xvf filename.tar ! unpack the file To transfer a zipped file, use the following steps from the ftp prompt: 1. ftp> binary ! set transfer mode to binary 2. ftp> get filename.zip ! transfer the file 3. ftp> bye ! leave ftp 4. $ unzip filename.zip ! uncompress your local copy ------------------------------------------------------------------- 9. REFERENCES For more publications see the Stable Isotope Lab web site at: http://instaar.colorado.edu/research/labs-groups/stable-isotope-laboratory/publications-detail/ Assonov, S.S., and Brenninkmeijer, C.A.M. 2003. On the 17O correction for CO2 mass spectrometric isotopic analysis. Rapid Communications in Mass Spectrometry 17(10): 1007-1016. Battle. M., M.L. Bender, P.P. Tans, J.W.C. White, J.T. Ellis, T. Conway, and R.J. Francey. 2000. Global carbon sinks and their variability inferred from the atmospheric O2 and d13C. Science 287: 2467-2470. Brand, W.A., S.S. Assonov, and T.B. Coplen. 2009. Correction for the 17O interference in d13C measurements when analyzing CO2 with stable isotope mass spectrometry (IUPAC Technical Report). Journal of Pure and Applied Chemistry 82(8): 1719-1733. Ciais, P., P.P. Tans, J.W.C. White, M. Trolier, R.J. Francey, J.A. Berry, D.R. Randall, R.J. Sellers, J.G. Collatz and D.S. Schimel. 1995. Partitioning of ocean and land uptake of CO2 as inferred by d13C measurements from the NOAA/CMDL global air sampling network, Journal of Geophysical Research 100: 5051-5070. Ciais, P., P.P. Tans, M. Trolier, J.W.C. White and R.J. Francey. 1995. A large northern hemisphere terrestrial CO2 sink indicated by the 13C/12C ratio of atmospheric CO2. Science 269: 1098-1102. Ciais, P., A.S. Denning, P.P. Tans, J.A. Berry, D.A. Randall, G.J. Collatz, P.J. Sellers, J.W.C. White, M. Trolier, H.A.J. Meyer, R.J. Francey, P. Monfray, and M. Heimann. 1997. A three-dimensional synthesis study of d18O in atmospheric CO2. 1. Surface fluxes. Journal of Geophysical Research 102: 5857-5872. Ciais, P., P.P. Tans, A.S. Denning, R.J. Francey, M. Trolier, H.A.J. Meyer, J.W.C. White, J.A. Berry, D.A. Randall, G.J. Collatz, P.J. Sellers, P. Monfray, and M. Heimann. 1997. A three-dimensional synthesis study of d18O in atmospheric CO2. 2. Simulations with the TM2 transport model. Journal of Geophysical Research 102: 5873-5883. Conway, T.J., P.P. Tans, L.S. Waterman, K.W. Thoning, D.R. Kitzis, K.A. Masarie, and N. Zhang. 1994. Evidence for interannual variability of the carbon cycle from the NOAA/CMDL global air sampling network. Journal of Geophysical Research 99: 22831- 22855. Cuntz, M., P. Ciais, G. Hoffmann, C.E. Allison, R.J. Francey, W. Knorr, P.P. Tans, J.W.C. White, I. Levin. 2003. A comprehensive global three- dimensional model of delta O-18 in atmospheric CO2: 2. Mapping the atmospheric signal. Journal of Geophysical, Research- Atmospheres 108(D17): article 4528. Evans, C.U. 2008. d18O of atmospheric carbon dioxide: Towards the development of an artifact free database from the NOAA/ESRL Carbon Cycle Cooperative Global Air Sampling Network. Masters Thesis, INSTAAR, University of Colorado, Boulder. Francey, R.J., P.P. Tans, C.E. Allison, I.G. Enting, J.W.C. White and M. Trolier. 1995. Changes in oceanic and terrestrial carbon uptake since 1982. Nature 373: 326-330. Gemery, P.A., M. Trolier, and J.W.C. White. 1996. Oxygen isotope exchange between carbon dioxide and water following atmospheric sampling using glass flasks. Journal of Geophysical Research - Atmospheres 101: 14415-14420. Ghosh, P., M. Patecki, M. Rothe and W.A. Brand. 2005. Calcite-CO2 mixed into CO2-free Air: A New CO2 in-Air Stable Isotope Reference Material for the VPDB Scale. Rapid Communications in Mass Spectrometry 19: 1097-1119. Miller, J.B., P.P. Tans, J.W.C. White, T.J. Conway, B.H. Vaughn. 2003. The atmospheric signal of terrestrial carbon isotopic discrimination and its implication for partitioning carbon fluxes. Tellus Series B-Chemical and Physical Meteorology 55(2): 197-206. Mook, W.G. and J. Jongsma. 1987. Measurement of the N2O correction for 13C/12C ratios of atmospheric CO2 by removal of N2O. Tellus 39B: 96-99. Tans, P.P., T.J. Conway, and T. Nakazawa. 1989a. Latitudinal distribution of the sources and sinks of atmospheric carbon dioxide from surface observations and an atmospheric transport model. Journal of Geophysical Research 94: 5151-5172. Tans, P.P, K.W. Thoning, W.P. Elliott, and T.J. Conway. 1989b. Background atmospheric CO2 patterns from weekly flask samples at Barrow, Alaska: Optimal signal recovery and error estimates, in NOAA Tech. Memo. (ERL ARL- 173). Environmental Research Laboratories, Boulder, CO, 131 pp. Tans, P.P., I.Y. Fung, and T. Takahashi. 1990. Observational constraints on the global atmospheric CO2 budget. Science 247: 1431-1438. Thoning, K.W., P.P. Tans, and W.D. Komhyr. 1989. Atmospheric carbon dioxide at Mauna Loa Observatory 2. Analysis of the NOAA GMCC data, 1974-1985, Journal of Geophysical Research 94: 8549-8565. Trolier, M., J.W.C. White, P.P. Tans, K.A. Masarie and P.A. Gemery. 1996. Monitoring the isotopic composition of atmospheric CO2: measurements from the NOAA Global Air Sampling Network. Journal of Geophysical Research 101: 25897-25916. Vaughn, B., Ferretti, D., Miller, J. & White, J. 2004: Stable isotope measurements of atmospheric CO2 and CH4. Handbook of Stable Isotope Analytical Techniques, vol 1, ch.14, Elsiever, 1248 p. -------------------------------------------------------------------