Skip to main content

Full text of "The use of computer intensive statistical modeling in estimating the variability of marine fouling communities."

See other formats


1111111 


mm 


^k 


HlB 


W0 


few 

9»lfil 


"■;  ■■:;■;' 


■BnHliil 

■HH1 

JRffiWw 

Hal 


awl 


HHBL. 


ffiHnlH 


Inn 


3Mft 


B— 

HP 


H 


Dudley  Knox  Library,  KPS 
Monterey,  CA  93943 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


THE  USE  OF  COMPUTER  INTENSIVE  STATISTICAL  MODELING  IN 
ESTIMATING  THE  VARIABILITY  OF  MARINE  FOULING 

COMMUNITIES 

by 

David  L.  Martin 

June  1983 


Thesis  Advisor: 


E.  C,  Haderlie 


Approved  for  public  release;  distribution  unlimited 

T209041 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


t.    REPORT  NUMBER 


2.  GOVT  ACCESSION  NO 


3.     RECIPIENT'S  CATALOG  NUMBER 


4.     TITLE  (and  Submit) 

The  Use  of  Computer  Intensive  Statistical 
Modeling  in  Estimating  the  Variability  of  Marine 
Fouling  Communities 


5.     TYPE  OF   REPORT  &   PERIOD  COVERED 

Master's  Thesis 
June  198.3 


S.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS 


8.  CONTRACT  OR  GRANT  NUMBERS 


David  L,  Martin 


i      PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


10.     PROGRAM  ELEMENT,  PROJECT,   TASK 
AREA  &   WORK  UNIT  NUMBERS 


II.     CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


12.  REPORT  DATE 

June  1983 


13.  NUMBER  OF  PAGES 


115 


14.     MONITORING  AGENCY  NAME  *   AOORESSf//  different  from  Controlling  Office) 


IS.     SECURITY  CLASS,  (ol  thla  report) 


15«.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


l«.     DISTRIBUTION  STATEMENT  (ol  thla  Report) 


Approved  for  public  release;  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (of  the  mbatract  entered  In  Block  20.  It  different  from  Report) 


It.     SUPPLEMENTARY  NOTES 


1».    KEY  WOROS  (Contlnua  on  reveree  aid*  II  nacaaamry  and  Idantlty  by  block  number) 

Biofouling,  Statistics,  Variability,  Bootstrap  Computer  Simulations, 
Maximum  Likelihood 


20.     ABSTRACT  'Continue  on  revetee  tide  II  naeaaamrr  and  Idantlty  by  block  number) 

The  variability  of  the  fouling  community  in  Monterey  Bay  was  investigated  by 
suspending  100  mild  steel  plates  in  Monterey  Harbor.   The  plates  were 
painted  with  either  a  non-toxic  control  paint  or  one  of  three  antifouling 
paints.   Following  the  monthly  retrieval  of  a  group  of  these  plates,  a 
census  of  the  fouling  organisms  was  conducted  and  initial  variability 
estimates  determined.   These  estimates  were  used  as  inputs  for  bootstrap 
simulations  of  theexperiment .   The  results  of  the  bootstrap  simulations 


do  ,; 


FORM 
AN  73 


1473  EDITION  OF   1  NOV  65  IS  OBSOLETE 

S/N  0102-  LF-  014-  6601 


1       SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Bnterec 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh«n  Dmtm  Bnffd) 


BLOCK  20:   ABSTRACT  (Continued) 


were  then  used  to  determine  an  appropriate  strategy  for  sampling  the 
fouling  community  in  Monterey  Bay.   The  results  indicate  that  twenty  to 
thirty  plates  are  required  to  resolve  ambiguities  concerning  the  mean 
percent  cover  of  a  group  of  plates  while  many  more  are  required  to  quantify 
the  variability  of  the  fouling  population. 


S    N  0102-  LF-014-  6601 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWTiwi  Datm  Enfrud) 


Approved  for  public  release;  distribution  unlimited 


The  Use  of  Computer  Intensive  Statistical  Modeling 
in  Estimating  the  Variability  of  Marine  Fouling  Communities 


by 


David  L.  Martin 
Lieutenant,  United  States  Navy 
B.S.,  University  of  Washington,  1976 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  OCEANOGRAPHY 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  1983 


/M3573 


ABSTRACT 

The  variability  of  the  fouling  community  in  Monterey 
Bay  was  investigated  by  suspending  100  mild  steel  plates  in 
Monterey  Harbor.  The  plates  were  painted  with  either  a  non- 
toxic control  paint  or  one  of  three  antifouling  paints. 
Following  the  monthly  retrieval  of  a  group  of  these  plates, 
a  census  of  the  fouling  organisms  was  conducted  and  initial 
variability  estimates  determined.  These  estimates  were  used 
as  inputs  for  bootstrap  computer  simulations  of  the 
experiment.  The  results  of  the  bootstrap  simulations  were 
then  used  to  determine  an  appropriate  strategy  for  sampling 
the  fouling  community  in  Monterey  Bay.  The  results  indicate 
that  twenty  to  thirty  plates  are  required  to  resolve 
ambiguities  concerning  the  mean  percent  cover  of  a  group  of 
plates  while  many  more  are  required  to  quantify  the 
variability  of  the  fouling  population. 


TABLE  OF  CONTENTS 

I .  INTRODUCTION 15 

A.  GENERAL  15 

1.  Sampling  Design  16 

2.  Previous  Research  on  Fouling  Community 
Variability  16 

B .  OBJECTIVE  18 

II.  METHODS  AND  MATERIALS  19 

A.  GENERAL , 19 

B.  PLATES  AND  PAINTS  19 

1.  Priming  Procedure  21 

2.  Painting  Procedure  23 

C.  DEPLOYMENT  PROCEDURE  24 

D.  FOULING  COMMUNITY  CENSUS  AND  IDENTIFICATION  . .  27 

1.  Sampling  Procedure  27 

2.  Identification  28 

III.  STATISTICS  29 

A.  EXPERIMENTAL  DATA 29 

1.  Percent  Cover  29 

2.  Similarity  29 

B.  STATISTICAL  MODELLING  31 

1.  Model  Alternatives  32 

2.  Procedure  34 

IV.  RESULTS  37 

A.   GENERAL  37 

5 


1.  Method  Verification  37 

2.  Explanation  of  Figures  43 

B.  RESULTS  FROM  MONTH   2  44 

1.  Experimental  Data  44 

2.  Computer  Simulations  45 

C.  RESULTS  FROM  MONTH   3  48 

Experimental  Data  48 

2 .   Computer  Simulations  50 

D.  RESULTS  FROM  MONTH   4  54 

1 .  Exper imental  Data  54 

2 .  Computer  Simulations  56 

E.  RESULTS  FROM  MONTH   5  56 

1 .  Exper  imental  Data  56 

2 .  Computer  Simulations  59 

F.  RESULTS  FROM  MONTH   6  63 

1.  Experimental  Data  63 

2.  Computer  Simulations  65 

G.  RESULTS  FROM  MONTH   7  68 

1 .  Exper  imental  Data  68 

2 .  Computer  Simulations  68 

H.   RESULTS  FROM  MONTH   8  70 

1.  Experimental  Data  70 

2 .  Computer  Simulations  73 

I .   RESULTS  FROM  MONTH   9  77 


1 .  Exper  imental  Data  77 

2.  Computer  Simulations  77 

J.   RESULTS  FROM  MONTH  10  79 

1.  Experimental  Data  79 

2 .  Computer  Simulations  79 

K.   RESULTS  FROM  MONTH  11  84 

1.  Experimental  Data  84 

2 .  Computer  Simulations  86 

V.   CONCLUSIONS  AND  RECOMMENDATIONS  90 

A.  DISCUSSION 90 

B.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH  92 

APPENDIX  A:  MICRON  22  ORGANO  METALLIC  POLYMER 

ANTIFOULING  PAINT  94 

APPENDIX  B:  NAVY  STANDARD  FORMULA  121  RED 

VINYL  ANTIFOULING  PAINT  95 

APPENDIX  C:  NAVY  STANDARD  FORMULA  170  BLACK 

CAMOFLAGE  ANTIFOULING  PAINT  96 

APPENDIX  D:  ZYNOLYTE  EPOXY  RUST  MATE  PAINT  97 

APPENDIX  E:  A  DISCUSSION  OF  THE  METHOD  OF  MAXIMUM 
LIKELIHOOD  ON  ITS  INCORPORATION  INTO  A 

MODEL  FOR  FOULING  COVER 98 

APPENDIX  F:  BOOTSTRAP  COMPUTER  SIMULATIONS  103 

APPENDIX  G:  TABULATED  MONTHLY  PERCENT  COVERAGE  VALUES 

FOR  NON-TOXIC  CONTROL  SURFACES  107 

APPENDIX  H:  LIST  OF  THE  SESSILE  SPECIES  IDENTIFIED 

BY  THE  RANDOM  POINT  CENSUS  AND  THE  MONTHS 
THEY  WERE  PRESENT  ON  THE  NON-TOXIC  CONTROL 

SURFACES  109 

LIST  OF  REFERENCES  Ill 

INITIAL  DISTRIBUTION  LIST  113 


LIST  OF  FIGURES 


1.  A  Diagram  of  Monterey  Bay  Showing  the  Deployment 
Site  at  the  Coast  Guard  Floating  Dock  20 

2.  The  Front  Side  of  One  of  the  Experimental  Plates  .   22 

3.  A  Drawing  Showing  the  Method  of  Attachment  of 

the  Suspending  Cable  and  Identification  Tag  25 

4.  A  Perspective  View  in  Cross  Section  Showing  the 
Deployment  of  the  Plates  at  the  Coast  Guard  Dock  .   26 

5.  Chart  Showing  the  Monthly  Mean  Similarity  Values 
for  the  Control  Surfaces  (Solid)  and  the  Anti- 
fouling  Coated  Surfaces  (Dashed)  39 

6.  A  Diagnostic  Plot  of  the  Model  With  the  Ranked, 
Normalized  Values  for  the  Individual  Epsilon 
Values  Plotted  on  the  Ordinate  (Abbreviated  as 
Z)  and  the  Theoretical  Order  Statistic  Plotted 
of  the  Abscissa.  The  Dashed  Line  Indicates 
Perfect  Correspondence  and  the  Dotted  Line  is 

the  Least  Squares  Best  Fit  for  the  Data  41 

7.  Chart  Showing  the  Arithmetic  Mean  of  the  Percent 
Fouling  Cover  From  the  Experimental  Data  (Solid) 
and  the  Bootstrap  Simulated  Mean  Percent  Cover 
(Dashed)  42 

8.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  2.   Dashed  Lines  Indicate  Mean 
Similarity  Values  46 

9.  Computer  Simulations  Using  Data  From  Month  2 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  47 


10.  Computer  Simulations  Using  Data  From  Month  2 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  49 

11.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  3.   Dashed  Lines  Indicate  Mean 
Similarity  Values  51 

12.  Computer  Simulations  Using  Data  From  Month  3 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  52 

13.  Computer  Simulations  Using  Data  From  Month  3 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  53 

14.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  4.   Dashed  Lines  Indicate  Mean 
Similarity  Values  55 

15.  Computer  Simulations  Using  Data  From  Month  4 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  57 


16.  Computer  Simulations  Using  Data  From  Month  4 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  58 

17.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  5.   Dashed  Lines  Indicate  Mean 
Similarity  Values  60 

18.  Computer  Simulations  Using  Data  From  Month  5 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  61 

19.  Computer  Simulations  Using  Data  From  Month  5 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  62 

20.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  6.   Dashed  Lines  Indicate  Mean 
Similarity  Values  64 

21.  Computer  Simulations  Using  Data  From  Month  6 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  66 


10 


22.  Computer  Simulations  Using  Data  From  Month  6 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  67 

23.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  7.   Dashed  Lines  Indicate  Mean 
Similarity  Values  69 

24.  Computer  Simulations  Using  Data  From  Month  7 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  71 

25.  Computer  Simulations  Using  Data  From  Month  7 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  72 

26.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  8.   Dashed  Lines  Indicate  Mean 
Similarity  Values  74 

27.  Computer  Simulations  Using  Data  From  Month  8 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  75 


11 


28.  Computer  Simulations  Using  Data  From  Month  8 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  76 

29.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  9.   Dashed  Lines  Indicate  Mean 
Similarity  Values  78 

30.  Computer  Simulations  Using  Data  From  Month  9 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  80 

31.  Computer  Simulations  Using  Data  From  Month  9 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  81 

32.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  10.   Dashed  Lines  Indicate  Mean 
Similarity  Values  82 

33.  Computer  Simulations  Using  Data  From  Month  10 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  83 


12 


34.  Computer  Simulations  Using  Data  From  Month  10 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  85 

35.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces 
for  Month  11.   Dashed  Lines  Indicate  Mean 
Similarity  Values  87 

36.  Computer  Simulations  Using  Data  From  Month  11 
Showing  the  Expected  Value  of  the  Mean  Percent 
Fouling  Cover  (the  Mean  of  the  200  Individual 
Group  Simulation  Percent  Fouling  Covers)  as  a 
Solid  Line  and  the  95%  Quantile  (Dashed)  of  the 
Expected  Mean  Percent  Fouling  Cover  as  a 

Function  of  Increasing  Sample  Size  88 

37.  Computer  Simulations  Using  Data  From  Month  11 
Showing:  (A)  the  Percent  Fouling  Coverage  That 
One  Can  Say  90%  of  the  Samples  Will  Have  Fouling 
Coverage  Less  Than  or  Equal  to  (With  95%  Confi- 
dence) ;  and  (B)  the  Standard  Deviation  of  Per- 
cent Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing 

Sample  Size  89 


13 


ACKNOWLEDGMENT 

I  express  my  sincere  gratitude  and  appreciation  to 
Distinguished  Professor  Eugene  C.  Haderlie  for  his  guidance 
and  many  excellent  suggestions  during  the  course  of  this 
study,  to  Professor  Donald  P.  Gaver  who  designed  the 
probabilistic  model  used  to  quantify  the  fouling  community 
variability,  to  Professor  Patricia  Jacobs  for  both  her 
critical  review  of  the  thesis  and  her  patience  in 
instructing  a  biologist  in  the  many  complexities  of  modern 
statistical  analysis,  and  finally  to  my  wife  whose  support 
and  understanding  made  the  successful  completion  of  this 
project  possible. 


14 


I.       INTRODUCTION 

A.       GENERAL 

The  term  biofouling  refers  to  the  settlement, 
attachment,  and  growth  of  marine  organisms  on  surfaces 
that  man  puts  into  the  ocean.  This  process  has  a  profound 
impact  on  naval  operations  due  to  the  fouling  of  ships 
hulls,  rudders,  salt  water  piping  systems,  sonar  domes, 
and  the  fouling  and  biodeter ioration  of  harbor  or  pier 
structures.  This  problem  has  been  estimated  to  cost  the 
Navy  several  hundred  million  dollars  a  year  (Woods 
Hole, 1952;  Fisher  et  al,1975)  due  in  part  to  increased  fuel 
requirements  caused  by  the  greater  frictional  drag  of 
fouled  ships  hulls,  the  increased  repair  or  replacement 
cost  of  piping  and  machinery  damaged  by  fouling  organisms, 
and  the  need  to  expend  funds  to  continually  remove  the 
fouling   organisms    from   vessels. 

The  principal  method  used  to  combat  the  problem  of 
biofouling  on  the  hulls  of  ships  has  been  the  application 
of  antifouling  paints.  In  general,  such  paints  contain  any 
of  several  metallic  compounds  which,  as  they  leach  out  of 
the  paint  matrix,    are   toxic   to   fouling   organisms. 

The  testing  of  antifouling  paints  prior  to  their 
release  for  general  use  involves  the  use  of  sampling  assay 
techniques.    In   simple    terms,    a  number   of   substrates      are 

15 


painted  with  the  paint  to  be  tested  and  are  deployed  in  the 
sea  for  some  arbitrary  length  of  time.  Following  the 
retrieval  of  these  substrates,  they  are  compared  with  non- 
toxic control  surfaces  that  have  been  exposed  for  the  same 
length  of  time  and  the  efficacy  of  the  antifouling  paint  in 
preventing  the  settlement  of  fouling  organisms  is 
ascertained. 

1.  Sampling  Design 

The  proper  number  of  plates  to  deploy  for  the 
testing  of  antifouling  paints,  or  for  that  matter,  the 
study  of  fouling  organisms  in  general,  has  always  been 
somewhat  arbitrary.  This  is  because  the  proper  number  of 
plates  to  deploy  to  sample  the  fouling  population  is  a 
direct  function  of  the  degree  of  variability  within  the 
population.  Despite  the  great  volume  of  fouling  research 
that  has  been  conducted,  the  study  and  quantification  of 
this  variability  has  only  very  recently  been  attempted. 

2.  Previous  Research  on  Fouling  Community  Variability 
Most  of  the  information  dealing  with  the 

variability  of  biofouling  communities  has  been  collected 
within  the  past  decade  and  is  often  somewhat  contradictory. 
This  contradiction  is  sometimes  caused  by  the  particular 
descriptors  of  the  fouling  community  (percent  cover, 
species  counts,  etc.)  the  researchers  used  in  the  analysis 
of  its  variability. 


16 


Research  concerning  the  variability  of  small 
numbers  of  panels  suspended  for  only  one  month  off  the 
Florida  coast  (Mook,1976)  showed  very  little  variability  in 
species  count.  These  results  should  be  interpreted 
cautiously  however   due    to   the   short   time  of   immersion. 

Similar  research  conducted  in  North  Carolina 
(Sutherland, 1974)  using  a  more  extensive  series  of  panels 
suggested  that  the  development  of  the  fouling  community  was 
extremely  variable.  This  research  supported  the  conclusion 
of  earlier  studies  in  California  (Boyd, 1972)  which  also 
found   significant   fouling   community  variability. 

Studies  conducted  in  Hawaii  (Schoener  et  al,1978) 
and  in  Massachusetts  (Osman,1977)  found  that  fouling 
community  variability  based  on  species  counts  was 
relatively  low.  This  conclusion  was  echoed  by  a  study  which 
analyzed  the  variability  of  identical  panels  in  terms  of 
total  percent  cover,  species  count,  and  inter-panel 
similarity  indices  (Schoener  and  Greene, 1980).  The  results 
of  the  study  indicated  that  approximately  ten  replicate 
panels  were  sufficient  to  resolve  to  a  high  degree  of 
confidence,     the    mean   value   of    these   descriptors. 

As  is  evident,  there  is  wide  disparity  between  the 
various  studies  conducted  to  date  on  what  the  variability 
of  the  biofouling  community  at  various  locations  truly  is. 
Before    meaningful    antifouling    paint    test   procedures    can    be 


17 


developed,  proper  sampling  techniques  based  on  the 
quantification  of  fouling  community  variability  must  be 
devised. 

B.      OBJECTIVE 

The  primary  objective  of  this  thesis  was  to  determine 
the  variability  of  the  biofouling  community  in  Monterey 
Bay.  Once  this  had  been  completed,  the  development  of  a 
appropriate  sampling  strategy  based  on  this  variability 
could  be  accomplished.  This  information  was  to  be  provided 
to  the  David  Taylor  Naval  Ship  Research  and  Development 
Center  so  that  modifications  to  present  antifouling  test 
procedures  could  be  undertaken  as  required. 


18 


II.    METHODS    AND    MATERIALS 

A.  GENERAL 

The  experiment  consisted  of  deploying  one-hundred  (100) 
mild  steel  plates,  each  painted  with  one  of  four  vinyl  or 
epoxy  based  paints,  in  Monterey  Harbor  (Figure  1)  for 
periods  of  up  to  11  months  and  determining  the  variability 
of  the  resulting  fouling  communities  that  settled.  Three  of 
the  paints  used  in  the  study  contained  antifouling 
compounds  while  the  fourth  (  the  control  surface  )  did  not. 
The  fouling  community  structure  on  each  plate  was  deter- 
mined destructively  (the  plates  were  not  redeployed  after 
study)  by  microscopic  analysis  and  the  fouling  population 
and  makeup  were  determined. 

B.  PLATES    AND    PAINTS 

The  one-hundred  plates  used  in  this  study  were 
fabricated  from  low  carbon,  mild  (structural)  steel  with 
dimensions  25.4cm  x  30.5cm  x  .16cm  .  A  small  (.64cm)  hole 
was  drilled  approximately  1.3  cm  down  from  the  midpoint  of 
the  top  edge  (one  of  the  edges  with  the  lesser  dimension) 
of  the  plates  to  allow  for  the  attachment  of  the  suspending 
line  for  deployment.  Additionally,  a  small  (.32cm  x  2cm) 
groove  was  milled  through  the  plate  approximately  6.5cm  to 
one    side    of    the    drilled    hole    for    the    attachment    of    the 

19 


Figure   1.  A  Diagram  of  Monterey  Bay  Showing  the  Deployment 
Site  at  the  Coast  Guard  Floating  Dock. 


20 


identification  tag.  For  consistent  reference,  the  front  of 
the  plate  was  chosen  as  that  side  which  would  face  the 
observer  when  the  plate  was  held  vertically  with  the  dril- 
led hole  at  the  top  and  the  identification  groove  to  its1 
left  (Figure  2). 

The  four  paints  used  in  the  experiment  were: 

1.  MICRON  22;  a  commercially  available  antifouling 
paint  containing  bis ( tr ibutyltin)  oxide  and  cuprous 
thiocyanate  as  the   antifouling   agents  (Appendix  A). 

2.  Navy  Standard  Formula  121  Red  Vinyl  Antifouling 
Paint;  the  discontinued  U.S.  Navy  antifouling  paint 
containing  cuprous  oxide  as  the  antifouling  agent 
(Appendix  B) . 

3.  Navy  Standard  Formula  170  Black  Camoflage  Vinyl 
Antifouling  Paint;  the  currently  used  standard 
antifouling  paint  of  the  U.S.  Navy  containing 
bis ( tr ibutytin)  oxide  and  tributyltin  fluoride  as  the 
antifouling  agents  (Appendix  C) . 

4.  Zynolyte  Epoxy  Rust  Mate;  a  commercially 
available  non-toxic  corrosion  resistant  epoxy  based  paint 
used  as  the  control  surface  (Appendix   D) . 

Allplates  were  sandblasted  then  primed  and  painted  in 
accordance  with  label  directions. 

1.   Priming  Procedure 

The  priming  procedure  consisted  of  first  applying 
one  coat  of  Navy  Standard  Formula  117  'Green  Wash  Primer' 
to  both  sides  of  all  plates  with  a  5cm  latex  rubber  brush 
and  allowing  this  to  dry  for  twenty-four  hours.  This  was 
followed  by  the  application  of  two  coats  of  Navy  Standard 
Formula  119  'Red  Lead'  to  both  sides  of  all  plates  using  a 
7.6cm  nylon  paint  roller.  The  first  coat  of  Formula  119  was 


21 


E 
u 

o 

6 
ro 


Figure 
Plates. 


2.   The  Front  Side  of  One  of  the  Experimental 


22 


allowed  to  dry  for  forty-eight  hours  before  the  second  coat 
was  applied,  and  an  additional  forty-eight  hours  drying 
time  was  allowed  prior  to  painting  with  the  vinyl  or  epoxy 
paints. 

2.   Painting  Procedure 

All  painting  was  accomplished  using  7.6cm  nylon 
paint  rollers.  The  plates  were  divided  into  four  groups 
(A,B,C,  and  D)  of  twenty-five  plates  each  and  labelled  and 
painted  as  follows: 

1.  Plates  Al  through  A25  were  painted  on  both  sides 
with  one  coat  of  MICRON  22  antifouling  paint. 

2.  Plates  Bl  through  B25  were  painted  on  both  sides 
with  one  coat  of  Navy  Standard  Formula  121  Antifouling 
Paint. 

3.  Plates  CI  through  C25  were  painted  on  both  sides 
with  one  coat  of  Navy  Standard  Formula  170  Antifouling 

paint. 

4.  Plates  Dl  through  D25  were  painted  on  both  sides 
with   one   coat   of   Zynolyte   Epoxy   Rust   Mate. 

The   paint   on   all   plates    was    then   allowed   to   cure 

for    96    hours    prior    to   deployment.    The    plates    were    labelled 

sequentially    in    each   group    (Al,A2,etc)    by    affixing,     through 

the    identification    groove,     an    embossed    DYMO    tape    label    to 

each    plate    using    monofilament    nylon    fishing    line.    This 

method  of   labelling    was    chosen    to   prevent    the   catalytic 

corrosion   problems   attendant   with   standard  bronze  or   copper 

tags    in   contact   with   the    steel   plates. 


23 


C.   DEPLOYMENT  PROCEDURE 

The  plates  were  randomly  divided  into  twelve  groups  of 
eight  plates  consisting  of  two  plates  from  each  of  the  four 
paint  groups.  To  each  plate  a  one  meter  length  of  .32  cm 
diameter  stainless  steel  cable  was  then  affixed  by  passing 
the  cable  through  the  drilled  hole  and  forming  a  loop  which 
was  closed  using  Nico-Press  crimp  fittings  (Figure  3).  A 
similar  loop  was  formed  at  the  distal  end  of  the  cable  to 
facilitate  attachment  at  the  deployment  site. 

On  22  and  23  May  1982,  the  plates  were  suspended  (by 
groups)  beneath  the  service  access  covers  that  extend  the 
length  of  the  floating  dock  at  the  Coast  Guard  Station 
Monterey.  Each  plate  was  individually  deployed  by  attaching 
the  distal  loop  of  the  cable  to  lOd  nails  driven  into  the 
dock  and  allowing  the  plate  to  hang  vertically  in  the  water 
(Figure  4).  The  plates  were  separated  by  a  minimum 
horizontal  distance  of  40cm  with  the  tops  of  the  plates 
approximately  one-half  meter  beneath  the  surface  of  the 
water.  Since  the  dock  rose  and  fell  with  the  tide,  the 
depth  of  immersion  remained  constant.  Water  depths  below 
the  plates  (at  MLLW)  ranged  from  approximately  three  meters 
depth  at  the  shallower  end  of  the  dock  to  more  than  ten 
meters  depth  at  the  seaward  end.  After  being  submerged  for 
one  month,  inspection  of  the  plates  revealed  that 
significant  galvanic  corrosion  had  occured  at  the  junction 


24 


a 


Monofilament  Nylon  line 
with  Identification  Tog 


Stainless  steel  cable 
with  Nico-  Press 
Crimp  fitting 


Figure   3.  A  Drawing  Showing  the  Method  of  Attachment  of 
the  Suspending  Cable  and  Identification  Tag. 


25 


SERVICE  ACCESS  COVER 


Figure   4.  A  Perspective  View  in  Cross  Section  Showing  the 
Deployment  of  the  Plates  at  the  Coast  Guard  Dock. 


26 


between  the  mild  steel  plates  and  the  stainless  steel 
cable.  As  a  result,  on  23  June  1982,  all  stainless  steel 
cables  were  removed  and  replaced  with  .95  cm  diameter  nylon 
line.  This  necessitated  the  exposure  of  each  plate  to  the 
atmosphere  for  approximately  thirty  seconds  during  the 
replacement  operation  but  ,  since  the  plates  remained 
moist,  it  was  felt  no  harm  was  done  to  the  fouling 
organisms    that   had    settled. 

D.       FOULING   COMMUNITY   CENSUS    AND    IDENTIFICATION 

Each  month  following  the  initial  deployment,  one  of  the 
groups  of  eight  plates  was  randomly  selected  for  retrieval 
and    study. 

1.      Sampling   Procedure 

For  each  group  of  plates,  a  sampling  grid  of  one- 
hundred  (100)  uniformly  distributed  random  points  was 
generated  and  graphically  plotted  by  computer  on  a  25.4  cm 
x  30.5  cm  output  sheet.  These  points  were  then  transferred 
manually  to  a  clear  plexiglass  cover.  The  fouling 
communities  on  the  plates  were  then  systematically  analyzed 
by  setting  the  plates  horizontally  in  a  shallow  container 
filled  with  seawater  and  positioning  the  plexiglass  cover 
over  the  top  coincident  with  the  edges  of  the  plate. 
Animals  beneath  plotted  points  were  then  censused  and 
identified  through  the  use  of  a  stereo  microscope.  Since 
the   plates   had   been   suspended   vertically    in    the    water,    none 

27 


of  the  sedimentation  problems  associated  with  horizontal 
deployment  strategies  developed.  Therefore  both  sides  of 
the  plates  were  analyzed  and  counted  as  separate 
substrates. 

2.      Identification 

To  eliminate  the  necessity  for  compound  microscope 
identification  of  settled  organisms  and  to  negate  the  ef- 
fects of  neuston  contamination  of  census  results,  only 
sessile,  attached  organisms  greater  than  .5mm  in  size  were 
counted  and  identified.  Identification  was  accomplished 
through  the  use  of  available  keys  and  literature; 
(Osburn,1952;  Knight  Jones, 1979;  Hader lie, 1974  ;  Morris 
et  al,1980;  Frazier,1937;  Smith  and  Carlton, 1975) .  In  all, 
some  thirty-two  taxa  were  identified  as  major  space  occu- 
piers   during    the    study    (see   Appendix   H   for    species    list). 


28 


III.    STATISTICS 

A.       EXPERIMENTAL    DATA 

The  data  obtained  from  the  census  and  identification  of 
the  organisms  on  the  plates  retrieved  each  month  was  used 
to  generate  the  statistics  that  described  the  fouling 
populations.  The  two  main  descriptive  quantities  used  to 
assess  the  variability  of  the  fouling  community  were  the 
total  percent  fouling  cover  on  each  plate  and  the 
similarity  of   the   fouling   organizations   between     plates. 

1.  Percent   Cover 

The  initial  estimate  of  the  total  percent  fouling 
cover  on  each  plate  was  calculated  by  dividing  the  number 
of  points  from  the  census  that  had  organisms  beneath  them 
(this  value  was  termed  St)  by  the  total  number  of  points 
censused  per  plate  (N  =  100  for  all  plates).  Since  this 
method  was  obviously  subject  to  some  unknown  degree  of 
uncertainty,  revised  estimates  for  the  total  percent 
fouling  cover  were  made  using  statistical  techniques 
described  below. 

2.  Similarity 

The  similarity  between  the  fouling  communities  on 
the  plates  retreived  each  month  was  determined  by  the 
calculation  of  the  Bray-Curtis  similarity  index, Ia 
(Whitaker,1952) .    This    index    is   defined    as: 

29 


Ia  =  ]jPmin(a,b) 
where  a  and  b  are  the  fractional  species  proportions 
present  on  plates  A  and  B.  These  fractional  species 
proportions  were  determined  from  the  initial  census  data  by 
dividing  the  total  number  of  instances  a  particular  species 
was  counted  during  the  plate  census  by  100.  In  this  study, 
empty  points  not  occupied  by  any  organism  were  treated  as  a 
separate  species.  As  an  example,  suppose  the  following  data 
were  collected: 


Species 

Proportion 
Plate  A   (  = 

on 
=a) 

Proportion 
Plate  B  {■■ 

on 
=b) 

min (a,b) 

Empty 
Species  #1 
Species  #2 
Species  #3 

.30 

0 
.46 
.24 

.24 
.11 
.42 
.23 

.24 

0 
.42 
.23 

then  the  Similarity  Index  =  Ia  =   T"*min(a,b)      =  .89 
The  similarity  index  is  thus  a  measure  of  the  degree  of 
variability  between  the  fouling  communities  on  the  two 
plates  in  terms  of  both  the  total  percent  cover  and  the 
species  composition. 

The  main  purpose  of  calculating  the  similarity 
index  between  plates  was  to  examine  whether  any  of  the 
antifouling  coated  plates  displayed  more  variability  than 
did  the  non-toxic  control  surfaces.  This  was  of  particular 
interest  since  a  previous  thesis  (Kelley,1981)  which  dealt 
with  marine  microfoulers  in  Monterey  Bay  had  shown  that 
surfaces  coated  with  an  organo-metallic  antifouling  paint 


30 


served  as  attractants  to  the  organisms  and,  hence,  were 
more  heavily  fouled,  showed  greater  species  diversity  ,  and 
were  more  variable  in  terms  of  their  fouling  communities 
than  were  the  control  surfaces.  Provided  that  the 
macrofouling  community  investigated  in  this  thesis  did  not 
behave  similarly,  that  is  if  the  fouling  communities  on  the 
non-toxic  control  surfaces  were  consistently  more  variable 
than  were  those  on  the  antifouling  coated  surfaces,  then 
any  sampling  strategy  which  could  discern  to  an  acceptable 
degree  of  error  the  amount  of  variability  of  the 
communities  on  the  control  surfaces  would  be  able  to  do  at 
least  as  well  concerning  the  antifouling  coated  surfaces. 
The  thrust  of  this  thesis  was  to  first  investigate  whether 
or  not  the  control  surfaces  for  each  month  exhibited  more 
variability  than  did  the  antifouling  coated  surfaces. 
Provided  this  criteria  was  met,  the  next  step  was  to  devise 
an  appropriate  sampling  strategy  for  the  control  surfaces 
using  sophisticated  statistical  modelling  to  extrapolate 
the  data  obtained  from  the  four  control  surfaces  retrieved 
each  month  to  any  desired  number  of  simulated  samples. 

B.   STATISTICAL  MODELLING 

Once  the  initial  percent  fouling  coverage  estimates  for 
the  four  non-toxic  control  surfaces  retreived  each  month 
were  determined,  the  data  had  to  be  manipulated  to  obtain 
better  estimates  of  the  variability  of  the  superpopulation 

31 


of  marine  fouling  organisms  the  plates  were  assumed  to 
be  sampling.  Since  there  were  only  four  substrates  examined 
per  month,  a  model  capable  of  extending  the  experimental 
data  to  simulated  samples  of  any  size  was  required.  The 
model  also  had  to  be  capable  of  allowing  and  quantifying 
inter-plate  variability  within  the  simulated  sample  set. 
The  techniques  used  to  successfully  meet  these  requirements 
have  only  recently  been  developed  and  this  thesis  is 
apparently  the  first  incorporation  of  these  techniques  into 
fouling  research. 

1.   Model  Alternatives 

One  of  the  methods  that  has  been  used  in  the  past 
(Schoener  and  Greene,  1980)  to  examine  the  degree  of 
variability  of  biofouling  coverage  values  has  been  to 
assume  that  the  fractional  coverage  estimates  have  a  normal 
distribution.  Using  this  simple  model,  the  upper  and  lower 
95%  confidence  limits  are  calculated  about  the  experimental 
mean  percent  fouling  cover  by  use  of  the  formula: 

the  95%  Confidence  Limits  =  x  [+/-]  1.96*  d  /  ^N~ 
where  x  is  the  mean  percent  fouling  cover,  d  is  the 
experimental  standard  deviation,  N  is  the  number  of  samples 
examined  ,  and  the  formula  is  derived  from  the  standard 
normal  distribution.  By  assuming  that  the  mean  (  x  )  and 
standard  deviation  (  a  )  do  not  vary  with  the  number  of 
plates  examined,  one  varies  N  in  the  formula  to  determine 


32 


the  effect  that  the  number  of  plates  has  on  the  approach  of 
the  upper  and  lower  confidence  limits  to  the  experimental 
mean.  It  must  be  pointed  out  however  that  this  model  has 
serious  drawbacks.  The  first  problem  with  this  approach  is 
the  assumption  that  fractional  fouling  coverage  values 
which  will  always  lie  between  zero  and  one  can  be  described 
by  a  normal  distribution  that  is  unbounded  in  range.  While 
this  assumption  permits  the  calculation  of  various  statis- 
tical parameters  using  well  known  analytical  formulas,  it 
is  obviously  weak  theoretically.  Secondly,  this  model 
requires  that  the  mean  and  standard  deviation  are  stable 
with  increasing  sample  size.  This  requirement  has  the 
effect  of  forbidding  inter-plate  variability.  Not  only  is 
such  a  restriction  biologically  untenable,  it  reduces  the 
formula  for  the  calculation  of  the  confidence  limits  to  : 

the  95%  Confidence  Limits  =[+/-]  constant/ -yN 
Note  that  the  above  equation  will  result  in  a  curve  similar 
to  [+/-]  1/JN  when  plotted.  Since  the  value  of  1/  J"n"~ 
decreases  by  nearly  70%  as  N  goes  from  one  to  ten,  use  of 
this  method  will  always  show  that  approximately  ten  plates 
are  sufficient  to  resolve  the  var  iability  of  the  mean 
percent  cover.  This  result  however  is  merely  an  artifact  of 
the  simplistic  model  used  in  its   calculation. 

Another  possible  method  that  could  be  used  to 
estimate  the  degree  of  variability  of  the  biofouling 


33 


community  would  be  to  invent  a  predictive  model  for  fouling 
populations.  Unfortunately,  such  a  model  is  considerably 
beyond  our  abilities  at  the  present  time.  Such  a  model 
would  require  advective/dif fusive  models  accurate  spatially 
to  microscale  resolution  .  It  would  also  require  the 
ability  to  parameterize  the  entire  range  of  biological 
forcing  functions  which  include  predation,  nutrient  supply, 
susceptibility  to  environmental  fluctuations,  behaviorism 
of  the  organisms  involved  and  planktonic  larva  survival 
rates  to  name  just  a  few.  Obviously,  such  sophistication  in 
a  model  is  not  likely  in  the  forseeable  future. 

Since  a  stochastic  model  which  assumed  the  fouling 
coverage  values  had  a  normal  distribution  was  insufficient 
due  to  its'  inherent  restrictions  which  forbid  interplate 
variability,  and  since  a  predictive  model  was  unattainable, 
this  study  used  a  probabilistic  model  to  explore  the 
variability  of  the  fouling  community.  In  such  a  model,  the 
variability  of  the  initial  experimental  data  is  determined 
and  the  statistical  descriptors  of  that  varibility  are  used 
as  inputs  for  computer  simulations  of  the  experiment. 
2.   Procedure 

Each   plate   (t)   was   assumed   to   have   a 

random, independent  proportion  of  fouling,  P  .  P.  was  assumed 

to  be  of  the  form:  e 

et 

P  =    e 


(1  +  e  r) 


34 


where  ^was  assumed  to  be  normally  distributed  with  unknown 
mean  {/j)  and  variance  (a  )  which  were  independent  of  the 
other  plates  censused  at  the  same  time.  Note  that  the  above 
equation  could  be  written  as: 
€t=  in(Pt  /(l  -  Pt  )) 
Each  month,  the  €  value  for  each  of  the  four  untreated 
plates  was  determined  as  were  the  mean  and  variance  of  the 
four  6t  values.  Revised  estimates  for  these  parameters  were 
then  determined  using  the  Method  of  Maximum  Likelihood  (see 
Appendix  E  for  mathematical  development) . 

Using  as  inputs  the  maximum  likelihood  values  for 
the  mean  (  /i  )  and  variance  (  a  )  of  the  monthly  epsilon 
values,  bootstrap  computer  simulations  (Efron,1979)  of  the 
experiment  were  conducted.  The  bootstrap  method  was  used  to 
simulate  200  groups  of  various  fixed  numbers  of  plates  (see 
Appendix  F  for  a  complete  discussion) . 

Using  the  results  from  the  bootstrap  computer 
simulations,  the  following  statistics  were  calculated: 


1.  The  expected  value  of  the  mean  percent  fouling 
coverage.  This  was  determined  by  calculating  the  mean  of 
the  200  simulated  group  means  for  each  of  the  various 
numbers  of  plates  per  group  simulated. 

2.  The  upper  95%  confidence  limit  about  the  expected 
value  of  the  mean  percent  fouling  cover  for  each  of  the 
various  numbers  of  plates  per  group  simulated.  This  was 
estimated  by  ordering  the  mean  percent  foulng  coverage 
values  for  each  of  the  200  simulated  groups  of  plates  in 
ascending  order  and  then  finding  the  95%  quantile  of  the 
expected  mean  percent  fouling  coverage  values. 


35 


3.  The  percent  fouling  cover  that  one  can  say  90%  of 
the  simulated  plates  will  have  fouling  coverage  less  than 
or  equal  to  (with  95%  confidence).  This  was  done  by 
finding  the  95%  quantile  (as  described  above)  of  the 
parameter   NEWCVR(I)    where: 

0  P 

NEWCVR(I)    =    e      /(l    +   e      ) 


A  A 


and     &  -     /!{!)+   1.285*     a  (I) 

Note   that       /i(I)    and      <j{I)      are    the    bootstrap   simulations 
of    the    mean   and    standard   deviation   of  for    each   of    the 

200    group    simulations.     The    value    1.285     is     the    90% 
quantile    for    the   standard   normal   distribution. 

The  purpose  of  calculating  the  parameter  NEWCVR  was 
to  estimate  how  many  plates  it  would  take  to  be  95% 
confident  of  capturing  90%  of  the  variance  of  the 
biofouling  population.  The  confidence  and  variance  values 
chosen  were  arbitrary,  and  the  model  can  be  readily 
modified  to  estimate  the  number  of  plates  required  to 
ascertain  any  degree  of  variance  to  any  desired 
confidence. 

4.  The  standard  deviation  about  the  mean  percent 
fouling  coverage  that  90%  of  the  simulated  plates  would 
have  was  also  estimated  .  This  was  done  by  finding  the 
mean  and  standard  deviation  of  the  200  NEWCVR(I)  values 
for  each  of  the  various  numbers  of  plates  per  group 
simulated   each   month. 


36 


IV.   RESULTS 

A.   GENERAL 

Experimental  data  were  collected  monthly  and  the 
initial  percent  cover  and  similarity  values  computed  (and 
the  coverage  estimates  tabulated)  for  the  second  through 
the  eleventh  months  of  the  experiment  (Appendix  G)  . 
Following  the  procedure  discussed  in  the  statistics 
chapter,  these  initial  estimates  were  used  as  starting 
values  for  Bootstrap/Maximum  Likelihood  computer 
simulations  of  the  experiment.  The  computer  simulation 
results  from  month  2  are  described  in  some  detail.  The 
simulation  using  the  data  from  months  3  through  11  followed 
the  same  procedure  and  are  each  presented  in  a  brief 
synopsis. 

1.   Method  Verification 

The  procedure  described  in  this  thesis  to  develop 
an  appropriate  sampling  strategy  for  antifouling  paint  test 
purposes  was  based  on  the  assumption  that  any  strategy 
which  could  ascertain  the  variance  of  the  biofouling 
community  on  non-toxic  control  surfaces  would  be  able  to  do 
at  least  as  well  regarding  antifouling  coated  surfaces. 
This  assumption  would  be  correct  provided  the  fouling  com- 
munities on  non-toxic  control  surfaces  were  more  variable 
than  were  those  on  the  antifouling  coated  surfaces.  By 

37 


analyzing  the  monthly  mean  similarity  values  for  the 
fouling  communities  on  these  two  types  of  surfaces  (Figure 
5),  it  is  clear  that  the  antifouling  coated  surfaces  dis- 
played consistently  less  variability  in  their  fouling 
structure  than  did  the  control  surfaces.  Therefore,  this 
assumption  is  considered  to  be  quite  strong. 

A  second  major  assumption  used  in  the  development 
of  this  procedure  was  that  the  monthly  epsilon  (6  )  values 
used  in  the  calculation  of  the  proportion  of  fouling  (P  ), 
were  independent  random  values  from  a  normal  distribution 
with  mean  (^  )    and  variance  (  a2  ) . 

If  the  epsilon  values  were  in  fact  distributed 
normally,  then  dividing  the  individual  epsilon  values  for 
each  month  by  the  monthly  standard  deviation  of  the  epsilon 
values  and  subtracting  the  mean  monthly  value  of  epsilon 
from  this  result,  would  transform  the  epsilon  distribution 
into  the  standard  normal  distribution.  By  plotting  the  (N) 
transformed  epsilon  values  ranked  in  ascending  order 
against  the  theoretical  order  statistic  obtained  by  finding 
the  inverse  function  of  the  standard  normal  distribution 
for  (j/N+1)  as  j  goes  from  1  to  N,  a  diagnostic  plot  of  the 
model  was  obtained.  If  the  model  was  perfect,  there  would 
be  a  one-to-one  correspondence  between  the  transformed 
epsilon  values  and  the  order  statistic.  Plotting  the  forty 
transformed  epsilon  values  obtained  in  this  study  against 


38 


o 
o 


•  — 

• •  -  - 

o 

\   < 

>- 

t— • 

SIMILRR 

0.50 
l 

LO 

•    — 

o 
o 

• 

o 

1               1 

1             II             1             11             I            ■ 

JUL  RUG  SEP  OCT  NOV  DEC  JflN  FEB  MflR  RPR 
1982  m   1983 

Figure   5.  Chart  Showing  the  Monthly  Mean  Similarity  Values 

for  the  Control  Surfaces  (Solid)  and  the  Antifouling  Coated 
Surfaces  (Dashed) . 


39 


the  order  statistic  (Figure  6)  shows  that  the  assumption  of 
a  normal  distribution  for  the  variable  epsilon  is  quite 
good.  The  least  squares  best  fit  line  for  the  data  had  a 
slope  of  1.07  (vice  1.00  for  the  theoretical  perfect 
correspondence)  and  the  correlation  coefficient  between  the 
transformed  epsilon  values  and  the  order  statistic  was 
0.98. 

Since  the  initial  experimental  data  for  percent 
fouling  cover  and  the  Bootstrap/Maximum  Likelihood  model 
were  not  independent  but  rather  were  coupled  by  the  esti- 
mates of  the  mean  and  variance  of  the  monthly  epsilon 
values,  one  would  reasonably  expect  the  final  mean  percent 
fouling  cover  predicted  by  the  Bootstrap  simulations  to 
approximate  the  actual  data  value  if  the  model  behaved 
reasonably.  A  chart  comparing  the  arithmetic  mean  of  the 
monthly  values  for  percent  cover  obtained  experimentally  to 
the  expected  value  of  the  mean  percent  cover  predicted  by 
the  bootstrap  simulations  (Figure  7) ,  shows  that  the  model 
agrees  quite  well  with  the  experimental  data.  In  those 
instances  where  there  was  a  significant  difference  between 
the  mean  percent  cover  obtained  from  the  data  and  that 
predicted  by  the  bootstrap  model,  analysis  of  the  actual 
data  suggested  that  simply  finding  the  arithmetic  mean  of 
the  four  monthly  data  percent  coverage  estimates  was  per- 
haps too  sensitive  to  unusually  high  or  low  values. 


40 


o 

CO 

« 

CM 

O 
OJ 

LO 


LO 

d 

o 

LO 

O 

I 


I 

in 


i 

LO 

CM 

I 


CO 

I 


.*  / 
f 


."  / 


■3.0-2.5-2.0-1.5-1.0-0.5     0.0     0.5     1.0     1.5     2.0     2.5     3.0 

PHI  INVERSE  OF  J  OVER  N+l 


Figure  6.  A  Diagnostic  Plot  of  the  Model  With  the  Ranked, 
Normalized  Values  for  the  Individual  Epsilon  Values  Plotted 
on  the  Ordinate  (Abbreviated  as  Z)  and  the  Theoretical 
Order  Statistic  Plotted  of  the  Abscissa.  The  Dashed  Line 
Indicates  Perfect  Correspondence  and  the  Dotted  Line  is  the 
Least  Squares  Best  Fit  for  the  Data. 


41 


o. 
o 


in 


0£ 

> 
O 

F-  o 

LJ 
CJ 
C£ 
LJ 
Q_ 


O 

■ 

o 


JUL  RUG  SEP  OCT  NOV  DEC  JRN  FEB  MflR  RPR 
1982  1983 

Figure  7.  Chart  Showing  the  Arithmetic  Mean  of  the  Percent 
Fouling  Cover  From  the  Experimental  Data  (Solid)  and  the 
Bootstrap  Simulated  Mean  Percent  Cover  (Dashed) . 


42 


In  summary,  it  is  felt  that  the  assumptions  made  in 
the  formulation  of  this  model  are  quite  good  and  that  the 
model  gives  very  reasonable  estimates  of  the  most  likely 
value  for  the  percent  fouling  cover. 
2.   Explanation  of  Figures 

The  similarity  indices  between  the  plates  retrieved 
each  month  were  plotted  on  two  separate  graphs  for  the 
control  surfaces  and  the  antifouling  coated  surfaces.  The 
range  of  the  similarity  values  from  zero  to  one  was  plot- 
ted on  the  ordinate  and  the  plate  designations  were  listed 
on  the  abscissa.  The  'F'  and  'B'  that  follow  the  plate 
group  number  refer  respectively  to  the  front  and  back  of 
the  plate.  To  obtain  the  similarity  value  between  any  two 
plates,  find  the  horizontal  line  above  the  abscissa  with 
arrowheads  which  terminate  at  points  above  the  desired 
plate  designations,  then  proceed  horizontally  to  the 
ordinate  to  find  the  similarity  index  value. 

The  expected  value  of  the  computer  simulated  mean 
percent  fouling  cover  and  the  upper  95%  confidence  limit 
about  the  mean  as  functions  of  increasing  sample  size  were 
plotted  for  each  month.  The  expected  value  of  the  mean 
percent  cover  was  dispayed  as  a  solid  line  and  the  upper 
95%  confidence  limit  was  dashed. 

Graphs  were  also  prepared  showing  the  results  of 
the  monthly  computer  simulations  used  to  determine  the  95% 


43 


confidence  limit  of  the  percent  fouling  coverage  that  90% 
of  the  simulated  plates  will  have  fouling  coverage  less 
than  or  equal  to.  The  ordinate  of  these  figures  was  gra- 
duated in  percent  cover  from  zero  to  one  hundred  percent. 
The  abscissa  was  labelled  with  the  number  of  plates  per 
group  used  in  the  simulations.  The  95%  confidence  limit  for 
the  percent  cover  of  90%  of  the  plates  was  then  plotted  as 
a  function  of  increasing  sample  size. 

The  same  units  for  the  ordinate  and  abscissa  were 
used  in  the  figures  showing  the  standard  deviation  about 
the  mean  percent  cover  that  90%  of  the  plates  would  have  as 
functions  of  increasing  sample  size.  For  these  figures, 
however,  the  percent  cover  label  on  the  ordinate  refers  to 
the  deviation  about  the  mean  percent  cover  and  not  the 
total  percent  fouling  cover  these  plates  will  have. 

B.   RESULTS  FROM  MONTH  2 
1.   Experimental  Data 

The  fouling  community  on  the  control  surfaces  after 
two  monhs  immersion  was  dominated  by  the  hydroid  Obelia 
spp.  Ectoproct  colonies  each  consisting  of  several  dozen 
zooids  had  also  settled  by  this  time.  These  initial 
bryozoans  were  primarily  Hippothoa  hyalina  and  Celloporaria 
b r_ u n n e a  although  the  ancestrula  stage  of  a  recently 
introduced  species,  ^a_ter_s__i£or_a  cucullata  was  also 
present.Watersipora    is    indigenous    to    the    Galapagos    Islands 

44 


and  has  not  been  reported  in  the  literature  farther  north 
than  the  Gulf  of  California  (Osburn,  1952) .  It  has  not  been 
noted  in  fouling  studies  conducted  in  Monterey  Bay  over  the 
last  twenty  years  and  its  appearance  on  the  plates  in  this 
study  is  probably  due  to  the  abnormally  warm  coastal  waters 
caused  by  the  recent  El  Nino  event. 

The  similarity  values  for  the  control  surfaces 
ranged  from  .74  to  .92  with  a  mean  of  .83  (Figure  8A).  The 
lack  of  any  settlement  on  the  antifouling  coated  surfaces 
resulted  in  similarity  values  of  1.0  for  all  of  those 
plates  (Figure  8B) . 

2.   Computer  Simulations 

The  bootstrap  computer  simulation  of  the  expected 
mean  percent  fouling  cover  (Figure  9)  resulted  in  a  final 
iterated  estimate  of  44.4%  for  the  simulated  mean  percent 
cover  (vice  44.5%  for  the  data  arithmetic  mean).  Note  that 
the  dashed  line  indicating  the  upper  95%  confidence  limit 
of  the  mean  percent  cover  does  show  some  inverse 
relationship  to  N  (the  number  of  samples)  but  does  not  show 
the  precipitous  approach  to  the  mean  as  N  goes  from  one  to 
ten  that  was  predicted  by  the  model  which  assumed  a  normal 
distribution  for  the  fractional  coverage  values. 

Perhaps  the  most  striking  result  of  the  computer 
simulations  was  that  of  the  effect  of  increasing  the 
number  of  simulated  plates  sampled  on  the  95%  confidence 


45 


o 

o 


[fl] 


>-• 

•—i 
Od  o 

-J  d 

•— • 

to 


t 


><- 


> 


B> 

> 


Dior  D10B 

PLATE  DESIGNATION 


D19F 


D19B 


o 

o 


(B) 


in 


>- 

0£  o 

az  w 


en 


L/l 


o 


£> 


fl21Ffl21BR25Ffl25BB13FB13BB15FB15BCHFC14BC24FC24B 
PLATE  DESIGNATION 


Figure  8.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  2. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


46 


or 

CjJ 

> 

o 
°c 

E-h  o- 

■z.  " 

CxJ 

o 

CiJ 
Q_ 


in- 
rg 


o 
o 


— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   Of  PLATES 


Figure  9.  Computer  Simulations  Using  Data  From  Month  2 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as    a   Function  of   Increasing    Sample    Size. 


47 


limit  of  the  percent  fouling  cover  that  90%  of  the  plates 
would  have  (Figure  10A).  Note  that  there  is  very  little 
dependence  on  N  on  the  percent  cover  that  one  can  say  90% 
of  the  plates  will  have  fouling  cover  less  than  or  equal 
to.  This  is  in  spite  of  the  fact  that  the  standard 
deviation  that  these  90%  of  the  plates  will  have  about  the 
mean  percent  cover  does  display  a  strong  inverse  dependence 
on  N  (Figure  10B).  While  these  two  results  might  seem 
dichotomous,  they  are  not.  The  standard  deviation  about  the 
mean  does  decrease  with  increasing  sample  size  just  as  one 
would  expect  from  the  Central  Limit  Theorem.  The  fact  that 
the  95%  confidence  limit  of  the  fouling  cover  of  90%  of  the 
plates  does  not  behave  similarly  is  simply  a  result  of  the 
fact  that  by  allowing  interplate  variability,  one  no  longer 
constrains  the  parameter  to  have  a  1/  N  dependence. 

C.   RESULTS  FROM  MONTH  3 
1.   Experimental  Data 

The  fouling  community  structure  on  the  control 
surfaces  after  three  months  immersion  remained  dominated  by 
the  hydroid  Obelia  spp.  The  four  species  of  spirorbid 
worms  that  live  in  Monterey  Bay  particularly  Circeis 
armor  icana,  were  also  present.  The  protozoan  Folliculina 
sp.  was  present  in  large  numbers  on  three  of  the  plates.  The 
bryozoans  were  represented  by  four  species  with 
Celloporaria  brunnea  dominating  this  group. 

48 


[fl] 


o 

8 


o 

on 
u 

> 
o 

Z  » 

LJ 
O 

Cl- 


in- 

CM 


■ 

o. 


0.0  5.0 


10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


o 


NO.   OP  PLATES 


> 
o 
o 

(B)| 

uj  a 
,     Q-  in 


1 1 1 1 1 1 \ 1 1 \ \ 1 1 1  I 

0.0  5.0    10.0   15.0  20.0  25.0  30.0.35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0 

NO.    OF  PLATES 


80.0 


Figure  10.  Computer  Simulations  Using  Data  From  Month  2 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have  As    a   Function   of    Increasing   Sample   Size. 


49 


The  similarity  values  ranged  from  .57  to  .86  with  a 
mean  of  .74  for  the  control  surfaces  (Figure  11A).  The 
mean  similarity  for  the  antifouling  coated  surfaces 
remained  at  1.0  (Figure  11B)  since  no  settlement  of 
organisms   had  occured. 

2.      Computer   Simulations 

The  computer  simulations  of  the  expected  value  of 
the  mean  percent  fouling  cover  (Figure  12)  resulted  in  a 
final  iterated  value  of  49.5%  vice  a  49.3%  arithmetic  mean 
of  the  initial  data.  The  upper  95%  confidence  limit  about 
the  mean  decreased  by  nearly  8%  as  the  number  of  simulated 
plates  per  group  increased  from  2  to  15  and  then  only 
decreased  an  additional  3%  as  the  number  of  plates  per 
group  was    increased    to   80. 

The  computer  simulations  of  the  95% confidence 
limit  of  the  fouling  cover  that  90%  of  the  plates  would  be 
fouled  to  an  equal  or  lesser  extent  showed  the  value 
remained  nearly  constant  at  approximately  70%  (Figure 
13A).The  standard  deviation  of  the  percent  cover  on  these 
plates  about  the  mean  percent  cover  showed  little  relative 
decrease  past  the  simulated  analysis  of  20  plates  per  group 
(Figure   13B) . 


50 


s 


~~ 

• 

u> 

.d 

< > 

f» 

& — 

o 

^1 

Jj^^l                      ^"S                      b»* 

>-> 

41 

Itfc 

ID 

SIMILAR] 

0.50 

^ 

w 

in 

CM 

o 

o 

-'■  -  ■  ■ 1-        ■■'■   ■ 

r                                    r                                    i 

D9P         D9B        D22F 
PLATE  DESIGNATION 


D22B 


o 
o 


(B) 


m 


C£  o 
CE  m 

-J  o 


zr 

CO 


a 
a 


<3" 


E> 


A6F    A6B  A12FA12B  B5F    B5B  B24F  B24B  C12F  C12B  C25F  C2SB 
PLATE  DESIGNATION 

Figure  11.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  3. 
Dashed  Lines   Indicate  Mean   Similarity  Values. 


51 


UJ 

Z> 

o 
°o 

-z.  » 

O 

ce 

CiJ 

a. 


CM 


o 

o . 


— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   Of  PLATES 


Figure  12.  Computer  Simulations  Using  Data  From  Month  3 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as   a  Function  of   Increasing   Sample   Size. 


52 


(A3 


o. 
o 


o 

an 
o 

°° 

£-<  d 

•^  in 

u 
o 


ir>- 


o 


1 1 1 1 1 1 1 1 1 1 1 1 ! 1 1 1 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.    OF  PLATES 


o 

QC  d- 

OJ  — 

> 

O 

o 

£B)  I 

a: 

Q_  in 


o 

d. 


— i 1 1 1 1 1 1 1 1 1 1 1 1         i         i         i 

0.0  5.0    10.0   15.0  20.0  25.0  30.0.  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 

NO.   OP  PLATES 


Figure  13.  Computer  Simulations  Using  Data  From  Month  3 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have  As    a   Function   of   Increasing   Sample   Size. 


53 


D.       RESULTS    FROM   MONTH    4 
1.      Experimental   Data 

The  fouling  community  on  the  control  surfaces  after 
four  months  immersion  was  dominated  by  the  bryozoan 
Watersipora  cuculatta.  Thehydroid  Obelia  spp.  had  been 
overgrown  to  a  large  extent  by  W^_  cuculatta  and  as  a  result 
its  relative  contribution  to  the  percent  cover  was 
diminished  considerably.  In  addition  to  W^_  cuculatta,  there 
were  six  other  species  of  either  upright  or  encrusting 
bryozoans  present.  Various  spirorbid  and  serpulid  worms 
were   present   in  limited  numbers. 

The  percent  coverage  estimates  varied  greatly  from 
plate  to  plate  with  a  maximum  of  84%  and  a  minimum  estimate 
on  one  plate  of  2%  cover.  The  wide  variability  was  perhaps 
partially  caused  by  the  ascendency  of  the  bryozoans  as  the 
dominant  organism  but  the  reason  for  the  nearly  total  lack 
of   fouling   on   two  of   the   surfaces    is   unknown. 

The  wide  range  in  the  percent  coverage  estimates 
was  mirrored  in  the  variability  of  the  similarity  indices 
between  the  control  surfaces.  The  values  ranged  from  .16  to 
.98  (Figure  14A).  The  antifouling  coated  plates  still  had 
no  fouling  so  again  the  mean  similarity  was  1.0  for  those 
plates. 


54 


8 


(R) 


in 


-J  o" 


en 


in 


o 
o. 


> 


> 


<T 


E> 


* 


> 


< 


D14F        D14B        D21P 
PLATE  DESIGNATION" 


D21B 


o 

a 


IB) 


in 


C£  o 
-J  o 

I— • 

to 


in 


o 
o 

o. 


B> 


A3F    R3B  fl22Ffl22BB14FB14BB2QFB20B  C4F    C4B    C6F    C6B 
PLATE  DESIGNATION 

Figure  14.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  4. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


55 


2.      Computer   Simulations 

The  final  iterated  value  for  the  simulated  mean 
percent  fouling  cover  was  31%  (Figure  15).  The  upper  95% 
confidence  limit  about  this  mean  was  quite  large  and 
decreased  by  nearly  40%  as  the  number  of  simulated  plates 
per   group  went   from   two   to   twenty    (Figure   15) . 

The  wide  variability  of  the  percent  coverage 
estimates  caused  the  95%  confidence  about  the  fouling 
coverage  estimate  of  90%  of  the  plates  to  be  quite  high 
(Figure  16A).  As  can  be  seen,  the  only  statement  that  can 
be  made  about  the  fouling  coverage  that  90%  of  the  plates 
willhave  is  that  the  coverage  will  be  something  less  than 
93%  (Figure  16A)  even  though  the  estimated  mean  percent 
coverage  was   31%. 

The  standard  deviation  that  these  90%  of  the  plates 
will  have  about  the  mean  was  also  quite  large  (Figure  16B) 
and  showed  significant  decline  out  to  approximately  the 
thirty  plates   per   group    simulation   point. 

E.       RESULTS    FROM   MONTH    5 
1.      Experimental   Data 

The  bryozoan  dominance  of  the  fouling  community  was 
firmly  established  by  the  fifth  month  of  immersion. 
Watersipora  cuculatta  was  the  primary  fouler  on  three  of 
the    four    control    surfaces    and    there    were    an    additional 


56 


> 
o 

°o 

t-  o 

2  m 

UJ 
O 

or 

CjlJ 


q 
d. 


0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 

NO.   OP  PLATES 


Figure  15.  Computer  Simulations  Using  Data  From  Month  4 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as   a  Function  of   Increasing   Sample   Size. 


57 


o 

8n 


in- 

U 

> 

o 
*-•  a- 

Z  in 

tflJa 

lj 


in- 


! 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0  5.0-10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   OF  PLATES 


o-i 


r> 

OH  o 


(B)< 


CJ   O 

OH  u^. 

o 


o. 


1 1 1 1 1 1 !  I  I  I  1  1  I  I  T  1 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 

NO.    OF  PLATES 


Figure  16.  Computer  Simulations  Using  Data  From  Month  4 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
oamples   Will  Have   As    a   Function   of    Increasing   Sample   Size. 


58 


eight    other    species    of    bryozoans    identified    during    the 
census. 

There  was  a  fairly  large  range  of  initial  percent 
cover  estimates  for  the  control  surfaces  with  a  minimum  of 
13%  and  a  maximum  of  96%.  The  mean  percent  coverage 
estimate  for  these  plates  was  63%  and  the  standard 
deviation   was    31%. 

The  similarity  indices  for  the  control  surface 
plates  ranged  from  .17  to  .67  with  a  mean  value  of  .45 
(Figure  17A).  The  s i m i 1  a r i t y  v a  1 u e s  for  the 
antifouling  coated  surfaces  ranged  from  .85  to  1.0  (Figure 
17B).  This  variability  was  caused  by  the  settlement  of 
Obelia  spp.  and  the  spirorbid  worm  Circeis  armor  icana  on 
several  of  the  antifouling  coated  plates. 
2.      Computer    Simulations 

The  final  iterated  estimate  for  the  mean  percent 
cover  was  69%  (Figure  18)  vice  the  arithemetic  mean  of  63% 
for  the  initial  data.  The  95%  coinfidence  limit  about  this 
mean  showed  an  appreciable  decrease  out  to  approximately 
twenty  plates   per   group   simulated. 

The  population  variability  again  resulted  in  the 
95%  confidence  limit  of  the  fouling  cover  of  90%  of  the 
plates  having  virtually  no  dependence  of  the  number  of 
samples  analyzed  (Figure  19A).  The  standard  deviation  about 
the    mean    that    these    90%    of    the    plates    would   have    showed 


59 


8 


V* 

Hi 

^ 

^ 

fe> 

az  o 

^ 

— -w 
— -t> 

2: 

^-     "- 

j^ 

^ 

*4 

^ 

in 

«1 

► 

^ 

0 
0 

P" " 

1 

1 

1 

D8F         D8B        Dllf 
PLRTE  DESIGNATION 


DUB 


o 
o 


(B) 


in 


CE  i" 
-3  o 


CD 


LT> 


Q 

O 


>"j i^i  1?^  ^  £j  ^j  i:^  fe: 


fllF    RIB  R16FR163  B8F    B3B  B10FB10B  C9F    C9B  C21FC21B 
PLRTE  DESIGNATION 

Figure  17.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  5. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


60 


o 

8-1 


a: 

LJ 
> 
O 

"o 

I—  O' 

2  « 

LiJ 

O 

tr 
a. 


04 


q 
o. 


0.0  5.0  10. 0  15. 0  20. 0  25. 0  30. 0  35. 0  40. 0  45. 0  50. 0  55. 0  60. 0  65. 0  70. 0  75. 0  80.0 

NO.   OP  PLATES 


Figure  18.  Computer  Simulations  Using  Data  From  Month  5 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as   a  Function  of   Increasing    Sample   Size. 


61 


(R) 


m 
K 

en 

> 
o 

■z  « 

LJ 

O 

CC 

LJ 

Cl- 


in- 

CNJ 


Q 
O. 


0.0  5.0 


- T 1 I 1 1 1 1 1 1 I I I 1 I 1 

10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   OF  PLATES 


O-i 


in- 


EB) 


o 
en  6- 

> 
O 

°° 
E— •  in 

"Z.  " 

LJ 

O 

S° 


o 
in' 


o 
d. 


! — 

0.0  5.0 


i        r 
10.0  15. 


- 1 1 r 

20.0  25.0  30 


1 1 \ 1 1 — 

.0  35.0  40.0  45.0  50.0  55.0 

NO.   OP  PLATES 


i         i         r 

60.0  65.0  70. 


1 1 

0  75.0  80.0 


Figure  19.  Computer  Simulations  Using  Data  From  Month  5 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As   a  Function   of   Increasing   Sample   Size. 


62 


little  decrease  past  approximately  20  plates  per  group 
simulated  (Figure  19B) . 

F.   RESULTS  FROM  MONTH  6 

1.   Experimental  Data 

The  fouling  community  structure  on  the  control 
surfaces  after  six  months  immersion  was  dominated  either  by 
bryozoans  or  the  hydroid  Obelia  spp.  depending  on  the  plate 
examined.  Two  species  of  solitary  tunicates,  Ascidia 
ceretodes  and  Styela  truncata  were  also  censused  for  the 
first  time  during  the  experiment. 

The  initial  percent  cover  estimates  for  the  control 
surfaces  ranged  from  16%  to  54%.  The  mean  percent  cover  for 
the  control  surfaces  was  37%  with  a  standard  deviation  of 
18%. 

The  antifouling  coated  surfaces  showed  fouling  on 
several  of  the  plates  with  the  barnacles  Megabalanus 
californicus  and  Balanus  crenatus  appearing  for  the  first 
time  on  these  surfaces.  Obelia  spp.  remained  the  dominant 
fouler  on  the  antifouling  coated  surfaces. 

The  similarity  values  for  the  control  surfaces 
ranged  from  .48  to  .88  with  a  mean  of  .70  (Figure  20A).  The 
mean  similarity  of  the  antifouling  coated  surfaces  was  .99 
(Figure  20B) . 


63 


[R] 


►—4 

-J  d 

I—* 

CO 


- 

<« — 

►<* 

> 

rfi 

fcs, 

^1. 

■%J- 

<~ 

► 

-4 

> 

^ 

► 

-■     ■-                 r,_ 

r  -       — 

T 

i 

D12F        D12B        D15F 
PLATE  DESIGNATION 


D15B 


o 
o 


(B) 


in 


>- 

E- 

CE  Lr>. 


CO 


o 

o 


fl8F 


A8B  R15FA15B 


B19FB19B 


B4F  B4B 
PLATE  DESIGNATION 


C3F 


C3B  C15FC15B 


Figure  20.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  6. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


64 


2.      Computer   Simulations 

The  computer  simulated  value  of  the  mean  percent 
cover  was  33%.  The  upper  95%  confidence  limit  about  the 
mean  showed  significant  decrease  out  to  approximately  20 
plates  per  group  simulated  (Figure  21).  Note  again  that 
there  was  a  significant  decrease  in  the  width  of  the  95% 
confidence  limit  about  the  mean  as  the  number  of  plates  per 
group   simulated  went   from   ten   to   twenty. 

The  rather  high  mean  similarity  value  for  the  sixth 
month  control  plates  and  by  inferrence  the  lessened 
interplate  variability  permitted  the  95%  confidence  limit 
of  the  fouling  cover  of  90%  of  the  plates  (Figure  22A)  to 
be  less  than  that  of  the  fourth  month  simulations  even 
though  the  mean  percent  coverage  value  for  the  sixth  month 
group  was  higher.  The  95%  confidence  limit  of  the  fouling 
cover  of  90%  of  the  plates  value  decreased  until 
approximately  the  twenty  plates  per  group  simulation  point 
and    then    remained   fairly  constant. 

The  standard  deviation  about  the  mean  percent  cover 
that  these  simulated  plates  displayed  decreased  only 
slightly  past  the  20  plates  per  group  simulated  point 
(Figure  22B).  Note  however  that  there  was  a  5%  decrease  in 
this  value  between  the  10  plates  per  group  simulated  point 
and   the   20   plates   per   group   simulated. 


65 


o 

o 

in- 
tx 

\ 

ERCENT  COVER 

50.0 

\ 
\ 

i  • 
\ 
\ 
\ 
\ 
\ 
\ 

Q_ 

V/"            " 

O 

LO  — 

CM 

O 

a 

l          ■!■-■,■     ■■    , 

1               1                1                1                1                1                1                1                !                1                1                1 

0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 

NO.   OF  PLATES 


Figure  21.  Computer  Simulations  Using  Data  From  Month  6 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as   a  Function  of   Increasing    Sample   Size. 


66 


IO- 
CS 

LJ 

> 
O 

°o 

t-H     O 

Z  in 

tfllB 

on 

LJ 


«- 


a 
d. 


— r^ — i 1 1 1 1 1 1 1 1 1 i 1 1 1 1 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   OP  PLATES 


to 


in- 

CM 


o 

CC  d- 

> 

O 


IB)  8 


o 
io" 


o 
d. 


I  I  I  1  I  I  I  I  1  I  I  I  I  I  l  i 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  SO.O  65.0  70.0  75.0  80.0 


NO.    OF  PLATES 


Figure  22.  Computer  Simulations  Using  Data  From  Month  6 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have  As    a   Function   of    Increasing   Sample   Size. 


67 


G.   RESULTS  FROM  MONTH  7 

1.  Experimental  Data 

The  fouling  community  structure  of  the  seventh 
month  conrol  plates  remained  dominated  by  bryozoans  with 
the  hydroid  Obelia  spp.  still  in  an  important  but  secondary 
role.  The  spirorbid  worm  Circeis  armoricana  was  dominant  on 
one  of  the  plates  and  was  common  on  all  the  control 
surfaces. 

The  primary  fouling  organisms  on  the  antifouling 
coated  surfaces  were  the  barnacle  Balanus  crenatus  ,  the 
bryozoan  Celloporaria  brunnea  ,and  the  hydroid  Obelia  spp. 
Several  of  these  plates  were  also  fouled  by  an  unknown  tube 
dwelling   amphipod,    possibly  of   the   genus   Ampithoe. 

The  percent  coverage  estimates  for  the  control 
surfaces  ranged  from  31%  to  62%  with  a  mean  of  49%.  The 
standard  deviation  of   these   values   was   11.5%. 

The  similarity  indices  between  the  control  surfaces 
ranged  from  .22  to  .88.  The  mean  similarity  value  for  the 
control    plates    was    .55     (Figure    23A). 

For  the  antifouling  coated  surfaces,  the  range  of 
similarity  values  was  from  .81  to  1.0.  The  mean  value  was 
.92    (Figure   23B). 

2.  Computer   Simulations 

The  final  iterated  value  for  the  simulated  mean 
percent    fouling    coverage    was    49%.    The    upper    95%    confidence 


68 


o 

o 


(R) 


^4 

M 

►        ' 

in 
d~ 

^ 

hi 

^^ 

t^ 

•— « 

^ 

N» 

-J  d" 
►—1 

^ 

^ 

^ 
"> 

in 
d~ 

^ 

► 

8 

o 

7 

1     '      -  '■ 

T 

i 

D23P        023B        D24F 
PLATE  DESIGNATION 


D24B 


o 
o 


(B) 


m 


f— 

or  o 
en  m 


CD 


in 
cm 


o 
o 


^#^^#i 


^     <^ =^4-^» — P> — |> — ► 


A7F    A7B  R13PR13B  B11FB11B  B16FB16BC10rC10BC13FC13B 
PLATE  DESIGNATION 

Figure  23.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  7. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


69 


limit    about    this    mean   remained    fairly   constant  past   the   20 
plates   per   group  simulated   point    (Figure   24) . 

The  simulated  percent  cover  that  one  could  say  90% 
of  the  plates  would  have  fouling  coverage  less  than  or 
equal  to  (Figure  25A),  decreased  slightly  until 
approximately  the  20  plates  per  group  simulation  point  and 
then  remained  fairly  constant  at  approximately  65% 
coverage.  The  standard  deviation  that  these  simulated 
plates  had  about  the  mean  percent  fouling  cover  showed 
little  decrease  past  the  30  plates  per  group  simulation 
point    (Figure   25B) . 

H.       RESULTS    FROM   MONTH    8 
1.      Experimental   Data 

The  fouling  community  structure  on  the  control 
surfaces  after  eight  months  immersion  was  dominated  nearly 
exclusively  by  bryozoans.  The  contribution  of  hydroids  to 
the  percentage  of  cover  had  been  reduced  on  three  of  the 
surfaces  to  less  than  5%.  Serpulid  worms,  solitary 
tunicates,  and  barnacles  were  beginning  to  emerge  as 
important  foulers  as  their  increasingly  large  size 
prohibited   their   overgrowth   by   encrusting   species. 

The  fouling  organisms  present  on  the  antifouling 
coated  surfaces  were  primarily  the  hydroid  Obelia  spp.  and 
barnacles.  The  spirorbid  worm  Circeis  armor icana  was  also 
present   on    several   of    the   plates. 

70 


a: 

CjJ 

> 
o 

E-«  o- 

QJ 
O 

a: 


CM 


O 
O. 


I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   OP  PLATES 


Figure  24.  Computer  Simulations  Using  Data  From  Month 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as   a   Function  of   Increasing    Sample   Size. 


71 


S-i 


(R) 


QJ 
> 
O 

"o 

t-  o- 
-^r  in 

LJ 
O 

LJ 
Q_ 


tn- 

CM 


O 
O. 


0.0  5.0 


10.0  15.0  20.0  25.0  30.0  35.0 


40.0  45 


1 1 l 1 1 1 1 

.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


o 
in-. 


NO.   OP  PLATES 


(B) 


o 
CC  a- 

> 
O 

o 


UJ 

o 
a: 
lj  a 

0-  in 


o 
d. 


1  I 1  l 1  l  l 1  l 1 1  l  i  i  i  i 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  50.0  65.0  70.0  75.0  30.0 

NO.   OF  PLATES 

Figure  25.  Computer  Simulations  Using  Data  From  Month  7 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have   As    a  Function   of   Increasing   Sample   Size. 


72 


The  estimates  for  the  percentage  of  fouling  cover 
ranged  from  48%  to  97%  for  the  control  surfaces.  The  mean 
of   these   values   was   76%    and   the   standard   deviation  was   19%. 

For  the  control  surfaces,  the  range  of  similarity 
values  was  from  .22  to  .67.  The  mean  value  was  .44  (Figure 
26A)  . 

The   similarity  values   for   the    antifouling  coated 

surfaces    ranged    from    .75    to    1.0.     The    mean    value    was    .94 
(Figure    263) . 

2.      Computer    Simulations 

The  computer  iterated  estimate  for  the  mean  percent 
fouling  coverage  for  the  control  surfaces  was  81%.  Note 
that  the  95%  confidence  limit  about  this  value  remained 
fairly  constant  past  approximately  the  25  plates  per  group 
simulation  point    (Figure   27) . 

The  high  estimate  for  the  simulated  mean  percent 
cover  and  the  variability  permitted  by  the  model  forced 
the  upper  95%  confidence  limit  about  the  percent  cover  that 
90%  of  the  plates  would  have  to  remain  greater  than  96%  for 
the  entire  range  of  the  simulations  (Figure  28A).  The 
standard  deviation  about  the  mean  percent  cover  that  these 
plates  would  have  showed  a  precipitous  decrease  out  to  the 
20  plates  per  group  simulated  point  and  then  decreased  only 
2%  out  to  the  80  plates  per  group  simulation  point  (Figure 
28B)  . 


73 


o 

o 


VH 

in 
d~ 

• 

^4 

^ 

hfc, 

> 

(R)  52- 

J^ 

S  — 

^ 

21 

a 

in 

•^ 

<l 

> 

C4 

^ 

► 

8 

o 

__.,.,.,.,..   .  „ 

T 

i 

1 

D18F        D18B        D25F 
PLATE  DESIGNATION 


D25B 


o 
o 


(B) 


E— 

en  o 
a:  "■> 
-J  d" 


en 


LO 

rsi 


o 

o 


Rl 9F  R 1 9B  R23F  R23B  B 18F  B 18B  B22F  B22B  C8F 
PLHTE  DESIGNATION 


C8B  C20FC20B 


Figure  26.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  8. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


74 


o 

81 


in 


CiJ 

> 
o 

i—  O 
•2L  » 
LJ 
O 
□C 

a. 


o 
o. 


0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   Of  PLATES 


Figure  27.  Computer  Simulations  Using  Data  From  Month  8 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as    a   Function   of    Increasing    Sample    Size. 


75 


[R3 


in- 


en 

CiJ 

> 
o 

°q 

(-<  d 

o 
en 

LJ 
Q_ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   OF  PLATES 


81 


(B) 


q 
in- 

CiJ 

> 
O 

°o 

E—  <=>■ 
Z  ~ 

GJ 
O 

cn 

UJ 

q 
in' 


0.0  5.0 


-I 1 "T 

10.0  15.0  20. 


25.0  30.0 

NO. 


i         r 

35.0  40 


r 

45. 


0  45.0 

OF  PLATES 


i         r 
50.0  55 


I  I  I  F  I 

,0  60.0  65.0  70.0  75.0  80. 


Figure  28.  Computer  Simulations  Using  Data  From  Month  8 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have   As    a   Function   of    Increasing   Sample   Size. 


76 


I.   RESULTS  FROM  MONTH  9 

1.  Experimental  Data 

The  bryozoan  Watersipora  cuculatta  dominated  the 
fouling  assemblages  on  three  of  the  four  control  surfaces 
while  the  hydroid  Obelia  spp.  dominated  on  the  fourth.  The 
total  number  of  species  represented  in  the  census  began  to 
stabilize  as  the  more  successful  species  excluded  others  in 
the  competition  for  the  rapidly  diminishing  space. 

The  dominant  fouling  organism  on  the  antifouling 
coated  surfaces  was  the  bryozoan,  Membranipora  membranacea. 
This  organism  is  usually  found  nearly  exclusively  on  the 
fronds  of  the  giant  kelp  Macrocystis  and  it  is  not  known 
what  chemical  or  other  stimulus  attracted  it  to  the 
antifouling  coated  surfaces. 

The  percent  coverage  estimates  for  the  control 
surfaces  ranged  from  35%  to  92%.  The  mean  percent  cover  was 
67%  and  the  standard  deviation  was  25%. 

The  similarity  values  for  the  plates  ranged  from 
.18  to  .74  for  the  control  surfaces  and  from  .60  to  1.0  for 
the  antifouling  coated  surfaces.  The  mean  similarity  value 
for  the  control  surfaces  was  .41  (Figure  29A)  and  was  .89 
for  the  antifouling  coated  surfaces  (Figure  29B) . 

2.  Computer  Simulations 

The  final  estimate  for  the  mean  percent  fouling 
cover  for  the  computer  simulated  plates  was  73%.  The  upper 


77 


s 


P4 

in 

•  *■ 

< — 

«ri 

t^ 

o 

^1 

> 

^ 

ID 

ILRRITY 

0.50 

31 

^1 

E% 

en 

^sT 

t> 

in 

^ 

< 

> 

_J^"d 

is 

CN 

j^*"< 

o 
o 

o 

_                        T                      __ .                      . 

D1F         DIB         D3F 
PLATE  DESIGNATION 


03B 


S 


a:  o 

(B)52- 

>— i 
CO 


o 

a 


> 


■  ij  ►- 


> 


nilFRHBR24FR24B  B1F    BIB  B21FB21B  C5F    C5B    C7F    C7B 
PLATE  DESIGNATION 

Figure  29.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  9. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


78 


95%  confidence  limit  about  this  mean  showed  little  decrease 
past  the  20  plates  per  group  simulation  point  (Figure  30) . 

The  95%  confidence  limit  on  what  one  can  say  90%  of 
the  simulated  plates  will  have  fouling  coverage  less  than 
or  equal  to  remained  above  95%  coverage  for  the  entire 
range  of  the  simulations  (Figure  31A).  The  standard 
deviation  that  these  plates  had  showed  little  change  past 
the  20  plates  per  group  simulation  point  (Figure  31B) . 

J.   RESULTS  FROM  MONTH  10 

1.  Experimental  Data 

Bryozoans  dominated  the  fouling  assemblages  on  all 
of  the  control  plates  with  Watersipora  cuculatta  occupying 
60%  of  the  space  on  one  of  the  surfaces.  Eleven  species  of 
bryozoans  were  identified  during  the  census. 

The  similarity  values  for  the  plates  ranged  from 
.20  to  .79  for  the  control  surfaces  and  from  .67  to  1.0  for 
the  antifouling  coated  surfaces.  The  mean  similarity  value 
for  the  control  surfaces  was  .47  (Figure  32A)  and  for  the 
antifouling  surfaces  was  .90  (Figure  32B). 

2.  Computer  Simulations 

The  final  iterated  value  for  the  computer 
simulation  of  the  mean  percent  cover  was  70%.  The  upper  95% 
confidence  limit  about  this  mean  showed  little  decrease 
past   the   20    plates    per    group    simulation   point    (Figure    33). 


79 


o 

8-1 


a: 

> 
o 

"o 

E-h  o 
■Z.  "» 

LJ 
O 

cc 

CjJ 


CM 


o 
d 


Vn 


— I 1 1 1 1 1 1 1 1 1 1 1 1 1        1        I 

0.0  5.0  10. 0  15. 0  20-.0  25".0  30. 0  35. 0  40. 0  45. 0  50. 0  55. 0  60. 0  65. 0  70. 0  75. 0  80.0 

NO.   OP  PLATES 


Figure  30.  Computer  Simulations  Using  Data  From  Month  9 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as    a   Function   of    Increasing    Sample    Size. 


80 


o 

S-i 


U1- 

LJ 

> 
o 

°° 

S-<  d- 

tmg 

LJ 
Q_ 


in- 


o 

CD. 


1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 

0.0  5.0    10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.    OF  PLATES 


O-i 

to 


xn- 


o 

en  d- 
LJ  oi 

> 

o 

f— '  ""> 

(B)|" 

£°. 

Q.  2 


o 
in' 


o 
d. 


~~ t 1 1 1 1 \ 1 1 1 1 i ; \ 1 1 1 

0.0  5.0    10.0  15.0  20.0  25.0  30.0.35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 

NO.   OP  PLATES 

Figure  31.  Computer  Simulations  Using  Data  From  Month  9 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have  As    a   Function   of    Increasing   Sample   Size. 


81 


8 


m 

«« 

► 

• 

• 

O 

1— « 

a:  o 

(FN  52- 

4£ 

M 

► 

*et 

I!* 

51 

^1 

fcb. 

p* 

en 

^ 

^ 

in 

CM 

^ 

— > 

o 
o 

a 

.       

r 

i 

1 

D5F         D5B        D16F 
PLATE  DESIGNRTION 


D16B 


o 

o 


in 


(B)  m- 


-J  a 


CO 


d 


o 

o 


fl9F    A9B  R17FR17B  B9F    B9B  B12F  B12B  CI  IF  CI  IB  C16F  C16B 
PLATE  DESIGNATION 

Figure  32.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  10. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


82 


in- 


OJ 

> 
O 

°o 
Z  « 

o 

LJ 
Cl- 


in' 

CM 


O 

d. 


— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


NO.   OF  PLATES 


Figure  33.  Computer  Simulatioms  Using  Data  From  iMonth  10 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as    a   Function  of    Increasing    Sample    Size. 


83 


The  percent  cover  that  one  can  say  90%  of  the 
plates  would  have  coverage  less  than  or  equal  to  (with  95% 
confidence)  remained  above  90%  cover  for  the  entire  range 
of  the  simulations  (Figure  34A).  The  standard  deviation 
these  plates  would  have  showed  little  decrease  past  the  20 
plates  per  group  simulation  point  (Figure  34B) . 

K.   RESULTS  FROM  MONTH  11 
1.  Experimental  Data 

The  fouling  communities  on  the  control  surfaces 
remained  dominated  by  bryozoans  with  Watersipora  cuculatta 
dominating  on  three  of  the  plates  and  the  upright  bryozoan 
Bugula  californica  dominating  on  the  fourth.  Of  the  seven- 
teen species  identified  on  these  control  surfaces,  thirteen 
were  bryozoans. 

The  predominant  fouling  organism  on  those 
antifouling  surfaces  that  showed  fouling  was  the  bryozoan 
Membranipora  membranacea.  The  hydroid  Obelia  spp.  and  the 
spirorbid  worm  Circeis  armor icana  were  also  in  evidence. 

The  initial  percent  coverage  estimates  for  the 
control  surfaces  ranged  from  64%  to  98%.  The  mean  value  was 
85%  and  the  standard  deviation  was  15%. 

The  similarity  values  for  the  plates  ranged  from 
.28  to  .79  for  the  control  surfaces  and  from  .64  to  1.0  for 
the  antifouling  coated  surfaces.  The  mean  similarity  value 


84 


o 

8-1 


(R) 


»■». 

LJ 

> 

o 

E-h  d- 

LJ 
O 

cc 

LJ 
Cl- 


in' 


o 

d. 


0.0  5.0 


10.0  15.0  20.0 


25.0  30.0  35.0  40. 


45.0  50.0  55.0 


60.0  65.0  70.0  75.0  80.0 


NO.   OP  PLATES 


to 


rg 


(B) 


o 

LJ  ™ 

> 
O 

t— i  in- 

•z.  — 

LJ 
O 


o 

t 

in' 


o 
d. 


T 


T 


0.0  5.0    10.0 


r 
15. 


r 
20. 


i         r 
25.0  30 


.0 

NO 


—\ i i i i i i i i i 

35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 


OF  PLATES 


Figure  34.  Computer  Simulations  Using  Data  From  Month  10 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples   Will  Have   As    a   Function   of    Increasing   Sample   Size. 


85 


for  the  control  surfaces  was  .52  (Figure  35A)  and  for  the 
antifouling  coated  surfaces  was  .87  (Figure  35B) . 
2.   Computer  Simulations 

The  final  iterated  value  for  the  computer 
simulations  of  the  mean  percent  cover  was  89%.  The  upper 
95%  confidence  limit  about  this  mean  decreased  by  only  7% 
over  the  entire  range  of  the  simulations  (Figure  36) . 

Quite  obviously,  with  such  a  large  valuefor  the 
mean  percent  cover,  the  percent  cover  that  90%  of  the 
plates  will  have  will  be  equally  large.  As  can  be  seen 
(Figure  37A)  ,  this  value  never  became  less  that  97%  cover 
over  the  entire  range  of  the  simulations.  The  standard 
deviation  about  the  mean  percent  cover  that  these  plates 
would  have  again  showed  little  decrease  past  the  20  plates 
per  group  simulation  point  (Figure  37B) . 


86 


s 


(fl) 


in 


t- 

OH   o 

a:  w 


z: 

I— I 

CO 


in 

CM 


o 
o 


< 


> 


> 


D18F        D18B        D25F 
PLATE  DESIGNATION 


D25B 


o 
o 


(B) 


1/1 


en  a 
CE  w 
-J  o" 

21 

CO 


in 


o 
a 


fe>^  ^^  |»<d  !► 


R4F    A4B  A20FA20B  B3F    B3B  B17FB17B  C2F    C2B  C23F  C233 
PLATE  DESIGNATION 


Figure  35.  Similarity  Graphs  for  (A)  Non-treated  Control 
Surfaces  and  (B)  Anti-fouling  Coated  Surfaces  for  Month  11. 
Dashed   Lines    Indicate   Mean   Similarity  Values. 


87 


o~ 

•■" 

\ 

* 

\  «, 

o 

in- 

r* 

Ql 

CjJ 

> 

O 

■ 

°o 

f_  o- 

Z  H 

LJ 

O 

cc 

tx3 

Q_ 

O 

ui- 

CM 

o 

o 

1         ]         1         1         1         1         1         1 

1          1          1          1          1          1          1          1 

0.05.0  10.0  15.0  20.025.0  30. 


NO.   OP  PLATES 


Figure  36.  Computer  Simulatioms  Using  Data  From  Month  11 
Showing  the  Expected  Value  of  the  Mean  Percent  Fouling 
Cover  (the  Mean  of  the  200  Individual  Group  Simulation 
Percent  Fouling  Covers)  as  a  Solid  Line  and  the  95% 
Quantile  (Dashed)  of  the  Expected  Mean  Percent  Fouling 
Cover    as   a   Function   of    Increasing    Sample    Size. 


88 


o 

Si 


tfl] 


in 

0£ 

UJ 

> 

o 

LJ 

o 

ce 
lj 

Cl- 


in- 

CM 


o 
d. 


0.0  5.0 


—i 1 1 1 — 

10.0  15.0  20.0  25.0 


o 
in-. 


- 1 1 1 1 1 1 1 1 1 1 1 

30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  70.0  75.0  80.0 

NO.   OF  PLATES 


(B) 


o 

CC  d 

Oj  -• 

> 

o 
o 


LJ 
O 

LJ  o 
Q_  in' 


o 
d. 


0.0  5.0    10.0 


i         i         r 
15.0  20.0  25. 


30.0  35.0  40.0 


50.0  55 


1 1 — 

0  60.0  65.0 


— I 1 1 

70.0  75.0  80.0 


NO.    OF  PLATES 


Figure  37.  Computer  Simulations  Using  Data  From  Month  11 
Showing  :  (A)  the  Percent  Fouling  Coverage  That  One  Can  Say 
90%  of  the  Samples  Will  Have  Fouling  Coverage  Less  Than  or 
Equal  to  (  With  95%  Confidence);  and  (B)  the  Standard 
Deviation  of  Percent  Coverage  (About  the  Mean)  That  These 
Samples  Will  Have  As  a  Function  of  Increasing  Sample  Size. 


89 


V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   DISCUSSION 

The  choice  of  an  appropriate  sampling  strategy  for  the 
study  of  fouling  organisms  in  Monterey  Bay  has  been  shown 
to  be  dependent  upon  the  type  of  information  desired. 

By  analyzing  the  width  of  the  upper  95%  confidence 

limit  about  the  expected  value  of  the  mean  percent  fouling 

cover,  estimates  were  made  concerning  the  optimum  number  of 

plates  to  be  deployed.  This  optimum  number  was  found  by 

determining  that  point  where  the  addition  of  additional 

plates   had  a  negligible  effect  upon  the  width  of  the 
confidence  interval.  In  nearly  all  the  cases  simulated,  the 

value  for  the  optimum  number  of  plates  to  deploy  appeared 

to  be  approximately  twenty. 

A  similar  procedure  was  followed  to  estimate  the 

optimum  number  of  plates  to  deploy  to  minimize  the  standard 

deviation  about  the  mean  percent  fouling  cover  that  90%  of 

the  plates  would  have.  The  optimum  number  of  plates  to 

deploy  to  satisfy  this  requirement  was  again  estimated  to 

be  approximately  twenty  for  most  of  the  simulations. 

However,  the  results  from  several  cases  suggested  that 

thirty  plates  would  be  a  more  appropriate  sample  size. 


90 


Based  on  these  results,  it  is  concluded  that  twenty 
plates  is  probably  the  minimum  number  of  plates  that  should 
be  deployed  to  obtain  accurate  estimates  of  the  mean 
percent  fouling  cover  of  a  group  of  plates  for  this 
locality.  Thirty  plates  per  group  is  probably  a  more 
appropriate  number  of  plates  to  deploy  to  insure  that 
accurate  results  are  obtained  for  groups  that  display 
heightened  variability  in  the  individual  plate  fouling 
coverage  estimates. 

The  experimental  results  showing  the  negligible  effect 
of  the  addition  of  more  plates  on  what  one  could  say  90%  of 
the  plates  would  have  fouling  coverage  less  than  or  equal 
to  (with  95%  confidence)  lead  to  the  conclusion  that  the 
inherent  variability  of  fouling  populations  is 
significantly  greater  than  previous  studies  indicated.  This 
means  that  while  twenty  to  thirty  plates  are  probably 
sufficient  to  resolve  ambiguities  concerning  the  mean  per- 
cent fouling  cover,  this  number  is  clearly  insufficient  to 
ascertain  with  any  high  degree  of  confidence  the  amount  of 
variability  of  a  large  segment  of  the  population.  For 
example,  if  it  is  desired  to  ascertain  with  95%  confidence 
the  variability  of  90%  of  the  fouling  population,  it  must 
be  understood  that  this  will  require  the  committment  of 
substantial  resources  to  the  study. 


91 


In  addition  to  the  development  of  an  appropriate 
fouling  community  sampling  strategy  for  Monterey  Bay,  more 
far  reaching  conclusions  can  be  drawn  concerning  the 
applicability  of  the  procedures  used  in  this  thesis  to 
other  locations.  Since  the  computer  modelling  procedure 
used  in  this  study  dealt  with  the  extension  of 
the  experimentally  observed  variability  and  contained  no 
site-specific  parameters,  it  is  believed  that  the  procedure 
can  be  directly  applied  to  the  study  of  fouling  community 
variability  at  any  desired  geographical  location.  This 
means  that  the  computer  programs  developed  for  this  study 
can  be  coupled  with  archived  fouling  coverage  data  from 
any  site,  depth,  season,  and  so  on  to  provide  information 
on   how  best   to   sample   the    fouling   populations. 

The  final  conclusion  drawn  from  this  study  is  that  the 
bootstrap  method  of  computer  intensive  statistical  analysis 
has  profound  implications  in  the  study  of  other  biological 
problems.  It  is  believed  that  this  method  could  be  applied 
to  a  wide  range  of  other  data  intensive  biological  areas 
including  fisheries  management,  larval  settlement  studies, 
the  pelagic  distributions  of  plankton  and  nekton,  and 
growth  rate  studies  to  name  just  a  few. 

B.       RECOMMENDATIONS    FOR   FURTHER    RESEARCH 

The  model  developed  in  this  thesis  should  be  applied  to 
a    number    of    other    geographical    locations    to    determine    the 

92 


differences  in  fouling  variability  in  terms  of  latitude, 
longitude,  depth,  environmental  stresses,  and  a  host  of 
other  biological  forcing  functions.  Since  the  computer 
program  already  exits  and  the  only  required  inputs  are 
archived  fouling  coverage  data,  this  should  present  no 
insurmountable  difficulties. 

Once  such  studies  have  been  completed,  the  development 
of  empirically  derived  sampling  strategies  for  any  location 
could  be  attempted. 

Finally,  the  use  of  the  bootstrap  or  other  computer 
intensive  statistical  techniques  should  be  vigorously 
investigated  in  other  areas  of  biological  interest.  It  is 
believed  that  the  use  of  these  techniques  might  well 
provide  answers  to  a  wide  range  of  biological  problems  that 
have  so  far  proved  intractable. 


93 


APPENDIX  A 


MICRON  22  ORGANO-METALLIC 
POLYMER  ANTIFOULING  PAINT 


INGREDIENTS  PERCENT  BY  WEIGHT 

Active: 

Bis  (tributyltin)  Oxide  11.7 

Cuprous  Thiocyanate  17.2 

Inert:  71.1 


100 


Elemental  Tin    4.4% 

Elemental  Copper  8.9% 

Paint  contains  1.1  lbs  of  Bis ( tr ibutyltin)  oxide  per 
gallon  and  1.6  lbs  of  Cuprous  Thiocyanate  per  gallon. 


Source:  Product  infromation  breakdown  on  label 


94 


APPENDIX  B 

NAVY  STANDARD  FORMULA  121 
RED  VINYL  ANTIFOULING  PAINT 


INGREDIENTS  PERCENT  BY  WEIGHT 

Cuprous  Thiocyanate  70.3 

Rosin  10.5 

Vinyl  Resin  2.7 

Tricreysl  Phosphate  2.4 

Methyl  Isobutyl  Ketone  8.1 

Xylene  5.6 

Antisettling  Agent  .4 


Source:  Department  of  the  Navy  Specification  MIL-P-15931C, 
Painty  Antifouling, Vinyl (Formula  Numbers  121  and  129) 


95 


APPENDIX  C 

NAVY  STANDARD  FORMULA  170  BLACK 
CAMOFLAGE  ANTIFOULING  PAINT 


INGREDIENTS  PERCENT  BY  WEIGHT 

Vinyl  Resin  17.5 

Bis  ( tributyltin)  oxide  4.2 

Tributyltin  Fluoride  18.1 

Carbon  Black  2.1 

Titanium  Dioxide  .8 

Ethylene  Glycolmonoethyl  3.0 
Ether  Acetate 

Normal  Propanol  11.1 

Normal  Butyl  Acetate  43.2 

Source:  Department  of  the  Navy  Military  Specification  DOD- 
P-245  8  8  ,  Paint ,  Ant  i  foul  ing,  Vinyl,  Camoflage(  Formula  numbers 

170, 171, 172, and  173) ,  2  May  1979. 


96 


APPENDIX  D 
ZYNOLYTE  EPOXY  RUST  MATE  PAINT 


INGREDIENTS  PERCENT  BY  WEIGHT 

Non-Volatile  (58.4%  of  total) 

Pigments  43.4 

Vehicle  56.6 

Epoxy  and  Menhaden 
Alkyd  Resins 

100. 


Volatile  (41.6%  of  total) 

Exempt  Mineral  Spirits  98.0 

Aromatic  Hydrocarbons  2.0 

100. 

Source:  Product  ingredient  breakdown  on  label 


97 


APPENDIX  E 

A  DISCUSSION  OF  THE  METHOD  OF  MAXIMUM  LINELIHOOD  AND  ITS 
INCORPORATION  INTO  A  MODEL  FOR  FOULING  COVER 


A.   GENERAL 

The  following  model  was  used  to  describe  the  properties 
of  the  proportion  of  a  plate  that  was  fouled.  Each  plate 
(t)  has  a  random  proportion  of  fouling, P  .  Given  P  the 
number  of  the  100  censused  sites  that  are  fouled  ,  (S  ),  on 
the  tth  plate  has  a  binomial  distribution: 

p{st  =  k)  =  c1^)  ptk  a  -  pt)100-k  (1> 

for  k  =  0,1,2. ,100  independent  of  the  other  plates.  The 
proportion  of   fouling, Pt/    is   assumed   to   have   the   form 

P     -         S  (2) 

(1  +  e  Z)      . 
where  e  has    a    normal    distribution    with    mean    fi        and 

variance  a2  independent  of  the  other  plates  censused  at 
the  same  time.  The  mean  fi  and  variance  <J  will  in  general 
be  a  function  of  the  amount  of  time  the  plate  is  submerged 
though  the  many  variables  involved  in  population  dynamics 
will   keep   the   function   from   being   linear. 

Assume  N  plates  are  inspected  at  a  time  with  the  result 
that  S      of   the  censused  sites  are  fouled  on  the  ttn  plate. 

98 


The   likelihood   function   for    the   model   is: 

N 

S.  100-S 

L--J\    <Ket;  u,c  }   C1™)  Ptt    (1-  Pt)  (3) 


t=l 


where 


^e^a2}  =  —  e  2(~       }       for  -  <  et  <  -  <4> 

Note  that  L  is  the  probability  of  observing  S„  ,..,SM 

1      N 

fouled  sites  on  each  of  the  plates. 

The  method  of  maximum  likelihood  is  to  find  those 

values  of  /i  ,  a   ,  and  €      ,  for  t  =  1,2, ...N  which  maximize 

L;  that  is,  those  values  which  maximize  the  probability  of 

observing  the  outcome.  These  values  also  maximize  the  log 

likelihood  function,  In  L,  where 
N  2 

X'   Jin  L  =  V1  ( £— j +  St2,nPt  -  £na  +  (100-StHn(l-Pt)}  (5) 

t=l  [+  constant] 

To  find  the  maximum  of  <X   the  partial  derivatives  of  J£ 
with  respect  to  fj.     ,  a     and  €t   are  set  equal  to  zero.  This 

results  in  the  equations: 

N 


t=l 


3y  ~  2-j  ' 


2 
a 


-   0 


solving  for  /i   results  in: 

N 


■*E 


£t 


(6) 


t=l 

99 


similarly  M 

.2 


3^  .V1     /£t-^"       1 


*  -E  «=* 


a3 
t=l 

solving    for    a2     results    in 

N 


-i}   =  0 


(7) 


a  =n2j  (£t - y) 

t=l 

and 

3y  Ce     -  y)       S  3P 

If  =--V~  +  ^^Pt  +  (100"Vn^y[-^i]  =  0    (8) 

t  a  t      t  t  t 

where  K>        9  £t  £t 

TSF1  =  4^-  C-V-3   =  —5- =  P^l-P.) 

tf£  d£  E.  £       „  t  t 

r  r     1+e  T         (1+e  V 

Equation    (8)    simplifies    to 


£ 

e 


2^  o     2  .  n  (9) 


e-     +  — (100  a  )  -  \i  -  S.a     =  0 

t  et  t 

1+e 

B.   SOLUTION  OF  THE  LOG  LIKELIHOOD  EQUATION 


A  recursive  method  was  used  to  solve  the  system  of 
equations  (6)  ,  (7)  ,  and  (9).  For  each  month,  the  100  sites 
were  censused  on  each  of  the  four  untreated  plates  and  S. 
the  number  of  fouled  sites  was  determined.  Initial  values 

A  A  «*  A  9 

for    p      ,  €.  /    ]i  ,    and     a      were  determined   as    follows 


Pt(0)    =   V100  (=PNOT, 


100 


£  (0)  =  HnCPt(0)/(l-Pt(0))]     (=EPSNOT) 

Jj(0)  =  i  2_   £t(0)  (=MUNOT) 

t=l 
4- 

S(0)  =  ?  H   (£t(0)  "  ^C°))2   (=SIGNOT) 
t=l 

The  values  /zfo)  and  (P(o)were  then  used  in  equation  (9)  and 
the  equation  solved  for  £.(1),  for  t  =  1,2,3,4.  Since 
equation  (9)  is  transcental,  Newton's  method  was  used  to 

A 

solve  for  £(1) .  Newton's  method  is  an  iterative  procedure 
in  which  the  (n+l)st  member  of  the  iterative  sequence  (Xn) 

is:  fCXn) 

Xn+1  =  Xn  +  FTXT 


n 


where  in  our  case 


X 
fCX  )  =  X„  +     e  y     (100  S2(0))  -  y(0)  -  StS2(0) 


n  n  X 

1+e 
and 


x~  Xru2.  -~~  "2, 


f'CX  )   =  1  +  [en/(l+en)Z]  100  a   (0) 
n 

Equation    (10)    was    iterated   until    iv        -  X  I    <  1  x  10~ 

•   n+1        n1 

The       e  (1),     t   =   1,2,3,4      were    then    used    to    compute    iterated 
values   of      \i   ,    G  ,    and     p       as    follows: 

e.(l)  £t(D 

Pt(l)   =  e  r       /(1+e  T       ) 


101 


^(1)  =  k  Yj    £t(1) 

t=l 


$1)   =i  ]P  (£t(1)  -  '^(1))2 


t=l 

and  equation  (9)  was  then  solved  with  M     =      /J.  (1)  and  a 
=  a   (1)  for   €f(2)  ,  t  =  1,2,3,4.  The  iterative  procedure 

was  continued  until  the  achievement  of  a  tolerance  value  of 

"6         1  , 

(1x10  )  calculated  by  |u(k)  -  u(k+l;  [  indicated  the  system  of 

equations  had  converged.  The  final  iterated  values  for 

Pt  (=PNEW),   €  (=EPSNEW),  a2   (=SIGNEW)  ,  and   fl   (=MUNEW) 

were  the  values  that  maximized  the  likelihood  function  for 

the  model.  SIGNEW  (a)     and  MUNEW  (£)  were  then  used  as 

initial  input  values  in  the  bootstrap  simulations  of  the 

experiment  (Appendix  F) . 


102 


APPENDIX    F 
BOOTSTRAP    COMPUTER   SIMULATIONS 

A.  AN  EXPLANATION  OF  THE  BOOTSTRAP  METHOD  OF  COMPUTER 
SIMULATIONS  OF  RANDOM  PROCESSES  AND  ITS  INCORPORATION 
INTO   AN   ASSUMED   STOCHASTIC   MODEL    FOR   FOULING    COVERAGE 

1.      Discussion 

The    common   statistical    tools    utilized    in    the    study 

of    biofouling    have    as    their     basis       the     simplifying 

assumption   that    the   data   collected    from    the   analysis   of 

such    communities    can    be    described    by    a    normal    or    Gaussian 

distribution.      That     is     to     say,      it     is     assumed     that 

fluctuations    in    the   values   of    some    experimentally   observed 

parameter    are    scattered    symmetrically   about    the    true   value 

of   the   parameter.    It   is   further    assumed    that   the   larger    the 

difference    between    the    the    true    value    and    the    observed 

value    of    the    parameter,     the    less    likely    it    is    that    the 

value    will    be    observed    experimentally    (Diaconis    and 

Efron,1983).      Many     years     of     experience     using     these 

assumptions     have     shown     that     even     if    the    dataare    only 

approximately  or    "pseudo"   normal,    the  Gaussian   theory   still 

works    quite    well.    If    however,    the    data   do   not    satisfy      the 

requirements    for    the    assumption   of    normality    or, if    the 

sample   size    is    such    that    the   various    tests    used    to  check 

for    normality   can    only   give    ambiguous    results,     it    is    clear 


103 


that  the  results  of  statistical  techniques  based  on  the 
assumption  of  normality  will  be  unreliable. 

Recent  developments  in  the  use  of  computer- 
intensive  techniques  for  statistical  analysis,  particularly 
the  invention  of  the  bootstrap  technique  (Ef ron,1977) ,  have 
enabled  the  computation  of  various  statistical  parameters 
without  the  necessity  of  assuming  a  Gaussian  distribution. 
This  technique  has  also  enabled  the  computation  of  those 
statistics  which  do  not  have  a  simple  analytical  formula. 
Prior  to  the  advent  of  large  main  frame  computers,  the 
difficulty  of  finding  numerical  solutions  to  non-linear 
problems  forced  statistical  methods  to  concentrate  on  those 
statistical  models  and  procedures  for  which  analytical 
results  could  be  obtained.  These  models  and  procedures  did 
give  useful  large  sample  size  results  for  the  common 
statistics  such  as  the  mean  and  variance  of  a  population. 
However,  they  ignored  other  important  statistical  questions 
that  did  not  have  analytical  formulas;  such  as,  the  degree 
of  variability  in  an  estimate  due  to  the  sample  being  of 
finite  size  or  the  confidence  limits  about  estimates  from  a 
finite  sample  (Diaconis  and  Efron,1983).  In  general  terms, 
the  bootstrap  method  consists  of  coupling  a  probabilistic 
model  with  the  data  gathered  experimentally  from  one  sample 
of  size  N  to  generate  a  large  number  of  simulated  samples 
of  size  N.  These  simulated  samples  are  then  analyzed  to 


104 


determine    the    variability   of    estimates    of    the    true    values 
of    the    statistics   of    the  population. 
2.      Procedure 

The  model  described  in  Appendix  E  was  used  in  a 
simulation  study  of  the  variability  of  the  estimate  of  the 
fraction  of  plate  coverage.  The  model  was  used  to  simulate 
200  groups  of  fixed  numbers  of  plates.  The  inputs  into  the 
simulation  were  the  values  for  the  maximum  likelihood 
estimates     of     the     mean      (MUNEW  11    )      and     variance 

(SIGNEW  a2)  for  the  monthly  epsilon  ( 6  )  values 
determined  from  the  analysis  of  the  fouling  cover  on  the 
four    untreated    plates. 

The  stochastic  model  used  in  this  case  assumed  a 
normal  distribution  for  the  epsilon  ( e.)  values  (  that  is 
€.  ~  N( n  f  a ) ) .  Using  available  computer  software,  the 
required  number  of  normally  distributed  random  numbers  with 
mean  equal  to  zero  and  variance  equal  to  unity  were 
generated   for   200    groups   of    the    following    number    of   plates: 


Number   of  Simulated  Required  Number (I) 

Plates  per   Group  of  Random  Numbers 


2 

I    =      400 

4 

I   =      800 

5 

I   =  1000 

10 

I    =   2000 

105 


Number  of  Simulated     Required  Number (I) 
Plates  per  Group        of  Random  Numbers 


15 

I  =  3000 

20 

I  -  4000 

30 

I  =  6000 

40 

I  =  8000 

80 

I  =16000 

Eacn  of  the  (I)  computer  generated  random  numbers 
from  the  standard  normal  distribution  were  then  transformed 
into  random  numbers(n^)  with  a  normal  distribution  with 
mean  fi  and  standard  deviation  o  by  multiplying  each 
invariable  by  a  (the  square  root  of  a  )  and  adding  //  . 
The  random  number  n^  was  transformed  to  give  PINIT(I) 
defined  as: 

PINIT(I)  =   eni/(l  +  eni)  ,  the  average  proportion  of 
the  ith  plate  fouled. 

The  initial  simulated  percent  cover  for  the  ifc" 
plate  was  determined  by  calculating  a  random  number  having 
the  properties  of  a  binomial  distribution  with  probability 
(P)  equal  to  PINIT(I)  and  N  equal  to  100.  The  resulting 
variable  was  termed  PNOT(I).  The  variables  EPSNOT(I) , MUNOT, 
SIGNOT,  PNEW(I)  ,EPSNEW(I)  ,  MUNEW,  and  SIGNEW  were  then 
calculated  using  the  method  of  maximum  likelihood  described 
in  Appendix  E. 


106 


APPENDIX  G 

TABULATED  MONTHLY  PERCENT  FOULING  COVERAGE  VALUES 
FOR  THE  NON-TOXIC  CONTROL  SURFACES 


MONTH 


PLATE  # 

1 
2 
3 
4 


FRACTIONAL  COVERAGE 

.50 
.50 
.24 
.54 


MONTH 


PLATE  £ 

1 
2 
3 

4 


FRACTIONAL  COVERAGE 

.33 
.49 
.69 
.46 


MONTH 


PLATE  # 

1 
2 
3 
4 


FRACTIONAL  COVERAGE 

.73 
.02 
.07 
.84 


MONTH 


PLATE  £ 

1 
2 
3 
4 


FRACTIONAL  COVERAGE 

.13 
.96 

.61 
.81 


MONTH 


PLATE  £ 

1 
2 
3 
4 


FRACTIONAL  COVERAGE 

.54 
.46 
.41 
.06 


MONTH 


PLATE  # 

1 
2 
3 
4 


FRACTIONAL  COVERAGE 

.54 
.47 
.62 
.31 


107 


MONTH  PLATE  £  FRACTIONAL  COVERAGE 

1 
8  2 

3 
4 


MONTH  PLATE  # 

1 

9  2 

3 

4 

MONTH  PLATE  £ 

1 

10  2 

3 
4 


MONTH  PLATE  £ 

1 
11  2 

3 
4 


.48 

.88 

.69 

.97 

FRACTIONAL  COVERAGE 

.91 

.49 

.92 

.35 

FRACTIONAL  COVERAGE 

.70 

.56 

.29 

.97 

FRACTIONAL  COVERAGE 

.83 

.95 

.64 

.98 

108 


APPENDIX  H 

LIST  OF  THE  SESSILE  SPECIES  IDENTIFIED  BY  THE  RANDOM  POINT 
CENSUS  AND  THE  MONTHS  THEY  WERE  PRESENT  ON  THE  NON-TOXIC 

CONTROL  SURFACES 


ORGANISMS 


MONTH    NUMBER 
23456789      10      11 


Protozoa: 

Folliculina  spp. 
Ephelota  gemmipara 

Coelentrata: 

Obelia  spp. 
Hydractinia  spp. 


Anthopleura  spp. 

Ectoprocta: (Bryozoans) 

Bugula  neretina 
Bugula  californica 
Watersipora  cucullata 
Hippothoa  hyalina 
Celloporaria  brunnea 
Cryptosula  pallianasa 
Schizoporella  unicornis 
Microporella  ciliata 
Microporella  californica 
Mernbranipora  membranacea 
Membranipora  serilamella 
Unknown  bryozoan  #1 
Unknown  bryozoan  #2 

Annelida: 

Circeis  armor icana 
Janua  nipponica 
Pileolaria  potswaldi 
Protolaeospira  exima 
Serpula  vermicular  is 
Anatides  groenlandica 

Arthropoda: 

Balanus  crenatus 
Megabalanus  californicus 
Amphipod  (unknown) 

Mollusca: 

Mytilus  edulis 


XXX 
X 


X   X 


X 


X   X 


X 


X 

X 

X 

X 

X 
X 

X 

X 
X 

X 

X 

X 
X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 
X 

X 

X 
X 
X 
X 

X 
X 

X 
X 
X 


XXX 
X 
X 


X 


X 


X 


X   X 


X 
X   X 

X 


X 


X 


X 


X   X 


109 


MONTH    NUMBER 
ORGANISMS  234567891011 

Echinodermata : 

Strogylocentrotus    spp.  X 

Chorda ta : 

Ascidia  ceretodes  X     XXX 

Styela  truncata  X 

Pyura  haustor  X 


110 


LIST    OF    REFERENCES 


Boyd, M., 1972.  Fouling  Community  Structure  and  Development 
in  Bodega  Harbor,  California.  Doctoral  Dissertation.  Unive- 
rsity   of    California, Davis. 

Diaconis,P.,  and  Efron,B.,1983.  Computer-Intensive  Methods 
in  Statistics.Scientif  ic  Amer  ican.,v.248,No.5;116-130. 

Ef ron,B.,1979.  Computers  and  the  Theory  of  Statistics: 
Thinking  the  Unthinkable.  I_n  SIAM  Rey_iew,  v. 21, No. 4. Octo- 
ber ,1979  :  460-480  . 

Fisher, E.C.,  Birbaum,L.S.  ,Depalma,  J. ,  Mur  acka,  J.S. ,  Dear  H. , 
and  Wood,F.G.  1975.  Survey  Report:  Navy  Biological  Fouling 
and  Biodeterioration.  Naval  Undersea  Center  Report  NUC-TP- 
456. 

Fraser  ,C.M. ,  1 9  3  7 .  Hy_dr_oids  o_f  the  Pacific  Coast  of  Canada 
and   the   United   States,    University   of   Toronto   Press. 297pp. 

Haderlie,E.C.  1974.  Growth  Rates,  Depth  Preference,  and 
Ecological  Succession  of  Some  Sessile  Marine  Invertebrates 
in  Monterey  Harbor.    Veliger. v.17    (supplement) :l-35. 

Kelley,P.R.,1981.  Scanning  Electron  Microscope  Observations 
of  Marine  Microorganisms  on  Surfaces  Coated  with 
Antifouling  Paint.  Master's  thesis,  Naval  Postgraduate 
School. Monterey,     California,     U.S.A. 

Knight-Jones,  P. , and  Knight-Jones  ,  E.  W.  197 9.  Spirorbidae 
(Polychaeta  Sedentaria)  from  Alaska  to  Panama.  Jj_  Zool. 
Lond.v. 189:419-458. 

Mook,  D. ,  1976.  Studies  of  Fouling  Invertebrates  in  the 
Indian   River.    Bulletin   of  Marine   Science. v. 26 : 610-615 . 

Morris,R.H.,Abbott,D.P.,and    Hader lie, E.C.,  1980.    Intertidal 
Inver teb rates     o  f     California. Stanford     University 
Press. 690pp. 

Osburn,R.C,  1952. Allen   Hancock   Pacific  Expeditions     v. 14, 
parts   1,2, and   3    ,    Univ.    of   Southern   California   Press. 841pp. 

Osman,R.W.  ,1977.  The  Establishment  and  Development  of  a 
Marine  Epifaunal  Community.  EcoljO£_ica_l  Monog_r_aphs_.  v. 
47:pp   37-63. 


Ill 


Schoener,A.,Long,E.R.  ,Depalma,J.R.,1978.      Geographic 
GVvariation    in    Artificial    Island    Colonization    Curves 
Ecology. v. 47 , No. 2:367-382. 

Schoener,A.,  and  Gr eene , C.  H. ,  19 8 0 .  Variability  Among 
Identical  Fouling  Panels  in  Puget  Sound,  Washington, 
U.S.A.In  Proceedings  of  the  Fifth  Internationl  Congress  on 
Marine    Fouling    and    Cor rosion. Barcelona    , Spain. May    1979. 

Smi  th,R.  I. ,  and  Car  1  ton  ,  J.T. ,  19  7  5  .  Eds. Light'  s  Manual, 
Intertidal  Invertebrates  of  the  Central  California  Coast, 
3rd   ed. , University  of   California   Press.    716pp. 

Sutherland, J. S., 1974.  Multiple  Stable  Points  in  Natural 
communities.    American  Naturalist. v.    108:859-873. 

Whitaker,R.H.,1952.  A  Study  of  Summer  Foliage  Insect 
Communities  in  the  Great  Smoky  Mountains. Ecological 
Monographs. v. 22, No. 1 :pp   1-42. 

Woods  Hole, 1952.  Marine  Fouling  and  its  Prevention.  U.S. 
Naval    Institute,    Annapolis, Md. 


112 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Professor  Robert  J.  Renard,  Code  63Rd  1 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California   93940 

4.  Professor  Christopher  N.  K.  Mooers,  1 
Code  68Mr 

Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Professor  Eugene  C.  Haderlie,  Code  68Hc  1 
Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  California   93940 

6.  Professor  Patricia  A.  Jacobs,  Code  55Jc  1 
Department  of  Meteorology 

Naval  Postgraduate  School 

7.  Professor  Donald  P.  Gaver,  Code  55Gv  1 
Department  of  Operations  Research 

Naval  Postgraduate  School 
Monterey,  California   93940 

8.  LT.  David  L.  Martin  1 
Route  3,  Box  667 

Yakima,  Washington   25314 

9.  Director  1 
Naval  Oceanography  Division 

Naval  Observatory 

34th  and  Massachusettes  Avenue  NW 

Washington,  D.C.   20390 


113 


10.  Commander 

Naval  Oceanography  Command  Central 

NSTL  Station 

Bay  St.  Louis,  Missouri   39522 

11.  Commanding  Officer 

Naval  Oceanographic  Office 

NSTL  Station 

Bay  St.  Louis,  Missouri   39522 

12.  Commanding  Officer 

Fleet  Numerical  Oceanography  Center 
Monterey,  California  93940 

13.  Commanding  Officer 

Naval  Ocean  Research  and  Development 

Activity 
NSTL  Station 
Bay  St.  Louis,  Missouri   39522 

14.  Commanding  Officer 

Naval  Environmental  Prediction  Research 

Facility 
Monterey,  California   93940 

15.  Chairman,  Oceanography  Department 
U.S.  Naval  Academy 

Annapolis,  Maryland   21402 

16.  Chief  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  Virginia   22217 

17.  Office  of  Naval  Research  (Code  480) 
Naval  Ocean  Research  and  Development 

Activity 
NSTL  Station 
Bay   St.  Louis,  Missouri   39522 

18.  Commander 

Oceanographic  Systems  Pacific 

Box  1390 

Pearl  Harbor,  Hawaii   96860 

19.  Mrs.  Anne  Harrington 
Hopkins  Marine  Station 
Pacific  Grove,  California   93950 


114 


20.  Dr.  Amy  Schoener 
Department  of  Oceanography 
University  of  Washington 
Seattle,  Washington  98195 

21.  Mrs.  Jean  Montemarino,  Code  2844 
David  W.  Taylor  Naval  Ship  Research 

and  Development  Center 
Annapolis,  Maryland   21402 


115 


2015^6 
20161*6 


Martin 

The  use  of  computer 
intensive  statistical 
modeling  in  estimating 
the  variability  of 
marine  fouling  communi- 
ties. 


«  -  JUt  M 

20  FEB 


2970* 

3  1175 


Thesis 

M3573 

c.l 


20161*6 


Martin 

The  use  of  computer 
intensive  statistical 
modeling  in  estimating 
the  variability  of 
marine  fouling  communi- 
ties. 


££ksssS££Sr 


*™2768  002  12812  6 

3DU0UY  KNOX  UBBARY 


m  ffffitaLdii 


■asfiS 


mSSm 


JilHI 

■■.■;■.•.'■.■:■: '■ 


•'■'V. 


9nS 


•■'■',-.■:»■< 


mm 


wt 


!  ■■'■  '■;••>.■: