United States Environmental Protection Agency
                           CBP/TRS 5/87

                            August 1987
A Comparison of Preservation
Techniques for Estuarine Water
       Samples for Analysis of
      Organic Carbon Fractions
                     Chesapeake
                       ^   Ba^
                       Program

-------
A COMPARASION OF PRESERVATION TECHNIQUES FOR ORGANIC
      CARBON ANAYSIS IN ESTUARINE WATER SAMPLES
                 ROBERT C. SIEGFRIED
               CHESAPEAKE BAY PROGRAM
            VIRGINIA WATER CONTROL BOARD
               2111 NORTH HAMILTON ST.
                RICHMOND, VA.  23230
                    JANUARY, 1987

-------
INTRODUCTION:

      In order  to determine  whether freezing  of  water  quality
samples was a acceptable method of preservation,  the Chesapeake
Bay  Program  Monitoring Subcommittee funded  a  comprehensive
comparison  study.   Virginia  Institute  of Marine Science
preformed a comparison of five preservation treatments for eight
water  quality parameters  (TP, TDP,  OP,  N02, N023,  TKN, Si,
TSS) .   These results were reported in Salley et al  1984.  The
State  of  Maryland conducted a comparison  of three preservation
treatments  for particulate  nitrogen, phosphorus  and  carbon
samples.   These  results were  reported  in Vaas  (1986) .  Old
Dominion  University provided the laboratory analyses  for total
organic carbon  (TOC)  and  dissolved  organic carbon  (DOC).  The
statistical analysis and graphical presentation of these  data  is
presented  in this  report prepared  by Chesapeake  Bay  Program
staff of the Virginia Water Control Board.

METHODS:

      Four  stations,  two  in the James  and two  in the York
rivers, were  sampled to give a range of salinities and nutrient
levels.  Five liter carboys  of water were filled at each  station
and  returned to  VIMS.   The  samples  for  carbon analysis were
transported to ODU.  Five different combinations of preservation
and  holding times  where employed.   Analyses for  the  'Day  0'
samples (DO)  were conducted 24 hours after collection.   Samples
where  frozen for 7 days  (D7F)  and 28 days  (D28F)  to  examine
freezing as a preservation  method.  Samples where also preserved
with  acid and  analyzed  within 48  hours  of collection  (Dl) and
after  a  28 day  (HT) holding time.   See Salley  et al (1986) for
further discussion of sampling.

      The  laboratory analyses  were  conducted  following EPA
Method Reference Number 415.1.   Inorganic carbon was purged from
the sample, which was then  digested to  convert organic carbon  to
C02 ,  which was  measured  with  an  infrared  detector.  Five
replicates   where  prepared  for  analysis   for  each
treatment-station  combination,  but  due  to the  loss of some
samples during  the  purging process,   the number of replicates
ranged  from three  to  five.   The quality assurance  and  quality
control measurements included  running standards, duplicates, and
spikes.   In  general,  for  each  treatment three  standards were
analyzed  during  the  run,  four samples where  spiked with known
concentrations,  and  30-40  duplicates  were  run.   The following
are the QA\QC results.

-------
Standard= 1.4  mg/1    Duplicates           Spiked Samples
Average % Diff.      Average % Diff.      Average % Recovery
                      from mean

HT    112              9.8                    129.9
Dl    122              6.8                    110.2
DO    102              7.0                    81.9
D7F   114              5.6                    89.0
D28F  108             14.4                    126.7


      Statistical methods similar to those used  by Salley et al
(1986)  were employed to test whether mean carbon concentrations
were different between the five treatments for each of the four
stations.   As with the other preservation comparison reports,
comparisons between the same treatments from different stations
were deemed not appropriate.  Comparisons between treatments from
the same  station where the major  interest  in  this  study.   \the
use of  four different stations was  to  obtain  a  range of sample
concentrations and matrixes.

      Due to  the limited  number of replicates  (3-5)  in each
treatment,   rigorous tests of the  data for  normality  and
homogenity  of variance  could not  be  performed.   Since these
tests  of the  assumptions of  parametric  statistics were  not
performed,  both parametric  and nonparametric  statistics were
applied to the data.

      For  each  of  the  four  stations,  a one  way ANOVA  was
performed to  test  the null hypothesis  that all  treatment means
were equal.   A parametric Tukey's multiple comparison test was
preformed  to identify which treatment means,  if any,  were
significantly  different.   Tukey's comparisons assume  equal
sample  size per treatment which was not always true within this
experiment.   A  Scheffe's comparison was conducted at both alpha
levels,  since  this  test  does not assume equal sample size,  but
it  is  not as  a powerful  test as  Tukey's if sample sizes are
equal.
      A nonparametric ANOVA similar to the Kruskal-Wallis Test
was performed by  assigning  ranks to  the  data for each station
then  performing a one-way ANOVA  on the  ranked data  (Vaas,
1986).   Tukey's  multiple comparisons were also conducted on the
ranked data to identify significantly different treatments.

-------
RESULTS:

      A  listing  of the  raw  data  with  treatment means  and
standard deviations for TOC  and  DOC are presented  in Tables  I
and II,  respectively.  Table III and IV contains the results of
the statistical analysis of the TOC  and DOC data, respectively.
Plots  of the  treatment means  for  each  station  for  TOC and DOC
are depicted in Figures 1-4  and Figures 5-8, respectively.
In  the discussion  of  the  results an alpha  level of 0.05 is
considered significant.  Results for  both  an alpha level of 0.05
and 0.01 are presented in the .tables.

TOTAL ORGANIC CARBON;

      For  all  stations  the  one-way  parametric ANOVA  was
significant, indicating  that not all  treatments were equal.  The
parametric  Tukey's multiple comparisons  identified the HT
treatment to be significant different from the DO treatment for
all stations.  For the  James River  samples, the HT treatment was
significantly different  from all other treatments, with no other
significant differences  observed in  the other treatment
comparisons.  The York 1 comparisons  identified the HT treatment
as  significantly different from the DO and D7F treatments.  The
D7F treatment was  also significantly  different from the D28F
treatment.   The York 2 comparisons identified the HT treatment
as significantly different  from the DO and D28F treatments.  The
D7F treatment was  also significantly  different from the D28F
treatment.

      The Scheffe's comparisons  identified  eleven of the 14
significantly  different  treatment  comparisons identified with
the Tukey' s comparisons.   The James  1  and York 1 comparisons
results  were the  same  as  with  Tukey's  test.   In  the James  2
comparisons, the HT treatment was significantly different from
the other treatments except the D28F treatment.   In the York  2
comparisons  only the  HT  treatment was significantly  different
from the DO treatment,  in contrast to  the Tukey's test  which
identified  three significant comparisons.

      The nonparametric ANOVA  (Kruskal-Wallis)  procedures also
indicated that, at  each station, not all  treatments were equal.
The nonparametric  Tukey's multiple  comparisons of the ranked
data identified a  few more  significantly  different comparisons
than the parametric  Tukey  comparisons.   The James 1 comparisons
were the same as  in the  parametric  Tukey's,  i.e.  HT  treatment
was significantly  different from all  other treatments.   The

-------
James  2  comparisons showed  all treatments were significantly
different from the  DO treatment except  the Dl treatment.   The  Dl
treatment  was significanlty  different from the  HT and D28F
treatments.   The  HT,  D7F  and  D28F treatments were not
significantly  different  from  each  other.    For  the York  1
comparisons  both the HT  and  D28F treatments were significantly
different from the  DO  treatment.   The D7F treatment was
significantly different from all  treatments  but  the  DO
treatment.   The York 2  comparisons indicated that  both  the  HT
and D28F treatments where significantly different  from  the  DO
treatment and the  HT treatment was  also  significantly different
from the  D7F  treatment.

DISSOLVED ORGANIC CARBON;

      The one-way parametric ANOVA was significant  for  only the
York 1 and York 2 stations.  The  treatment means for the James 1
and James  2  stations were  not  significantly  different.  The
Tukey's  multiple comparisons  for the York 1 station identified
the D28F treatment as significantly different  from  the  DO,  HT
and D7F  treatments.  The  York 2 comparisons identified  the D28F
treatment   as  significant  different from  the DO  and  HT
treatments.   The  Scheffe's   comparisons identified the same
significant differences as the  Tukey's comparisons.

      The nonparametric ANOVA (Kruskal-Wallis)  procedure
indicated the  same  results as the parametric  ANOVA.  The
nonparametric  Tukey's multiple  comparisons  of  the ranked data
identified a few more significantly different comparisons than
the parametric comparisons.   For the York 1 comparisons  the  Dl
treatment was  significantly  different from the DO and  HT
treatments.    As  in  the parametric comparisons,  the D28F
treatment was significantly  different  from all treatments but
the Dl treatment.   York 2 comparisons  indicated that  the D28F
treatment was  significantly  different from the DO and  HT
treatments, which is the same as  in the parametric comparisons.

DISCUSSION:

      The basic  question  of this research is if the preservation
methods employed preserve the original  organic carbon content  of
the samples.   The DO treatment was the  control against which the
other  treatments  where  compared.   The HT  treatment is the
currently EPA  approved and  accepted  method.   The freezing
treatments were compared  to both  the control treatment  (DO) and
the EPA  approved treatment (HT)  to determine  if freezing  could
be adopted as comparable to the EPA approved method.

-------
      Examination of the plots of the TOC results,  (Figures  1  -
4) ,  shows a general  trend  for  increasing  TOC  with increased
holding time, whether frozen or acidified.   While a decrease in
organic  carbon could be explained by poorly preserved samples,
an  increase in  not  easily  explained.   Day  to  day variation in
the  accuracy of  the  analytical  method may represent  as  much
variation  as that introduced  with  the  use  of  different
preservation methods.   Based on the Scheffe's  comparisons,  the
HT  treatment was  significantly different from the  control
treatment in three  out of  four  stations.  The D7F treatment
appears  to  not  be  significantly different from  the DO
treatment.   Examination  of the  plots  and  raw  data shows an
increase in variation and mean concentration when comparing  the
D28F  treatment to  the control.    The  D28F treatment  is  not
significanlly different  from the  DO  treatment in  3 out of  4
stations  based on  the Scheffe's  comparisons.The EPA approved
holding time of 28 days with acidification (HT)  is significantly
different from the  control  for three  out  of four stations,  and
significantly different from the other treatments in  six out of
the  other twelve  possible comparisons.   Of the  five treatments,
the  EPA approved holding time appears to be the least favorable
for the preservation of TOC  samples  based on comparison with  the
control treatment.  When freezing is compared to HT,  especially
the  short term  freezing  treatment,  this preservation method
appears  comparable or  better.    When compared  to  the  control
treatment,  freezing is  somewhat guestionable,  especially when
held for 28  days.

      The plots  of the  DOC results (Figure 5-8)  do  not depict
any  general pattern to the data.   The  statistical analyses
indicate  that there are no significant  differences  between  the
treatments  from  the James  1  and  2  stations.   The raw data
indicates a larger amount  of variance in the D28F treatment  than
the other treatments.   The HT treatment  for  the James 2  station
was  also  highly variable  (range 1.8 to  4.5).   The EPA approved
28 day with acidification  treatment  (D28F) does  not appear to be
significantly different  from the  control  treatment  (DO) ,  in
contrast with what was  found in  the TOC  analyses.   For the  York
1 and 2  stations,  the  D28F  treatments  were significantly
different from the  'control1 treatment  (DO), as  well  as  the HT
treatment,  yet  not  from  the Dl  treatment.    All  seven  day
freezing  treatments (D7F)  where not significanly different  from
the control.  Out of 8  possible comparisons between the freezing
treatments and the control  treatment, only 2 where significantly
different.   When  compared  to the control  or the HT treatment,
freezing  appears to be an acceptable method of preserving  DOC
samples  despite  a  few significant differences.   Since these
diffenences  where detected  in the long term freezing treatment,
the length of storage time  for frozen samples should be as short
as possible.

-------
                              REFERENCES
Salley,  B.  A., J.  G. Bradshaw  and B. J.  Neilson. 1986.  Results of
Comparative Studies of Preservation Techniques for Nutrient Analysis on
Water Samples. Gloucester Piont, Virginia Institute of Marine Science.

Vaas,  P.  A.  1986.   Freezing  of  Estuarine Nutrient Samples  as a
Preservation  Technique:  The Analysis of Particulate Nitrogen, Carbon,
and  Phosphorus Fractions.   Maryland Office of Environmental  Programs,
Ecological Modeling and Analysis Division.   Technical Report No. xxx.

-------
 TABLE I
  TOTAL

JAMES 1
ORGANIC  CARBON  ANALYSES

     HT       Dl       DO
  MEAN
STD.DEV.
               5.8
               4.9
               5.6
               3.7
    5.0
    0.8
             2.7
             2.5
             2.6
             3.0
2.7
0.2
2.2
2.4
2.4
2.9
3.3

2.6
0.4
D7F

2.9
2.7
2.5
2.6
2.5

2.6
0.1
D28F

 2.8
 3.1
 2.8
 3.3
 3.0

 3.0
 0.2
JAMES 2
YORK 1
     HT
     HT
 Dl
 Dl
 DO
 DO
D7F
D7F
D28F





MEAN
STD.DEV.
4.9
4.6
5.1
4.8
5.0
4.9
0.2
3.2
2.9
2.9
3.1
3.3
3.1
0.2
3.1
3.0
2.8
2.8
3.2
3.0
0.2
3.3
3.3
3.2
3.4

3.3
0.1
3.1
3.3
5.8
3.5
3.6
3.9
1.0
D28F





MEAN
STD.DEV.
5.4
3.9
3.8
4.3

4.4
0.6
3.4
3.6
3.6
3.1
4.2
3.6
0.4
3.4
3.5
3.2
3.1
3.1
3.3
0.2
2.8
2.8
2.9
3.2
3.3
3.0
0.2
3.6
3.3
3.9
4.5
4.7
4.0
0.5
YORK 2
  MEAN
STD.DEV.
     HT

    3.6
    3.6
    3.7
    3.6
    0.0
 Dl

3.6
3.0
3.3
2.8
3.2
0.3
 DO

2.5
2.9
2.5
3.0
2.5
2.7
0.2
D7F

2.8
2.6
3.1
2.7
2.9

2.8
0.2
D28F

 3.6
 3.3
 2.7
 4.3
 4.2

 3.6
 0.6

-------
       TABLE  II
DISSOLVED ORGANIC  CARBON  ANALYSES
JAMES  1
 HT
 Dl
 DO
D7F
D28F





MEAN
STD.DEV.
2.6
2.0
2.5
2.6

2.4
0.2
2.3
2.2
2.2
2.5

2.3
0.1
2.1
2.2
2.2
2.5
2.6
2.3
0.2
2.8
2.5
2.3
2.3
. 2.1
2.4
0.2
2.4
2.6
2.2
3.9
2.5
2.7
0.6
JAMES 2
  MEAN
STD.DEV.
 HT

2.5
1.8
4.3
2.4
2.7
0.8
 Dl

3.5
2.9
2.8
3.0
3.6

3.2
0.3
 DO

2.6
2.4
2.5
2.5
2.6

2.5
0.1
D7F

2.7
2.6
2.5
2.8
2.7
0.1
D28F

 2.2
 2.5
 3.4
 2.4
 3.2

 2.7
 0.5
YORK 1
 HT
 Dl
 DO
D7F
D28F





MEAN
STD.DEV.
2
2
2
1

2
0
.1
.3
.2
.7

.1
.2
2
2
2
2
3
2
0
.5
.5
.7
.6
.2
.7
.3
2.
2.
2.
2.
2.
2.
0.
3
4
2
2
2
3
1
2.
2.
2.
2.
2.
2.
0.
2
1
3
4
6
3
2
2
2
3
3
4
3
0
.8
.3
.3
.8
.6
.4
.8
YORK 2
  MEAN
STD.DEV.
 HT

2.9
3.0
3.1
3.0
0.1
 Dl

3. 1
2.7
2.9
2.5
2.8
0.2
 DO

3.2
2.6
3.1
3.0
2.5
2.9
0.3
D7F

2.5
2.2
3.0
2.3
2.5

2.5
0.3
D28F

 2.6
 2.4
 2.0
 2.4
 1.9

 2.3
 0.3

-------
TABLE III.


                             TOTAL ORGANIC CARBON


                  James 1        James 2         York 1           York 2

One-way ANOVA      0.0001        0.0001          0.0013           0.0038

PARAMETRIC
TUKEY'S Mult.  DO Dl HT D7F     DO Dl HT D7F     DO Dl HT D7F     DO Dl HT D7F
Comparisons  Dl  .             Dl .             Dl .             Dl .
             HT  # #           HT # #           HT*.          HT*.
            D7F  . .  #       D7F^, .  #       D7F'^  .  #      D7F .  .  .
           D28F  . .  #  .   D28F ',.' .  *  .   D28F ...  *  D28F $  .  .  *


Scheffe's     DO Dl HT D7F     DO Dl HT D7F     DO Dl HT D7F     DO Dl HT D7F
Comparisons Dl .             Dl .             Dl .             Dl .
            HT # #           HT # #           HT*.          HT..
           D7F . .  #       D7F . .  #       D7F .  .  #      D7F .  .  .
          D28F . .  #  .   D28F ....   D28F ...  *  D28F *  ...
                  James 1      James 2       York 1         York 2

KRUSKAL-WALLIS     0.0018      0.0001        0.0001         0.0041
NONPARAMETRIC
ANOVA

NONPARAMETRIC   ''
TUKEY'S Mult.  DO Dl HT D7F    DO Dl HT D7F     DO Dl HT D7F    DO Dl HT D7F
Comparisons  Dl  .            Dl .             Dl .             Dl .
             HT  # *          HT # #           HT #  .          HT*.
            D7F  . .  #      D7F * .  .       D7F .  *  #      D7F .  .  *
           D28F  . .  *  .  D28F # *  .  .   D28F *  .  .  #  D28F *  ...
# = Significant difference between means at alpha=0.01.
* = Significant difference between means at alpha=0.05.
. = No significant difference between means at alpha=0.05.

D0= DAY 0  (CONTROL)
Dl = DAY 1
HT = EPA APPROVED HOLDING TIME, 28 DAYS WITH ACIDIFTCATION
D7F = SAMPLE FROZEN FOR SEVEN DAYS
D28F = SAMPLE FROZEN FOR 28 DAYS

-------
TABLE IV;

                           DISSOLVED ORGANIC CARBON
James 1
0.4568
James 2
0.4136
York 1
0.0026
York 2
0.0072
One-way ANOVA

PARAMETRIC
TUKEY'S Mult.  DO Dl HT D7F    DO Dl HT D7F    DO Dl HT D7F     DO Dl HT D7F
Comparisons  Dl  .            Dl .            Dl .             Dl .
             riJL  • •          riP • •          rxP •  *          Hi •  •
            D7F  ...      D7F . .  .      D7F ...      D7F . ..  .
           D28F  ....  D28F ....  D28F #  .  #  *  D28F *  .  *  .
Scheffe's      DO Dl HT D7F    DO Dl HT D7F    DO Dl HT D7F     DO Dl HT D7F
Comparisons  Dl .            Dl .            Dl .             Dl .
             HT . .          HT . .          HT .  .          HT .  .
            D7F ...      D7F . .  .      D7F ...      D7F .  .  .
           D28F ....  D28F ....  D28F *  .  *  *  D28F *  .  *  .
James 1
0.6715
James 2
0.1034
York 1
0.0002
York 2
0.0075
KRUSKAL-WALLIS
NONPARAMETRIC
ANOVA

NONPARAMETRIC
TUKEY'S Mult.  DO Dl HT D7F    DO Dl HT D7F    DO Dl HT D7F     DO Dl HT D7F
Comparisons  Dl  .            Dl  .            Dl *'   ^         Dl .
             HT  . .          HT  . .          HT .  | '         HT .  .
            D7F  . .  .      D7F  . .  .      D7F ...      D7F .  .   .
           D28F  ....  D28F  ....  D28F #  .  #  *  D28F *  .   *   .


* = Significant difference between means at alpha=0.05
. = No significant difference between means at alpha=0.05

D0= CAY 0  (CONTROL)
Dl = DAY 1
HT = EPA APPROVED HOLDING TIME, 28 DAYS WITH ACIDIFICATION
D7F = SAMPLE FROZEN FOR SEVEN DAYS
D28F = SAMPLE FROZEN FOR 28 DAYS

-------
 FIGURE 1
COMPARISON OF TOTAL ORGANIC CARBON
CONCENTRATION BY TREATMENT
e.o
3.3
a.o
2: 4.0
O
^ "— ^
o
0 3.3
3.O-
2.3-
2.O-
AT JAMES 1
\
\ '.
\\
\ \\
\\v
-tl^^^^
'--.. ..-•''
HT D1 DO D7F D28F
LEGEND
	 1- SD
	 - SD



TREATMENT
FIGURE 2
  COMPARISON  OF  TOTAL  ORGANIC  CARBON

         CONCENTRATION  BY TREATMENT
                         AT JAMES 2
 UJ
 o
                                                   LEGEND
                                                       mean


                                                       -e SD


                                                       - SD
             HT
                    D1
                           DO
                                  D7F     D28F

-------
FIGURE 3
cor
a.o-
o.o-

a.o-
CONCENTRATION
(MG/L)
V f .ť
a o a
3.O-
a.o-
a.o-
\/IPARISON OF TOTAL ORGANIC CARBON
CONCENTRATION BY TREATMENT
AT YORK 1


\ •
/
X \ 7
X Vx \ //
'"'"^^^-"r^'''

HT D1 DO D7F D28F
TREATMENT
LEGEND
	 mean
	 -t- SD
	 - SD



FIGURE 4
  COMPARISON  OF TOTAL  ORGANIC CARBON

         CONCENTRATION  BY  TREATMENT

                          AT YORK  2
 CJ
 z
 o
             HT
                    Dl      DO     D7F

                        TREATMENT
                                          D2SF
                                                     LEGEND
                                                         mean


                                                         -t- SD


                                                         - SD

-------
  FIGURE 5
cc
<.oo-
J.BŤ
J.38
CONCENTRATION
(MG/L)
-MM
ť H a
00-
1.13-
.00-
0-
DMPARISON OF DISSOLVED ORGANIC CARBON
CONCENTRATION BY TREATMENT
AT YORK 1
/
/ /
; /^^^^/>
.-'


HT D1 DO D7F D28F
TREATMENT
LEGEND
	 +• SO
. 	 — SD



FIGURE 6
      COMPARISON OF  DISSOLVED.  ORGANIC  CARBON
               CONCENTRATION  BY TREATMENT
                             AT YORK 2
     4.3O
 O
 (=
 UJ
 O
I o
                                                           LEGEND
                                                               moan


                                                               + SD

                                                               - SD
               HT
                       D1      DO      D7F

                           TREATMENT
D28F

-------
 FIGURE 7
COMPARISON OF DISSOLVED ORGANIC CARBON
CONCENTRATION BY TREATMENT
4.SO
3.84-
3.38-

CONCENTRATION
(MG/L)
— M M
'Ť u a
ID 0 -
1.13-
.sŤ-
o-
AT JAMES 1

x
X
r____^r^:^




HT D1 DO D7F D2BF
LEGEND
maon
	 -t- SD
	 - SD




TREATMENT
FIGURE  8
      COMPARISON OF  DISSOLVED  ORGANIC CARBON
               CONCENTRATION BY  TREATMENT
                             AT JAMES 2
     3.94-
     3.38- •
 O
 (=
   o
 O
            LEGEND

            	 moon

            	  -t- SD

            	  - SD
               HT      D1       DO      D7F
                           TREATMENT
D28F

-------