Data

Investigation of evaporation and biodegradation of fuel spills in Antarctica - a chemical approach using GC-FID

Australian Antarctic Data Centre
SNAPE, IAN ; HARVEY, PAUL ; FERGUSON, SUSAN ; RAYNER, JOHN LENDON ; REVILL, ANDREW
Viewed: [[ro.stat.viewed]] Cited: [[ro.stat.cited]] Accessed: [[ro.stat.accessed]]
ctx_ver=Z39.88-2004&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Adc&rfr_id=info%3Asid%2FANDS&rft_id=info:doi10.4225/15/5afe71576cb33&rft.title=Investigation of evaporation and biodegradation of fuel spills in Antarctica - a chemical approach using GC-FID&rft.identifier=10.4225/15/5afe71576cb33&rft.publisher=Australian Antarctic Data Centre&rft.description=Evaporation Model for hydrocarbon spills. Developed by the Australian Antarctic Division to simulate fractionation of Special Antarctic Blend (SAB) and other diesel range fuels during evaporation. Text version of notes for excel model, please read me. This Package of files includes a pdf of the scientific paper, a readme word document and 3 excel files. The purpose of the 3 excel files are as follows: Excel File 1. Evap Model_V1_single temperature.xls. Evaporation predictions with a single sample at one temperature. This file requires input of initial composition, composition after weathering and a single temperature. Excel File 2. Evap Model_V1_five temperatures.xls. As above but with the input of 5 different temperatures. Excel File 3. Evap Model_development version, derivation of data and AAD examples.xls This is an earlier development version of the numerical model. This file is intended to allow others to see how available Antoines equation parameters were originally fitted for use with other R+UCM regions. This file includes a worksheet where the fitting parameters for Antoines eqn were calculated, a necessary task to be able to extrapolate to compounds for which explicit vapour pressure data are not available. The data input sheet has some differences to the singe temperature and five temperature model as the R+UCM regions are split up differently. Also the assumptions about how the different classes of compound (i.e. the Aliphatic and aromatic classes) behave can be altered. Some raw data from the AAD evaporation experiment is included and plotted. The Temperature comparisons worksheet summarises how ratios of interest change and allows plotting of one ratio at different temperatures (see worksheet Figs 2 and 3). Read me about background to this excel file Background Notes This excel file estimates the relative evaporation rates of different hydrocarbons from a hydrocarbon mixture (i.e. a fuel). In this excel file, a single temperature is considered for an evaporating hydrocarbon mixture. The desired temperature and details about the hydrocarbon mixture are entered in the Main Input page. The model was developed for use with the bulk fuels used by the Australian Antarctic Division at Casey, Davis, Mawson and Macquarie Island. This fuel, Special Antarctic Blend (SAB) starts at C9. As some of the spill sites undergoing remediation are mixed with heavier lube range hydrocarbons the options for inputs go to C36 - a typical maximum when Total Petroleum Hydrocarbon analysis is undertaken. The model uses thermodynamic data where such data are available. Estimated thermodynamic data are used for components when specific data are unavailable. Other important information can be found in the sections listed below. References GC-FID data used in the model Temperature and Vapour pressure estimation Other corrections How the calculations are done Experimental data (that supports the approach taken in this excel file) Read me about References Notes about the scientific publication that this model is reported in. Paper can be obtained as a pdf file from the Australian Antarctic Data Centre. The paper contains details on the evaporation experiments at +20 degrees C and -20 degrees C. The results from the experiments agree well with the calculated fractionation rates that this Excel file produces. Read me about GC-FID data for this model Notes about GC-FID data used in this excel model. This model was set up to directly use GC-FID data outputs for each of the compounds of GC regions listed. Consistent units are required for the areas of each measured region (and GC bias needs to be low across the range of fuel components). Note that the summation of all the regions and compounds = total area identified in the chromatogram from C9 up to C36. This range covers the observed range of components in SAB, Arctic Blend Diesel and lubes that may have been spilled in the same area. Read me about temperature / vapour pressure estimation Notes about Temperature correction method Raoults law and the Antoine equation are used to calculate the composition of the evaporating portion of the fuel mixture. Vapour pressure data were obtained for a range of available hydrocarbons C9 and larger in the temperature ranges that covered the site temperature ranges, approx -20 degrees C to +20 degrees C. The available data were mostly limited to n-alkanes. This situation was exacerbated because other compounds of interest are solids at these temperatures when pure. Consequently vapour pressure data are not available for these components in the liquid form at these temperatures. The available n-alkane vapour pressure data were combined and a best fit of these data were determined as a function of effective Carbon Number (ECN). This allowed the estimation of the vapour pressure of fuel components with ECN's between the n-alkane ECN's. As a refinement each region of the chromatogram was split into 5 classes - 3 x aliphatic fractions and 2 x aromatic fractions. The ECN used to estimate Vapour pressure of each class was slightly modified from the average retention time relative to an n-alkane. This modification was carried out to correct for specific GC-column / compound interactions - interactions that increase from acyclic to cyclic to polycyclic to hindered aromatics [i.e. multiply alkylated] to unhindered aromatics. Examination of the evaporation behaviour of n-alkanes vs. the nearest ECN regions confirms the need for the correction with the calculated evaporation profiles better matching the observed profiles. When an n-alkane's evaporation rate is compared to a GC region with distinctly different vapour pressures these corrections to ECN make little or no difference to the predicted selectivity during evaporation. Read me about other corrections Notes about other corrections used in this model. This model does not include systems that are under diffusional control (i.e. limited by diffusion of hydrocarbons within the soil / fuel mixture). The assumptions in the model are for an evenly mixed liquid that is evaporating. The experiment was set up to avoid this problem (by rotating the flasks). Soil with evaporating fuels may well be affected by this and other problems. Read me about how calculations are done Notes about how the calculation is done The fuel is divided into a number of fuel classes and specific compounds are identified. These classes are listed in the Main Input Page ready for use with the example data or other fuel data in the appropriate fuel range. R+UCM stands for Resolved + Unresolved Complex Mix but needs to be calculated excluding the specifically identified compounds like n-alkanes. Specifically identified compounds are excluded from the R+UCM so they are not double counted. Mol fraction of each component is estimated. For specific compounds an exact molecular mass is known. For other R+UCM classes molecular mass is estimated from known compounds in that region. Vapour pressure of the pure component or region is estimated with Antoines equation for the temperature required. This pure vapour pressure estimate is combined with liquid phase mol fraction to calculate the gas phase composition at each evaporation step. The evaporating portion (i.e. the gas phase portion) is removed from the liquid phase portion. The Excel sheet is setup such that this subtraction accounts for approximately 1% of the initial fuel amount. The need to calculated mol fractions then back-calculated mass remaining is the reason it is not exactly 1% by mass at each step. To avoid numerical errors, division by zero errors and rounding errors many of the calculations contain and the IF formula. When a fuel component is greater than 0 the mol fractions are calculated, otherwise a value of 0 is returned. Read me about AAD experimental data Notes about AAD experiments that back up the model (a full description is in the scientific publication) Small portions of Special Antarctic Blend (SAB) fuel were placed into vials. Each vial was placed into a +20 degrees C or -20 degrees C chamber. A slow stream of nitrogen passed into the top of each vial to remove the evaporating portion of the fuel. The vials were slowly rotated to ensure even mixing of the residual fuel during rotation. Periodically a vial was removed, weighed to calculate mass fuel evaporated, and analysed with GC-FID apparatus. After a range of vials were analysed at different levels of evaporation at the 2 temperatures a data set was obtained to validate the numerical model.&rft.creator=SNAPE, IAN &rft.creator=HARVEY, PAUL &rft.creator=FERGUSON, SUSAN &rft.creator=RAYNER, JOHN LENDON &rft.creator=REVILL, ANDREW &rft.date=2006&rft.coverage=northlimit=-53.0; southlimit=-70.0; westlimit=60.0; eastLimit=160.0; projection=WGS84&rft.coverage=northlimit=-53.0; southlimit=-70.0; westlimit=60.0; eastLimit=160.0; projection=WGS84&rft_rights=This data set conforms to the CCBY Attribution License (http://creativecommons.org/licenses/by/4.0/). Please follow instructions listed in the citation reference provided at http://data.aad.gov.au/aadc/metadata/citation.cfm?entry_id=fuel_evaporation when using these data.&rft_subject=environment&rft_subject=CONTAMINANT LEVELS/SPILLS&rft_subject=EARTH SCIENCE&rft_subject=HUMAN DIMENSIONS&rft_subject=ENVIRONMENTAL IMPACTS&rft_subject=OIL SPILLS&rft_subject=RECLAMATION/REVEGETATION/RESTORATION&rft_subject=LAND SURFACE&rft_subject=SOILS&rft_subject=PETROLEUM&rft_subject=SOLID EARTH&rft_subject=EARTH GASES/LIQUIDS&rft_subject=ENVIRONMENTAL ASSESSMENTS&rft_subject=ENVIRONMENTAL GOVERNANCE/MANAGEMENT&rft_subject=PERMAFROST&rft_subject=FROZEN GROUND&rft_subject=SOIL TEMPERATURE&rft_subject=SAB&rft_subject=Special Antarctic Blend&rft_subject=Fuel&rft_subject=Fuel evaporation&rft_subject=Fuel Spill&rft_subject=Contamination&rft_subject=Hydrocarbons&rft_subject=Biodegradation&rft_subject=Soil&rft_subject=Bioremediation&rft_subject=GC-FID&rft_subject=Model&rft_subject=GC-FID > Gas Chromatography - Flame Ionization Detector&rft_subject=LABORATORY&rft_subject=FIELD INVESTIGATION&rft_subject=OCEAN > ARCTIC OCEAN&rft_subject=OCEAN > SOUTHERN OCEAN > MACQUARIE ISLAND&rft_subject=CONTINENT > ANTARCTICA&rft_subject=GEOGRAPHIC REGION > POLAR&rft_subject=OCEAN > SOUTHERN OCEAN > HEARD AND MCDONALD ISLANDS&rft_place=Hobart&rft.type=dataset&rft.language=English Access the data

Licence & Rights:

view details

This data set conforms to the CCBY Attribution License (http://creativecommons.org/licenses/by/4.0/). Please follow instructions listed in the citation reference provided at http://data.aad.gov.au/aadc/metadata/citation.cfm?entry_id=fuel_evaporation when using these data.

Access:

Open view details

These data are publicly available for download from the provided URLs.

Brief description

Evaporation Model for hydrocarbon spills.

Developed by the Australian Antarctic Division to simulate fractionation of Special Antarctic Blend (SAB) and other diesel range fuels during evaporation.

Text version of notes for excel model, please read me.

This Package of files includes a pdf of the scientific paper, a readme word document and 3 excel files. The purpose of the 3 excel files are as follows:

Excel File 1. Evap Model_V1_single temperature.xls. Evaporation predictions with a single sample at one temperature. This file requires input of initial composition, composition after weathering and a single temperature.

Excel File 2. Evap Model_V1_five temperatures.xls. As above but with the input of 5 different temperatures.

Excel File 3. Evap Model_development version, derivation of data and AAD examples.xls This is an earlier development version of the numerical model. This file is intended to allow others to see how available Antoines equation parameters were originally fitted for use with other R+UCM regions. This file includes a worksheet where the fitting parameters for Antoines eqn were calculated, a necessary task to be able to extrapolate to compounds for which explicit vapour pressure data are not available. The data input sheet has some differences to the singe temperature and five temperature model as the R+UCM regions are split up differently. Also the assumptions about how the different classes of compound (i.e. the Aliphatic and aromatic classes) behave can be altered. Some raw data from the AAD evaporation experiment is included and plotted. The Temperature comparisons worksheet summarises how ratios of interest change and allows plotting of one ratio at different temperatures (see worksheet Figs 2 and 3).

Read me about background to this excel file
Background Notes

This excel file estimates the relative evaporation rates of different hydrocarbons from a hydrocarbon mixture (i.e. a fuel).

In this excel file, a single temperature is considered for an evaporating hydrocarbon mixture. The desired temperature and details about the hydrocarbon mixture are entered in the Main Input page.

The model was developed for use with the bulk fuels used by the Australian Antarctic Division at Casey, Davis, Mawson and Macquarie Island. This fuel, Special Antarctic Blend (SAB) starts at C9. As some of the spill sites undergoing remediation are mixed with heavier lube range hydrocarbons the options for inputs go to C36 - a typical maximum when Total Petroleum Hydrocarbon analysis is undertaken.

The model uses thermodynamic data where such data are available. Estimated thermodynamic data are used for components when specific data are unavailable.

Other important information can be found in the sections listed below.
References
GC-FID data used in the model
Temperature and Vapour pressure estimation
Other corrections
How the calculations are done
Experimental data (that supports the approach taken in this excel file)

Read me about References
Notes about the scientific publication that this model is reported in.

Paper can be obtained as a pdf file from the Australian Antarctic Data Centre. The paper contains details on the evaporation experiments at +20 degrees C and -20 degrees C. The results from the experiments agree well with the calculated fractionation rates that this Excel file produces.

Read me about GC-FID data for this model
Notes about GC-FID data used in this excel model.

This model was set up to directly use GC-FID data outputs for each of the compounds of GC regions listed. Consistent units are required for the areas of each measured region (and GC bias needs to be low across the range of fuel components). Note that the summation of all the regions and compounds = total area identified in the chromatogram from C9 up to C36. This range covers the observed range of components in SAB, Arctic Blend Diesel and lubes that may have been spilled in the same area.

Read me about temperature / vapour pressure estimation
Notes about Temperature correction method

Raoults law and the Antoine equation are used to calculate the composition of the evaporating portion of the fuel mixture.

Vapour pressure data were obtained for a range of available hydrocarbons C9 and larger in the temperature ranges that covered the site temperature ranges, approx -20 degrees C to +20 degrees C. The available data were mostly limited to n-alkanes. This situation was exacerbated because other compounds of interest are solids at these temperatures when pure. Consequently vapour pressure data are not available for these components in the liquid form at these temperatures.

The available n-alkane vapour pressure data were combined and a best fit of these data were determined as a function of effective Carbon Number (ECN). This allowed the estimation of the vapour pressure of fuel components with ECN's between the n-alkane ECN's.

As a refinement each region of the chromatogram was split into 5 classes - 3 x aliphatic fractions and 2 x aromatic fractions. The ECN used to estimate Vapour pressure of each class was slightly modified from the average retention time relative to an n-alkane. This modification was carried out to correct for specific GC-column / compound interactions - interactions that increase from acyclic to cyclic to polycyclic to hindered aromatics [i.e. multiply alkylated] to unhindered aromatics. Examination of the evaporation behaviour of n-alkanes vs. the nearest ECN regions confirms the need for the correction with the calculated evaporation profiles better matching the observed profiles.

When an n-alkane's evaporation rate is compared to a GC region with distinctly different vapour pressures these corrections to ECN make little or no difference to the predicted selectivity during evaporation.

Read me about other corrections
Notes about other corrections used in this model.

This model does not include systems that are under diffusional control (i.e. limited by diffusion of hydrocarbons within the soil / fuel mixture). The assumptions in the model are for an evenly mixed liquid that is evaporating. The experiment was set up to avoid this problem (by rotating the flasks).

Soil with evaporating fuels may well be affected by this and other problems.

Read me about how calculations are done
Notes about how the calculation is done

The fuel is divided into a number of fuel classes and specific compounds are identified. These classes are listed in the Main Input Page ready for use with the example data or other fuel data in the appropriate fuel range.

R+UCM stands for Resolved + Unresolved Complex Mix but needs to be calculated excluding the specifically identified compounds like n-alkanes. Specifically identified compounds are excluded from the R+UCM so they are not double counted.

Mol fraction of each component is estimated. For specific compounds an exact molecular mass is known. For other R+UCM classes molecular mass is estimated from known compounds in that region.

Vapour pressure of the pure component or region is estimated with Antoines equation for the temperature required. This pure vapour pressure estimate is combined with liquid phase mol fraction to calculate the gas phase composition at each evaporation step.

The evaporating portion (i.e. the gas phase portion) is removed from the liquid phase portion. The Excel sheet is setup such that this subtraction accounts for approximately 1% of the initial fuel amount. The need to calculated mol fractions then back-calculated mass remaining is the reason it is not exactly 1% by mass at each step.

To avoid numerical errors, division by zero errors and rounding errors many of the calculations contain and the IF formula. When a fuel component is greater than 0 the mol fractions are calculated, otherwise a value of 0 is returned.

Read me about AAD experimental data
Notes about AAD experiments that back up the model (a full description is in the scientific publication)

Small portions of Special Antarctic Blend (SAB) fuel were placed into vials. Each vial was placed into a +20 degrees C or -20 degrees C chamber. A slow stream of nitrogen passed into the top of each vial to remove the evaporating portion of the fuel. The vials were slowly rotated to ensure even mixing of the residual fuel during rotation.

Periodically a vial was removed, weighed to calculate mass fuel evaporated, and analysed with GC-FID apparatus. After a range of vials were analysed at different levels of evaporation at the 2 temperatures a data set was obtained to validate the numerical model.

Issued: 2006-08-01

Data time period: 2003-11-01 to 2003-11-30

This dataset is part of a larger collection

Click to explore relationships graph

160,-53 160,-70 60,-70 60,-53 160,-53

110,-61.5

text: northlimit=-53.0; southlimit=-70.0; westlimit=60.0; eastLimit=160.0; projection=WGS84

Other Information
Identifiers