Difference between revisions of "Field Data Analysis Guide"

From Jimenez Group Wiki
Jump to: navigation, search
m (Make campaign specific adjustments to fragmentation table)
m (Make campaign specific adjustments to fragmentation table)
Line 321: Line 321:
 
The purpose of the fragmentation table is to apportion the ToF-AMS signal among the different chemical species present in the aerosol (ammonium, chloride, nitrate, organics, and sulfate).  The following discussion assumes that the duty cycle correction has been applied to the data (this routinely occurs during data pre-processing in Squirrel unless that function has specifically been disabled by the user) and that the data has been processed to the point of calculating the air beam correction.<br/>   
 
The purpose of the fragmentation table is to apportion the ToF-AMS signal among the different chemical species present in the aerosol (ammonium, chloride, nitrate, organics, and sulfate).  The following discussion assumes that the duty cycle correction has been applied to the data (this routinely occurs during data pre-processing in Squirrel unless that function has specifically been disabled by the user) and that the data has been processed to the point of calculating the air beam correction.<br/>   
 
<br/>
 
<br/>
Details of the signal apportionment are presented in [http://cires.colorado.edu/jimenez/Papers/Allan_Frag_Table_Published.pdf Allan et al., J. Aerosol Sci., 2004].<br/>
+
 
 +
The fundamentals and an example (for NH4+) of the frag. table are presented in [http://cires.colorado.edu/jimenez/Papers/Allan_Frag_Table_Published.pdf Allan et al., J. Aerosol Sci., 2004].<br/>
 
<br/>
 
<br/>
 
Below are detailed discussions and suggestions for modifying the ToF-AMS fragmentation table.<br/>
 
Below are detailed discussions and suggestions for modifying the ToF-AMS fragmentation table.<br/>

Revision as of 13:02, 15 February 2011

Contents


Introduction

The ToF-AMS Field Data Analysis Guide is dedicated to the analysis of complex ToF-AMS field data sets and is intended to establish standard practices in this area on a community-wide basis. It is assumed that users will have some familiarity with Igor and the community based AMS analysis tool, Squirrel. Preliminary AMS data can be generated while a project is underway without worrying about some details; this guide is intended for identifying all post-project analysis steps involved in generating final data.

A Message for Contributors

We want to encourage active participation by all ToF-AMS users in the evolution of the information contained within this wiki and welcome the addition of content that is beneficial to the community as a whole. However, please DO NOT delete any content from this page!! Significant time, effort, and deliberation has gone into the information contained in this page. Rather than deleting content, please feel free to voice your concerns by posting a comment to the discussion page where others can contribute (please be sure to include a topic to be referenced in responses).

ToF-AMS Process Flowchart

ToF-AMS Analysis Flowchart.png

Process Calibrations

Determine Flow Rate Calibration

Basic Principle: AMS flow rate is used to convert MS signal intensity to units of mass concentration. Flow rate is derived from lens pressure (analog input #3) and converted using a conversion factor. Post-correction of flow rate MUST be performed if there is:
  1. any deviation between the derived and calibrated flow rates;
  2. errors in the logged flow rate value;
  3. electronic noise, or;
  4. drift in the zero of the baratron.
In the absence of errors, flow correction is made by multiplying the logged flow rate by the ratio of derived to calibrated flow rates. To correct for errors, the calibrated flow rate must be used.
NOTE: Save backup of flow rate wave (located in root:diagnostics) and construct new, corrected flow rate wave in Squirrel data analysis experiment. This new wave MUST be named “flow rate”.
Data Recording: A flow calibration MUST be performed at least once in the field using a Gilibrator.
How Calculated:
Flow ratecorrected = flow ratederived x (flow ratecalibrated/flow ratederived)
How Applied: Flow rate is used in the calculation of Squirrel component mass concentrations:
Squirrel Code:
Function: SQ_MS_Concentration
ugfac=ug_op?((aw/Navo)*1e12/(ioneff_wave[p1]*flowrate[p1]*rie*ce)):(1)
Recommendations:
Additional Resource: See Q-AMS manual for more information on correcting flow rate.

Determine Time Line of Single Ion Values

Basic Principle: The single ion (SI) value is used to convert MS signal intensity to units of mass conc. SI values logged by the data acquisition software may not be the true values. Unlike threshold values, single ion values can (and in many cases must) be post-corrected by applying a new set of values prior to pre-processing the data.
Data recording: It is recommended that the threshold analysis be performed daily to have adequate data for timeline of SI values.
How Modified: Single ion values should be obtained in the field using the threshold analysis window within the data acquisition software. If post-correction of the single ion values is required, this change can be made by clicking the Modify SI button in the “Optional Checks & Advanced Features” (HDF Index tab).
NOTE: Any change in single ion values MUST be made prior to pre-processing.
How Applied: SI value is applied in Squirrel to convert bit*ns to ions/s (Hz)
Squirrel code:
Function: sq_Hz(index_pos,datawav)
wave ion_wave=root:diagnostics:ionSingleStr_logged
datawav*=(1/ion_wave[series_index[index_pos[p]]]) * ToFPulserInHz[series_index[index_pos[p]]]
Recommendations:
  • Add original and modified SI time series to the MDP.

Process all IE/RIE_NH4 Calibrations

Basic Principle: The ionization efficiency (IE) calibration determines the ionization efficiency of ammonium nitrate. Quantification of all non-refractory AMS components are based on the linearity of this IE (i.e., the IE for all species can be related to the measured IE of a single species, in our case NO3).
Data Recording: It is recommended that all IE calibrations be recorded periodically during sampling with the data acquisition software running in GenAlt Mode.
How Derived: ToF-AMS IE calibrations are typically performed in the field using BFSP mode. Software is available to analyze these IE calibrations at the ToF-AMS Resource Page. A tutorial regarding BFSP IE calibrations is also available here. This software should be used to determine IENO3 and RIE_NH4 values for each field IE calibration. RIE_NH4 should be double checked using MS mode data.
How Applied: Calibrated ion efficiency is applied in the calculation of Squirrel component mass concentrations.
NOTE: Airbeam correction also influences the IE and will be discussed later (add link).
Recommendations:
  • Add IE calibration data (IE, AB, and flow rate) to the MDP.
Additional Resource: See Q-AMS manual for more information on IE calibrations.

Evaluate SI, IE, RIE_NH4 to use, based on SI, IE, AB, flow, and IE/AB for the whole campaign

Basic Principle: Based on results of threshold analyses and IE calibrations performed during the campaign select single ion (SI), ionization efficiency (IE), and ammonium relative ionization efficiency (RIE_NH4) values for the entire campaign.
How Derived:
  • Ideally, these values are simply the averages calculated from all evaluations (i.e., threshold analyses and IE calibrations) performed during the campaign.
  • However, it is likely that the time lines of these values will need custom corrections. This is particularly true of the time series of the time line of the IE value.
  • A critical parameter used to finalize the IE time series is the (IE/AB)Correction value.
How Applied: After deciding on correct values:
  • SI values can be modified using the Modify SI button in the “Optional Checks & Advanced Features” (found in the "HDF Index" tab)
  • IE values:
  • if no correction needed, the average IE will be used during the AB correction.
  • When doing AB correction, “ioneff” is written to single constant value.
  • If autoset is not used when performing an airbeam correction, then “ioneff_logged” is never used.
  • If correction needed, see Discussion
  • for Squirrel versions earlier than 1.38, use a post-correction function
  • For Squirrel versions 1.38 and later, use “ioneff_to_AB_ratio_rel” and “user_corr_fac”.
  • RIE_NH4 value should be used to modify the batch table.
Recommendation:
  • Add final SI to the MDP along with logged and user corrected IE values.
NOTE: When modifying SI and IE values, “*_logged” waves should never be changed as these changes will be overwritten whenever you perform a squirrel diagnostics.
IE/AB Correction
Basic Principle: The ToF-AMS ionization efficiency (IE) measures the sensitivity of the instrument to a measured quantity of a specific compound, typically NH4NO3, during instrument calibration. Throughout sampling, this sensitivity is reflected by changes in the instrument air beam (AB). In most cases, changes in sensitivity are observed in the AB and will be accounted for when performing an AB correction. However, there are instances where this may not the case and a custom correction is needed. The parameter that can be used to evaluate the need for a custom correction is the IE/AB ratio ((IE/AB)calibration) obtained from each IE calibration. If the (IE/AB)calibration for calibrations performed throughout the sampling period do not change significantly an IE/AB correction is not needed. If these values are not constant, a correction is required.
In Squirrel versions 1.38 and greater, an IE/AB specific correction wave “ioneff_to_ab_rel_ratio” has been implemented to aid the user. However, in Squirrel versions prior to 1.38, this wave must be created and the correction must be implemented manually.
Application of IE/AB Correction
Basic Principle: The intended use of the “ioneff_to_ab_rel_ratio” wave is to account for circumstances where the (IE/AB)calibration changes over time. In Squirrel versions 1.38 and higher, this wave is present by default with values set to 1. For earlier Squirrel versions, this wave will need to be created. In either case, “ioneff_to_ab_rel_ratio” MUST be modified to account for changes in (IE/AB)calibration values.
How Calculated: The values in the “ioneff_to_ab_rel_ratio” wave can be calculated by first establishing a correct, initial IE/AB value (IE/AB)t=0 (typically from an early IE calibration). The value of the required correction factor at each subsequent timestep, i, is then calculated by dividing (IE/AB)t=0 by the IE/AB ratio at this time (IE/AB)i:


ioneff_to_ab_rel_ratio[i] = (IE/AB)t=0 / (IE/AB)[i]

A time series of (IE/AB)[i] can be constructed by interpolating linearly between the (IE/AB)calibration values from each IE calibration.
How Applied:
  • In Squirrel versions 1.38 and above “ioneff_to_ab_rel_ratio” is automatically applied to adjust the species mass concentrations through the SQ_MS_Concentration function.
  • In Squirrel versions prior to 1.38 a post-correction function should be used to apply a corrected IE time series (that has had the IE/AB correction factor applied), through the following convention:
ioneff_corr = ioneff * ioneff_to_ab_rel_ratio
Recommendations:
  • Add the IE/AB corrections wave, if created, to the MDP.

Process size calibrations

Basic Principle: The parameters contained in the data acquisition software to describe particle time of flight may not be correct and need to be re-evaluated using the polystyrene latex (PSL) particle size calibrations performed in the field.
Data Recording: It is recommended that size calibrations be performed during the field deployment using 100, 300, 500 , 700, and 900 nm PSL particles. Sizes at or below 100 nm can be obtained with the classifier after it is calibrated at 100 nm.
How Derived: The size calibration is calculated by fitting a plot of particle size (Dva,PSL) vs. average particle velocity (Vpart, in m/s) using the following fit function:
Vpart = f(Dva,PSL) = Vl + (Vg - Vl) / (1 + (Dva,PSL/D*) ^ b)
where:
  • Dva = vacuum aerodynamic diameter;
  • Vl = lens velocity, and;
  • Vg = gas velocity
NOTES:
1. The time component of the Vpart calculation is obtained using the left (leading) edge of the ToF distribution. For choppers with a duty cycle >=2% and for calibrations with PSL, one-half of the chopper width should be added to the edge time. One-half chopper width is calculated as (0.5*duty cycle)/chopper frequency (e.g., 71.4 us for a 2% chopper at 140 Hz).
2. The distance component of the Vpart calculation is the length of particle flight chamber. For the AMS, these are:
  • 0.293 m (short chamber);
  • 0.39 m (long chamber).
3. To convert nominal diameter to vacuum aerodynamic diameter:
Dva = dp * p * S
where:
  • Dp = nominal diameter
  • S = Jayne shape factor (S = 1 for PSL), and
  • p = density (p = 1.05 for PSL)
4. Other fit functions can be used as long as that fit interpolates the data points.
5. For the default fit equation, initial guesses for each coefficient are as follows:
Vl = (10-25)
Vg = (450-650)
D* = 20
b = 0.50
In some cases, the lens velocity resulting from the fit can be negative. Under these circumstances, either Vg or Vl should be held constant within the range of values suggested above. It should be noted, however, that constaining Vg within the suggested range can still result in a negative Vl.
6. Periods with different flowrates (changes in orifice, clogging) need to have different calibrations.
How Applied: Parameters obtained by fit will be used to recalculate aerodynamic diameters when performing corrections in Squirrel.

Process Unit Resolution Data

Load data into Squirrel, Inspect Diagnostics

NOTE: Prior to loading data in squirrel, it is strongly recommended that the user review the General and Technical FAQs listed on the Squirrel main page as well as the several tutorials that are available on the Squirrel Help page

Squirrel requires AMS data to be saved in .HDF format. Under normal operating conditions, data is saved automatically by the data acquisition software in both .ITX and .HDF format. If .HDF files were not saved during sampling, the user can use the HieDI file converter to convert the .ITX files.

To load .HDF data into Squirrel:

  • Press “Get Index” button;
  • Browse to folder containing .HDF data files.

After indexing the .HDF data, Squirrel creates a diagnostics plot that can be used to diagnose the data set.

Implement Corrected Flow Rate and Single Ion Values, if any
Basic Principle: As discussed previously, incorrect flow rate and single ion values must be corrected prior to Squirrel data pre-processing because these values are used to convert AMS signal intensity to mass concentration during this step. Once data has been pre-processed, the single ion values are “set” and correcting any errors requires repeating data pre-processing. Therefore, it is strongly suggested that the user devote significant effort to ensure that the single ion values are correct and final before moving on. If no changes are made, pre-processing will use the flow rate and single ion values that were captured by the data acquisition software. These values are shown in the Squirrel diagnostics plot, which can be used as a reference.
Data Recording: N/A
How Derived: N/A
How Applied: Recommended ways of modifying the flow rate and single ion values prior to data pre-processing were presented previously. See previous flow rate and single ion value discussions for more details.
Recommendations:
Notes:
  1. Generally, it is recommended that the AMS flow calibration be expressed in terms of mass rather than volumetric flow.
  2. Reported mass concentrations are dependent on the calibrated instrument flow rate. Therefore, the conditions under which the instrument flow rate was calibrated in the field should be reported in all manuscripts (e.g., at STP or ambient and, if ambient, under what conditions).
Discard Bad Runs (from notebook, large noise, or low sampling time)

Bad runs can sometimes produce huge concentrations that are not real. Often new users forget to remove those runs, and their averages, scatter plots, etc. get hugely distorted. Reasons for bad runs include the inlet valve that was closed during part or all the run, cell phone noise, some electronic switch (e.g. heater bias from manual to computer) being changed during the run, and many others.

To discard runs, just type the run numbers in the appropriate box in Squirrel and hit "discard." Try to be thorough, and remove all runs identified as bad in the notebook. If you are not sure whether a run is bad, look at its mass spectra and size distributions, and it generally becomes clear whether the data is good or whether it has some major problems.

Perform m/z Calibration

Basic Principle: Under normal operating conditions, the user calibrates the instrument via the data acquisition software during instrument start-up. Afterwards, the ToF-AMS data acquisition software automatically calibrates m/z during autosave operation to account for drift in the calibration over time. Squirrel m/z calibration SHOULD be performed as a first step in ToF-AMS data analysis to ensure correct assignment of peak location in order to convert the raw data to stick format. The m/z calibration has to be performed prior to data pre-processing in order for the Squirrel fitting parameter values to be available as an option during pre-processing.
Data Recording: Initial m/z calibrations are required after restarting EITHER the instrument OR the acquisition software.
How Derived: N/A
How Applied: The m/z fitting procedure is initiated by clicking the “Check m/z calibration” button and then the “Begin mz fitting” button after the SQ_MzCalibration_Panel window appears.
Recommendations:
  1. For most fragments, use closed spectrum for m/z fitting (exceptions are air peaks (N+, O+, and Ar) for which open spectra can be used);
  2. C+, O+, and the Tungsten series have been found to be problematic for some data sets. In this case, these fragments should be removed when performing the m/z calibration;
  3. It is recommended that the user have at least one peak greater than m/z 90;
  4. At present, the calibration as performed cannot reliably be extrapolated to masses larger than m/z 200. This may be fixed in future Squirrel releases. As a result, if the exclusive focus of the analysis is large mass fragments, the user should remove all low mass fragments from the mass table prior to performing the mass calibration.
Mzcalibpanel.png

1. Using the inset table, the user can modify the fragments used in the m/z calibration prior to beginning m/z fitting. To remove a fragment from the fitting procedure, place a zero in the corresponding cell in the SpecType column. The table can be reset to default values by clicking the “Revert to Default” button. Peaks from the closed spectra are typically better with the exception of the air peaks, for which open spectra can be used.
2. After selecting the fragments to include in the fitting procedure, the procedure is started by clicking “Begin mz fitting”.
3. The accuracy of the m/z calibration and mass calibration for each run are displayed during m/z fitting. The accuracy parameter evaluates the “goodness of fit” of the Squirrel calibration relative to that of the acquisition program.
4. These plots display the results of m/z fitting from each run. The raw data (black dots) and the Squirrel fit (red line) are plotted along with the designated peak center determined from the data acquisition (black vertical line) and Squirrel (yellow vertical line) m/z calibrations.
5. The overall goodness of the Squirrel m/z calibration, relative to the calibration from the data acquisition program, is reported (ratios (data acquisition/Squirrel) for intercept, slope, and power). Deviation from unity represents a “bad” fit by the Squirrel m/z calibration.
6. “Bad” fits can be “filtered” after the m/z fitting procedure, if needed. The threshold criteria for filtering is the deviation of the ratios found in (4) from unity. The value of this threshold is at the discretion of the user as is the method used to fill in the bad fit value.

Perform Baseline Subtraction

Basic Principle: ToF-AMS mass concentrations are derived from the difference spectra. Both raw open and closed spectra can have a significant non-zero baseline that is transferred to the raw difference spectrum as well. The “check baseline” routine allows the user to examine and make corrections to the baseline of the difference spectrum. This correction ensures that any baseline contribution to the stick value for each m/z fragment is removed prior to calculating the difference stick spectra.
Data Recording: N/A
How Derived: N/A
How Applied: The check baseline procedure is started by clicking the “Check baseline” button, which will open the SQ_Baseline_Panel.
Recommendations:
  1. Use mass defect wave as peak center and remove additional regions from stick complement (“Advanced Options” tab).
  2. Is there an objective way to say whether the baseline correction is good enough? (numerical estimate of corrected vs uncorrected)

Pre-process data

Perform air beam correction

Basic Principle: As discussed above, the single ion strength (SI) and ionization efficiency (IE) is a fundamental quantities that need to be known in order report AMS mass concentrations. While it is not possible to continuously measure the SI and IE throughout the course of a sampling period, the air beam (AB) provides an indirect measure because of it's dependence on both. For this reason, it is important to note the AB as well as the IE when performing calibrations.
The AB correction is used to account for changes in SI and/or IE throughout the course of sampling relative to a designated reference period. The AB correction factor (rer. corr_fact in Squirrel) compensates for changes in the AB over time by adjusting AMS mass concentrations as follows:
component conc.[i,ABcorrected] = component[i] * Corr_Fact

The AB correction factor is calculated for each run in the sampling period as the ratio of the AB at the reference period to that measured at each step, i. Typically, reference periods are chosen immediately FOLLOWING an IE calibration (where the precise IE and its corresponding AB are well known).
If the IE/AB ratio does not remain constant throughout the sampling period, the AB correction will NOT produce the correct answers. In addition to the AB correction, a seperate correction is required in this case (see this wiki, section 1.4.1).
How Calculated:
Corr_Fact = AB[reference]/AB[i]
Squirrel code:
Function: SQ_MS_Airbeam()
corr_fact=(j==detector_num[p])?(airbeamCorrToFTypeValue[1][j]/airbeam[p]):(corr_fact[p])
How Applied:
component conc.[i,ABcorrected] = component[i] * Corr_Fact
Squirrel Code:
Function: SQ_MS_Concentration()
ugfac_wav_tmp*=corr_fact[p]
Performing AB correction in Squirrel: The AB correction can be found under the Corrections:Airbeam tab. Two options are available in Squirrel when calculating the AB correction: auto-set and non-auto-set.
Active auto-set
When the auto-set option is active, the reference period is selected by the user with the Squirrel marquee function (ref. Squirrel tutorials). Upon selecting an AB reference period, that range is automatically entered into the squirrel AB corrections panel and average values for the MS AB, IE, and flow rate from that range are calculated when performing the AB correction. It is important to use average values from the data set immediately following IE calibration instead of values from that calibration for consistency with data (e.g., the pressure in the inlet during calibration may be different from that observed during sampling which will cause the IE calibration to have a different AB value).
Inactive auto-set
If the auto-set option is not active, the MS AB, IE, and flow rate values from the reference period MUST be manually entered by the user (for the high resolution ToF-AMS, values for BOTH V and W modes must be entered paying special attention to the ToF/ion type).
Recommendations:
  • Add corr_fact and corr_fact_flow to the MDP.
NOTE: It is IMPORTANT to note that, when the AB correction is perfomed, the IE value from the reference period REPLACES reference IE (wave name = ioneff) values in the Squirrel experiment if they exist. Therefore, if any custom modifications have been made to the IE time series prior to performing the AB correction, the user MUST reapply these modifications via a post correction function (add link to future section in this wiki).
Correct for Flow
Basic Principle: The objective of the AB correction is to compensate for all changes in instrument sensitivity relative to the reference time period. While the majority of these AB changes can be attributed to changes in the IE and/or SI, the AB can also change with changes in flow rate. If the AB correction is applied before removing the flow rate changes errors are introduced into the data. Removing flow rate changes from the AB correction is accomplished by applying the "Correct AB for flowrate" option within the Corrections:Airbeam panel. (add screen shot of the flow rate box in the Squirrel panel).
NOTE: The flow rate is ALWAYS applied to calculate mass concentrations regardless of whether this correction is applied.
How Calculated:
Corr_fact_flow =  (Flowrate[i] + Flowrate offset))/(Flowrate[reference] + Flowrate offset)

Squirrel Code:
Function: SQ_MS_Airbeam()
If no smoothed flowrate:
corr_fact_flow=(flowrate+airbeam_ref_flowOffset_sv)/(airbeam_ref_flow_sv+airbeam_ref_flowOffset_sv)
If smoothed flowrate:
corr_fact_flow=(flowrate_smoothed+airbeam_ref_flowOffset_sv)/(airbeam_ref_flow_sv+airbeam_ref_flowOffset_sv)
How Applied:
Corr_fact = corr_fact*corr_fact_flow
Flowrate offset
Basic Principle: In order to properly apply the flow rate correction, the flow rate offset must be determined. The flow rate offset is that flow rate at which the AB is zero. The flow rate offset is instrument specific and a range of offsets are possible between instruments. However, instrument specific flow rate offsets should not change drastically. Therefore, it is recommended that this offset be determined experimentally at least once for each instrument.

NOTE: If this offset is not applied, you will have an error in your determined mass concentrations as large as the percent change in the flow rate.

How Calculated:
If a significant change in flow rate exists in a data set:
1. find region where flow rate changes over time (e.g., region where orifice is clogging);
2. for only that region, make a scatter plot of flow rate vs. AB;
3. do a linear regression from this scatter plot with a non-zero intercept (alternatively, a line can be drawn by eye);
4. the flow rate offset is equal to the value of the y-intercept.
Once the flowrate offset has been determined, enter this value in the flow rate offset cell.
NOTES/RECOMMENDATIONS:
1. If the flow rate offset cannot be determined from data existing in the data set, this value must be determined experimentally (use 3 or 4 different orifice sizes or a needle valve on the front of the inlet and only choose those values that lie within the range of values that are relevant to the data set).

Make campaign specific adjustments to fragmentation table

The purpose of the fragmentation table is to apportion the ToF-AMS signal among the different chemical species present in the aerosol (ammonium, chloride, nitrate, organics, and sulfate). The following discussion assumes that the duty cycle correction has been applied to the data (this routinely occurs during data pre-processing in Squirrel unless that function has specifically been disabled by the user) and that the data has been processed to the point of calculating the air beam correction.

The fundamentals and an example (for NH4+) of the frag. table are presented in Allan et al., J. Aerosol Sci., 2004.

Below are detailed discussions and suggestions for modifying the ToF-AMS fragmentation table.

Air Fragmentation Patterns

The order of these steps is important, you should get the air corrections right first (especially frag_O16[16] before doing the water corrections in the next step.

Historically, only frag_O16[16] has been modified to adjust the air fragmentation pattern. Currently, modifications to the air fragmentation patterns at m/z 29 (at frag_air[29]) and m/z 44 (at frag_co2[44]) have also been required for ToF-AMS data. As with the modification to frag_O16[16], these adjustments require the user to look at the DIFFERENCE spectra during either filter periods (preferable) or otherwise periods with very low signal.

It is important to note that the frag_CO2[44] coefficient is a combination of a number of factors including:

1. the mixing ratio of CO2 in air (0.00037);
2. the RIE of CO2 from the literature (1.36);
3. the reciprocal of the fraction nitrogen in air (1.28), and;
4. a correction for mz15 fragmentation of nitrogen (1.14).

If a change is required in the frag_CO2[44] term, it is recommended that the RIE term be modified since the other terms are likely to be constant.

Modifications to the frag_O16[16], frag_air[29], and frag_CO2[44] coefficients can be made in one of two ways. In the first method, the user can examine the filter period average spectrum and manually modify each respective coefficient until the contribution of the corresponding fragment or component ([NH4] in the case of frag_O16[16]) in the mass spectrum is zero.
Alternatively, ratios for mz16:mz14, mz29:mz28, and mz44:mz28 can either be calculated by linear regression or by calculating the ratio for the entire set of filter periods and then calculating the average ratio across all filter periods. After calculating these ratios, the default coefficients are modified as follows:

  • The measured mz16:mz14 ratio is used to modify the default coefficient at frag_O16[16];
  • The measured mz29:mz28 ratio is used to modify the default coefficient at frag_air[29];
  • The measured mz44:mz28 ratio is used to modify the default coefficient at frag_CO2[44] (in order to modify only the RIE term of the frag_CO2[44] coefficient, multiply 1.36*(measured mz44:mz28 ratio/0.000734) and use this number to replace the default RIE value).


NOTES:

1. As seen in the water fragmentation pattern, water can contribute to m/z 16 if the air sampled into the instrument during filter sampling is not dry. If this is found to be the case, the water contribution should be subtracted using the water fragmentation pattern along with m/z 18.
Water Fragmentation Pattern

Water typically dominates the closed spectrum so the signals at m/z's 16, 17, and 18 can be used to determine the fragmentation pattern for water.

Begin by looking at the CLOSED signals for these ions, plot m/z 16 and m/z 17 versus m/z 18, do a linear regression through zero for each plot to find the slopes, and compare those ratios with the default frag table coefficients. If significantly different from default values, the slopes of the mz16:mz18 and mz17:mz18 plots are used as user-defined coefficients at frag_RH-, frag_SO3-, frag_organic-, and frag_water[16] and [17], respectively.

NOTES:

1. mz17:mz18 is fairly stable over a wide range of conditions and is independent of heater temperature or pumping time (slight change may be due to varying ToF-AMS conditions);
2. In data where the m/z 16 signal is low due to low water in the background (long pumping times), the signal at m/z 16 is constant relative to m/z 18 and the slope becomes flat. For those cases, the current recommendation is to use the default value;
3. If possible, look at a high resolution spectrum at m/z's 16, 17, and 18 to ensure that there is no contribution from NH4 to these masses.
4. The water fragmention pattern is used in four places in the frag table and the coefficients MUST be changed at each fragment when altering the default coefficient values .
Other Fragmentation Table Adjustments
Modifying frag_NH4_17[17] and frag_NH4[17]

The ratio of NH4_17:NH4_16 has been measured in laboratory studies using a Q-AMS to be ~1.1:1. However, different NH4 souces in the aerosol (such as amines) can have different ratios and different AMS conditions could also change the NH4 fragmentation pattern. To check this ratio plot NH4_17 vs NH4_16 and do a linear regression through zero. Under certain circumstances, these fragmentation patterns can be modified, however this is NOT recommended.

In extreme circumstances ONLY (where NH4_17 to NH4_16 ratio is very far off from 1.1:1 and is noisy) the frag_NH4_17[17] and frag_NH4[17] can modified to force a 1.1:1 ratio, but again this is NOT recommended. To force this 1.1:1 ratio, the following changes can be made (after making these changes, make sure NH4 loadings during filter periods do not change drastically):

  • Change frag_NH4_17[17] to blank ("")
  • Change frag_NH4[17] to "1.1*frag_NH4[16],1*frag_NH4_17[17]"
Modifying Nitrate Fragment (frag_org[30])

The current value for Frag_Org[30] is based on isotopes of Frag_Org[29]. Some non-nitrogen-containing organic material could also be present at 30, making the isotope attribution too low for Frag_Org[30]. In some cases, the ratio between NO3_30 and NO3_46 can be forced to be the same as for ammonium nitrate. This would be reasonable if there are no other types of nitrate present, such as sodium, calcium, or organic nitrate which tend to have a higher 30/46. In the absence of high resolution data, the attribution of 30 should be something for the data analysis person to decide.

Modifying Sulfate Fragments

The fragmentation pattern for sulfate is not relevant for most analyses where Frag_H2SO4 and Frag_SO3 are not examined and only the combined Frag_sulfate is used. However, there are some instrumental differences for the amount of Frag_Sulfate[18].

The diagnostics ("frag checks") plot for SO4 should be plotted. Each scatter plot between the sulfate fragments (mz64:mz48, mz80:mz48, mz81:mz48, and mz98:mz48) should have linear relationships. Although the linearity of data is important, the exact slopes of lines are unique to each instrument (these slopes can also be compared with (NH4)2SO4 calibration if required). To check for the presence of organic sulfates, the correlation of any deviation from linearity with organics can be investigated.

Modification of the sulfate fragments can be made, but are not often necessary (the possible exception being when sulfate is low and organics are high). Any adjustments can be made at frag_org[80] and [81]. If there is evidence of MSA, organo-sulfites or organo-sulfates in the aerosol, the sulfate fragments will need to be adjusted.

Adjust Frag_Org[39]
  • The Squirrel default doesn't have any organic signal at m/z 39. This is because there is an interference from K+ at that m/z, which can be very large and variable in some datasets. However in many datasets, after checking the HR data in PIKA and verifying that K+ is negligible for your dataset, one can set Frag_Org[39] = 39.


Time Dependent RH
Time Dependent Gas Phase CO2


The chart below presents a quick reference guide to the frag table modifications discussed above.

ToF-AMS Fragmentation Table Modification Guide
Step Ultimate Target Adjustment required at Spectrum to use Diagnostic Default coefficient Comments
1 Water frag_RH[16,17],
frag_SO3[16,17],
frag_org[16,17],
frag_water[16,17]
Closed mz16:mz18 ratio,
mz17:mz18 ratio
0.04,
0.25
1. User defined coefficients are measured mz16:mz18 and mz17:mz18 ratios (if significantly different from default values).
2. When making a change to these coefficients, changes MUST be made at each of the four listed locations in the frag table.
3. HR spectrum can also be consulted to make sure there is no NH4 contributions at mz16,17.
2 NH4(I) frag_NH4[17] Difference (NH4_17:NH4_16) ratio 1*frag_NH4_17[17] 1. Changes in these fragmentation patterns is generally NOT recommended.
2. Under certain circumstances, such as when there is significant deviation from 1.1:1 ratio and if the signal is noisy, this ratio can be “forced” to 1.1:1 by:
- Changing frag_NH4_17[17] to blank (""), and;
- Changing frag_NH4[17] to "1.1*frag_NH4[16],1*frag_NH4_17[17]“.
3. Make sure NH4 loadings during filter periods do not change drastically.
3 NH4(II) frag_O[16] Difference from filter period [NH4] ~ 0 0.353 1. User defined coefficient should be evaluated such that [NH4] ~ 0.
2. An alternative method for calculating this coefficient is to calculate the mz16:mz14 ratio (after subtracting any water contribution using he water fragmentation pattern).
4 Organic (I) frag_air[29] Difference from filter period org_29 ~ 0 0.00736 1. User defined coefficient should be evaluated such that org_29 ~ 0.
2. An alternative method for calculating this coefficient is to calculate the mz29:mz28 ratio.
3. The default coefficient is based on isotopic ratios and needs to be changed for ToFs.
4. Purpose of this adjustment is to account for residual signal that is not captured in the baseline subtraction.
5 Organic (II) frag_CO2[44] Difference from filter period org_44 ~ 0 1.36 1. Looking at difference spectrum during filtered period, adjust this coefficient until org_44 is ~ 0.
2. An alternative method for calculating this coefficient is to calculate the mz44:mz28 ratio.
3. Value of 1.36 is RIE of CO2 from literature.
4. Also look at both diagnostics plots ("frag checks"):
- Look at intercept of m/z 44 to organics plot (should go through zero);
- Look at intercept of m/z 43 to m/z 44 plot (supposed to go through zero, but not as crucial as 44 vs org).
6 Nitrate frag_organic[30] All difference NO3_46:NO3_30 (ambient), NO3_46:NO3_30 (IE calibration) 0.022 1. Look at diagnostics plot ("frag checks") for NO3_46 vs NO3_30.
2. Check ambient ratio compared with ratio determined during NH4NO3 calibration.
3. If difference between ambient and cal 46 to 30 ratios, frag_organic[30] coefficient can be adjusted to match that from IE calibration.
4. However, this assumes that only NH4NO3 is present in the aerosol.
5. HR spectra can be consulted to make a better determination of this ratio.
6. In the absence of this HR data, the attribution of mz30 should be something for the data analysis person to decide based on available information.
7 Sulfate frag_organic[80],
frag_organic[81]
All difference SO4 “frag checks” See Comments 1. Each ratio (m/z’s 64, 80, 81, 98) should have linear relationships with m/z 48 (linearity of this data is important).
2. Exact slopes of lines are unique to each instrument.
3. Check for deviations from linearity that correlate with organics.
4. Sulfate fragments rarely need adjustment (perhaps only when organics are high and sulfate is low, see #7 below).
5. Ambient sulfate fragmentation patterns can also be compared with sulfate calibration.
6. If you have MSA, organo-sulfites or organo-sulfates, will need to adjust SO4 frags.
7. When organics >> SO4, organics will dominate the signal at many fragments typically associated with SO4 (e.g., m/z's 64, 80, 81, 98, etc.) and SO4 concentrations can be negative as a result. In this case, the signal intensity of these fragments should be calculated using both the signal intensity at m/z 48 (SO+, contribution of organics at this mass is rare) and the ratio of other characteristic SO4 fragments to m/z 48 obtained from an ammonium sulfate standard to calculate the intensity of larger fragments.
8 Water frag_RH[18],
frag_water[18]
      1. No consensus on whether or not to change frag_organic[18] (default is 1*frag_organic[44]).
2. May be better to use value < 1 for coefficient here.
3. There is no consensus yet on this.
4. SO4 frags should have informed user as to whether frag_sulfate[18] needed to be changed (most likely no).
If important, can be determined with lab calibration of AS in Argon.
5. These two things very important, but need more lab work to determine these better.
6. After organic and sulfate water accounted for, rest is then attributing water vapor between gas and particle phase.
7. Can inform this decision with measurement of RH from another instrument.
8. Can use pToF to see partitioning of gas phase and particulate phase water.



ToF-AMS Fragmentation Table, Selected User Modifications
Frag Default TEXAQS C-ToF MIRAGE V-ToF INTEX-B V-ToF Whistler, HR-ToF SOAR, HR-ToF SOAR, C-ToF
frag_RH[16],
frag_SO3[16],frag_org[16],
frag_water[16]
0.04 0.05       default default
frag_RH[17],
frag_SO3[17],frag_org[17],
frag_water[17]
0.25 0.28       default 0.20549
frag_NH4_17[17],
frag_NH4[17]
  no change,
no change
      no change,
no change
no change,
no change
frag_O[16] 0.353 0.361 0.353 0.3614 0.313 0.3535 0.3565
frag_air[29] 0.00736 0.0049 0.0057 0.0062   0.00670 0.00156
frag_CO2[44] 1.36 1.45 1.15 1.60   2.32 0.36
frag_organic[30] 0.022 no change   0.022   default default
frag_RH[18] 0.01 0.008   0.01   default default

Evaluate the fragmentation patterns of the main species ("Alice's Diagnostics")

Evaluate ammonium balance

Generate time series and average spectra of unit-resolution data

Choose preliminary CE from observed composition and accumulated AMS experience

The data processing steps described above produce a concentration of each species (as measured by the AMS) from first principles. E.g. 5.2 ug m-3 of sulfate at 4:00-4:05 PM. In the first studies using the AMS for ambient data (Jimenez et al., 2003; Drewnick et al., 2004; Alfarra et al., 2004) it was observed that the AMS concentration had the same structure in time as for other collocated instruments (PILS species, TEOM mass, SMPS volume), but that it was low. Comparing individual species such as sulfate it was found at the AMS concentrations needed to be multiplied by ~2 to match the concentrations measured by other quantitative instruments. This is normally expressed as a collection efficiency CE = 0.5.

Early on it was thought that CE < 1 arose from excessive broadening of the particle beam at the exit of the aerodynamic lens inside the AMS, leading to some particles not hitting the vaporizer and thus not being detected. This led to the development of an automated Beam Width Probe (BWP, Huffman et al. 2005). However the use of the BWP in many ambient studies (summarized in Salcedo et al., 2007) has shown that this is not a problem.

Tim Onasch et al. during the Ron Brown 2004 study showed convincingly (by comparing PToF and light scattering particle pulses) that the cause of CE < 1 was particle bounce at the vaporizer. Huffman et al. (AS&T 2005) expressed CE = E_b * E_L * E_s to clearly separate the different effects that may reduce the signal reaching the AMS vaporizer. Only E_b is due to particle bounce.

The current status (as of June 2007) of the E_b issue is as follows:

(a) Aerodyne has tried several new vaporizer designs in an attempt to reduce or eliminate bounce. However these have not been successful so far.

(b) For that reason, at present CE is estimated based on the observed composition and phase of the particles, and verified based on intercomparisons with other quantitative instruments.

Fortunately CE ~ 0.5 has been found to hold in most field studies, including those dominated by sulfate and/or oxygenated organics (OOA). Laboratory studies can be more variable depending on the particle phase and volatility used. Known deviations from CE ~ 0.5 include:

(a) Liquid organic particles or particles with a large liquid coating, which are collected with CE = 1. Note that OOA and solid organic particles do bounce (e.g. SOA measurements by Bahreini et al., 2005).

(b) Sulfate-dominated particles can be collected with CE > 0.5 when the particles are acidic, i.e. there is not enough NH4+ to neutralize all the sulfate. This has been seen in many field studies. Based on the Onasch et al. measurements in 2004, CE = 0.5 until the sulfate is in the form of bisulfate, at which point CE increases linearly and reaches CE = 1 for H2SO4.

(c) Particles dominated by ammonium nitrate (NH4NO3 > 50% of the mass). In this case CE also increases, which can be parameterized in the same way as (b), but with the X-axis being the NH4NO3 mass fraction.

(d) Particles with high water content that retain significant water in the AMS lens. This was pointed out by James Allan in his 2004 JGR paper, and also shown by Matthews et al. with laboratory particles (AST, 2008). However organic-dominated ambient particles do not seem to pick up enough water to drive their CE to 1 using this effect (e.g. per experiments from Tim Onasch and Tim Bates in the Ron Brown during ICARTT 2004). For this reason this method has not been pursued, and the use of a dryer in front of the AMS is highly recommended to remove this potential composition-dependent effect on E_b.

(e) (a)-(d) above assume internally mixed particles, which is very often the case for ambient data. However if you see a particle size distribution mode that is clearly externally mixed, it will have its own CE. We have seen this sometimes with NH4NO3 condensation into small particles in the early morning in the Eastern US, when those particles can have E_b ~1 even though the accumulation mode has E_b ~ 0.5. It can also happen when a fresh power plant or volcanic plume with lots of acidic sulfate is embedded in a more aged aerosol, and this has been seen in multiple aircraft experiments.

IMPORTANT NOTE: disagreements between concentrations from the AMS and other instruments can arise for many other reasons than variations in Eb. For example the AMS user may have made mistakes in its inlet setup or lens alignment (causing large particle losses), or the threshold setting, the IE calibrations or IE/AB corrections, or in determining the single ion, or in other AMS data processing procedures. And even more often in our experience, other instruments have problems (of functioning, software, or inlets) that may not be noticed, especially standard instruments that don't receive much attention during a campaign (such as an SMPS or OPC). We have seen every type of quantitative instrument produce wrong concentrations at least once in the field (including PILS, TEOM, etc.). If you think you are seeing a CE that is changing all over the place, almost certainly the problem is somewhere else in the AMS or other instruments. The CE / E_b SHOULD NOT BE THE WASTEBASKET OF ALL DISAGREEMENTS. Rather CE should be chosen based on the composition-based rules above, and remaining disagreements should be investigated by looking at all aspects of the AMS and of other instruments.

Look a MS and size distributions of very large outliers, blacklist if appropriate

Work up comparisons of total mass and volume with other instruments

Make size distributions for a handful of periods with good signal

Work up comparisons of size distributions for those same periods with other instruments

Reevaluate CE choice based on intercomparisons

Start looking at the scientific questions of interest

Mission Diagnostics Plot

The mission diagnostics plot (MDP)is useful to diagnose both instrument performance during a field campaign and also to more easily determine where various corrections to the data are required. Throughout the above discussions, the user is prompted to add data to this plot at various steps within the analysis process. It is STRONGLY suggested that a MDP be produced for every data set. In Squirrel versions 1.39 (?) and above, this plot is created using various user interfaces within the data analysis software. However, for analysis conducted in earlier versions of the software, this plot must be produced manually.

Sample image of mission diagnostic plot (produced from SOAR1 data). In this plot, calibration data are represented by circles and Squirrel data are represented by lines.

MissionDiagnosticsPlot.png

Determining Precision Directly from Field Data

See the following procedure adopted by NASA as part of the Tropospheric Airborne Measurement Evaluation Panel (TabMEP) process.

Reporting your data

At a gathering of many European AMS users in May, 2009, it was suggested that users routinely report their data with supplemental information about the measurement configuration and analysis. The list of suggested supplemental fields is:

AMS_Type: C, V, or W, and ionization type if not EI
AMS_SamplingTime: length of time of typical run, and any other info (2 minute data saving alternating V and W mode)
AMS_LensType: standard, high throughput, etc
AMS_LensPressure: typical lens pressure
AMS_FlowRate: typical flow rate
AMS_PrsTempFlowCal: P & T of room the AMS was in during the flowrate calibration, i.e. 1 atm, 30 C, 20 cc/min
AMS_STPConversionFactor: conversion factor for STP (273K 1 atm)
AMS_IEOverAB: description of IE/AB ratio used for data processing (can be a single value or several) If using a single value, report the IE/AB ratio as calculated from the IE and AB values used in the airbeam correction.
AMS_RIENH4: description of the RIE used in data processing for the Ammonium calibration
AMS_CE: description of the CE value used (e.g. single value or based on acidity etc.)
AMS_PMCut: PM cut on inlet e.g. Size cut for inlet line at 2.5 micrometers
AMS_InletRH: inlet dried or ambient (how?) Inlet RH < 40% using nafion drier
AMS_LocalToUTCTime: for ground-based station, time relative to UTC e.g UTC + 2h
AMS_RelatedInstrumenationInfo: i.e. other instruments that shared the inlet

Users who are required to submit data using the ICARTT text format can use one of two Igor tools to help with this. The first tool was written by James Allen and uses one ipf: ICARTT_AMS_FormattingTool_V2_0.ipf. There is a second tool written mostly by Donna Sueper that uses a pair of ipfs: 2010GeneralMacros.ipf and 2010IcarttFormat.ipf. This second tool has a more detailed interface and is useful for those wanting more guidance. Information about the ICARTT file format can be found at http://esrl.noaa.gov/csd/tropchem/ICARTT_data_mgt15.pdf.