EPA-660/3-75-013
MAY 1975
Ecological Research Series
An Analysis of the Dynamics of DDT
in Marine Sediments
National Environmental Research Center
Office of Research and Development
U.S. Environmental Protection Agency
Corvallis, Oregon 97330
-------
RESEARCH REPORTING SERIES
Research reports of the Office of Research and Development,
U.S. Environmental Protection Agency, have been grouped into
five series. These five broad categories were established to
facilitate further development and application of environmental
technology. Elimination of traditional grouping was consciously
planned to foster technology transfer and a maximum interface in
related fields. The five series are:
1. Environmental Health Effects Research
2. Environmental Protection Technology
3. Ecological Research
4. Environmental Monitoring
5. Socioeconomic Environmental Studies
This report has been assigned to the ECOLOGICAL RESEARCH STUDIES
series. This series describes research on the effects of pollution
on humans, plant and animal species, and materials. Problems are
assessed for their long- and short-term influences. Investigations
include formation, transport, and pathway studies to determine the
fate of pollutants and their effects. This work provides the technical
basis for setting standards to minimize undesirable changes in living
organisms in the aquatic, terrestrial and atmospheric environments.
EPA REVIEW NOTICE
This report has been reviewed by the Office of Research and
Development, EPA, and approved for publication. Approval does
not signify that the contents necessarily reflect the views and
policies of the Environmental Protection Agency, nor does mention
of trade names or commercial products constitute endorsement or
recommendation for use.
-------
EPA-660/3-75-013
MAY L9'75
AN ANALYSIS OF THE DYNAMICS OF DDT IN MARINE SEDIMENTS
By
John H. Phillips
Hopkins Marine Station
Pacific Grove, California 93950
Eugene E. Haderlie
Naval Postgraduate School
Monterey, California 93940
Welton L. Lee
California Academy of Sciences
San Francisco, California 94118
Grant No. R800365
Program Element IBA025
21AIS, Task No. 08
Project Officer
Dr. Milton H. Feldman
Coastal Pollution Branch
Pacific Northwest Environmental Research Laboratoi
National Environmental Research Center
Corvallis, Oregon 97330
NATIONAL ENVIRONMENTAL RESEARCH CENTER
OFFICE OF RESEARCH AND DEVELOPMENT
U.S. ENVIRONMENTAL PROTECTION AGENCY
CORVALLIS, OREGON 97330
For Sale by the National Technical Information Service
U.S. Department of Commerce, Springfield, VA 22151
-------
ABSTRACT
The concentration of the three chlorinated hydrocarbons, DDT, DDD, and DDE, were
measured in sediments at 57 stations in Monterey Bay on the Central California coast
during 1970 and 1971. Mean concentration in parts per billion was DDT 3.1, DDD 2.3,
and DDE 5.4. Maximum concentrations were DDT 19.3, DDD 8.7, DDE, 20.5 parts per
billion. The distribution of the three compounds within South Monterey Bay was
charted. During 1973 nineteen of the original stations, representing locations that were
low, intermediate, and high concentrations in the original survey, were resampled. The
mean concentration approximately three years later were DDT 15.5, DDD 2.3, and DDE
5.4 parts per billion with maximum levels of DDT 83.1, DDD 11.4, and DDE 17.5 parts
per billion. A chart of the concentrations in South Monterey Bay revealed essentially
the same distribution of chlorinated hydrocarbons.
Two approaches to the estimation of annual system rates for input, I, output, O, decay,
D, and internal translocation, Tj and TQ, expressed as decimal fractions of the existing
concentration were developed, and Fortran programs that permit rapid estimations were
written. The mean annual system rates obtained were for DDT, 1+1.30, O-.059, D-.036,
TI and TO t -80 with a residence time of 11 years and life time of 29 years. An I of 1.30
means the amount of input is 130% of the existing concentration per year. The mean
annual rates obtained for DDD were, I + 0.25, O - 0.11, D - 0.025, Tj and TQ ± 0.20 with
residence time of 7 years and life time of 44 years. The rates for DDE were I + 0.28,
O - 0.10, D - 0.027, TO and Tj + 0.22 with residence time of 8 years and life time of 39
years. The approaches to these estimates are dependent upon variability in net rates of
change at the various stations and an approach to evaluation of the standard deviation
of the estimated rates relative to distributions of net rates with minimal variance is pre-
sented.
Laboratory assays were developed to determine the relative rate of decomposition in
sediment placed under conditions selective for various physiologically different kinds
of microorganisms. l^C rjng labelled substrates were used in all assays. Decay of the
three chlorinated hydrocarbons under aerobic conditions without additional nutrients
was greater than decay under anaerobic conditions. The addition of accessory energy
and carbon sources such as sodium acetate did not increase the rate of decay under
anaerobic conditions. There was some decay under anaerobic conditions suggesting
mechanisms of ring cleavage not involving incorporation or oxygen prior to ring split.
Nitrate as an accessory electron acceptor increased the rate of decomposition under
anaerobic conditions. Degradation products formed from the parent compounds in-
cluded water soluble intermediates as well as carbon dioxide.
The QIQ for the decay process as determined by laboratory assays incubated at 10° and
20° C. is 2. 5.
This report was submitted in fulfillment of Grant No. R 800365 by Hopkins Marine Station.
Work was completed under sponsorship of the Environmental Protection Agency as of 1974.
A.GEHOY
-------
CONTENTS
Page
Abstract ii
List of Figures iv
List of Tables v
Acknowledgements vii
Sections
I Conclusions 1
II Recommendations 2
III Introduction 3
IV Methods 9
V Results and Discussion 13
VI References 58
VII Appendices 59
iii
-------
FIGURES
No. Page
1 The study area, Monterey Bay. Sampling stations are indicated 5
by number.
2 DDT as a percent of the total concentration of DDT, ODD, and 19
DDE plotted for data obtained in 1970 and 1971. Circled numbers
indicate actual percents in excess of 50%.
3 ODD as a percent of the total concentration of DDT, DDD, and 20
DDE plotted for data obtained in 1970 and 1971. Circled numbers
indicate actual percents in excess of 50%.
4 DDE as a percent of the total concentration of DDT, DDD, and 21
DDE plotted for data obtained in 1970 and 1971. Circled numbers
indicate actual percents in excess of 50%.
5 Total concentration in parts per billion of DDT, DDD, and DDE 22
from data obtained in 1970 and 1971. Circled numbers indicate
actual concentrations in excess of 50 ppb.
6 Total concentration in parts per billion of DDT, DDD, and DDE 23
from data obtained in 1973. The blank portions of the area were
not sampled. Circled numbers indicate actual concentrations in
excess of 50 ppb.
7 Model of the system of sediment compartments and this system's 29
relation to other systems.
8 Composite chart of the translocation of DDT compounds based 49
upon the rates of change, K, at individual stations in the
southern portion of Monterey Bay.
IV
-------
TABLES
No. Page
1 Concentration of DDT, DDD, and DDE in sediments 6
of the Monterey area land drainage system in 1972
(State of Calif., 1974).
2 Typical decay assay protocol. 11
3 Concentrations of DDT, DDD, and DDE in marine sediment 14
samples from Monterey Bay.
4 Levels of DDT, DDD, and DDE as percent of total residues 16
in marine sediment samples from Monterey Bay.
5 Variance of sampling measured at Station 38. 18
6 First page of computer output showing concentration of 25
pollutant compounds in sediment from sample stations at
first sampling time. Cj identifies as concentrations at time one.
7 Second page of computer output showing concentration of 26
pollutant compounds in sediment from sample stations at the
second sampling time. G£ identifies as concentrations at time two.
8 Third page of computer output showing percent of total of each 27
of the three compounds in sediments from sample stations at the
first sampling time. C^ identifies as data for time one.
9 Fourth page of computer output showing percent of total of each 28
of the three compounds in sediments from sample stations at the
second sampling time. €2 identifies as data for time two.
10 Fifth page of computer output showing the rate of change, K, 31
for DDT in each sediment compartment.
11 Sixth page of computer output showing the rate of change, K, 32
for DDD in each sediment compartment.
12 Seventh page of computer output showing the rate of change, K, 33
for DDE in each sediment compartment.
13 Eighth page of computer output showing a summary of the 3 5
annual system rates expressed as decimal fractions of the mean
concentration of DDT present in the system.
-------
No. Page
14 Ninth page of computer output showing a summary of the 36
annual system rates expressed as decimal fractions of the mean
concentration of ODD present in the system.
1 5 Tenth page of computer output showing a summary of the 37
annual system rates expressed as decimal fractions of the mean
concentration of DDE present in the system.
16 Comparison of estimates obtained from the 16 and 19 station 38
data sets and using actual paired sample analyses. Standard
deviations and coefficients of variation are included.
17 Comparison of estimates obtained from 49 and 57 station 40
data sets and using sample analyses paired with mean concentration
levels. Standard deviations and coefficients of variation are included.
18 Standard deviations and standard errors of distributions of K 43
with minimal variance for given values of I, Tj and TQ, (O+D),
and n.
19 Comparison of uncorrected and corrected standard deviations of 44
system estimates.
20 Mean of the estimates for the South Monterey Bay system and 47
associated descriptive statistics.
21 Total amounts of DDT, DDD, and DDE in the South Monterey 48
Bay study area based on the mean concentrations at the two
sample times, and expected amounts affected by the mean of the
estimates of system rates.
22 Results of a laboratory assay of annual rate of decay of DDT to CO2, 5 1
^COo' exPressed as a decimal fraction of the initial concentration
of DDT maintained at 10° C under aerobic conditions.
23 Results of laboratory assays of the annual rates of decay to CO2, 52
> ar>d the effect of environmental variables on the process.
24 Rates of decay to water soluble compounds and CC>2 determined 55
by laboratory assays.
VI
-------
ACKNOWLEDGEMENTS
The financial support of the U.S. Environmental Protection Agency and interest of
Dr. Milton H. Feldman are gratefully acknowledged.
A gift from the Forest Park Foundation permitted the acquisition of equipment es-
sential to this project and other studies on chlorinated hydrocarbons in our laboratories.
This support and interest in environmental problems has been greatly appreciated.
Finally, the authors would like to acknowledge Philip Murphy, Will McCarthy, Barbara
Cunningham, Anne Edwards, and Charles Bates for their technical assistance in this
project.
Vll
-------
SECTION I
CONCLUSIONS
Chlorinated hydrocarbons associated with sediment particles tend to concentrate in
sedimentation basins which may be at some distance from the input source.
Although the use of chlorinated hydrocarbon pesticides has declined sharply the levels
of three materials has continued to increase in marine sediments. The principal source
of this additional pollutant load in this instance appears to be more related to translo-
cation of these materials absorbed to sediments of adjacent land drainage systems.
The dynamics of chlorinated hydrocarbons in the coastal marine environment, although
complex, are susceptible to study. Approaches to the estimation of rates of input, decay,
and translocation can be developed and assessed by continued analysis of environmental
samples.
The measurement of decay rate by laboratory assay appears to have its greatest utility
in the determination of the effect of environmental conditions on the process of decay.
Duplication of conditions existing in situ in the laboratory can only be approximated
and then only for a limited time. The laboratory work, short term in its execution, serves
only as a guide to what is happening in the environment.
-------
SECTION II
RECOMMENDATIONS
The complexities of the dynamics of coastal pollution by chlorinated hydrocarbons
necessitates an initial survey of the concentration of these environmental contaminants at
a large number of stations. Once basins of accumulation are established and principal
translocation paths established a much smaller number of stations require surveillance at
later points in time. It doesn't appear to be essential to monitor exactly the same stations
in any surveillance program as long as the set of surveillance stations includes established
basins and positions along translocation pathways.
It is recommended that initial intensive surveys be carried out in the coastal marine envi-
ronment adjacent to major agricultural and industrial areas which are known to produce
or utilize poorly degraded environmental contaminants such as the chlorinated hydro-
carbons.
Monterey Bay is a very useful model coastal marine environment for the establishment
and testing of approaches to system rate estimation. Continued surveillance of this area is
recommended.
It is also recommended that work be done on extending the approach to estimation of
system rates explored with respect to sediments to other environmental systems including
populations of organisms. It would appear desirable to concentrate initially upon abun-
dant and useful indicator organisms rather than commercially desirable or affected species.
Finally, it is recommended that additional effort be expended on the study of laboratory
assays of decay not only as approximations of the environment but as useful preparations
for elucidating the conditions inhibitory and stimulatory to the decay process.
-------
SECTION III
INTRODUCTION
Although the accumulation of chlorinated hydrocarbons in the marine ecosystem has
been a matter of concern for some time, methods for assessing the rates of accumula-
tion, decay, and translocation have been lacking. The problem is not unique to the ma-
rine environment, and methods for assessment of the dynamics of chemical pollutants
in general are needed for meaningful analysis of the residue measurements tabulated in
most investigations. Without an assessment of rates such tabulations generally permit
only the detection of some general trend of increase or decrease in concentration during
the period of study. In many cases, however, the amount of variability is so great that
the number of samples required to show such general trends is prohibitive. Yet we have
both the data available and a need to use these data for meaningful assessment. In addi-
tion, before any feasible monitoring activity geared to control and regulatory strategies
are designed and implemented, a means of assessing any new tabulations is required as
a determinant in the design of such activities. Whatever systems of assessment may be
developed in the future it cannot be expected that they will overcome the variability
that plagues environmental sampling. Rather, such systems should be expected to pro-
vide an estimate of this variability and a confidence interval for any derived parameter
of environmental change.
Several models stressing one or another aspect of the dynamics of pesticides in the en-
vironment have been presented (Hamaker 1966, Robinson 1967, Woodwell 1967, Har-
rison et al. 1970, and Eberhardt et al. 1971), but there still appears to be a need for a
general approach that provides a means of estimating rates of input, decay, and trans-
location from some minimal number of analyses. The study presented here is an attempt
to fill this need.
The data used here for these estimations consists of analyses of marine sediment samples
for l,l,2-trichloro-2,2-bis (p-chlorophenyl) ethane, DDT; l,l-dichloro-2,2-bis (p-chloro-
phenyl) ethane, DDD; and l,l-dichloro-2,2-bis (p-chlorophenyl) ethylene, DDE. The
rates of decay at a sampling site and translocation away from a sampling site are difficult
to separate through the approach to estimation presented. Laboratory measurements of
the rate of C ring labelled DDT in marine sediments held under a variety of conditions
are also presented. These measurements reflect decay to the point of CO 2 release rather
than conversion to any one of a variety of other metabolites including DDD and DDE,
but are useful in assessing the method of estimation based upon environmental samples
alone.
The analysis of DDT residue levels in marine sediments reported herein is only a part of
a larger study correlating the levels of pollutants with density and composition of benthic
populations. Other results of this study will be reported elsewhere.
3
-------
THE STUDY AREA
This study was carried out in Monterey Bay located in the central coastal region of
California. Figure 1 shows the study area and the location of the forty-nine Stations
from which sediment samples were obtained. The figure also shows several geographical
features pertinent to this investigation. The bottom of Monterey Bay is divided by a
major submarine canyon over 3800 meters in depth at its deepest point. The sampling
effort was concentrated in the southern portion of the bay with no sampling beyond
the 200 fathom, 365 meter, line. Residue levels of DDT, ODD, and DDE were first
measured in samples from this southern portion of the bay during 1970 and nineteen
of these stations were resampled in 1973. A small number of stations were sampled in
the northern part of the bay during 1971.
Monterey Bay is the recipient of drainage from a major agricultural area, the Salinas
Valley, where DDT was used in large amounts for a period of twenty years. Usage of
this pesticide and DDD has decreased sharply since 1969. A tabulation of use was
started in 1970 when 33,931 pounds was applied to 19,387 acres in Monterey County.
This input level was further reduced in 1971 to 4,697 pounds, and in 1972 to
10 pounds on 20 acres (Calif. Dept. of Agriculture 1970, 1971, 1972). Final tabulations
for 1973 will probably show levels of input similar to those of 1972. Although the use
of DDT in the area adjacent to Monterey Bay has declined sharply since 1970, the level
of DDT in marine sediments appears to be increasing as more of this pesticide finds its
way to the sea via the drainage system of the neighboring agricultural area. The decrease
in usage on adjacent land and apparent increase in concentration in the marine sediments
of the area suggests that continued study of the Monterey area is of particular interest
in determining the time lag between terrestrial input and marine accumulation of persis-
tent chemical pollutants.
Although in the past, when DDT was being regularly applied on the adjacent lands, the
atmosphere was an important source of input to the bay; at the present time the major
source of input appears to be the Salinas River which drains the inland agricultural areas.
This river flows directly into the bay only intermittently. Most of the time the mouth
of the river is blocked by a bar of sand that is removed only at times of heavy rainfall
to prevent flooding. During this investigation this event occurred Jan. 13, 1970, Nov. 30,
1970, Dec. 29, 1971, Nov. 16, 1972, Nov. 17, 1972, and Nov. 20, 1973. Input directly
by the river has, therefore, not been continuous.
Analyses of the sediment samples from the river bed along its course in 1972 (State of
California, 1974) showed considerable variation in the relative abundance and concen-
tration of the three compounds. Table 1 gives the results of these analyses and the ap-
proximate location of the samples relative to the mouth of the river.
During the periods when the mouth of the river is blocked, there is a sluggish flow north
to Elkhorn Slough which served as the mouth of the river until 1908. This flow is joined
by drainage from Trembladero Slough which receives water and sediments from the Re-
clamation Canal that flows through the City of Salinas to the east and beyond the right-
hand margin of the figures. The Reclamation Canal receives effluents from food proces-
sing plants and other industries, and analyses of its sediment in 1972 (State of Calif.,
1974) revealed the levels also listed in Table 1.
-------
SANTA
CRUZ
TREMBLADERO
SLOUGH
Figure 1. The study area, Monterey Bay. Sampling stations are indicated by number.
5
-------
Table 1. CONCENTRATION OF DDT, ODD, AND DDE IN SEDIMENTS OF THE
MONTEREY AREA LAND DRAINAGE SYSTEM IN 1972 (STATE OF
CALIF., 1974)
Salinas River
distance from mouth
(kilometers)
42
25
8
3
Reclamation Canal
distance from mouth
of Elkhorn Slough
(kilometers)
20
(ppb)
DDT
1.0
120.
150.
16.
0.12
7,000.
21,000.
ODD
1.3
1000.
620.
30.
45,000.
150,000.
DDE
- 20.
360.
10,000.
-------
RATIONALE OF DESCRIBED WORK
Selection of Study Site and Source of Marine Sediments for Decay Assays—For the
estimation of rates governing the dynamics of a chlorinated hydrocarbon pollutant
in marine sediments an area with the following characteristics appeared most desirable.
(1) The marine area should be adjacent to a land area for which there exists an account-
ing of input to the environment through normal use. The use of DDT and DDD within
the State of California has been subject to such accounting on a square mile section
basis since 1970 (Calif. Dept. of Agriculture 1970). Such accounting is available only
for normal agricultural and related uses. Therefore, areas which receive or have received
less well determined inputs from chlorinated hydrocarbon manufacture, such as the ocean
adjacent to Los Angeles, are less desirable for this type of study. (2) In order to assess
translocation within the study area it would appear desirable to select a marine area
with a limited number of point sources of input rather than one subject to diffuse in-
put by way of the atmosphere. (3) The area should be one open to general oceanic in-
fluence rather than a closed system so that translocation of the pollutant out of the
system by dilution or dissemination can be assessed. (4) As a source of materials for
laboratory assays of decay the area should be one which has had a long exposure to
the pollutant, thus insuring the establishment of microbial systems with the capacity
for decomposition of the pollutant. (5) The area should be known to be contaminated
with the pollutant. (6) The area should be accessible to sampling and close to the re-
quired analytical capability.
Monterey Bay, and in particular the southern portion of Monterey Bay, has these char-
acteristics and was selected as the study site and source of materials for the development
of laboratory assays for the rate of decay of DDT, DDD, and DDE.
Survey of Residue Levels in Monterey Bay Sediments—In order to assess the variability
in concentration and distribution of the three compounds in the sediments of Monterey
Bay thirty-seven sample sites were selected for analysis in the southern portion of the
bay which receives water and sediments from the agricultural area of Monterey County
by way of the Salinas River. An additional eleven sample sites in the northern portion
of the bay were selected in order to assess any augmenting effect of additional river
input sources such as the San Lorenzo and Pajaro Rivers that drain areas of Santa Cruz
and San Benitio Counties lying adjacent to Monterey County and Monterey Bay.
Determination of the Amount of Change in Residue Levels with Time—In order to as-
sess the magnitude of change in the concentration of DDT and related compounds a
subset of the original survey sampling stations was resampled and analyzed after ap-
proximately three years. Nineteen of the original sample stations were selected as this
subset. The selection was made on a basis of accessibility and representations of stations
showing a broad range of residue concentrations as determined in the original survey.
Determination of the Variance of Sampling—One additional sample station, number 38,
which had never before been sampled was added to the resampled subset and sampled
three times on the same day. Three aliquots from each of these samples were analyzed
for the three compounds to provide an estimate of the variability of sampling.
-------
Approaches to the estimation of rates and Dynamics of the compounds in Sediments-
Using the tabulated data obtained from the sampling programs various approaches to
the estimation of the rates of input, translocation, and decay were developed for the
system of sample sites. Considerable attention was directed to estimation of variance
of these derived rates.
Development of Laboratory Assay Methods for the Determination of Decay Rate-
Measurement of decay rate based on changes in residue level observed by repeated
sampling from the environment are subject to error due to translocation to or away
from the sample site. Therefore, a means of estimating decay rate in a closed system
not susceptible to such error would be desirable. A variety of preparations using ^C
ring labelled compounds were established for such estimations.
Effect of Environmental Variables on Decay Rate—Any closed system preparation is
by its very nature selective for one or another metabolic type of microorganism. The
initial conditions and conditions which subsequently develop may have a marked ef-
fect upon the observed rate of decomposition through the election of particular micro-
bial populations. Therefore, it was necessary to study the process of decay as influenced
by a number of environmental variables chosen to encourage one or another of the ma-
jor metabolic types of microorganisms.
-------
SECTION IV
METHODS
ANALYSIS OF SEDIMENT SAMPLES
Samples of sediment were collected by Shipek grab or shallow dredge. Between 50
and 70 grams of wet sediment were placed in a 250 ml bottle and mixed with 30-50
grams of granular anhydrous sodium sulfate. The sediment was extracted with 50 ml
of acetone:hexane, 1:1, by shaking for four hours. The acetone, hexane was decanted
and filtered through a fritted glass filter or silicon-treated phase separation paper into
a separatory funnel. Three additional 50 ml portions of hexane were used to wash the
sediment and added to the original extractant.
The extract was washed with three 200 ml portions of water followed by dehydration
of the extract by passage through a 2x5 cm column of anhydrous sodium sulfate and
concentration in a Kuderna-Danish concentrator to less than 10 ml. The extract was
then cleaned by shaking first with 1 ml of concentrated sulfuric acid and finally with
approximately 0.1 ml of mercury. The analysis was performed in a Beckman GC-4 Gas
Chromatograph with electron capture detector, using a mixed bed column of Chromo-
sorb W, 80-100 mesh, DMCS treated, and acid washed, containing 5% DC-200 and 5%
QF1.
Although the efficiency of extraction is difficult to assess, the effect of concentration
and clean-up procedures can be measured by the use of ^C labelled materials added
just prior to extraction with acetone, hexane. Recovery was 73.9% for DDT, 94.4%
for DDD, and 84.8% for DDE, and these figures were used to correct the results of
analyses.
LABORATORY DECAY ASSAYS
A variety of preparations have been investigated for their applicability to decay assay
preparations. These preparations have included sealed stationary aliquots of sediment
and l^c labelled substrate as well as ones in which the sediment with labelled substrate
was subjected to continuous percolation or periodic gas flow. Maintenance of percolat-
ing systems for the length of time required to measure the very slow rates of decay is
not feasible, and it is difficult to maintain a large number of preparations under condi-
tions whereby they may be subjected to periodic gas flow and trapping of metabolic CO2-
Therefore, sealed stationary preparations have proved to be the only feasible type of
preparation so far developed. The most convenient container for such preparations has
been 125 ml Hypovials, Pierce, Rockford, Illinois, No. 12995, fitted with Teflon liners.
The preparation of decay assays is as follows. Sediment is collected as for samples for
residue analysis, packed in ice, and brought to the laboratory within a few hours. The
sediment is rinsed through screen with 16 mesh to the inch to remove macroscopic in-
fauna and refrigerated. Aliquots of the slurried sediment are removed for dry weight
-------
determination. A volume of the slurried sediment equivalent to 24 grams dry weight
is delivered to a sterile Hypovial and seawater, with or without additional nutrients,
is added to give a volume of 98 ml total. One ml each of ^C and *4C substrate ad-
sorbed to sterile sediment is added giving a final volume of 100 ml. The preparation
may be gassed with nitrogen to produce an anaerobic environment prior to sealing.
All incubators are in the dark for periods of generally twelve weeks. All preparations
are set up in quintuplicate. A typical protocol is presented in Table 2.
i >\
C Substrate Preparation—2.4 grams of either, l,l-bis-(p-chlorophenyl)-2,2,2-tri-
chloro ethane, p-p'DDT 99+% No. 10, 002-1; 2,2-bis-(p-chlorophenyl)-l,l-dichloro
ethylene, No. 12, 289-7 (B 3964); or 2,2-bis-(p-chlorophenyl)-l,l-dichloroethane,
puriss B 3959 Aldrich Chemical Co. Inc., Milwaukee, Wisconsin, were dissolved in
10 ml of acetone. To 10 grams of dried sterile sediment 1 ml of acetone solution
was added and the sediment wet with an additional 3 ml of acetone. The acetone
was evaporated off at room temperature and 96 ml of distilled water added to slurry
the sediment and its adsorbed substrate. One ml contains 2.4 x 10* ug of substrate
on 0.1 gram of sediment per ml. Similar preparations were made giving 2.4 x 10 ug
and 21.6 ug of substrate on 0.1 gram of sediment per ml.
14C-DDT Substrate Preparation-Uniformly ring labelled DDT, Amersham/Searle Corp.,
63.9 u Ci/mg in benzene was used for preparation of the substrate. The original 250 u Ci
preparation was diluted with acetone and 240 ug in 4 ml was added to 10 grams of dried
sterile sediment. The acetone was removed by evaporation at room temperature and 96
ml of distilled water added to give 2.4 ug *4C-DDT and 0.1 gram of sediment per ml.
A similar preparation was made giving 0.24 ug ^C-DDT and 0.1 gram of sediment per ml.
14C-DDD Substrate Preparation-14C-DDT was converted to 14C-DDD by the method
of Murphy (1970) and purity of the product confirmed by gas chromatography. The
resulting material was used to prepare substrate as described above for *4C-DDT.
14C-DDE Substrate Preparation-14C-DDT was converted to 14C-DDE by the method
of Gunther and Blinn (1950) and purity of the product confirmed by gas chromatography.
The resulting material was used to prepare substrate as described above for *4C-DDT.
Analysis of Decay Assays—After incubation for generally 12 weeks CO2 was trapped
by the addition of 1.5 ml of 5 N NaOH to the Hypovial. The base was introduced by
syringe and the ampoule resealed with tape. Syringe delivered 5 ml aliquots of the
basic slurried sediment were transfered to 25 ml Erlenmeyer flasks containing magnetic
stirring bars. The flasks were stoppered with Top stoppers, K-882310, fitted with plastic
center wells, K-882320, both from Kontes Glass Co., Vineland, N.J. The center wells
contained an accordian pleated Whatman No. 1 filter paper wick, 2.5x5 cm. (3-phenyl-
ethylamine, 0.15 ml, was delivered to the well and wick by syringe through the stopper.
While the sediment in the flask was gently stirred on a magnetic stirrer 0.25 ml of 5 N
F^SO^. was added to the sediment. The flasks were then held for 24 hours at room
temperature after which time the wicks were removed to scintillation vials to which
was added 15 ml of Toluene-omnifluror. Appropriate preparations for background
10
-------
Table 2. TYPICAL DECAY ASSAY PROTOCOL.
Hypovial
No.
1-5
6-10
11-15
16-20
21-25
Slurried
(grams)
24
24
24
24
24
Sediment
(ml)
59
59
59
59
59
Seawater
plus
nutrients
(ml)
39
39
39
39
39
12C
Substrate
(ug) (ml)
2400
240
21.6
0
0
1
1
1
0
0
14c
Substrate
(ug) (ml)
2.4
2.4
2.4
2.4
0.24
1
1
1
1
1
Total
Substrate
(ppm)
100
10
1
0.1
0.01
Total
volume
(ml)
100
100
100
100
100
11
-------
measurement were also made. The amount of CC>2 was determined in a Nuclear
Chicago Corp. Unilux II. Diffusion time and trapping volume of 0-phenylethylamine
were established through tests using a standard preparation of Na ^COj.
DECAY AS AFFECTED BY ENVIRONMENTAL VARIABLES
The effect of temperature was determined by comparing the amount of decomposition
at 10° and 20°C, and the effect of oxygen, nitrate, and sulfate as terminal electron ac-
ceptors in the presence and absence of cometabolizable sodium acetate and ethanol was
determined by appropriate additions to the Hypovials.
12
-------
SECTION V
RESULTS AND DISCUSSION
SURVEYS OF RESIDUE LEVELS IN MONTEREY BAY SEDIMENTS
The concentration in parts per billion of the three compounds, DDT, ODD, and DDE
in sediment samples collected during the three sampling periods are presented in Table 3.
Table 4 presents the same set of analyses in terms of the percent of total residues for
each of the three compounds.
The variance of sampling at Station 38 can be assessed from the data presented in
Table 5. The greatest variation in results can be observed with respect to DDT, the
compound also showing the greatest loss during the extraction, concentration, and
cleanup procedures as mentioned in the section on methods.
The data obtained in the 1970 and 1971 samplings is presented in Figures 2, 3, and 4,
where the distribution of DDT and its two derivatives is displayed in terms of percent
of the concentration of total DDT derivatives. Figures 5 and 6 show the distribution in
terms of the total concentration of DDT and its two derivatives in parts per billion.
Figure 5 shows the distribution in 1970 and 1971, and Figure 6 shows the distribution
as indicated by the analyses of the smaller number of samples obtained in 1973.
The small number of sample stations in the northern portion of the bay did not reveal
any unusual augmentation in concentrations of the three compounds due to input from the
San Lorenzo and Pajaro Rivers although the percent composition of DDT derivatives
does indicate differences between the northern and southern portions of the bay.
If particular attention is paid to the southern portion of the bay for which there is the
greatest information, the distributions suggest a number of characteristics of the system.
After input with sediments from the Salinas River, and perhaps also through Elkhorn
Slough, these materials are subjected to considerable translocation due to the currents
operating within the south bay. The highest concentration of DDT derivatives is to be
found at a considerable distance from the mouth of the river. Close to the mouth of
the river, however, the sediments show a high percentage of DDT which is characteristic
of some of the sediments within the drainage system. These high DDT percentages are
also found at the more distant points where the highest concentrations of derivatives
are found as well. Over much of the area in terms of percent, however, DDE represents
the major compound.
These plots of distribution reflect input over a considerable period of time. During this
time the major routes of input may have changed considerably as has the relative con-
centrations of the three derivatives in these input sources. Nevertheless, the apparent
constancy of location of major basins of deposition is remarkable. Areas with high
concentrations in 1970 have become even more heavily contaminated in 1973.
13
-------
Table 3. CONCENTRATIONS OF DDT, ODD, AND DDE IN MARINE SEDIMENT
SAMPLES FROM MONTEREY BAY.
Station
1
2
3
4
5
6
7
8
g
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
LOCATION
Latitude Longitude
3647.2512148.90
3646.85 121 53.50
3646.35 121 49.00
3646.05121 51.00
3646.00121 57.00
3645.45121 50.00
3645.30121 54.00
3645.20121 54.00
3645.10121 52.00
3645.10121 50.00
3645.0012149.00
3644.60121 50.50
3644.25121 50.35
3644.20121 52.25
3644.00121 50.00
3644.00121 49.50
3643.75 121 54.45
3643.50121 51.80
3643.35 121 56.25
3643.18121 57.00
3643.00121 51.00
3642.90121 58.00
3642.55 121 53.30
3642.50121 50.30
3641.70121 55.00
3641.55 121 55.50
3641.50121 52.00
3641.00121 51.00
3640.90121 56.40
3640.50121 53.50
3640.08121 54.05
3639.80121 54.50
3639.80121 5T.50
3639.10121 53.08
3639.10121 53.08
3637.95 121 52.50
3637.77 121 51.83
Date
8-23-70
11-15-70
2-20-70
11-15-70
5-29-70
11-15-70
5-29-70
5-29-70
5-29-70
5-29-70
2-20-70
2-20-70
11-15-70
8-23-70
5-29-70
2-20-70
1M5-70
2-20-70
8-23-70
2- 8-70
5-29-70
2-20-70
8-23-70
8-23-70
2-20-70
2- 8-70
5-29-70
11-15-70
2-20-70
5-29-70
2- 8-70
5-29-70
2- 9-70
2- 8-70
2- 8-70
2-20-70
2- 8-70
DDT
(ppb)
8.36
1.63
5.71
4.28
2.14
2.04
0.0
2.65
4.48
6.42
3.67
0.52
0.26
5.20
0.0
0.69
1.02
1.73
1.12
0.0
6.12
0.0
13.20
19.30
1.22
0.0
2.85
0.0
1.32
2.55
0.0
2.04
0.0
2.44
8.67
2.65
0.49
ODD
(ppb)
3.67
6.76
0.71
6.61
0.93
1.17
2.50
4.26
5.14
8.67
0.40
0.18
0.19
7.50
0.19
0.14
0.38
2.64
0.25
5.00
1.30
0.35
5.73
0.65
0.53
2.35
2.50
1.61
1.61
1.76
0.82
1.91
1.42
0.66
0.66
2.79
0.21
DDE
(ppb}
5.76
14.70
1.02
10.70
4.00
1.80
4.51
6.51
4.51
7.01
0.45
0.28
0.45
15.50
0.35
2.75
0.70
2.40
0.65
20.50
6.01
1.92
13.00
2.75
2.40
7.01
8.01
4.26
9.02
6.76
3.25
5.26
8.52
2.40
3.00
10.00
0.50
TOTAL
17.79
23.09
7.44
21.59
7.07
5.01
7.01
13.42
14.13
22.10
4.52
0.98
0.90
28.20
0.54
3.58
2.10
6.77
2.02
25.50
13.43
2.27
31.93
22.70
4.15
9.36
13.36
5.87
11.95
11.07
4.07
9.21
9.94
5.50
12.33
15.44
1.20
14
-------
Table 3. (continued) CONCENTRATIONS OF DDT, ODD, AND DDE IN MARINE
SEDIMENT SAMPLES FROM MONTEREY BAY.
Station
39
40
41
42
43
44
45
46
47
48
49
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
38
LOCATION
Latitude Longitude
3654.8012201.00
3657.10121 56.20
3656.70121 59.20
3655.50121 52.60
3655.10121 56.70
3653.60121 57.50
3653.00121 55.00
3652.30121 59.80
3651.00121 49.80
3650.80121 53.60
3650.20121 50.20
3647.25 121 48.90
3646.85 121 53.50
3646.35 121 49.00
3646.05 121 51.00
3645.10121 50.00
3645.00121 49.00
3644.20121 52.25
3644.00121 49.50
3643.75 121 54.45
3643.35 121 56.25
3643.18 121 57.00
3642.90121 58.00
3642.55 121 53.30
3641.70121 55.00
3641.55121 55.50
3640.90121 56.40
3639.10121 53.08
3637.95 121 52.50
3637.77 121 51.83
3638.47 121 51.68
Date
11-24-71
11-10-71
11-24-71
11-10-71
11-10-71
11-24-71
11-10-71
11-24-71
11-10-71
11-24-71
11-10-71
7- 9-73
7- 9-73
7- 9-73
7- 9-73
7- 2-73
7- 2-73
6-21-73
7- 2-73
8- 9-73
8- 9-73
6-21-73
8- 9-73
6-21-73
7-16-73
7-16-73
7-16-73
8- 9-73
6-21-73
7-16-73
9-21-73
DDT
(ppb)
0.60
1.62
0.93
0.85
0.81
1.13
1.21
1.27
1.16
1.62
0.78
1.06
9.50
1.10
3.63
0.92
2.18
30.60
0.96
5.41
72.70
63.10
0.93
29.90
1.14
0.68
0.70
1.18
83.10
0.54
0.62
ODD
(ppb)
1.90
8.15
2.75
1.58
3.07
2.54
2.01
3.81
1.27
5.61
1.48
0.53
11.40
0.53
5.43
0.39
0.72
6.07
0.06
4.54
3.19
0.79
0.90
4.32
2.74
2.20
1.11
0.42
0.95
0.20
0.38
DDE
(ppb)
2.00
5.54
4.48
0.66
2.59
2.47
1.88
5.06
1.13
6.72
1.29
0.56
17.50
0.63
6.91
0.52
0.83
11.20
0.23
17.30
12.00
3.48
6.06
12.20
10.49
8.67
5.67
2.44
3.34
0.40
2.72
TOTAL
4.50
15.31
8.16
3.09
6.47
6.14
5.10
10.14
3.56
13.95
3.55
2.15
38.40
2.26
15.97
1.83
3.73
47.87
1.25
27.25
87.89
67.37
7.89
46.42
14.37
11.55
7.48
4.04
87.39
1.14
3.72
15
-------
Table 4. LEVELS OF DDT, ODD, AND DDE AS PERCENT OF TOTAL RESIDUES IN
MARINE SEDIMENT SAMPLES FROM MONTEREY BAY.
Station
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
LOCATION
Latitude Longitude
3647.25 121 48.90
3646.85 121 53.50
3646.35 121 49.00
3646.05 121 51.00
3646.00121 57.00
3645.45 121 50.00
3645.30121 54.00
3645.20121 54.00
36 45. TO 121 52.00
3645.20121 50.00
3645.00121 49.00
3644.60121 50.50
3644.25 121 50.35
3644.20121 52.25
3644.00121 50.00
3644.00121 49.50
3643.75 121 54.45
3643.50121 51.80
3643.35 121 56.25
3643.18 121 57.00
3643.00121 51.00
3642.90121 58.00
3642.55 121 53.30
3642.50121 50.30
3641.70121 55.00
3641.55 121 55.50
3641.50121 52.00
3641.00121 51.00
3640.90121 56.40
3640.50121 53.50
3640.08 121 54.05
3639.80121 54.50
3639.80121 51.50
3639.10121 53.08
3639.10121 53.08
3637.95 121 52.50
3637.77121 51.83
Date
8-23-70
11-15-70
2-20-70
11-15-70
5-29-70
11-15-70
5-29-70
5-29-70
5-29-70
5-29-70
2-20-70
2-20-70
11-15-70
8-23-70
5-29-70
2-20-70
11-15-70
2-20-70
8-23-70
2- 8-70
5-29-70
2-20-70
8-23-70
8-23-70
2-20-70
2- 8-70
5-29-70
11-15-70
2-20-70
5-29-70
2- 8-70
5-29-70
2- 9-70
2- 8-70
2- 8-70
2-20-70
2- 8-70
DDT
(%)
46.99
7.06
76.75
19.82
30.27
40.72
0.0
19.75
31.71
29.05
81.19
53.06
28.89
18.44
0.0
19.27
48.57
25.55
55.45
0.0
45.57
0.0
41.34
85.02
29.40
0.0
21.33
0.0
11.05
23.04
0.0
22.15
0.0
44.36
70.32
17.16
40.83
ODD
(%)
20.63
29.28
9.54
30.62
13.15
23.35
35.66
31.74
36.38
39.23
8.85
18.37
21.11
26.60
35.19
3.91
18.10
39.00
12.38
19.61
9.68
15.42
17.95
2.86
12.77
25.11
18.71
27.43
13.47
15.90
20.15
20.74
14.29
12.00
5.35
18.07
17.50
DDE
(%)
32.38
63.66
13.71
49.56
56.58
35.93
64.34
48.51 .
31.92
31.72
9.96
28.57
50.00
54.96
64.81
76.82
33.33
35.45
32.18
80.39
44.75
84.58
40.71
12.11
57.83
• 74.89
59.96
72.57
75.48
61.07
79.85
57.11
85.71
43.64
24.33
64.77
41.67
16
-------
Table 4. (continued)
LEVELS OF DDT, ODD, AND DDE AS PERCENT OF TOTAL
RESIDUES IN MARINE SEDIMENT SAMPLES FROM MONTEREY
BAY.
Station
39
40
41
42
43
44
45
46
47
48
49
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
38
LOCATION
Latitude Longitude
3654.8012201.00
3657.10121 56.20
3656.70121 59.20
3655.50121 52.60
3655.10121 56.70
3653.60121 57.50
3653.00121 55.00
3652.30121 59.80
3651.0012149.80
3650.80121 53.60
3650.20121 50.20
3647.2512148.90
3646.85 121 53.50
3646.35 121 49.00
3646.05 121 51.00
3645.10121 50.00
3645.00121 49.00
3644.20121 52.25
3644.00121 49.50
3643.75 121 54.45
3643.35121 56.25
3643.18 121 57.00
3642.90121 58.00
3642.55 121 53.30
3641.70121 55.00
3641.55 121 55.50
3640.90121 56.40
3639.10121 53.08
3637.95 121 52.50
3637.77 121 51.83
3638.47 121 51.68
Date
11-24-71
DDT
(%)
13.33 .
11-10-71 : 10.58
11-24-71
11-10-71
11-10-71
11-24-71
11-10-71
11-24-71
11-10-71
11-24-71
11-10-71
7- 9-73
7- 9-73
7- 9-73
7- 9-73
7- 2-73
7- 2-73
6-21-73
7- 2-73
8- 9-73
8- 9-73
6-21-73
8- 9-73
6-21-73
7-16-73
7-16-73
7-16-73
8- 9-73
6-21-73
7-16-73
9-21-73
11.40
27.51
12.52
18.40
23.73
12.52
32.58
11.61
21.97
49.30
24.74
48.67
22.73
50.27
58.45
63.92
76.80
19.85
82.72
93.66
11.79
64.41
7.93
5.89
9.36
29.21
95.09
47.37
16.67
ODD
(%)
42.22
53.23
33.70
51.13
47.45
41.37
39.41
37.57
35.67
40.22
41.69
24.65
29.69
23.45
34.00
21.31
19.30
12.68
4.80
16.66
3.63
1.17
11.41
9.31
19.07
19.05
14.84
10.40
1.09
17.54
10.22
DDE
(%)
44.44
36.19
54.90
21.36
40.03
40.23
36.86
49.90
31.74
48.17
36.34
26.05
45.57
27.88
43.27
28.42
22.25
23.40
18.40
63.49
13.65
5.17
76.81
26.28
73.00
75.06
75.80
60.40
3.82
35.09
73.12
17
-------
Table 5. VARIANCE OF SAMPLING MEASURED AT STATION 38.
Sample
1
2
3
Subsample
1
2
3
1
2
3
1
2
3
DDT
(ppb)
.687
.772
.550
.561
.706
.801
.663
.398
.405
ODD
(ppb)
.430
.470
.370
.345
.333
.280
.439
.315
.418
DDE
(PPb)
3.01
2.90
2.85
2.89
2.38
2.57
2.63
2.96
2.32
TOTAL
(PPb)
4.13
4.14
3.77
3.80
3.42
3.65
3.73
3.67
3.14
Mean .6159 .3778 2.7233 3.7167
Variance .02167 .00416 .06574 .09841
Standard Deviation + .1472 + .0645 + .2564 + .3137
Standard Error + .0491 + .0215 + .0855 + .1046
95% Confidence Limits t.1131 +.0495 + .1971 + .2411
18
-------
%DDT 1970 8 1971
Figure 2. DDT as a percent of the total concentration of DDT, ODD, and DDE plotted for data
obtained in 1970 and 1971. Circled numbers indicate actual percents in excess of 50%.
19
-------
% ODD 1970 a 1971
Figure 3. ODD as a percent of the total concentration of DDT, ODD, and DDE plotted for data
obtained in 1970 and 1971. Circled numbers indicate actual percents in excess of 50%.
20
-------
%DDE 1970 8 1971
:z: 0-10%
Figure 4. DDE as a percent of the total concentration of DDT, ODD, and DDE plotted for data
obtained in 1970 and 1971. Circled numbers indicate actual percents in excess of 5QP/
21
-------
TOTAL DOT DERIVATIVES
/ PPB 1970 a 1971
Figure 5. Total concentration in parts per billion of DDT, ODD, and DDE from data obtained in
1970 and 1971. Circled numbers indicate actual concentrations in excess of 50 ppb.
22
-------
TOTAL DDT DERIVATIVES
1973
Figure 6. Total concentration in parts per billion of DDT, ODD, and DDE from data obtained in
1973. The blank portions of the area were not sampled. Circled numbers indicate actual
concentrations in excess of 50 ppb.
23
-------
ANALYSIS OF DYNAMICS
An approach to the analysis of the dynamics of sediment systems has been developed
and has led to the development of Fortran programs permitting the rapid evaluation of
data. The discussion of the approach to analysis will refer to output from these programs.
The programs themselves with explanatory documentation are to be found in an appen-
dix at the end of this report.
The first program requires sampling at the same set of stations at two points in time.
The residue levels measured in sediments from the 19 stations sampled in both 1970
and 1973 constitute the data set used by this program. These data are presented as the
first two pages of output, see Tables 6 and 7, followed by two pages showing the per-
cent composition of total derivatives, see Tables 8 and 9. From the sums and means in
Tables 6 and 7 it would appear that while DDT has shown an increase of several-fold
the concentrations of DDD and DDE have changed very little. With respect to these
latter two compounds input must be rather closely balanced with respect to output
and decay. The changes in levels detected at individual stations must be a reflection
of the rates of input of new material, output or removal both geographically and into
other parts of the ecosystem, decay or decomposition within the sediment, and finally
a shifting about of the material from sampling station to sampling station due primarily
to the action of currents. The obvious complexity of the effect of these various rates
has made the analysis of such a system extremely difficult. The approach presented
here has necessitated the making of several simplifying assumptions. The utility of the
method and the validity of the assumptions must await further evaluation, and the
approach is intended more as a beginning than a final answer to the needs for methods
of data analysis.
Figure 7 presents a diagram of the essential features of the system as it is envisaged.
The individual stations where sediment samples were obtained are considered as com-
partments within the system of sediments in the southern portion of Monterey Bay.
The diagram indicates that this system has a relationship to all other systems both
geographical and of other kinds where the three compounds occur. Systems of dif-
ferent kinds would include the water above the sediment, the atmosphere above the
water, organisms, etc. The effect of the rate of input, I, the rate of output, O, the rate
of decay, D, and the rates of internal translocation, T] and TQ, on the concentration
within the system and within compartments is indicated.
A comparison of Figures 5 and 6 suggests that with continued input areas with the
higher concentrations tend to increase in concentration due to the movement of the
compounds within the system to these sinks or basins. Therefore, the amount of in-
crease within any sediment compartment would appear to be related to the concen-
tration already existing in that compartment. A similar relationship between the
amount of decrease and concentration is less easily deduced from these Figures.
However, the results of laboratory assays to be discussed in a later section have not
revealed either a saturation of the decay process nor a stimulation by induction and
selection of microbial populations that can be related to the concentration of these
compounds. Instead the amount of decomposition appears to be a function of con-
centration. That the amount of translocation would be similarly related to concen-
tration seems apparent.
24
-------
Table 6. Fl RST PAGE OF COMPUTER OUTPUT SHOWING CONCENTRATION OF
POLLUTANT COMPOUNDS IN SEDIMENT FROM SAMPLE STATIONS
AT FIRST SAMPLING TIME. Cj IDENTIFIES AS CONCENTRATIONS AT
TIME ONE.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
LOCATION
Latitude Longitude
3647.25 121 48.90
3646.85 121 53.50
3646.35 121 49.00
3646.05 121 51.00
3645.10121 50.00
3645.00121 49.00
3644.20121 52.25
3644.00121 49.50
3643.75 121 54.45
3643.35121 56.25
3643.18 121 57.00
3642.90121 58.00
3642.55121 53.30
3641.70121 55.00
3641.55 121 55.50
3640.90121 56.40
3639.10121 53.08
3637.95 121 52.50
3637.77 121 51.83
TOTALS
Mean
Standard Deviation
Standard Error
95% Confidence Limits
Date
8-23-70
11-15-70
2-20-70
11-15-70
5-29-70
2-20-70
8-23-70
2-20-70
11-15-70
8-23-70
2- 8-70
2-20-70
8-23-70
2-20-70
2- 8-70
2-20-70
2- 8-70
2-20-70
2- 8-70
DDT
(ppb)
8.36
1.63
5.71
4.28
6.42
3.67
5.20
0.69
1.02
1.12
0.0
0.0
13.20
1.22
0.0
1.32
2.44
2.65
0.49
59.4199
3.1274
+ 3.4385
t 0.7889
t 1.6574
ODD
(ppb)
3.67
6.76
0.71
6.61
8.67
0.40
7.50
0.14
0.38
0.25
5.00
0.35
5.73
0.53
2.35
1.61
0.66
2.79
0.21
54.3199
2.8589
* 2.9296
+ 0.6721
t 1.4121
DDE
(ppb)
5.76
14.70
1.02
10.70
7.01
0.45
15.50
2.75
0.70
0.65
20.50
1.92
13.00
2.40
7.01
9.02
2.40
10.00
0.50
125.9899
6.6310
+ 6.0673
t 1.3919
t 2.9245
TOTAL
17.79
23.09
7.44
21.59
22.10
4.52
28.20
3.58
2.10
2.02
25.50
2.27
31.93
4.15
9.36
11.95
5.50
15.44
1.20
239.7298
12.6174
+ 10.1773
t 2.3348
t 4.9055
1
25
-------
Table 7. SECOND PAGE OF COMPUTER OUTPUT SHOWING CONCENTRATION OF
POLLUTANT COMPOUNDS IN SEDIMENT FROM SAMPLE STATIONS AT
THE SECOND SAMPLING TIME. C2 IDENTIFIES AS CONCENTRATIONS AT
TIME TWO.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
LOCATION
Latitude Longitude
3647.25 121 48.90
3646.85 121 53.50
3646.35 121 49.00
3646.05121 51.00
3645.10121 50.00
3645.00121 49.00
3644.20121 52.25
3644.00121 49.50
3643.75 121 54.45
3643.35 121 56.25
3643.18 121 57.00
3642.90121 58.00
3642.55 121 53.30
3641.70121 55.00
3641.55 121 55.50
3640.90121 56.40
3639.10121 53.08
3637.95 121 52.50
Date
7- 9-73
7- 9-73
7- 9-73
7- 9-73
7- 2-73
7- 2-73
6-21-73
7- 2-73
8- 9-73
8- 9-73
6-21-73
8- 9-73
6-21-73
7-16-73
7-16-73
7-16-73
8- 9-73
6-21-73
DDT
(ppb)
1.06
9.50
1.10
3.63
0.92
2.18
30.60
0.96
5.41
72.70
63.10
0.93
29.90
1.14
0.68
0.70
1.18
83.10
3637.7712151.83 7-16-73 ; 0.54
ODD
(ppb)
0.53
11.40
0.53
5.43
0.39
0.72
6.07
0.06
4.54
3.19
0.79
0.90
4.32
2.74
2.20
1.11
0.42
0.95
0.20
DDE
(ppb)
0.56
17.50
0.63
6.91
0.52
0.83
11.20
0.23
17.30
12.00
3.48
6.06
12.20
10.49
8.67
5.67
2.44
3.34
0.40
TOTAL
2.15
38.40
2.26
15.97
1.83
3.73
47.87
1.25
27.25
87.89
67.37
7.89
46.42
14.37
11.55
7.48
4.04
87.39
1.14
TOTALS 309.3296 46.4899 120.4299 476.2488
Mean 16.2805 2.4468 6.3384 25.0657
Standard Deviation t 26.9909 t 2.8805 t 5.7417 t 29.2362
Standard Error t 6.1921 t 0.6608 t 1.3172 + 6.7072
95% Confidence Limits t 13.0097 +1.3884 + 2.7675 t 14.0919
26
-------
Table 8. THIRD PAGE OF COMPUTER OUTPUT SHOWING PERCENT OF TOTAL OF
EACH OF THE THREE COMPOUNDS IN SEDIMENTS FROM SAMPLE STATIONS
AT THE FIRST SAMPLING TIME. Cj IDENTIFIES AS DATA FOR TIME ONE.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
LOCATION
Latitude Longitude
3647.25 121 48.90
3646.85 121 53.50
3646.3512149.00
3646.05 121 51.00
3645.10121 50.00
3645.00121 49.00
3644.20121 52.25
3644.00121 49.50
3643.75 121 54.45
3643.35 121 56.25
3643.18 121 57.00
3642.90121 58.00
3642.55 121 53.30
3641.70121 55.00
3641.55 121 55.50
3640.90121 56.40
3639.10121 53.08
3637.95 121 52.50
3637.77 121 51.83
Date
8-23-70
11-15-70
2-20-70
11-15-70
5-29-70
2-20-70
8-23-70
2-20-70
11-15-70
8-23-70
2- 8-70
2-20-70
8-23-70
2-20-70
2- 8-70
2-20-70
2- 8-70
2-20-70
2- 8-70
DDT
(%)
46.99
7.06
76.75
19.82
29.05
81.19
18.44
19.27
48.57
55.45
0.0
0.0
41.34
29.40
0.0
11.05
44.36
17.16
40.83
ODD
(%)
20.63
29.28
9.54
30.62
39.23
8.85
26.60
3.91
18.10
12.38
19.61
15.42
17.95
12.77
25.11
13.47
12.00
18.07
17.50
DDE
(%)
32.38
63.66
13.71
49.56
31.72
9.96
54.96
76.82
33.33
32.18
80.39
84.58
40.71
57.83
74.89
75.48
43.64
64.77
41.67
TOTALS 586.7412 351.0149 962.2397
Mean 30.8811 18.4745 50.6442
Standard Deviation t 24.2998 t 8.6373 t 22.2953
Standard Error t 5.5748 t 1.9815 t 5.1149
95% Confidence Limits t 11.7126 t 4.1632 t 10.7464
27
-------
Table 9. FOURTH PAGE OF COMPUTER OUTPUT SHOWING PERCENT OF TOTAL
OF EACH OF THE THREE COMPOUNDS IN SEDIMENT FROM SAMPLE
STATIONS AT THE SECOND SAMPLING TIME. C2 IDENTIFIES AS DATA
FOR TIME TWO.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
LOCATION
Latitude Longitude
3647.25 121 48.90
3646.85 121 53.50
3646.35 121 49.00
3646.05121 51.00
3645.10121 50.00
3645.00121 49.00
3644.20121 52.25
3644.00121 49.50
3643.75121 54.45
3643.35 121 56.25
3643.18121 57.00
3642.90121 58.00
3642.55121 53.30
3641.70121 55.00
3641.55 121 55.50
3640.90121 56.40
3639.10121 53.08
3637.95 121 52.50
3637.77 121 51.83
Date
7- 9-73
7- 9-73
7- 9-73
7- 9-73
7- 2-73
7- 2-73
6-21-73
7- 2-73
8- 9-73
8- 9-73
6-21-73
8- 9-73
6-21-73
7-16-73
7-16-73
7-16-73
8- 9-73
6-21-73
7-16-73
DDT
(%)
49.30
24.74
48.67
22.73
50.27
58.45
63.92
76.80
19.85
82.72
93.66
11.79
64.41
7.93
5.89
9.36
29.21
95.09
47.37
ODD
(%)
24.65
29.69
23.45
34.00
21.31
19.30
12.68
4.80
16.66
3.63
1.17
11.41
9.31
19.07
19.05
14.84
10.40
1.09
17.54
DDE
(%)
26.05
45.57
27.88
43.27
28.42
22.25
23.40
18.40
63.49
13.65
5.17
76.81
26.28
73.00
75.06
75.80
60.40
3.82
35.09
TOTALS 862.1616 294.0427 743.7920
Mean 45.3769 15.4759 39.1469
Standard Deviation +29.2068 + 9.2122 +24.6220
Standard Error + 6.7005 + 2.1134 + 5.6487
95% Confidence Limits t 14.0777 + 4.4403 +11.8679
28
-------
OTHER SYSTEMS
n _. . .
SYSTEM
fc. 1 T
COMPARTMENT
C
-TO +TJ
C
COMPARTMENT
0.
ZC-DC
SYSTEMS
C2=C,(I + I-0-D)N
COMPARTMENTS
N
VC.CI+I+WO-D) „ _
C, = CONCENTRATION OF RESIDUE AT TIME I
C2=CONCENTRATION OF RESIDUE AT TIME 2
I = RATE OF fNPUT OF RESIDUE
0=RATE OF OUTPUT OF RESIDUE
D = RATE OF DECAY
T0 = RATE OF TRANSLOCATION OUT OF A COMPARTMENT
OF TRANSLOCATION INTO A COMPARTMENT
Figure 7. Model of the system of sediment compartments and this system's relation to other
systems.
29
-------
Therefore, for the estimation of the overall rate of change in a compartment, i.e., the
resultant of the various rates affecting concentration, the following expression was
solved for K,
GI and €2 are the concentrations within the compartment at time one and time two,
N is the length of the time interval in years, and e is the natural logarithm base. K is a
nominal percentage rate in the form of a decimal fraction resulting in continuous com-
pounding, and is converted to an annual rate for the expression,
C2 = C1(1+K)N (2)
The results of these calculations for the three compounds are presented as the fifth, sixth,
and seventh pages of computer output in Tables 10, 11, and 12. In these tables the values
of K are sorted into positive and negative values for purposes discussed below. Compart-
ments which showed a zero concentration at time one were adjusted by substitution of
0.004 ppb, a value generally just below the level of detection in the analyses.
The standard deviation of these estimates was approximated through the use of the ex-
pression for the standard deviation of a function of two random variables (Papoulis, 1965),
^ (IK)2 Q2 + (3K)2 a2 + 2^ ^ (T (3)
- ac/ c: X C2 'dcl dC2 Cic2
For ease in computation only two variables at a time were used in developing this ap-
proximation to the standard deviation.
If we assume that the rate of change within the system can be approximated by the mean
rate of change of its separate compartments, the mean of the K values becomes an esti-
mate of the rate of net change of the system.
Net rate of change = I - (O+D) (4)
This net rate of change is unaffected by the rates of internal translocation, Tj and TQ,
which are equal in magnitude and opposite in sign. The net rate of change is the sum of
two other mean rates. One is the rate of input, I, which can be estimated by the mean
of the positive K's, and the other is obtained as the mean of the negative K's and may
be taken as an estimate of (O+D) in equation 4.
The mean of the differences between each K and the net rate of change, that is the mean
deviation from the mean of K, becomes an estimate of TQ and Tj. The results of these
calculations are included in Tables 10, 11, and 12.
The separation of the rate O and D is more difficult and several approaches have been
attempted. The decimal fraction of the input rate that is translocated within the system,
Tj/I, differs from compound to compound: DDT, 0.665; ODD, 0.882; and DDE, 0.860.
One explanation for this difference is that they reflect differences in the rates of decom-
position within the sediments. Based upon this assumption the rate O and D have been
estimated by the following equations,
30
-------
TaBle 10. FIFTH PAGE OF COMPUTER OUTPUT SHOWING THE RATE OF CHANGE, K,
FOR DDT IN EACH SEDIMENT COMPARTMENT.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
C2DDT
1.06
9.50
1.10
3.63
0.92
2.18
30.60
0.96
5.41
72.70
63.10
0.93
29.90
1.14
0.68
0.70
1.18
83.10
0.54
CT DDT
8.36
1.63
5.71
4.28
6.42
3.67
5.20
0.69
1.02
1.12
0.0
0.0
13.20
1.22
0.0
1.32
2.44
2.65
0.49
N
2.8795
2.6493
3.3836
2.6493
3.0959
3.3644
2.8301
3.3644
2.7342
2.9644
3.3671
3.4685
2.8301
3.4027
3.4356
3.4027
3.5014
3.3342
3.4356
+K
0.0
0.9452
0.0
0.0
0.0
0.0
0.8706
0.1031
0.8408
3.0866
16.6503
3.8113
0.3350
0.0
3.4588
0.0
0.0
1.8105
0.0287
-K
-0.5119
0.0
-0.3854
-0.0603
-0.4661
-0.1434
0.0
0.0
0.0
0.0
0.0
0.0
0.0
-0.0197
0.0
-0.1701
-0.1874
0.0
0.0
+K + -K
-0.5119
0.9452
-0.3854
-0.0603
-0.4661
-0.1434
0.8706
0.1031
0.8408
3.0866
16.6503
3.8113
0.3350
-0.0197
3.4588
-0.1701
-0.1874
1.8105
0.0287
+K - Net
R
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
1.5079
15.0715
2.2325
0.0
0.0
1.8800
0.0
0.0
0.2317
0.0
-K - Net
R
-2.0906
-0.6336
-1.9641
-1.6390
-2.0449
-1.7222
-0.7082
-1.4756
-0.7380
0.0
0.0
0.0
-1.2438
-1.5985
0.0
-1.7488
-1.7661
0.0
-1.5501
Totals 309.3296 59.4199 60.0930 31.9407 -1.9442 29.9964 20.9236 -20.9235
Mean 16.2805 3.1274 3.1628 1.6811 -0.1023 1.5788 1.1012 -1.1012
S.D. +26.9909 1 3.4385 1 0.3 100 t 0.9016 t 0.0984 + 1.0000 + 0.8738 + 0.1262
S.E. t 6.1921 1 0.7889 +0.0711 + 0.2068 t 0.0226 t 0.2294 + 0.2005 t 0.0289
95% C.L. + 13.0097 1 1.6574 1 0.1494 t 0.4346 + 0.0474 t 0.4820 t 0.4212 + 0.0608
31
-------
Table 11. SIXTH PAGE OF COMPUTER OUTPUT SHOWING THE RATE OF CHANGE, K,
FOR ODD IN EACH SEDIMENT COMPARTMENT.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
C9 ODD
e.
0.53
11.40
0.53
5.43
0.39
0.72
6.07
0.06
4.54
3.19
0.79
0.90
4.32
2.74
2.20
1.11
0.42
0.95
0.20
CT ODD
3.67
6.76
0.71
6.61
8.67
0.40
7.50
0.14
0.38
0.25
5.00
0.35
5.73
0.53
2.35
1.61
0.66
2.79
0.21
N
2.8795
2.6493
3.3836
2.6493
3.0959
3.3644
2.8301
3.3644
2.7342
2.9644
3.3671
3.4685
2.8301
3.4027
3.4356
3.4027
3.5014
3.3342
3.4356
+K
0.0
0.2181
0.0
0.0
0.0
0.1909
0.0
0.0
1.4774
1.3607
0.0
0.3130
0.0
0.6206
0.0
0.0
0.0
0.0
0.0
-K
-0.4893
0.0
-0.0828
-0.0715
-0.6328
0.0
-0.0720
-0.2226
0.0
0.0
-0.4219
0.0
-0.0950
0.0
-0.0190
-0.1035
-0.1211
-0.2761
-0.0141
+K + -K
-0.4893
0.2181
-0.0828
-0.0715
-0.6328
0.1909
-0.0720
-0.2226
1.4774
1.3607
-0.4219
0.3130
-0.0950
0.6206
-0.0190
-0.1035
-0.1211
-0.2761
-0.0141
+K - Net
R
0.0
0.1360
0.0
0.0
0.0
0.1089
0.0
0.0
1.3953
1.2787
0.0
0.2309
0.0
0.5386
0.0
0.0
0.0
0.0
0.0
-K - Net
R
-0.5714
0.0
-0.1648
-0.1536
-0.7148
0.0
-0.1541
-0.3047
0.0
0.0
-0.5039
0.0
-0.1770
0.0
-0.1011
-0.1856
-0.2031
-0.3582
-0.0961
Totals 46.4899 54.3199 60.0930 4.1806 -2.6218 1.5588 3.6884 -3.6884
Mean 2.4468 2.8589 3.1628 0.2200 -0.1380 0.0820 0.1941 -0.1941
S.D. 1 2.8805 +2.9296 +0.3100 t 0.7233 1 0.2767 +1.0000 +0.7233 +0.2767
S.E. +0.6608+0.6721+0.0711 1 0.1659 +0.0635 +0.2294 +0.1659 +0.0635
95% C. L. + 1.3884 + 1.4121 +0.1494 +0.3486 +0.1334 t 0.4820 +0.3486 +0.1334
32
-------
Table 12. SEVENTH PAGE OF COMPUTER OUTPUT SHOWING THE RATE OF CHANGE, K,
FOR DDE IN EACH SEDIMENT COMPARTMENT.
Station
1
2
3
4
10
11
14
16
17
19
20
22
23
25
26
29
34
36
37
C2DDE
0.56
17.50
0.63
6.91
0.52
0.83
11.20
0.23
17.30
12.00
3.48
6.06
12.20
10.49
8.67
5.67
2.44
3.34
0.40
G! DDE
5.76
14.70
1.02
10.70
7.01
0.45
15.50
2.75
0.70
0.65
20.50
1.92
13.00
2.40
7.01
9.02
2.40
10.00
0.50
N
2.8795
2.6493
3.3836
2.6493
3.0959
3.3644
2.8301
3.3644
2.7342
2.9644
3.3671
3.4685
2.8301
3.4027
3.4356
3.4027
3.5014
3.3342
3.4356
+K
0.0
0.0680
0.0
0.0
0.0
0.1996
0.0
0.0
2.2318
1.6740
0.0
0.3929
0.0
0.5426
0.0638
0.0
0.0047
0.0
0.0
-K
-0.5549
0.0
-0.1327
-0.1522
-0.5684
0.0
-0.1085
-0.5217
0.0
0.0
-0.4094
0.0
-0.0222
0.0
0.0
-0.1275
0.0
-0.2803
-0.0629
+K + -K
-0.5549
0.0680
-0.1327
-0.1522
-0.5684
0.1996
-0.1085
-0.5217
2.2318
1.6740
-0.4094
0.3929
-0.0222
0.5426
0.0638
-0.1275
0.0047
-0.2803
-0.0629
+K - Net
R
0.0
0.0
0.0
0.0
0.0
0.0818
0.0
0.0
2.1141
1.5563
0.0
0.2752
0.0
0.4249
0.0
0.0
0.0
0.0
0.0
-K - Net
R
-0.6726
-0.0497
-0.2505
-0.2699
-0.6861
0.0
-0.2262
-0.6394
0.0
0.0
-0.5272
0.0
-0.1399
0.0
-0.0539
-0.2453
-0.1130
-0.3980
-0.1806
Totals 120.4299125.9899 60.0930 5.1774 -2.9407 2.2367 4.4522 -4.4522
Mean 6.3384 6.6310 3.1628 0.2725 -0.1548 0.1177 0.2343 -0.2343
S.D. +5.7417 +6.0673 f 0.3100 +0.7781. +0.2243 + 1.0024 +0.7761 +0.2262
S.E. + 1.3172 t 1.3919 +0.0711 + 0.1785 +0.0515 + 0.2300 +0.1781 f 0.0519
95% C.L +2.7675 +2.9245 tO. 1494 +0.3750 1 0.1081 +0.4831 +0.3741 +0.1091
33
-------
O = T!(O+D) (5)
i
D = (1.0-T!)(O+D)orD = (O+D)-O (6)
I
The residence time, TR, and lifetime, TL, in years, are calculated as the corresponding
reciprocals.
TR=1.0/(0+D) (7)
TL=1.0/D (8)
The last three pages of computer output present a summary of these estimations and
are presented in Tables 13, 14, and 15.
The effect of substitution of'a minimal value for zero concentrations was investigated
by reducing the set of sample stations to sixteen and elimination of all stations showing
a zero concentration of DDT at time one. While there was some effect upon the esti-
mates of rates as the system was reduced in size, only the estimates of TQ for DDT
were significantly different when tested by the "test of equality of the means of two
samples whose variances are assumed to be unequal" (Sokol and Rohlf, 1969). The
difference between the other estimates was very small compared to the standard
deviation of these estimates. Table 16 presents for comparison the set of rates for
the nineteen and sixteen station data sets.
The approach to analysis of the data which provided these estimates of system rates
requires sampling at the same stations at two different times. However, as presented
in Table 3, there is additional data available with respect to the south bay system at
time one. This additional data can not be used by the approach to analysis presented
so far. More stations were sampled in the first sampling period than were sampled in
the second, and the approach requires pairs of samples identical except for time of
sampling. An additional program was written to permit analysis of a system where
sampling does not meet the requirements of the first approach. This second program
treats all samples as unpaired and evaluates the rate of change, K, at the different
sample locations by comparison of the actual measurement at that station at time
one or time two with the mean concentrations of the system at either time one or
time two. That is, a measurement at time one is paired with the mean concentration
at time two and vice versa for the evaluation of K. Further the time interval, N, is
evaluated as the interval between the time of actual sample of one sampling time
and the mean time of the other sampling period. Equation 1 becomes,
C2 = G! eKN 9.
with N = T2 - TI
34
-------
Table 13. EIGHTH PAGE OF COMPUTER OUTPUT SHOWING A SUMMARY OF THE ANNUAL
SYSTEM RATES EXPRESSED AS DECIMAL FRACTIONS OF THE MEAN CONCEN-
TRATION OF DDT PRESENT IN THE SYSTEM.
System of Rates for DDT
Net rate of change =
Translocation into compartments =
Translocation out of compart-
ments
Input
Output and Decay =
Output from System =
Decay
Lifetime in years =
Residence time in years =
Summary Equation for the System-
DDT Mean C2 Mean C-j
16.2805 = 3.1274
Net = + 1.5788
T| = + 1.1012
T0 = - 1.1012
I = + 1.6811
0+D = - 0.1023
0 = - 0.0670
D = - 0.0353
TL = 28.3322
TR = 9.7724
' Tl
(1.0+ 1.6811 + 1.1012
S.D.
* 1.0000 t
+ 0.8738 +
+ 0.1262 "*"
f 0.9016 t
t 0.0984 t
+ 0.0644 t
+- 0.0339 t
+ 27.2386 t
* 9.3952 t
T0 0
- 1.1012-0.0670-
S.E.
0.2294
0.2005
0.0289
0.2068
0.0226
0.0148
0.0078
6.2490
2.1554
D
0.0353)
95%
Limit
+ 0.4820
* 0.4212
* 0.0608
+ 0.4346
+ 0.0474
"*" 0.0311
* 0.0164
t 13.1291
•t 4.5285
N
3.1628
35
-------
Table 14. NINTH PAGE OF COMPUTER OUTPUT SHOWING A SUMMARY OF THE ANNUAL
SYSTEM RATES EXPRESSED AS DECIMAL FRACTIONS OF THE MEAN CONCEN-
TRATION OF ODD PRESENT IN THE SYSTEM.
System of Rates for ODD
Net rate of change
Translocation into
compartments
Translocation out of
compartments
Input
Output and Decay
Output from System
Decay
Lifetime in years
Residence time in years
Summary Equation for the
ODD Mean C2
2.4468 =
= Net = + 0.0820 *
= T| = + 0.1941 t
= T0 = - 0.1941 t
= I = + 0.2200 t
= 0+D = - 0.1380 ^
= 0 = - 0.1217 t
= D = - 0.0162 t
= TL = 61.5459 *
= TR = 7.2469 +
System -
Mean GI 1 T|
2.8589 (1.0 + 0.2200 + 0.1941-
S.D.
1.0000 t
0.7233 f
0.2767 t
0.7233 t
0.2767 +-
0.2441 t
0.0326 t
123.4241 *
14.5330 +
T0 0
0.1941 -0.1217
S.E.
0.2294
0.1659
0.0635
0.1659
0.0635
0.0560
0.0075
28.3154
3.3341
D
- 0.0162)
'
95%
Limit
+- 0.4820
* 0.3486
"*" 0.1334
* 0.3486
t 0.1334
* 0.1177
* 0.0157
? 59.4907
+ 7.0049
N
3.1628
36
-------
Table 15. TENTH PAGE OF COMPUTER OUTPUT SHOWING A SUMMARY OF THE ANNUAL
SYSTEM RATES EXPRESSED AS DECIMAL FRACTIONS OF THE MEAN CONCEN-
TRATION OF DDE PRESENT IN THE SYSTEM.
95%
System of Rates for DDE S.D. S.E. Limit
Net rate of change = Net = + 0.1177 - 1.0024 - 0.2300 - 0.4831
Translocation into
compartments = T( = + 0.2343 + 0.7761 + 0.1781 + 0.3741
Translocation out of
compartments = TQ = - 0.2343 * 0.2262 * 0.0519 + 0.1091
Input =l =+ 0.2725 t 0.7781 t 0.1785 t 0.3750
Output and Decay = 0+D = - 0.1548 + 0.2243 t 0.0515 t 0.1081
Output from System =0 =- 0.1331 + 0.1929 "*" 0.0442 t 0.0930
Decay =D =- 0.0217 + 0.0314 + 0.0072 t 0.0151
Lifetime in years = TL = 46.1286 t 66.8453 t 15.3354 t 32.2196
Residence time in years = TR =~ 6.4611 t 9.3629 t 2.1480 t 4.5129
Summary Equation for the System-
DDE Mean C2 Mean GI I T( TQ 0 OD N
6.3384 = 6.6310 (1.0 + 0.2725 + 0.2343 - 0.2343 - 0.1331 - 0.0217) 3-1628
37
-------
Table 16. COMPARISON OF ESTIMATES OBTAINED FROM THE 16 AND 19 STATION DATA SETS AND
USING ACTUAL PAIRED SAMPLE ANALYSES STANDARD DEVIATIONS [S.D.] AND COEF-
FICIENTS OF VARIATION [C.V.] ARE INCLUDED.
DDT
C1 (ppb)
C2 (ppb)
Net
I
0+D
T0
T,
0
D
TL (years)
TR (years)
ODD
CT (ppb)
C2 (ppb)
Net
I
0+D
T0
T|
0
D
T|_ (years)
TR (years)
DDE
G! (ppb)
C2 (ppb)
Net
I
0+ D
T0
T,
0
D
TL (years)
TR (years)
16 STATION DATA SET
Estimate
3.7137
15.2887
+ 0.3798
+ 0.5013
- 0.1215
- 0.3534
+ 0.3534
- 0.0857
- 0.0358
27.9014
8.2294
2.9137
2.6625
+ 0.1054
+ 0.2417
- 0.1363
- 0.2088
+ 0.2088
- 0.1177
- 0.0186
53.8306
7.3364
6.0350
6.3887
+ 0.1368
+ 0.2950
- 0.1582
- 0.2563
+ 0.2563
- 0.1374
- 0.0208
48.1189
6.3211
S.D.
+ 3.4446
± 26.3645
+ 1.0000
+ 0.7556
+ 0.2444
± 0.2591
± 0.7409
+ 0.1723
+ 0.0721
+ 56.1105
+ 16.5496
± 3.0908
+ 3.0921
+ 1 .0000
i 0.7279
+ 0.2721
i 0.2721
+ 0.7279
+ 0.2350
+ 0.0371
+ 107.4660
+ 14.6462
+ 5.4299
+ 6.2166
± 1.0030
+ 0.7843
+ 0.2186
+ 0.2211
+ 0.7818
± 0.1899
+ 0.0287
+ 66.4924
+ 8.7347
C.V.
(%)
92.8
172.4
263.3
150.7
138.2
73.3
209.6
201.1
201.4
201.1
201.1
106.1
116.1
948.8
301.2
199.6
130.3
348.6
199.7
199.5
199.6
199.6
90.0
97.3
733.2
265.9
138.2
86.3
305.0
138.2
138.0
138.2
138.2
19 STATION DATA SET
Estimate
3.1274
16.2805
+ 1.5788
+ 1.6811
- 0.1023
- 1.1012
+ 1.1012
- 0.0670
- 0.0353
28.3322
9.7724
2.8589
2.4468
+ 0.0820
+ 0.2200
- 0.1380
- 0.1941
+ 0.1941
+ 0.1217
- 0.0162
61.5459
7.2469
6.6310
6.3384
+ 0.1177
+ 0.2725
- 0.1548
- 0.2343
+ 0.2343
- 0.1331
- 0.0217
« 46.1286
6.4611
S.D.
+ 3.4385
+ 26.9909
+ 1 .0000
+ 0.9016
+ 0.0984
+ 0.1262
+ 0.8738
+ 0.0644
+ 0.0339
i 24.2386
+ 9.3952
+ 2.9296
i 2.8805
t 1 .0000
+ 0.7233
+ 0.2767
+ 0.2767
+ 0.7233
+ 0.2441
+ 0.0326
±123.4241
± 14.5330
+ 6.0673
+ 5.7417
± 1 .0024
+ 0.7781
+ 0.2243
+ 0.2262
+ 0.7761
+ 0.1929
+ 0.0314
+ 66.8543
+ 9.3629
C.V.
(%)
109.9
165.8
63.3
53.6
96.2
11.5
79.3
96.1
96.0
96.1
96.1
102.5
117.7
1219.5
328.8
200.5
142.6
5009.4
200.6
201.2
200.5
200.5
91.5
90.6
851.7
285.5
144.9
96.5
332.1
144.9
144.7
144.9
144.9
38
-------
T2 = mean time of second sampling period
TI = time of actual sampling in first sampling period
andC2 = C1eKN (10)
withN = T2-T;i
T2 = time of actual sampling in second sampling period
TI = mean time of first sampling period.
Table 17 presents the estimates of the system obtained using this pairing with means
approach. Once again the effect of substitution of a minimal value for zero concentra-
tions was explored by eliminating stations with zero concentration thus providing the
subset of 49 samples from the complete set of 57. Except for the estimates of TQ for
DDT, there was no significant difference between the two sets of estimates once again,
nor are these estimates significantly different from either of the sets of estimates based
on the 16 and 19 station data sets. The principal effect of inclusion or exclusion of the
zero level values with substitution of a minimal value is upon the estimates of the rates
of input, I, translocation, Tj and TQ, and the net rate. The stations showing a zero
concentration of DDT at time one show high positive rates of change, and therefore,
have a particularly marked effect on the positive rate estimates as well as those based
to at least some extent upon these positive rate estimates.
The second approach which uses sample values paired to mean values should find use
in the analysis of systems where real paired values are impossible to obtain. Animals
which are sacrificed at the time of sampling obviously can not be resampled at another
point in time. The use of sample values at one sample time paired to the mean value
of another permits estimation of system rates for the population. The comparison be-
tween the two approaches to these estimates that is presented here indicates that the
use of mean values in pairing gives a close approximation of rate estimates obtained
with real paired values.
Both of these approaches to the estimation of system rates are dependent upon vari-
ability in concentration level and rate of c.iange within compartments. It is essential
to these methods of analysis that individual compartments show the effect of the
various processes to different degrees. If all the concentration levels and rates of
change within compartments were the same, it would be possible to gain an estimate
of net rate of change only. Therefore, these approaches to estimation of system rates
are dependent upon variability in environmental samples of the system and make use
of this variability for estimating the rates of the various processes.
39
-------
Table 17. COMPARISON OF ESTIMATES OBTAINED FROM THE 49 AND 57 STATION DATA SETS AND
USING SAMPLE ANALYSES PAIRED WITH MEAN CONCENTRATION LEVELS. STANDARD
DEVIATIONS [S.D.] AND COEFFICIENTS OF VARIATION [C.V.] ARE INCLUDED.
DDT
G! (ppb)
C2 (ppb)
.Net
I
0+D
TO
T|
0
D
TL (years)
TR (years)
ODD
CT (ppb)
C2 (ppb)
Net
I
0+D
TO
T|
0
D
T|_ (years)
TR (years)
DDE
G! (ppb)
C2 (ppb)
Net
I
0+D
TO
T|
0
D
TL (years)
TR (years)
49 SAMPLE DATA SET
Estimate
3.9576
15.4975
+ 0.5905
+ 0.6819
- 0.0913
- 0.3234
+ 0.3234
- 0.0433
- 0.0480
20.8292
10.9502
2.4107
2.3435
+ 0.1283
+ 0.2703
- 0.1420
- 0.2095
+ 0.2095
- 0.1101
• 0.0319
31.3031
7.0424
5.1138
6.1575
+ 0.1748
+ 0.2802
- 0.1054
- 0.1946
+ 0.1946
- 0.0732
• 0.0322
31.0400
9.4853
S.D.
+ 4.1746
+ 26.5034
± 1 .0000
± 0.6374
± 0.3626
+ 0.3966
i 0.6034
+ 0.1720
+ 0.1906
+ 82.7111
i 43.4823
+ 2.5354
+ 2.8415
+ 1.0000
± 0.6357
+ 0.3643
+ 0.3653
+ 0.6347
+ 0.2823
+ 0.0820
+ 80.3119
i 18.0682
± 4.4111
± 5.6469
± 1.0010
± 0.6628
± 0.3382
+ 0.3466
+ 0.6544
+ 0.2348
± 0.1033
± 99.5728
± 30.4277
C.V.
(%)
105.4
171.0
169.3
93.5
397.2
122.6
186.6
397.2
397.1
397.1
397.1
105.2
121.3
779.4
235.2
256.5
174.4
303.0
256.4
257.1
256.6
256.6
86.3
91.7
572.7
236.5
320.9
178.1
336.3
320.8
320.8
320.8
320.8
57 SAMPLE DATA SET
Estimate
3.1019
15.4975
+ 2.2567
+ 2.3233
- 0.0667
- 1.4256
+ 1.4256
- 0.0409
- 0.0258
38.8090
14.9951
2.2743
2.3435
+ 0.1587
+ 0.2813
- 0.1226
- 0.2039
+ 0.2039
- 0.0889
- 0.0337
29.6518
8.1558
5.3681
6.1575
+ 0.1793
+ 0.2785
- 0.0993
- 0.1906
+ 0.1906
- 0.0679
- 0.0314
31.8905
10.0735
S.D.
+ 4.0336
+ 26.5034
± 1.0000
± 0.9204
+ 0.0796
+ 0.1513
+ 0.8487
+ 0.0488
+ 0.0307
± 46.2947
± 17.8875
± 2.3532
± 2.8415
+ 1.0000
+ 0.6329
± 0.3671
+ 0.3698
+ 0.6311
t 0.2662
± 0.1010
± 88.7883
± 24.4216
± 4.8069
+ 5.6469
± 1.0009
i 0.6787
± 0.3222
± 0.3311
± 0.6697
+ 0.2204
± 0.1018
± 103.4957
± 32.6922
C.V.
(%)
130.0
171.0
44.3
39.6
119.3
10.6
59.5
119.3
119.0
119.3
119.3
103.5
121.3
630.1
225.0
299.4
180.9
309.5
299.4
299.7
299.4
299.4
89.5
91.7
558.2
243.7
324.5
173.7
351.4
324.6
324.2
324.5
324.5
40
-------
For any set of estimates of I, (O+D), Tj and TQ, based on a number of samples, n,
there is a distribution of K's with a minimal variance. The members of the distribu-
tion can be determined through one of the following sets of equations:
Where the net rate of change, I + (O+D), is positive,
j = nl - nTj and j is an integer obtained without rounding. (11)
I + (O+D)
I + (O+D) + njj = Kj, K2 . . . K: (12)
If
K Lnl
1
K + l (13)
n(0+D) = Kj + 2,Kj + 3...Kn (14)
n-j-1
If
K = nl
1
n(0+D) = KJ + l,Kj + 2...Kn (15)
n-j
Where the net rate of change, I + (O+D), is zero,
— and j is an integer obtained without rounding. (16)
nTi_= K1,K2..'.Kj (17)
j
njo = Kj + l.Kj + 2... K2j (18)
If 2j L n,
Kn = 0.0 (19)
Where the net rate of change, I + (O+D), is negative,
j = n(O+D) - nTQ (20)
I + (O+D)
I + (0+D) + nTn = Ki K7 . . . K: (21)
5=L -1! *•» J
j
41
-------
J
If £_ K L n (O+
1
n(0+D)-jK1 = Kj+1 (22)
_nl_= Kj+2, Kj+3 . . . Kn (23)
n-j-1 .
J
If
1
..Kn (24)
n-j
The variances of such distributions are the minimal variances that will permit the estima-
tions of I, Tj and TQ, and (O+D) with a given number of samples. This variance is less af-
fected by the number of samples than it is by the difference between the values of I, Tj
and TQ, and (O+D) as can be seen in Table 18. The lowest standard deviations are observed
where Tj is low. Where I is increased relative to Tj, the standard deviation is reduced as well
but not to the same extent. For example, I = 2.0, Tj = 1.2 has a ratio of 0.6 as does I = 1.5,
Tj = 0.9, however, the latter has the lower standard deviation. The unavoidable variance
related to any series of values of I, Tj and TQ, O+D, and n has significance to survey design.
The greater the amount of internal translocation due to Tj and TQ the greater the unavoid-
able variance of the estimation of K. Increasing the number of sampling points has only a
minor effect upon the variance although it has a marked effect upon the standard error and
95% confidence limits of the estimates.
The corrected standard deviations with associated standard errors and 95% confidence
limits can be calculated using Subroutine FACTOR which will be found in the Appendix.
The correction is imposed following the calculation of the standard deviation of K using
equation 3, but only with respect to first moment as is true for the other estimations of
standard deviations.
The variance is corrected as follows,
2 2 \2 2 2
s - s \ s = s (25)
K calc. Mm. \ K K corr.
2
SK calc.
2 2
Where s is the variance calculated by equation 3, smjn is the variance of the distribution
of K's with minimal variance, s,, . is the variance of the distribution of K's calculated
by equation 3, and Sj, is the corrected variance of K. This correction appears to be
justifified because the variance of interest is that which is related to the variance of a sys-
tem with particular characteristics as compared to a similar system with minimal unavoid-
able variance. Table 19 presents a comparison of uncorrected standard deviations from
Tables 16 and 17 and the corresponding corrected values. The system estimates for
42
-------
Table 18. STANDARD DEVIATIONS AND STANDARD ERRORS OF DISTRIBUTIONS OF K WITH MINIMAL
VARIANCE FOR GIVEN VALUES OF I, T, AND TQ, (0+D) AND n.
1
2.00
1.75
1.50
1.50
1.50
1.50
1.50
Tl
1.20
1.20
1.20
1.20
1.20
0.90
0.60
0+D
-0.15
-0.15
-0.15
-0.30
-0.60
-0.15
-0.15
Net
1.85
1.60
1.35
1.20
0.90
1.35
1.35
n = 5
S.D. S.E.
+ 2.7524
* 3.4084
+- 3.3586
+- 3.3719
t 3.4249
+- 2.0724
"t 1.4335
+ 1.2309
+ 1.5243
+ 1.5020
+ 1.5080
t 1.5317
* 0.9268
+ 0.6411
n= 10
S.D. S.E.
t 2.5965
t 2.7758
+ 3.1663
+ 3.1785
* 3.2267
+ 1.9558
t 1.3528
+ 0.8211
"!" 0.8778
* 1.0013
t 1.0051
* 1.0204
+- 0.6185
t 0.4278
n = 20
S.D. S.E.
t 2.5338
^ 2.7107
t 3.0831
^ 2.8433
+ 2.7077
* 1.9124
t 1.3063
"!" 0.5666
* 0.6061
^ 0.6894
* 0.6358
* 0.6055
+- 0.4276
"*" 0.2921
43
-------
Table 19. COMPARISON OF UNCORRECTED AND CORRECTED STANDARD DEVIATIONS OF
SYSTEM ESTIMATES
DDT
Net
I
0+D
TO
T|
0
D
TL
TR
ODD
Net
I
0+D
TO
T|
0
D
TL
TR
DDE
Net
I
0+D
TO
T|
0
D
TL
TR
16Sarm
Uncorrected
± 1.0000
+ 0.7556
± 0.2444
± 0.2591
± 0.7409
± 0.1723
± 0.0721
± 56.1105
± 16.5496
± 1 .0000
± 0.7279
± 0.2721
± 0.2721
+ 0.7279
± 0.2350
± 0.0371
+107.4660
± 14.6462
± 1.0030
+ 0.7843
± 0.2186
+ 0.2211
± 0.7818
± 0.1899
± 0.0287
± 66.4924
± 8.7347
Die Set
Corrected
± 0.2751
+ 0.2806
± 0.0907
± 0.0962
± 0.2751
± 0.0640
± 0.0268
±20.8365
± 6.1457
± 0.3604
± 0.2623
± 0.0981
± 0.0981
+ 0.2623
± 0.0847
± 0.0134
+38.7279
± 5.2781
± 0.3602
± 0.2817
± 0.0785
± 0.0794
± 0.2808
± 0.0682
± 0.0103
±23.8815
± 3.1372
19 Sample Set
Uncorrected
± 1 .0000
+ 0.9016
± 0.0984
± 0.1262
± 0.8738
± 0.0644
± 0.0339
± 27.2386
+ 9.3952
+ 1.0000
± 0.7233
± 0.2767
± 0.2767
± 0.7233
± 0.2441
± 0.0326
±123.4241
± 14.5330
± 1.0024
+ 0.7781
± 0.2243
+ 0.2262
+ 0.7761
± 0.1929
+ 0.0314
+ 66.8543
± 9.3629
Corrected
± 0.5986
± 0.5397
i 0.0589
± 0.0755
± 0.5231
+ 0.0386
+ 0.0203
±16.3047
± 5.6239
± 0.3860
± 0.2792
± 0.1068
± 0.1068
+ 0.2792
± 0.0942
+ 0.0126
±47.6463
± 5.6103
+ 0.4716
+ 0.3661
+ 0.1055
+ 0.1064
+ 0.3651
± 0.0907
± 0.0148
±31 .4484
± 4.4049
49 Sam
Uncorrected
± 1.0000
± 0.6374
± 0.3626
± 0.3966
± 0.6034
± 0.1720
± 0.1906
+82.7111
+43.4823
± 1.0000
+ 0.6357
+ 0.3643
± 0.3653
± 0.6347
± 0.2823
+ 0.0820
±80.3119
±18.0682
± 1.0010
± 0.6628
± 0.3382
± 0.3466
± 0.6544
+ 0.2348
± 0.1033
±99.5728
±30.4277
pie Set
Corrected
+ 0.3366
± 0.2145
± 0.1221
± 0.1335
+ 0.2031
± 0.0579
± 0.0642
+27.8380
+ 14.6348
± 0.3419
+ 0.2174
± 0.1246
+ 0.1249
± 0.2170
± 0.0965
± 0.0280
±27.4603
± 6.1779
+ 0.4379
± 0.2900
±.0.1479
± 0.1516
± 0.2863
± 0.1027
± 0.0452
±43.5593
±13.3109
57 Sam
Uncorrected
+ 1.0000
+ 0.9204
± 0.0796
± 0.1513
± 0.8487
± 0.0488
± 0.0307
± 46.2947
+ 17.8875
± 1 .0000
± 0.6329
± 0.3671
± 0.3689
± 0.6311
± 0.2662
± 0.1010
± 88.7883
+ 24.4216
± 1 .0009
+ 0.6787
± 0.3222
± 0.3311
± 0.6697
± 0.2204
± 0.1018
±103.4957
± 32.6922
pie Set
Corrected
+ 0.5379
+ 0.4951
± 0.0428
± 0.0814
± 0.4565
± 0.0263
± 0.0165
±24.9713
± 9.6215
± 0.3521
± 0.2228
± 0.1293
± 0.1299
± 0.2222
± 0.0937
± 0.0356
±31.2619
+ 8.5987
± 0.4545
± 0.3082
± 0.1463
± 0.1504
+ 0.3041
± 0.1001
+ 0.0462
+46.9942
+14.8445
44
-------
DDT obtained from the four data sets did show some significant differences when
compared using these corrected estimates of the standard deviation. The estimates
obtained with the 49 and 57 sample sets were significantly different at the .05 level
for Net, I, TQ, and Tj. The estimates obtained with the 16 and 57 sample sets were
significantly different for Net, I, and TQ, and the estimates of TQ for the 19 and 57
data sets were also significantly different. These differences would appear to be
primarily the result of inclusion or exclusion from the system of sites where there are
major increases in the concentration of DDT rather than the effect of substitution of
a minimal value for the concentration at time one. The estimation of TQ in systems
showing a positive Net rate of change are particularly sensitive to significance testing
due to their relatively low standard deviations that result from the distribution of
variance between Tj and TQ.
If we keep in mind the limitations imposed by the variability of the data, the estimates
can be used to gain a picture of the flux of these pollutants in the study area. The area
of south Monterey Bay is approximately 280 square kilometers, or 69,190 acres in size.
The density of the sediments on a dry weight basis averages 1.32 grams per cm'. Table
20 gives the mean of the estimates for system concentrations and rates that were ob-
tained by the two approaches to analysis and the four data sets. Standard deviations,
standard errors, 95% confidence limits, and coefficients of variation for these means are
included. These latter descriptive statistics refer only to the variation of the estimates
and do not include the effect of compartment variability discussed above.
Table 21 uses the mean of the estimates and gives the total amounts of these chlorinated
hydrocarbons in the area and the concentration in pounds per acre based upon the mean
concentrations at the two times of sampling. These total amounts are estimated as being
present in the top 10 cm of sediment, a depth generally sampled with the collecting gear
used. Considering that the usual level of application on land is 2 pounds to the acre the
total level of these compounds per acre has reached somewhat more than 1/100 of the
land applications level.
The estimated annual rates of input, I, as seen in Table 20, average 130% for DDT, 25%
for DDD, and 28% for DDE. The corresponding amounts of these materials expected in
the next year are indicated in Table 21. Expected loss due to translocation, output, and
decay based on the estimated annual rates, O+D, 10% for DDT, 13% for DDD, and 13%
for DDE, are also shown. The resulting net effect for the year period following the last
sample time in 1973 gives the expected values shown, Table 21. The expected change in
the amount of the total chlorinated hydrocarbons derived from DDT amounts to an in-
crease of 182%. The amounts translocated within the system are presented in Table 21
along with a separation of the expected loss into that expected from output and decay.
All of the projections, of course, assume that the estimated rates reflecting flux of these
materials in the past three years will persist for the next year period.
The K values for the individual compartments can also be used to present a composite
view of the translocation of the three compounds within the system and principal points
of geographical exit. The stations at their geographical location are connected with arrows
45
-------
pointing from more negative to less negative K values and ending in basins with positive
K values. The result is a kinematic graph representing the movement of these materials
within the system. It is composite with respect to the time interval under consideration
and would appear to represent the result of several events of translocation. Figure 8
presents such a graph developed for the 19 station data set. The large double arrows in-
dicate the main offshore forces that drive the inshore circulation and correlated with the
kinematic expression of circulation within the system.
46
-------
Table 20. MEAN OF THE ESTIMATES FOR THE SOUTH MONTEREY BAY SYSTEM AND
ASSOCIATED DESCRIPTIVE STATISTICS.
C1 DDT (ppb)
DDD (ppb)
DDE (ppb)
C2 DDT (ppb)
DDD (ppb)
DDE (ppb)
Net DDT
DDD
DDE
I DDT,
DDD
DDE
0 + D DDT
DDD
DDE
TQ & T! DDT
DDD
DDE
0 DDT
DDD
DDE
D DDT
DDD
DDE
TL DDT (years)
DDD (years)
DDE (years)
TR DDT (years)
DDD (years)
DDE (years)
Mean
3.4752
2.6144
5.7870
15.6411
2.4491
6.2605
+ 1.2015
+ 0.1186
+ 0.1522
+ 1.2969
+ 0.2533
+ 0.2816
- 0.0955
- 0.1347
- 0.1294
+- 0.8009
+ 0.2041
* 0.2190
- 0.0592
- 0.1096
- 0.1029
- 0.0362
- 0.0251
- 0.0265
28.9680
44.0829
39.2945
10.9868
7.4454
8.0853
S.D.
+ 0.4281
+- 0.3196
t 0.6837
S.E.
* 0.2141
*• 0.1598
+ 0.3419
+ 0.4375 t 0.2188
* 0.1504 t 0.0752
+ 0.1207 t 0.0604
+ 0.8764
+ 0.0327
+ 0.0298
t 0.8587
t 0.0278
+ 0.0096
+ 0.0096
+ 0.0229
+ 0.0084
+ 0.5504
+ 0.0071
+ 0.0318
+ 0.0212
+ 0.0146
+ 0.0375
+- 0.0091
+ 0.0090
+ 0.0061
* 7.4078
1 16.0370
t 9.0835
* 2.8952
* 0.4893
* 1.9717
t 0.4382
t 0.0164
^ 0.0149
+ 0.4294
* 0.0139
* 0.0048
* 0.0114
* 0.0042
"!" 0.0157
* 0.2752
* 0.0036
+ 0.0159
* 0.0106
+ 0.0073
* 0.0187
^ 0.0045
+ 0.0143
"^ 0.0031
+ 3.7039
* 8.0185
t 4.5418
+ 1.4476
t 0.2447
+ 0.9859
95%C.L.
+- 0.6812
+- 0.5086
* 1.0878
* 0.6961
* 0.2393
* 0.1921
t 1.3944
+ 0.0521
* 0.0475
+ 1.3663
"!" 0.0442
+ 0.0152
* 0.0364
* 0.0134
+ 0.0499
* 0.8756
+ 0.0113
* 0.0505
"!" 0.0338
+ 0.0233
* 0.0596
+ 0.0145
* 0.0143
* 0.0097
"hi. 7858
"^25.5148
* 14.4519
+ 4.6062
* 0.7785
* 3.1371
C.V.
12.3
12.2
11.8
2.8
6.1
1.9
72.9
27.6
19.6
66.2
11.0
3.4
24.0
6.2
27.3
68.7
3.5
14.5
35.8
13.3
36.4
25.1
35.9
23.0
25.6
36.4
23.1
26.4
6.6
24.4
47
-------
Table 21. TOTAL AMOUNTS OF DDT, ODD, AND DDE IN THE SOUTH MONTEREY BAY STUDY
AREA BASED ON THE MEAN CONCENTRATIONS AT THE TWO SAMPLE TIMES, AND
EXPECTED AMOUNTS AFFECTED BY THE MEAN OF THE ESTIMATES OF SYSTEM RATES
Amount at Sample Time 1
Amount at Sample Time 2,
3 years later
Expected input for next
year interval
Expected loss for next
year interval
Expected amounts due to Net
change for next year interval
Expected amount translocated
within the system in next
year interval
Expected amount Output to
other systems in next year
interval
Expected amount Decayed
in next time interval
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
DDT
ODD
DDE
TOTAL
Kilograms
128
97
214
439
579
91
232
932
753
23
65
841
58
12
30
100
1274
102
267
1643
463
18
51
532
35
7
23
65
21
2
6
29
f
Pounds
284
213
472
969
1276
200
511
1987
1659
50
143
1852
128
26
66
220
2807
224
588
3619
1020
40
112
1172
77
22
51
150
46
5
14
65
Pounds/ Acre
0.004
0.003
0.007 .
0.014
0.018
0.003
0.007
0.028
0.024
0.001
0.002
0.027
0.0018
0.0004
0.0010
0.0032
0.041
0.003
0.008
0.052
0.015
0.001
0.002
0.018
0.0011
0.0003
0.0007
0.0021
0.0007
0.0001
0.0002
0.0010
48
-------
CIRCULATION OF DDT DERIVATIVES
Figure 8. Composite chart of the translocation of DDT compounds based upon the rates
of change, K, at individual stations in the southern portion of Monterey Bay.
49
-------
DEVELOPMENT OF LABORATORY ASSAY METHODS FOR DETERMINATION OF
DECAY RATE
Of the various preparations tested for the assay of decay rate, the sealed hypovial prepar-
ations described in the Methods section have best met the following desired criteria.
(1) Preparations must be capable of being sealed to prevent loss of the chlorinated hydro-
carbon and its degradation products including CO2- (2) The containers must be readily
sterilized and of materials that prevent contamination by other chlorinated hydrocarbons.
(3) The preparations must be easily manipulated with respect to the establishment of
aerobic and anaerobic conditions. (4) The preparation must be susceptible to replication
both in terms of individual preparations and aliquots from the same preparation.
The most convenient estimate of decay cari be obtained by measurement of the amount
of ^CO2 produced from ring labelled substrate after an interval of time. Knowing the
initial concentrations of substrate the decay to carbon dioxide can be expressed as a deci-
mal fraction of this initial concentration. The decimal fraction is the DCQ... Table 22
presents the results of an assay of DDT to CO2 under aerobic conditions at 10°C. Two
aliquots from each of five preparations at four concentrations of DDT were analysed for
their ^CO2 content. There is no significant difference between the DCQ^ measurements
at the four concentrations of DDT. Therefore, over the range from 100 parts per billion
to 100 parts per million there was neither a stimulation of the decay process nor a satura-
tion of the decay process by substrate. Table 23 presents the results of assays for DCQ-
of DDT, ODD, and DDE. This Table also includes the results of assays in which the effect
of environmental variables on the D was determined.
The Q10 for DCO2 of DDT calculated from the aerobic 10° and 20° assays is 2.50. The
remaining assays where DDT is the substrate were designed to determine the participation
of various physiologically different microbiol populations in the decay process. Aerobic
conditions without additional nutrients gave the maximum DCQ«. The decay process was
inhibited by anaerobiosis, but a rate 27% of the aerobic rate remained. The addition of
nitrate as an additional electron acceptor under anaerobic conditions permitted an in-
crease in the anaerobic rate. The three highest concentrations of nitrate, 5 X 10~1% to
5 X 10~3% were inhibitory but below these concentrations the anaerobic rate becomes
68% of the aerobic rate at 5 X 10"^% sodium nitrate.
The addition of a possible cometabolite, sodium acetate, somewhat removes the inhibi-
tory effect of 5 X 10~1% sodium nitrate probably by its lowering of the nitrate level
through denitrification. However, at none of the levels of sodium acetate tested did the
anaerobic rate reach the level with 5 X 10"5% sodium nitrate alone. The effect of the
addition of cometabolites on decay in the presence of nitrate reducing systems must be
tested at lower concentrations of nitrate.
Sulfate, present in the seawater, was available as an electron acceptor under anaerobic
conditions. Attempts to stimulate sulfate reduction systems by the addition of ethanol
under anaerobic conditions were successful. However, the anaerobic decay of DDT was
not increased over the rate observed with optimum nitrate concentrations and in the
absence of added electron donors such as sodium acetate.
50
-------
Table 22. RESULTS OF A LABORATORY ASSAY OF ANNUAL RATE OF DECAY OF DDT
TO C02, DCQ2, EXPRESSED AS A DECIMAL FRACTION OF THE INITIAL CONCEN-
TRATION OF DDT MAINTAINED AT 10°C. UNDER AEROBIC CONDITIONS.
DDT
100ppm
10 ppm
1 ppm
100 ppb
Prepar
1
1
2
2
3
3
4
4
5
5
1
1
2
2
3
3
4
4
5
5
1
1
2
2
3
3
4
4
5
5
1
1
2
2
3
3
4
4
5
5
DC02
.0046
.0045
.0048
.0042
.0059
.0056
.0050
.0045
.0046
.0053
.0050
.0052
.0058
.0048
.0045
.0056
.0056
.0057
.0051
.0056
.0050
.0059
.0045
.0046
.0062
.0057
.0058
.0052
.0063
.0058
.0057
.0057
.0045
.0047
.0051
.0051
.0055
.0058
.0063
.0053
Means
.00455
.00450
.00575
.00475
.00495
.00510
.00530
.00505
.00565
.00535
.00545
.00455
.00595
.00550
.00605
.00570
.00460
.00510
.00565
.00580
S.D.
t .000071
t .000424
t .000212
t .000354
t .000495
t .000141
t .000707
+ .000778
t .000071
t .000354
t .000636
t .000071
t .000354
t .000424
+ .000354
t .0000
t .000141
t .0000
t .000212
t .000707
Means
.00490
.00529
.00550
.00537
S.D.
t .000544
t .000436
t .000638
t .000542
Mean
.00527
S.D.
t .000570 .
51
-------
Table 23. RESULTS OF LABORATORY ASSAYS OF THE ANNUAL RATES OF DECAY TO CO2,
Dc02- AND THE EFFECT OF ENVIRONMENTAL VARIABLES ON THE PROCESS.
Conditions
Aerobic, 10°C
Aerobic, 20°C
Anaerobic, 10°C
Anaerobic, 10°C
5x lO'HNaNC-s
5x 10'2%NaN03
5x 10'3%NaN03
5 x 10'4% NaN03
5x10-5%NaN03
-
5x 10-6%NaN03
5x 10'7%NaN03
Substrate
DDT 100 ppm
10 ppm
1 ppm
100 ppb
DDT 100 ppm
10 ppm
1 ppm
100 ppb
DDT 100 ppm
10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DC02
.0050
.0053
.0055
.0054
.0100
.0111
.0167
.0154
.0012
.0013
.0015
.0018
.0013
.0016
.0016
.0017
.0018
.0020
.0024
.0024
.0027
.0030
.0036
.0036
.0037
.0034
.0037
.0036
.0025
.0032
.0031
.0032
.0031
Mean
.00529
.01320
Q10 2.50
.00145
.00150
.00183
.00250
.00340
.00360
.00310
.00313
S.D.
t .00023
t .00335
t .00027
t .00017
t .00015
t .00017
t .00035
t .00017
t .00056
t .00006
52
-------
Table 23 CONTINUED. RESULTS OF LABORATORY ASSAYS OF THE ANNUAL RATES OF
DECAY TO C02, Dc02- AND THE EFFECT OF ENVIRONMENTAL
VARIABLES ON THE PROCESS.
Conditions
Anaerobic, 10°C,
5x 1Q-1%NaN03
5 x 10"1% Na Acetate
5 x 10"2% Na Acetate
5 x 10'3% Na Acetate
5x 10'4%Na Acetate
5x 1Q-5%Na Acetate
5x 10'6%Na Acetate
5 x 10'7%Na Acetate
Aerobic, 10°C
5x 10'1%Na Acetate
5 x 10"2% Na Acetate
5x 1 Q-3% N a Acetate
Concentration
DDT 10ppm
1 ppm
100ppb
DDT 10 ppm
1 ppm
100 ppb
. DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
- - - - -
DC02
.0011
.0008
.0008
.0008
.0008
.0010
.0022
.0022
.0023
.0022
.0025
.0024
.0022
.0023
.0023
.0019
.0022
.0023
.0024
.0024
.0024
.0031
.0033
.0031
.0034
.0031
.0027
.0025
.0023
.0023
Mean
.00090
.00087
.00223
.00237
.00227
.00213
.00240
.00317
.00307
.00237
S.D.
+ .00017
t .00012
+ .00006
+ .00015
+ .00006
+ .00021
t .00000
+ .00012
t .00035
t .00012
53
-------
Table 23 CONTINUED. RESULTS OF LABORATORY ASSAYS OF THE ANNUAL RATES OF
DECAY TO C02, DcQ2> AND THE EFFECT OF ENVIRONMENTAL
VARIABLES ON THE PROCESS.
Conditions
5 x 10'4% Na Acetate
5x 10'5%Na Acetate
5x 10'6%Na Acetate
5x 10'7%Na Acetate
Anaerobic, 10°C
5x 10'1%Ethanol
5x 10'2%Ethanol
5x 10'3%Ethanol
5x 10'4%Ethanol
5x 10'5%Ethanol
5x 10'6%Ethanol
5x 10'7%Ethanol
Aerobic, 10°C
Aerobic, 10°C
Concentration
DDT 10ppm
1 ppm
100ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
1 00 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
1 00 ppb
DDT 10 ppm
1 ppm
100 ppb
DDT 10 ppm
1 ppm
100 ppb
ODD 100 ppm
10 ppm
1 ppm
100 ppb
DDE 100 ppm
10 ppm
1 ppm
100 ppb
DC02
.0028
.0030
.0031
.0027
.0030
.0029
.0025
.0027
.0029
.0028
.0030
.0025
.0007
.0005
.0001
.0027
.0028
.0027
.0034
.0031
.0027
.0029
.0030
.0030
.0034
.0032
.0030
.0022
.0023
.0024
.0023
.0022
.0025
.0016
.0015
.0015
.0023
.0030
.0028
.0031
.0041
Mean
.00297
.00287
.00270
.00277
.00043
.00273
.00307
.00297
.00320
.00230
.00233
.00173
.00325
S.D.
t .00015
t .00015
t .00020
+ .00025
t .00031
t .00006
t .00035
t .00006
t .00020
+ .00010
+ .00015
t .00000
t .00058
54
-------
Table 24. RATES OF DECAY TO WATER SOLUBLE COMPOUNDS AND C02 DETERMINED
BY LABORATORY ASSAYS.
Laboratory Assays
DC02
DWS
Laboratory Assays
Corrected by Q^Q
DC02
DWS
•*
Estimations from
Field Data
D
DDT
.00529
.01539
.00600
.01746
.0362
S.D.
t .00023
t .000817
t .00026
+ .00093
ODD
.00173
.00309
:00196
.00351
.0251
S.D.
t .00036
t .00052
t .00041
t .00059
DDE
.00325
.00459
.00369
.00521
.0265
S.D.
+ .00058
+ .00074
+ .00066
t .00084
55
-------
The addition of sodium acetate as an extra electron donor under aerobic conditions was
inhibitory to the aerobic decay process. However, since there was hydrogen sulfate pro-
duced in these preparations the inhibition may have been due to the competition for the
available oxygen and the production of anaerobic conditions.
In summary, decay to CO2 appears to be primarily due to the activity of aerobic micro-
organisms. The process attains the greatest rate where there is no unusual competition for
oxygen. Since the known mechanisms for splitting aromatic rings involve the addition of
oxygen to the aromatic nucleus prior to splitting, these observations are not unexpected.
However, some considerable activity remains under anaerobic conditions even where an
additional oxidizable substrate such as sodium acetate or ethanol is present to remove
any traces of residual oxygen. The results also indicate that nitrate and sulfate may be
acceptable electron acceptors in the oxidation of aromatic compounds under anaerobic
conditions. The mechanisms for anaerobic ring split have not been elucidated. Finally,
The QIQ for the decay process under aerobic conditions presents no surprise as to its
magnitude.
A comparison of the DCQ? f°r DDT, DDD, and DDE reveals a similar relationship to the
total decay rates, D, estimated for South Monterey Bay in that
DDDD,C02 Just ** DDDT > DDDE >DDDD- See Table 24-
For purposes of analysis the process of decay can be divided into a series of steps as follows,
DDT —* LS 0^ WS 1^2 CO2
DDD -^ LS ^S* WS -^£2 CO2
DDE -^ LS ^»> WS
where LS represents lipid soluble degradation products of the starting compound and WS
represents water soluble degradation products of the starting compound.
Water soluble degradation products were measured as water soluble C after high speed
centrifugation of samples from the initial preparations followed by acidification to remove
14C02.
values presented in Table 24 are based on the sum of the C present in water solu-
ble form plus that present as 14CO2. Attempts at determining the amount of lipid soluble
degradation products were unsuccessful. The high levels of the starting compound still
present in the preparations made quantification by gas chromatography difficult. Thin
layer chromatography was more successful but revealed that the sodium hydroxide added
to stop further biological breakdown and to absorb ^4CO2 from the gas phase caused
conversion of a considerable amount of the DDT to DDD.
While laboratory assays of decay rate have revealed rates compatible with the field esti-
mation, it has not been possible to use this approach for full appraisal of the method of
56
-------
estimation of field rates. If we take the difference between the values of Dyys obtained
from laboratory assays and D obtained from field estimations the rates of decay of the
parent compounds to lipid soluble breakdown products, DL§, are .0187 for DDT, .0216
for ODD, and .0213 for DDE under aerobic conditions at 11°C, the mean temperature
of the sediments. It should be noted that although every precaution was taken to ensure
purity of starting materials in laboratory assays, the amounts of decomposition in three
month periods is extremely small and trace contaminants containing labell could have a
large effect upon the results. In addition it must be emphasized that conditions in labora-
tory preparations poorly approximate conditions in the field. Therefore, their value is
more in terms of results obtained by comparisons between preparations rather than com-
parisons between laboratory preparation and field observation.
57
-------
SECTION VI
REFERENCES
Calif. Dept. of Agriculture, 1970, 1971, 1972, 1973. Pesticide Use Report. Data Proces-
sing Center, Calif. Dept. of Agri., Sacramento.
Eberhardt, L. L., R. L. Meeks, and T. J. Peterle. Food chain model for DDT kinetics in a
freshwater marsh. Nature. 230:60-62. 1971.
Gunther, A., and R. C. Blinn. The DDT-type compound as source material in organic syn-
thesis. J. Chem. Educ. 27:654-658. 1950.
Hamaker, J. W. Mathematical prediction of cumulative levels of pesticides in soil. Adv. in
Chem. Sci. 60. Amer. Chem. Soc., Wash., D.C. 1966.
Harrison, H. L., O. L. Louchs, J. W. Mitchell, D. F. Parkhurst, C. R. Tracy, D. G. Watts,
and V. J. Yannacone, Jr. System studies on DDT transport. Sci. 170:503-508. 1970.
Murphy, P. G. Effects of salinity on uptake of DDT, DDE and ODD by fish. Bull. Envir.
Cont. andTox. 5:404-407. 1970.
Papoulis, A. Probability, random variables, and stochastic processes. McGraw-Hill Book
Co.,N. Y. 1965.
Robinson, J. Dynamics of organochlorine insecticides in vertebrates and ecosystems.
Nature. 215:33-35. 1967.
Sokal, R. R. and F. J. Rolf. Biometry. W. H. Freeman and Co., San Francisco. 1969. .
State of California. Monterey basin pilot monitoring project. To be released in 1974.
Woodwell, G. M. Toxic substances and ecological cycles. Sci. Amer. 216:24-31. 1967.
58
-------
APPENDIX A
PROGRAM FOR ESTIMATING SYSTEM RATES BASED ON REAL PAIRED SAMPLE VALUES.
This program for calculation of estimates of rates of input, output, translocation, and decay
was written in Fortran IV level G, and was run on an IBM 360/67. In our experience 112k
was used and the program required approximately 40 seconds per run. A maximum of 60
stations, 7 chemical compounds, and 2 sample times is permitted with the program as written.
The time interval is calculated in the subroutine, LEAPYR, through use of a calendar table
described below. K values are calculated using double precision, and confidence intervals are
estimated through use of a table of "t values."
There are eight cards which precede the data deck. Their formats and content are as follows:
First three cards, FORMAT (1X,13F6.3/13F6.3/4F6.3), contain the table of t values.
The following numbers are punched using the indicated format:
First card, 12.706 4.303 3.182 2.776 2.571 2.447 2.365 2.306 2.262 2.228 2.201
2.1792.160
Second card, 2.145 2.131 2.120 2.110 2.101 2.093 2.086 2.080 2.074 2.069 2.064
2.0602.056
Third card, 2.052 2.048 2.045 2.042.
Fourth card, FORMAT (1214), contains numbers for calculation of time intervals.
The following numbers are punched using the indicated format:
0 31 59 90 120 151 181 212 243 273 304 334.
Fifth card, FORMAT (215), contains the number of stations followed by the number
of chemical compounds in the data set.
Sixth through eighth cards, FORMAT (10A8), contain the names of the chemical
compounds entered, left justified, followed by the word TOTAL, followed by the
concentration level repeated once for each chemical compound. Any remaining
portion of the three cards is left blank. The set of name cards used with the data
analyzed in the present case was as follows:
First Card
DDT DDD DDE TOTAL PPB PPB PPB PPB PERCENT PERCENT
Second Card
PERCENT
The third card was left blank.
The data is organized using FORMAT (1X,I2,2(A4,A2),I2,2(1X,I2),7F7.2). The first variable
is the station number. The next six fields store the location in terms of.latitude and longitude.
The next three variables store the month, day, and year, and the remaining fields store the
measured concentrations of each chemical compound.
An optional subroutine FACTOR may be called by placing a card before the END card with CALL
FACTOR.
59
-------
0001
0002
0003
0004
0005
0006
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
PROGRAM FOR ESTIMATING SYSTEM RATES BASED ON
REAL PAIRED SAMPLE VALUES.
DIMENSION TABLE(30),MONTH(12),ALOC(2,60)6))TOT(10)8),STD(23,8)
1,STE(23,8),CL95(23,8))VAR1(7),VAR2(7),VAR3(7))SUM1(7),SUM2(
27),SUM3(7),SUM4(7),COV1(7),COV2(7)
REAL *4MEAN,MR(7),M(17)8)
REAL*8X(10,60,7),V2(60,7),NAME(23),V1(7)
INTEGER CST(2,60),CDATE(2)60)3)
COMMON X, TABLE,IA,I,K,KD,ID
COMMON/BLK1/NAME,TOT,M,STD)STE,CL95,ALOC,YR,CST,CDATE,MONTH,L1,L
COMMON/BLK2/MR
READ (5 ,45) TABLE
READ (5, 46) MONTH
READ(5,47)IA,ID
CALCULATE INDEXES.
AI NUMBER OF STATIONS CONVERTED TO A REAL NUMBER
IP1 1D+1
IP2 ID + 2
I2TP2 2 * ID + 2
I2TP3 2 * ID + 3
I3TP2 3 * ID + 2
AI=IA
IP1=ID+1
IP2=ID+2
I2TP2=2*ID+2
I2TP3=2*ID+3
I3TP2=3*ID+2
CLEAR X ARRAY.
DO 1 1=1,10
DO 1 J=1,IA
DO 1 K=1,IP1
1 X(IJ,K)=0.0
0020
0021
C
C
C
C
C
WRITE (6,50)
READ IN DERIVATIVE NAMES AND CONCENTRATION LEVEL ON UP TO 3 CAR
READ (5,48) NAME
READ IN DATA.
60
-------
0022 DO 2 1=1,2
C
0023 DO2J=1,IA
0024 2 READ (5,49) CST(I1J)1(ALOC(I,J,L))L=1)6))(CDATE(I,J)L)1L=1,3),(X(I
1,J,K),K=1,ID)
C
C
C COMPUTE TOTAL OF EACH STATION.
C
0025 DO 3 1=1,2
C
0026 DO 3 J = 1,IA
C
0027 DO 3 L=1,ID
0028 3 X(I,J,IP1)=OC(I,J,L)+X(I,J,IP1)
C
C
C WRITE HEADING OF FIRST TWO PAGES.
C
0029 DO 5 1=1,2
0030 L=I
0031 WRITE (6,51)1
0032 WRITE (6,53) (NAME(N),N=1,IP1)
0033 WRITE (6,52) (NAME(N),N=IP2,I2TP2)
0034 WRITE (6,54)
C
0035 DO4K=1,IP1
0036 CALL STDEV (TOTAL,MEAN,SD,SE,CL)
0037 TOT(I,K)=TOTAL
0038 M(I,K)=MEAN
0039 STD(I,K)=SD
0040 STE(I,K)=SE
0041 4 CL95(I,K)=CL
C
C L1=NUMBER OF SETS COMPUTED.
C
C WRITE FIRST TWO PAGES.
0042 L1=IP1
0043 CALL PRINT
0044 WRITE (6,53) (NAME(N),N=1,IP1)
0045 WRITE (6,52) (NAME(N),N=IP2,I2TP2)
0046 WRITE (6,54)
0047 CALL PRINT2
0048 5 CONTINUE
C
C
C COMPUTE PERCENTS.
C
61
-------
0049 DO 8 1=3,4
0050 L1=ID
0051 L=I-2
0052 WRITE (6,51) L
0053 WRITE (6,53) (NAME(N),N=1,ID)
0054 WRITE (6,52) (NAME(N),N=I2TP3,I3TP2)
0055 WRITE (6,54)
C
0056 DO 6 K=1,ID
C
0057 DO6J=1,IA
0058 6X(I,J,K)=X(L)J,K)/X(L)J,IP1)*100.
C
C
0059 DO7K=1,3
0060 CALL STDEV (TOTAL,MEAN,SD,SE,CL)
0061 TOT(I,K)=TOTAL
0062 M(I,K)=MEAN
0063 STD(I,K)=SD
0064 STE(I,K)=SE
0065 7 CL95(I,K)=CL
C
0066 CALL PRINT
0067 WRITE (6,53) (NAME(N),N=1,ID)
0068 WRITE (6,52) (NAME(N),N=I2TP3,I3TP2)
0069 WRITE (6,54)
0070 CALL PRINT2
0071 8 CONTINUE
C
C
0072 DO10J=1,IA
C
0073 DO10L=1,IA
0074 IF (CST(1,J).EQ.CST(2,L)) GO TO 9
0075 GO TO 10
0076 9CALLLEAPYR(J)
C
0077 DO10K=1,ID
0078 X(5,J,K)=YR
0079 10 CONTINUE
C
C
C CALCULATE TOTAL AND MEAN OF N.
C
62
-------
OQ80
0081
0082
0083
0084
0099
0100
0101
0102
0103
0104
0105
DO 12K=1,ID
TOT(5,K)=0.
DO 11 J=1,IA
11 TOT(5,K)=TOT(5)K)+X(5,J,K)
12M(5,K)=TOT(5,K)/AI
0085
0086
0087
0088
0089
0090
0091
0092
0093
0094
0095
0096
0097
0098
C
C
C
C
C
C
C
C
C
C
C
C
C
DO 14K=1,ID
V=0.0
DO 13J=1,IA
13 V=(M(5,K)-X(5,J,K))**2+V
STD( 5 ,K)=SQRT( V/(AI-1 .0))
14 CALL STDEV2 (STD(5,K),STE(5,K))CL95(5,K))
CALCULATE K VALUES.
DO 15K=1,ID
SUM1(K)=0.0
DO 15 J=1,IA
IF (X(1J,K).EQ.O) X(1,J,K)=.004
IF (X(2J,K).EQ.O) X(2,J,K)=.004
V=(DLOG10(X(2J,K))-DLOG10(X(1J,K)))/(X(5J,K))
V2(J,K)=10.**V-1.0
15 SUM1(K)=SUM1(K)+V2(J,K)
SORT RVALUES.
C
C
C
C
DO 17 K=1,ID
DO 17J = 1,IA
IF(V2(J,K).GT.O)GOTO 16
X(7,J,K)=V2(J,K)
GO TO 17
16X(6,J,K)+V2(J,K)
17X(8,J,K)=X(7)J)K)+X(6IJ,K)
CALCULATE K-NET.
63
-------
0106
DO19K=1,ID
0107
0108
0109
0110
0111
0112
0113
0114
C
C
C
C
C
C
DO 19J=1,IA
V=X(8,J,K)-SUM1(K)/AI
IF(V.GT.O)GOTO18
X(10,J,K)=V
GO TO 19
18 X(9,J,K)=V
19 CONTINUE
COMPUTE SUM AND MEAN FOR K VALUES.
DO 21 K=1,ID
0115
0116
0117
0118
0119
0120
0121
0122
0123
0124
0125
0127
0128
0129
0130
C
C
C
C
C
C
C
C
DO 21 1=6,10
V=0.0
DO20J=1,IA
20 V=V+X(IJ,K)
TOT(I,K)=V
21 M(I,K)=V/AI
CALCUALTE STANDARD DEVIATION, STANDARD ERROR, AND 95% CONFIDE:
LIMITS OF K VALUES.
DO22K=1,ID
SUM 1 (K)=0.0
SUM 2 (K)=0.0
SUM 3 (K)=0.0
22 SUM 4 (K)=0.0
DO 23 J=1,IA
V2(J,K)=DLOG(X(2J,K))-DLOG(X(1,J,K))
SUM2(K)=V2(J,K)+SUM2(K)
SUM 3 (K)=(DLOG(X(1J,K))-ALOG(M(1,K)))**2+SUM3(K)
64
-------
0131 23 SUM 4 (K)=(DLOG(X(2,J,K))-ALOG(M(2,K)))**2+SUM4(K)
C
C
0132 DO24K=1,ID
0133 VARl(K)<43429/M(l,K))**2*SUM3(K)/(AI-1.0)+(-.43429/M(2,K))**2
1*SUM4(K)/(AI-1.0)
0134 24 V1(K)=SUM2(K)/AI
0135 DO25K=1,ID
0136 VAR2(K)=((1.0/M(5,K))**2*VARl(K))+(-Vl(K)/(M(5,K)*»2))**2*STD(5,K)
1**2
0137 VAR2(K)=10.0**VAR2(K)
0138 STD(8,K)=SQRT(VAR2(K))
0139 25 CALL STDEV2 (STD(8,K),STE(8>K),CL95(8,K))
C
C
C
C CALCULATE THE DISTRIBUTION OF VARIANCE BETWEEN +K AND -K
C
0142 DO30K=1,ID
0143 V=0.0
C
0144 DO27J=1,IA
0145 IF (X(6,J,K)) 27,27,26
0146 26 V=(X(6,J,K)-M(8,K))**2+V
0147 27 CONTINUE
C
C
0148 V=V/(AI-1.0)
0149 W=0.0
C
C
0150 DO29J=1,IA
0151 IF (X(7J,K» 28,29,29
0152 28W=(X(7,J,K)-M(8,K))**2+W
0153 29 CONTINUE
C
0154 W=W/(AI-1.0)
0155 U=V+W
0156 V=STD(8,K)**2*(V/U)**2
65
-------
0157 STD(6,K)=SQRT(V)
0158 W=STD(8,K)**2*(W/U)**2
0159 STD(7,K)=SQRT(W)
0160 CALL STDEV2(STD(6,K),STE(6,K),CL95(6,K))
0160 30 CALL STDEV2 (STD(7,K),STE(7>K),CL95(7,K))
C
C
C CALCULATION OF STANDARD DEVIATION K-NET AND ITS DISTRIBUTION.
C
C
C
0161 DO35K=1,ID
0161 V=0.0
0162 - W=0.0
0163 DO34J=1,IA
0164 IF(X(9,J,K)) 32,32,31
0165 31 V=V+(X(9,J,K)**2)
0166 32 IF(X(10,J,K)) 33,34,34
0167 33 W=W+(X(10,J,K)**2)
0168 34 CONTINUE
C
C
0169 V=V/(AI-1.0)
0170 W=W/(AI-1.0)
0171 STD(9,K)=SQRT((V/(V+W))**2*(STD(8,K)**2))
0172 CALL STDEV2(STD(9,K),STE(9,K),CL95(9,K))
0173 STD(10,K)=SQRT((W/(V+W))**2*(STD(8,K)**2))
0174 35 CALLSTDEV2(STD(10,K),STE(10,K),CL95(10,K))
C
C
C
0175 CALLPRINT3
C
CALCULATE 0 AND ITS STANDARD DEVIATION
C
0176 DO41K=1,ID
0177 M(l 1,K)=(M(9,K)/M(6,K))*M(7,K)
0178 STD(11,K)=SQRT(STD(7,K)**2*((M(9,K)/M(6,K))**2))
0179 CALL STDEV2(STD(11,K),STE(11,K),CL95(11,K))
C
C
C
C
C CALCULATION OF D
C
66
-------
c
c
0192 M(12,K)=M(7,K)-M(11,K)
C
C
C CALCULATION OF STANDARD DEVIATION OF D
C
0193 STD(12,K)=SQRT(STD(7,K)**2*(1.-M(9,K)/M(6,K))**2)
0194 CALL STDEV2 (STD(12,K))STE(12,K),CL95(12,K))
C
C
C CALCULATE TL.
C
C
0196 M(13,K)=-1.0*(1.0/M(12,K))
0197 STD(13,K)=SQRT(STD(12,K)**2*(1.0/M(12,K)**2)**2)
0198 41 CALL STDEV2 (STD(13,K),STE(13,K),CL95(13)K))
C
C
C CALCULATE TR.
C
0199 DO42K=1,ID
0200 M(14,K)=-1.0*(1.0/M(7,K))
0201 STD(14,K)=SQRT(STD(7,K)**2*(1.0/M(7,K)**2)**2)
0202 42 CALL STDEV2 (STD(14,K),STE(14,K))CL95(14,K))
C
C
0203 DO44K=1,ID
0204 WRITE (6,55) NAME (K)
0205 WRITE (6,56) NAME (K),M(8,K),STD(8,K),STE(8,K),CL95(8,K)
0206 WRITE (6,57) NAME (K),M(9,K),STD(9,K),STE(9,K),CL95(9,K)
0207 WRITE (6,58)
0208 WRITE (6,59) NAME (K),M(10,K),STD(10,K),STE(10,K),CL95(10,K)
0209 WRITE (6,60)
0210 WRITE (6,61) NAME (K),M(6,K),STD(6,K))STE(6,K),CL95(6,K)
0211 WRITE (6,62) NAME (K),M(7,K),STD(7,K),STE(7,K),CL95(7,K)
0212 WRITE (6,63) NAME (K),M(11,K),STD(11,K),STE(11,K),CL95(11,K)
0213 WRITE (6,64) NAME (K),M(12)K),STD(12,K),STE(12,K),CL95(12,K)
0214 WRITE (6,65) NAME (K),M(13,K),STD(13,K),STE(13,K),CL95(13,K)
0215 WRITE (6,66)
0216 WRITE (6,65) NAME (K),M(14,K),STD(14,K),STE(14,K),CL95(14,K)
0217 WRITE (6,67)
67
-------
0218 DO43L=1,3
0219 43 WRITE (6,54)
C
0220 WRITE (6,68)
0221 WRITE (6,69) M(5,K)
0222 WRITE (6,70) NAME(K),M(2,K)>l(l,K),M(6,K)Al(9,K)JVl(10,K),M(ll,K),M
K12.K)
0223 44 CONTINUE
C
0224 CALL FACTOR
0225 STOP
C
0226 45 FORMAT (1X,13F6.3/13F6.3/4F6.3)
0227 46 FORMAT (1214)
0228 47 FORMAT (215)
0229 48 FORMAT (10A8)
0230 49 FORMAT (1X,I2,2(2A4,A2),I2,2(1X,I2),7F7.2)
0231 50 FORMAT('l')
0232 51 FORMAT('1>,/2X,I1)/3X,'STATION',3X,'LATITUDE',3X),13X,'=NET'13X,A81'=>13X,4F11.4/)
0238 57 FORMAT(2X,'MEANOF + (K-NET) = T',5X,A8,'=>)3X14F11.4)
0239 58 FORMAT(27X,T/)
0240 59 FORMAT(2X,'MEAN OF - ( K - NET ) = T',5X,A8,'=',3X,4F11.4)
0241 60 FORMAT(27X,'OV)
0242 61 FORMAT(2X,'MEANOF + K')llX,'=r,5X,A8,'=',3X)4F11.4//)
0243 62 FORMAT(2X,,3X,4F11.4//)
0244 63 FORMAT(26X,'O',5X,A8,'=',3X,4F11.4/)
0245 64 FORMAT(26X,'D',5X>A8,'=',3X,4F11.4/)
0246 65 FORMAT(26X>'T',5X1A8,'=',3X,4F11.4)
0247 66 FORMAT(27X,'LV)
0248 67 FORMAT(27X,'R')
0249 68 FORMAT(13X,'MEAN C',6X,'MEAN C>,16X,T,10X,>6XI<-1>4X,T,6X>
-<-l,5X,'O',5X,'->,5Xf€D>,9X,'NVl9X,<2MlX,
-------
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
C
C
SUBROUTINE PRINT
DIMENSION TABLE(30),MONTH(12),ALOC(2,60,6),TOT(10,8),STD(23,8)
1,STE(23)8),CL95(23,8)
REAL*4MEAN,M(17,8)
REAL *8X(10,60,7),NAME(23)
INTEGER CST(2)60))CDATE(2,60,3)
COMMON X,TABLE,IA,I)K,KD,ID
COMMON/BLKl/NAME)TOT)M)STD)STE)CL95)ALOC)YR)CST,CDATE,MONTH)Ll,L
DO 1 J = 1,IA
1 WRITE (6,3) CST(L,J),(ALOC(L,J,K),K=1,6),(CDATE(L,J,K),K=1,3),(X(I
1J,K),K=1,L1)
SKIP TO BOTTOM OF PAGE
N=(68-(IA+6))/2
0011
0012
0013
0014
0015
0016
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
C
C
C
C
DO2J=1,N
2 WRITE (6,4)
RETURN
3 FORMAT (5X,I2,5X,2A4,A2)2X,2A4,A2)2X,I2,2('-',I2),8F11.2)
4 FORMAT (/)
END
SUBROUTINE PRINT2
DIMENSION TABLE(30),MONTH(12),ALOC(2,60,6),TOT(10,8),STD(23,8)
1,STE(23,8),CL95(23,8)
REAL *4MEAN,M(17,8)
REAL *8X(10,60,7),NAME(23)
INTEGER CST(2)60),CDATE(2,60,3)
COMMON X,TABLE,IA,I,K>KD,ID
COMMON/BLK1/NAME)TOT,M,STD)STE,CL95,ALOC,YR,CST,CDATE,MONTH,L1(L
WRITE (6,1) (TOT(IJ)J=1,L1)
WRITE (6,2) (M(IJ),J=1,L1)
WRITE (6,3) (STD(I,J),J=1,L1)
WRITE (6,4) (STE(IJ)J=1,L1)
WRITE (6,5) (CL95(I,J),J = 1,L1)
RETURN
0014
0015
0016
0017
0018
0019
C
C
1 FORMAT (34X,'TOTALS',6X,7F 10.4)
2 FORMAT (/34X,'MEAN',8X,7F10.4)
3 FORMAT (/34X,'S.D.',8X,7F10.4)
4 FORMAT (/34X,'S.E.',8X,7F10.4)
5 FORMAT (/34X,'95% CL',6X,7F10.4)
END
69
-------
0001 SUBROUTINE PRINTS
0002 DIMENSION TABLE(30),MONTH(12)ALOC(2,60,6),TOT(10,8),STD(23,8)
1,STE(23,8),CL95(23,8)
0003 REAL *8X(10,60,7)
0004 REAL * 8NAME(2 3)
0005 REAL *4MEAN,M(17,8)
0006 INTEGER CST(2,60),CDATE(2,60)3)
0007 COMMON X,TABLE,IA,I,K,KD,ID
0008 COMMON /BLK1/ NAME,TOT,M)STD)STE,CL95,ALOC,YR,CST>CDATE(MONTH,L1,L
C
0009 DO 2 K=1,ID
0010 WRITE (6,3)
0011 WRITE (6,4)
0012 WRITE (6,5) NAME(K),NAME(K)
0013 WRITE (6,6)
0014 WRITE (6,8)
C
0015 DO1J = 1,IA
0016 1 WRITE (6,7) CST(1J),X(2,J,K),X(1 J,K),X(5,J,K),(X(IXJ,K),IX=6,10
1)
C
0017 WRITE (6,8)
0018 WRITE (6,17) TOT(2,K),TOT(1,K),TOT(5,1),(TOT(L,K),L=6,10)
0019 WRITE (6,16)
0020 WRITE (6,14) NAME(K)
0021 WRITE (6,17) M(2,K),M(1,K),M(5,1),(M(L,K),L=6,10)
0022 WRITE (6,9)
0023 WRITE (6,14) NAME(K)
0024 WRITE (6,17) STD(2,K),STD(1,K),STD(5,1),(STD(L,K))L=6,10)
0025 WRITE (6,10)
0026 WRITE(6,13)NAME(K)
0027 WRITE (6,17) STE(2,K),STE(1,K),STE(5,1),(STE(L,K))L=6,10)
0028 WRITE (6,11)
0029 WRITE (6,13) NAME(K)
0030 WRITE (6,17) CL95(2,K),CL95(1,K),CL95(5,1),(CL95(L,K),L=6110)
0031 WRITE (6,12)
0032 WRITE(6,15)NAME(K)
0033 2 CONTINUE
C
0034 RETURN
C
0035 3 FORMAT (T.IX.'STATION')
0036 4 FORMAT (12X,'C',9X,,9X,I1',52X,,10X,
-------
0041
0042
0043
0044
0045
0046
0047
0048
0049
0050
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
C
C
9 FORMAT (V,94X,'MEANS')
10 FORMAT (V,94X,'S.D.')
11 FORMAT (<+>,94X,)
12 FORMAT (V,94X,'95% CONFIDENCE LIMITS')
13 FORMAT ('+', 99X.A8)
14 FORMAT (V,102X,A8)
15 FORMAT (V,116X,A8)
16 FORMAT ('+'>94X,'TOTALS')
17 FORMAT(/9X,5F10.4,3F11.4)
END
SUBROUTINE LEAPYR (J)
DIMENSION TABLE(30), MONTH(12), ALOC(2,60,6), TOT(10,8), STD(23,8)
1,STE(23)8),CL95(23,8)
REAL *4MEAN,M(17,8)
REAL *8X( 10,60,7)
REAL *8NAME(23)
INTEGER TOT,YR1,YR2,DA1,DA2,DAYS
INTEGER CST(2,60),CDATE(2,60,3)
COMMON/BLK1/NAME,TOT,M,STD,STE,CL95,ALOC,YR,CST,CDATE,MONTH,L1,L
COMMON X,TABLE,IA,I,K,KD,ID
DAYS=0
NT=0
MO1=CDATE(1,J,1)
DA1=CDATE(1,J,2)
YR1=CDATE(1J,3)
DA2=CDATE(2,L,2)
YR2=CDATE(2,L,3)
MO2=CDATE(2,L,1)
AMO=MO1
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
DO4I=YR1,YR2
A=I
LEAP=0
IZ=A/4.
Z=IZ
Z=Z*4.
IF(I.EQ.YR1)GOTO 1
GO TO 2
1 DAYS=365-(MONTH(MO1)+DA1)
IF (Z.EQ.A.AND.AMO.LT.3.) LEAP=1
GOTO 3
2 IF (Z.EQ.A) LEAP=1
3 NT=DAYS+LEAP+NT
4DAYS=365
71
-------
0033 IF(LEAP.EQ.1)GOTO5
0034 GO TO 6
0035 5 IF (MO2.LT.3) NT=NT-1
0036 6 YR=NT-365+MONTH(MO2)+DA2
0037 YR=YR/365.
0038 RETURN
0039 END
C
C
0001 SUBROUTINE TDIST (T)
0002 REAL *8X(10,60,7)
0003 DIMENSION TABLE(30)
0004 COMMON X.TABLE.IA.I.K.KD.ID
0005 11=IA-1
0006 AI=I1
0007 IF (II) 1,1,2
0008 1 WRITE (6,11) I
0009 GO TO 10
0010 2IF(I1.LT.31)GOTO9
0011 IF(I1.LT.41)GOTO3
0012 GO TO 4
0013 3 TINT=((2.042-2.021)/10.)*(AI-30.)
0014 T=TINT+2.042
0015 GO TO 10
0016 4IF(I1.LT.61)GOTO5
0017 GO TO 6
0018 5 TINT=((2.021-2.000)/20.)*(AMO.)
0019 T=TINT+2.021
0020 GO TO 10
0021 6IF(I1.LT.121)GOTO7
0022 GO TO 8
0023 7 TINT=((2.000-1.980)/40.)*(AI-60.)
0024 T=TINT+2.000
0025 GO TO 10
0026 8 T=1.960
0027 GO TO 10
0028 9T=TABLE(I1)
0029 10 RETURN
C
0030 11 FORMAT ('I'.'I IN T TABLE =',I3)
0031 END
C
C
72
-------
0001
0002
0003
0004
0005
0006
SUBROUTINE STDEV (SUMX,XBAR,STD,STE,CL$)
REAL *8X(10,60,7)
DIMENSION TABLE(30)
COMMON X,TABLE,IA,I,K,KD,ID
DEV=0.
SUMX=0.
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
C
C
DO 1 J=1,IA
1 SUMX=SUMX+X(I,J,K)
AI=IA
XBAR=SUMX/AI
DO 2 J=1,IA
DEV=(XBAR-X(I ,J ,K)) * * 2+DEV
2 CONTINUE
STD=SQRT(DEV/(AI-1.))
STE=STD/SQRT(AI)
CALL TDIST (T)
CL$=T*STE
END
SUBROUTINE STDEV2 (STD,STE,CL$)
REAL *8X( 10,60,7)
DIMENSION TABLE(30)
COMMON X,TABLE,IA,I,K,KD,ID
AI=IA
STE=STD/SQRT(AI)
CALL TDIST (T)
CL$=T*STE
RETURN
END
73
-------
c
c
0001
0002
0003
0004
0005
0006
0007
C
c
c
c
0008
SUBROUTINE FOR PROGRAM FOR ESTIMATING SYSTEM RATES
BASED ON REAL PAIRED SAMPLE VALUES.
SUBROUTINE FACTOR
DIMENSION TABLE (30), MONTH(12), ALOC(2,60,6), TOT(10,8), STD(23,8)
1,STE(23,8), CL95(23,8), VAR1(7), VAR2(7), VAR3(7), SUM1(7), SUM2(
27), SUM3(7), SUM4(7), COV1(7), COV2(7)
REAL *4MEAN,M(17,8),MR(7)
REAL*8X(10,60,7),V2(60,7);NAME(23),V1(7)
INTEGER CST(2,60),CDATE(2,60,3)
COMMON X,TABLE,IA,I,K,KD,ID
COMMON/BLK1/NAME,TOT,M,STD,STE,CL95,ALOC,YR,CST,CDATE,MONTH,L1,L
COMMON /BLK2/ MR
CALCULATE CORRECTION FACTOR FOR STANDARD DEVIATION
AI=IA
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
C
C
DO 14 K=1,ID
IF (M(8,K)) 1,4,5
1 JX=(AI*M(7,K)-AI*M(10,K))/M(8,K)
VJ=JX
V=(((AI*M(10,K))/VJ)**2)*VJ
IF ((M(8,K)+((AI*M(10,K))/VJ))*VJ-(AI*M(7,K))) 3,3,2
2 V=V+(AI*M(7,K)-VJ*(M(8,K)+((AI*M(10,K))/VJ))-M(8,K))**2
V=V+(((AI*M(6,K))/(AI-VJ-1.0))-M(8,K))**2*(AI-VJ-1.0)
GO TO 8
3V=V+(((AI*M(6,K))/(AI-VJ))-M(8,K))**2*(AI-VJ)
GO TO 8
4 JX=AI/2.0
VJ=JX
V=((AI*M(6,K)/VJ)**2)*VJ
V=V+((AI*M(7,K)/VJ)**2*VJ
GO TO 8
5 JX=(AI*M(6,K)-AI*M(9,K))/M(8,K)
VJ=JX -
V=(((AI*M(9,K))/VJ)**2)*VJ
IF ((M(8,K)+((AI*M(9,K))/VJ))*VJ-(AI*M(6,K))) 6,7,7
6 V=V+(AI*M(6,K)-VJ*(M(8,K)+((AI*M(9,K))/VJ))-M(8,K))**2
V=V+(((AI*M(7,K))/(AI-VJ-1.0))-M(8,K))**2*(AI-VJ-1.0)
GO TO 8
7 V=V+(((AI*M(7,K))/(AI-VJ))-M(8,K))**2*(AI-VJ)
8 V=V/(AI-1.0)
W=0.0
74
-------
f)035 DO9J=1,IA
0036 9 W=W+(X(8J,K)-M(8,K))**2
C
c
0037 W=W/(AI-1.0)
0038 C=((W-V)/W)**2
C
C CALCULATE CORRECTED STD,6,7,AND 8
C STD(15,K) IS CORRECTED STD(6,K)
C
0039 STD(15,K)=SQRT(C*STD(6,K)**2)
0040 CALL STDEV2 (STD(15,K))STE(15,K))CL95(15,K))
C
C STD(16,K)IS CORRECTED STE(7,K)
C
0041 STD(16,K)=SQRT(C*STD(7,K)**2)
0042 CALL STDEV2 (STD(16,K),STE(16,K))CL95(16,K))
C
C STD(17,K)IS CORRECTED STD(8,K)
C
0043 STD(17,K)=SQRT(C*STD(8,K)**2)
0044 CALL STDEV2 (STD(17,K),STE(17,K),CL95(17,K))
C
C CALCULATE CORRECTED STD(9,K) AND STD(10,K)
C
0045 V=0.0
C
C
0046 DO11J=1,IA
C
0047 IF (X(9J,K)) 11,11,10
0048 10 V=V+X(9J,K)**2
0049 11 CONTINUE
C
C
C
0050 V=V/(AI-1.0)
0051 W=0.0
C
C
0052 DO13J=1,IA
C
0053 IF(X(10J,K))12,13,13
0054 12 W=W+(X(10J,K))**2
0055 13 CONTINUE
C
C
C 75
-------
0056 W=W/(AI-1.0)
C
C STD(18,K) IS CORRECTED STD(9,K)
C
0057 STD(18,K)=SQRT((V/(V+W))**2*C*(STD(8,K)**2))
0058 CALL STDEV2 (STD(18,K))STE(18,K),CL95(18,K»
C
C
C STDU9.K) IS CORRECTED STD(IO.K)
C
0059 STD(19,K)=SQRT((W/(V+W))**2*C*(STD(8,K)**2))
0060 CALL STDEV2 (STD(19,K),STE(19,K),CL95(19IK))
C
C CALCULATE CORRECTED STD(20,K) CORRECTED STD(11,K)
C
0061 STD(20,K)=SQRT(STD(16,K)**2*(M(9,K)/M(6,K))**2)
C
C
0062 CALL STDEV2 (STD(20)K))STE(20,K),CL95(20,K))
C CALCULATE STD(21,K) CORRECTED STD(12,K)
0063 STD(21,K)=SQRT(STD(16,K)**2*(1.0-M(9,K)/M(6,K))**2)
0064 CALL STDEV2 (STD(21,K),STE(21)K),CL95(21,K))
C
C CALCULATE STD(22,K) CORRECTED STD(13,K)
C
0065 STD(22,K)=SQRT(STD(21,K)**2*(1.0/M(12,K)**2)**2)
0066 CALL STDEV2 (STD(22,K),STE(22,K))CL95(22>K))
C
C CALCULATE STD(23,K) CORRECTED STD(14,K)
C
0067 STD(23,K)=SQRT(STD(16,K)**2'(1.0/M(7,K)**2)**2)
0068 CALL STDEV2 (STD(23,K),STE(23,K))CL95(23)K))
C
C
0069 14 CONTINUE
C
C
0070 DO 15 K=1,ID
C
0071 WRITE (6,16) NAME(K)
0072 WRITE (6,17) NAME(K),M(8,K),STD(17,K),STE(17,K)>CL95(17,K)
0073 WRITE (6,18) NAME(K),M(9,K),STD(18,K),STE(18,K),CL95(18,K)
0074 WRITE (6,19)
0075 WRITE (6,20) NAME(K),M(10,K),STD(19,K),STE(19,K),CL95(19,K)
76
-------
007$ WRITE (6,21)
Ot)77 WRITE (6,22) NAME(K),M(6,K),STD(15,K),STE(15,K),CL95(15,K)
0078 WRITE (6,23) NAME(K),M(7IK),STD(16IK),STE(16,K),CL95(16IK)
0079 WRITE (6,24) NAME(K),M(11,K),STD(20,K))STE(20,K),CL95(20,K)
0080 WRITE (6,25) NAME(K),M(12,K))STD(21)K),STE(21,K),CL95(21,K)
0081 WRITE (6,26) NAME(K),M(13,K))STD(22,K),STE(22)K),CL95(22,K)
0082 WRITE (6,27)
0083 WRITE (6,26) NAME(K))M(14)K),STD(23,K),STE(23,K),CL95(23,K)
0084 WRITE (6,28)
0085 15 CONTINUE
C
C
0086 16 FORMAT ('l',40X,'CORRECTED STANDARD DEVIATIONSV/,2X,'RATES OF CHA
INGE FOR ',A8,30X,'S.D.',7X,'S.E.',4X,'95% LIMIT'//)
0087 17 FORMAT (2X,'MEAN OF K',13X,'= NET'.SX.AS.^'.S
0088 18 FORMAT (2X/MEAN OF + ( K - NET ) = T'.SX.AS.^
0089 19 FORMAT (27X.T/)
0090 20 FORMAT (2X,'MEAN OF - ( K - NET ) = T'.SX.AS.
0091 21 FORMAT (27X/OV)
0092 22 FORMAT (2X,'MEAN OF + K',1 1X,'= I',5X,A8>'=>,3X,4F11.4//)
0093 23 FORMAT (2X,'MEAN OF- K',11X,'= O + D'.lX.AS.
0094 24 FORMAT (26X,'O',5X,A8,'=>,3X,4F11.4/)
0095 25 FORMAT (26X,'D>,5X)A8,'=')3X,4F11.4/>
0096 26 FORMAT (26X,'T',5X,A8,'=',3X)4F11.4)
0097 27 FORMAT (27X/LY)
0098 28 FORMAT (27X,'R()
C
0099 RETURN
0100 END
77
-------
APPENDIX B
PROGRAM FOR ESTIMATING SYSTEM RATES BASED ON SAMPLE VALUES PAIRED TO MEANWALUES
This program for calculation of estimates of input, output, translocation, and decay was written in
Fortran IV level G, and was run on an IBM 360/67. In our experience 112k was used and the pro-
gram required approximately 40 seconds per run. A maximum of 60 stations, 7 chemical compounds,
and 2 sample times is permitted with the program as written.
The time interval is calculated in the subroutine, NCOMP, which calls the subroutine, LEAPYR. K
values are calculated using double precision, and confidence intervals are estimated through use of
a table of "t values."
There are eight cards which precede the data deck. Their formats and content are as follows:
First four cards, as in preceding program.
Fifth card, Format (315), contains the number of stations at time one, followed by the
number of stations at time two, followed by the number of chemical com-
pounds in the data set.
Sixth through eighth cards, Format (10A8), as in preceding program.
The data is organized as in the preceding program but is sorted chronologically.
An optional subroutine FACTOR may be called by placing a card before the END card with CALL
FACTOR.
C PROGRAM FOR ESTIMATING SYSTEM RATES BASED ON
C SAMPLE VALUES PAIRED TO MEAN VALUES.
C
C
0001 DIMENSION TABLE(30),MONTH(12),ALOC(2,60,6),TOT(10,8),STD(23,8).
1, STE(23,8),CL95(23,8),VAR1(7),VAR2(7),VAR3(7),SUM1(7),SUM 2(
27), SUM3(7), SUM4(7), COV1(7), COV2(7), IA(2), AI(2)
0002 REAL *4MEAN,M(17,8),MR(7)
0003 REAL *8X(10,60,7),V2(60,7),NAME(23),V1(7).
0004 INTEGER CST(2,60),CDATE(2,60,3)
0005 COMMON X,TABLE,IA,IB,I,K,KD,ID
0006 COMMON /BLK2/MR
0006 1 FORMAT (1X,13F6.3/13F6.3/4F6.3)
0007 READ (5,1) TABLE
0008 2 FORMAT (1214)
0009 READ (5,2) MONTH
0010 READ (5,3) IA(1)'IA(2),ID
0011 3 FORMAT (315)
0012 C
78
-------
C CALCULATE INDEXES.
C AI NUMBER OF STATIONS CONVERTED TO A REAL NUMBER.
C AI3 AI(1) + AI(2)
C IA3 IA(1) + IA(2)
C IP1 ID + 1
C IP2 ID + 2
C I2TP2 2 * ID + 2
C I2TP3 2 * ID + 3
C I3TP2 3 * ID + 2
C J2T IA(1) + IA(2)
0013 AI(1)=IA(1)
0014 AI(2)=IA(2)
0015 AI3=AI(1)+AI(2)
0016 IA3=IA(1)+IA(2)
0017 IP1=ID+1
0018 IP2=ID+2
0019 I2TP2=2*ID+2
0020 I2TP3=2*ID+3
0021 I3TP2=3*ID+2
0022 J2T=IA(1)+IA(2)
C
C CLEAR X ARRAY.
C
0023 DO 4 1=1,10
C
0024 DO4J = 1J2T
C
0025 DO4K=1,IP1
0026 4X(IJ,K)=0.0
C
0027 WRITE (6,9)
C
C READ IN DERIVATIVE NAMES AND CONCENTRATION LEVEL ON UP TO 3 CARDS.
0028 READ (5,5) NAME
0029 5 FORMAT (10A8)
C
C READ IN DATA.
0030 6 FORMAT (1X)I2,2(2A4,A2),I2,2(1X,I2),7F7.2)
C
0031 DO 7 1=1,2
0032 IB=IA(I)
C
0033 DO7J=1,IB
0034 7 READ (5,6) CST(I,J),(ALOC(I J,L),L=1,6),(CDATE(I J,L),L=1,3),(X(II
1J,K),K=1,ID)
C
C
C COMPUTE TOTAL OF EACH STATION.
79
-------
0035 DO 8 1=1,2
0036 IB=IA(I)
C
0037 DO8J = 1,IB
C
0038 DO 8 L=1,ID
0039 8 X(IJ,IP1)=X(IJ,L)+X(IJ,IP1)
C
C
C WRITE HEADING OF FIRST TWO PAGES.
C
0040 DO 15 1=1,2
0041 IB=IA(I)
0042 L=I
0043 9 FORMAT ('!')
0044 10 FORMAT (<1','C'/2X,I1,/3X,'STATION',3X,'LATITUDE')3X,'LONGITUDE>,
1 5X/DATE')
0045 WRITE (6,10) I
0046 11 FORMAT (48X,8(3X.A8))
0047 12 FORMAT (V,47X,8(3X,A8))
0048 WRITE (6,12) (NAME(N),N=1,IP1)
0049 13 FORMAT (/)
0050 WRITE (6,11) (NAME(N),N=IP2,I2TP2)
0051 WRITE (6,13)
C
0052 DO14K=1,IP1
0053 CALL STDEV (TOTAL.MEAN.SD.SE.CL)
0054 TOT(I,K)=TOTAL
0055 M(I,K)=MEAN
0056 STD(I,K)=SD
0057 STE(I,K)=SE
0058 14CL95(I,K)=CL
C
C L1=NUMBER OF SETS COMPUTED.
C
C WRITE FIRST TWO PAGES.
0059 L1=IP1
0060 CALL PRINT
0061 WRITE (6,12) (NAME(N),N=1,IP1)
0062 WRITE (6,11) (NAME(N),N=IP2,I2TP2)
0063 WRITE (6,13)
0064 CALL PR1NT2
0065 15 CONTINUE
C
C
C COMPUTE PERCENTS.
80
-------
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079
0080
0081
0082
0083
0084
0085
0086
0087
0088
0089
C
c
C
c
DO 181=3,4
IB=IA(I-2)
L1=ID
L=I-2
WRITE (6,10) L
WRITE (6,12) (NAME(N),N=1,ID)
WRITE (6,11) (NAME(N),N=I2TP3,I3TP2)
WRITE (6,13)
DO 16K=1,ID
DO 16J=1,IB
16X(IJ,K)=X(L,J,K)/X(LJ,IP1)*100.
DO 17K=1,3
CALL STDEV (TOTAL,MEAN,SD,SE,CL)
TOT(I,K)=TOTAL
M(I,K)=MEAN
STD(I,K)=SD
STE(I.K)=SE
17 CL95(I,K)=CL
CALL PRINT
WRITE (6,12) (NAME(N),N=1,ID)
WRITE (6,11) (NAME(N),N=I2TP3,I3TP2)
WRITE (6,13)
CALL PRINT2
18 CONTINUE
0090
0091
0092
0093
0094
0095
0096
0097
C
C
C
C
C
CALL NCOMP
CALCULATE TOTAL AND MEAN OF N.
DO20K=1,ID
TOT(5,K)=0.
DO 19J = 1,IA3
19 TOT(5,K)=TOT(5,K)+X(5,J,K)
20M(5,K)=TOT(5,K)/AI3
DO22K=1,ID
V=0.0
81
-------
0098
0099
0100
0101
0102
0103
0104
0105
0106
0107
0108
0109
0110
0111
0112
0113
0114
0115
0116
0117
0118
0119
0120
0121
0122
0123
C
c
C
c
c
c
c
c
c
c
c
c
c
DO21J=1,IA3
21 V=(M(5,K)-X(5J,K))**2+V
STD( 5 ,K)=SQRT( V/( AI3-1.0))
22 CALL STDEV2 (STD(5,K))STE(5,K),CL95(5,K))
CALCULATE K VALUES.
DATA IN TWO SETS ARRANGED CHRONOLOGICALLY
CALCULATE K VALUES
DO24K=1,ID
SUM1(K)=0.0
IB=IA(2)
DO 23 J=1,IB
IF (X(2,J,K).EQ.O.) X(2,J,K)=.004
V=(DLOG 10(X(2, J ,K))-ALOG 10(M( 1 ,K)))/(X( 5 J ,K»
V2(J,K)=10.**V-1.0
23 SUM1(K)=SUM1(K)+V2(J,K)
IB=IA(1)
DO24J = 1,IB
IF (X(1,J,K).EQ.O.) X(l J,K)=.004
V=ALOG10(M(2,K))-DLOG10(X(1,J,K)))/X(5,J+IA(2),K)*.43429)
V2(J+IA(2),K)=10.**(V*.43429)-1.Q
24 SUM1(K)=SUM1(K)+V2(J+IA(2),K)
SORT VALUES
DO 26K=1,ID
D026J = 1J2T
IF (V2(J,K).GT.O.) GO TO 25
X(7J,K)=V2(J,K)
GO TO 26
25 X(6J,K)=V2(J,K)
26X(8J,K)=X(7J,K)+X(6J,K)
CALCULATE K-NET
DO28K=1,ID
82
-------
0124* DO28J=1,J2T
01-25 V=X(8J,K)-SUM1(K)/AI3
0126 IF (V.GT.O) GO TO 27
0127 X(10J,K)=V
0128 GO TO 28
0129 27 X(9J,K)=V
0130 28 CONTINUE
C
C COMPUTE SUM & MEAN FOR K VALUES
C
0131 DO30K=1,ID
C
0132 DO 30 1=6,10
0133 V=0.0
C
0134 DO29J=1,J2T
0135 29 V=V+X(I,J,K)
C
0136 TOT(I,K)=V
0137 30M(I,K)=V/AI3
C
0138 DO 31 1=6,10
C
0139 DO31K=1,7
0140 STD(I,K)=0.0
0141 STE(I,K)=0.0
0142 31 CL95(I,K)=0.0
C
C CALCULATE STANDARD DEVIATION, STANDARD ERROR, AND 95% CONFIDENCE
C LIMITS OF K VALUES.
0144 DO 32 K=1,ID
0145 SUM1(K)=0.0
0146 SUM2(K)=0.0
0147 SUM3(K)=0.0
0148 32 SUM4(K)=0.0
0149 IB=IA(2)
C
0150 DO33J=1,IB
0151 V2(J,K)=DLOG(X(2,J ,K))-DLOG(X( 1J ,K))
0152 SUM2(K)=V2(J,K)+SUM2(K)
0153 SUM3(K)=(DLOG(X(1,J,K))-ALOG(M(1,K)))**2+SUM3(K)
0154 33 SUM4(KHDLOG(X(2J,K))-ALOG(M(2,K)))'*2+SUM4(K)
C
0155 DO34K=1,ID
C
83
-------
0156
0157
0158
0159
0160
0166
0167
0168
0169
C
C
C
VARl(K)<43429/M(l,K))**2*SUM3(K)/(AI3-1.0)+(-.43429/M(2,K))*»2
1*SUM4(K)/(AI 3-1.0)
34V1(K)=SUM2(K)/AI3
DO36K=1,ID
VAR2(K)=((1.0/M(5,K))**2*VARl(K))+(-Vl(K)/M(5,K)**2))**2
-------
0197
0198
0199
0200
0201
0202
0203
0203
0204
0205
0206
0207
0208
0209
0222
0223
0224
0225
0226
0227
0228
0229
0230
0231
0232
0233
0234
0235
0236
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
V=V/(AI3-1.0)
W=W/(AI3-1.0)
STD(9,K)=SQRT(((V/(V+W))* *2*(STD(8,K)* *2))
CALLSTDEV2(STD(9,K),STE(9,K),CL95(9,K))
STD(10,K)=SQRT(((W/(V+W))**2*(STD(8,K)**2))
46 CALL STDEV2(STD(10,K),STE(10,K),CL95(10)K))
CALL PRINT3
CALCULATE O AND ITS STANDARD DEVIATION
DO 52 K=1,ID
M(11,KHM(9,K)/M(6,K))*M(7,K)
STD(11,K)=SQRT(STD(7,K)**2*((M(9,K)/M(6,K))**2))
CALL STDEV2 (STD(11,K),STE(11,K),CL95(11,K))
CALCULATION OF D
M(12,K)=M(7,K),-M(11,K)
CALCULATION OF STANDARD DEVIATION OF D
STD(12,K)=SQRT(STD(7,K)**2*(1.-M(9,K)/M(6,K))**2)
CALL STDEV2 (STD(12)K),STE(12,K),CL95(12,K))
CALCULATE TL.
M(13,K)=-1.0*(1.0/M(12,K))
STD(13,K)=DSQRT(STD(12,K)**2*(1.0/M(12,K)**2)**2)
52 CALL STDEV2(STD(13,K),STE(13,K),CL95(13,K))
CALCULATE TR.
DO 53 K=1,ID
M(14,K)=-1.0*(1.0/M(7,K))
STD(14,K)=SQRT(STD(7,K)**2*(1.0/M(7,K)**2)**2)
53 CALL STDEV2 (STD(14)K),STE(14,K),CL95(14,K))
DO 71 K=1,ID
WRITE (6,54) NAME(K)
54 FORMAT ('1',1X,'RATES OF CHANGE FOR))A8,30X,'S.D.')7X)'S.E.',4X,
1 '95%LIMIT'//)
WRITE (6,55) NAME (K),M(8,K),STD(8,K),STE(8,K),CL95(8,K)
55 FORMAT (2X,'MEAN OF KM3X,'=NET',3X,A8,'=',3X,4F11.4/)
WRITE (6,56) NAME(K),M(9,K)STD(9,K),STE(9,K),CL95(9)K)
56 FORMAT (2X, 'MEAN OF + ( K-NET ) = T' ,5X,A8,'=>f3X>4F11.4)
WRITE (6,57)
85
-------
0237 57 FORMAT (27X.T/)
0238 WRITE (6,58) NAME(K)rM(10,K),STD(10)K))STE(10)K))CL95(10,K)
0239 58 FORMAT (2X,'MEAN OF - ( K - NET ) = T',5X,A8,'=',3X,4F11.4)
0240 WRITE (6,59)
0241 59FORMAT(27X,'OV)
0242 WRITE (6,60) NAME(K),M(6,K),STD(6,K),STE(6,K),CL95(6,K)
0243 60 FORMAT (2X/MEAN OF + K',11X,'= r,5X,A8,'=',3X,4Fl 1.4//)
0244 WRITE (6,61) NAME(K),M(7,K),STD(7,K),STE(7,K),CL95(7,K)
0245 61 FORMAT (2X,'MEAN OF - K',11X,'= O + D',1X,A8,'=',3X,4F11.4//)
0246 ' WRITE (6,62) NAME(K),M(11,K),STD(11,K),STE(11,K),CL95(11,K)
0247 62 FORMAT (26X,'O')5X,A8,'=',3X,4F11.4/)
0248 WRITE (6,63) NAME(K),M(12,K),STD(12,K),STE(12tK),CL95(12,K)
0249 63 FORMAT (26X,'D')5X,A8,<=',3X,4F11.4/)
0250 WRITE (6,64) NAME(K),M(13,K),STD(13,K),STE(13,K),CL95(13,K)
0251 64 FORMAT (26X,'T>,5X,A8,'=',3X,4F11.4)
0252 WRITE(6,65)
0253 65 FORMAT(27X,'LV)
0254 WRITE(6,64) NAME(K),M(14,K),STD(14,K),STE(14,K),CL95(14,K)
0255 WRITE(6,66)
0256 66 FORMAT(27X,'R')
0257 DO 67 L=l,3
0258 67WRITE(6,13)
0259 WRITE(6,68)
0260 68 FORMAT(13X,'MEAN C',6X,'MEAN C'.ieX.T.lOX.'T'^X.'-'^X.'T'^X,
-<-')5X,'OI,5X,'-',5X,'D>)9X,'N>/19X,'2>,llX,'l',27X,'r,llX,'OV)
0261 69FORMAT(/97X,F11.4)
0262 WRITE(6,69) M(5,K)
0263 WRITE(6,70) NAME(K),M(2,K),M(1,K)M(6,K),M(9,K)>1(10,K)JVI(11(K))
-M(12,K)
0264 70 FORMAT(2X,A8,F10.4,' =',F10.4,' ( 1.0 +',F10.4,' +',F10.4,3(F12.4)
V )')
0265 71 CONTINUE
0266 CALL FACTOR
0267 STOP
0268 END
C
C
0001 SUBROUTINE PRINT
0002 DIMENSION TABLE(30),MONTH(12),ALOC(2,60,6),TOT(10,8),STD(23,8).
1, STE (23,8),CL95(23,8),IA(2),AI(2)
0003 REAL *4MEAN,M(17,8)
0004 REAL *8X(10,60,7),V2(60,7),NAME(23),V1(7)
0005 INTEGER CST(2,60),CDATE(2,60,3)
0006 COMMON X,TABLE,IA,IB,I,K,KD,ID
0007 COMMON /BLK1/ NAME/TOT,M)STD)STE,CL95,ALOC,YR,CST,CDATE,MONTH,L1,I
C
86
-------
0008 DO 1 J=1,IB
0009. 1 WRITE (6,3) CST(L,J),(ALOC(LJ,K),K=1,6),(CDATE(LJ,K),K=1,3),(X(I
1J,K),K=1,L1)
C
C SKIP TO BOTTOM OF PAGE
0010 N=(68-(IB+6))/2
C
0011 DO2J=1,N
0012 2 WRITE (6,4)
C
0013 RETURN
C
0014 3 FORMAT (5X)I2)5X,2A4,A2,2X)2A4,A2,2X,I2,2(<-')I2))8F11.2)
0015 4 FORMAT (/)
0016 END
C
C
0001 SUBROUTINE PRINT2
0002 DIMENSION TABLE(30))MONTH(12),ALOC(2)60,6))TOT(10,8))STD(23,8).
1, STE(23,8),CL95(23,8), IA(2),AI(2)
0003 REAL *MEAN,M(17,8)
0004 REAL *8X(10,60,7),V2(60,7),NAME(23),V1(7)
0005 INTEGER CST(2,60),CDATE(2)60,3)
0006 COMMON X.TABLE.IA.IB.I.K.KD.ID
0007 COMMON /BLK1/ NAME,TOT,M,STD,STE,CL95,ALOC,YR,CST,CDATE,MONTH,L1,L
0008 WRITE(6,1)(TOT(IJ),J = 1,L1)
0009 WRITE (6,2) (M(I,J),J=11L1)
0010 WRITE (6,3) (STD(IJ)J=1,L1)
0011 WRITE (6,4) (STE(IJ)J=1,L1)
0012 WRITE (6,5) (CL95(I,J),J = 1,L1)
0013 RETURN
C
0014 1 FORMAT (34X,'TOTALS',6X,7F 10.4)
0015 2 FORMAT (/34X,'MEAN',8X,7F10.4)
0016 3 FORMAT (/34X,'S.D.',8X,7F10.4)
0017 4 FORMAT (/34X)'S.E.',8X,7F10.4)
0018 5 FORMAT (/34X,'95% CL',6X,7F10.4)
0019 END
C
C
0001 SUBROUTINE PRINT3
0002 DIMENSION TABLE(30),MONTH(12))ALOC(2,60,6),TOT(10,8),STD(23,8)
1,STE(23,8),CL95(23,8),IA(2),AI(2)
0003 REAL •4MEAN,M(17,8)
0004 REAL *8X(10,60,7),V2(60,7),NAME(23),V1(7)
0005 INTEGER CST(2,60),CDATE(2,60,3)
0006 COMMON X.TABLE.IA.IB.I.K.KD.ID
87
-------
0007 COMMON /BLK1/ NAME)TOTrM)STD,STE,CL95,ALOC,YR)CST,CDATE,MONTH)Ll,L
0008 DO 19K=1,ID
0009 WRITE (6,3)
0010 WRITE (6,4)
0011 WRITE (6,5) NAME(K),NAME(K)
0012 WRITE (6,6)
0013 WRITE (6,7)
0014 IB=IA(2)
C
0015 DO1J=1,IB
0016 1 WRITE (6,17) CST(2,J),X(2)J,K),X(5,J,K),(X(IXJ,K),IX=6,10)
C
0017 JPIA=IA(2)
0018 IB=IA(1)
C
0019 DO2J=1,IB
0020 JPIA=JPIA+1
0021 2 WRITE (6,18) CST(1,J),X(1,J,K))X(5,JPIA,K),(X(IX>JPIA,K),IX=6,10)
C
0022 WRITE (6,7)
0023 WRITE (6,16) TOT(2,K),TOT(1,K),TOT(5,1),(TOT(N,K),N=6,10)
0024 WRITE (6,15)
0025 WRITE (6,13) NAME(K)
0026 WRITE(6,16)M(2,K),M(1,K),M(5,1),(M(N,K),N=6,10)
0027 WRITE (6,8)
0028 WRITE (6,13) NAME(K)
0029 WRITE (6,16) STD(2,K),STD(1,K),STD(5,1),(STD(N,K),N=6,10)
0030 WRITE (6,9)
0031 WRITE (6,12) NAME(K)
0032 WRITE (6,16) STE(2,K),STE(1,K),STE(5,1),(STE(N,K),N=6,10)
0033 WRITE (6,10)
0034 WRITE(6,12)NAME(K)
0035 WRITE (6,16) CL95(2,K),CL95(1,K),CL95(5,1),(CL95(N,K),N=6,10)
0036 WRITE (6,11)
0037 19 WRITE (6,14) NAME(K)
0038 RETURN
C
0039 3 FORMAT ('1 MX/STATION')
0040 4 FORMAT (12X,'C',9X,
-------
0047
0048
0049
0050
0051
0052
0053
0054
0055
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
C
C
11 FORMAT (V,94X,'95% CONFIDENCE LIMITS')
12 FORMAT ('+', 99X.A8)
13 FORMAT (V,102X,A8)
14 FORMAT (V,116X,A8)
15 FORMAT ('+')94X)'TOTALS')
16 FORMAT (/9X,5F10.4,3F11.4)
17 FORMAT (4X,I2,1X)F10.2,12X,3F10.4,3F11.4)
END
SUBROUTINE TDIST (KA,T)
REAL *8X(10,60,7)
DIMENSION TABLE(30), IA(2), AI(2)
COMMON X)TABLE,IA,IB,I)K,KD)ID
I1=KA-1
AK=I1
IF (II) 1,1,2
1 WRITE (6,11) I
GO TO 10
2IF(I1.LT.31)GOTO9
IF(I1.LT.41)GOTO3
GO TO 4
3TINT=((2.042-2.021)/10.)*(AK-30.)
T=TINT+2.042
GO TO 10
4IF(Il.Lt.61)GOTO5
GO TO 6
5TINT=((2.021-2.000)/20.)*(AK-40.)
T=TINT+2.021
GO TO 10
6IF(I1.LT.121)GOTO7
GO TO 8
7TINT=((2.000-1.980)/40.)*(AK-60.)
T=TINT+2.000
GO TO 10
8T=1.960
GO TO 10
9T=TABLE(I1)
10 RETURN
11 FORMAT ('1VI INT TABLE =',I3)
END
C
C
89
-------
0001
0002
0003
0004
0005
0006
SUBROUTINE STDEV (SUMX,XBAR,STD,STE,CL$)
REAL *8X( 10,60,7)
DIMENSION TABLE(30), IA(2), AI(2)
COMMON X,TABLE,IA,IB,I,K,KD,ID
DEV=0.
SUMX=0.
0007
0008
0009
0010
DO 1 J = 1,IB
1SUMX=SUMX+X(I,J,K)
AIUHA(I)
XBAR=SUMX/AI(I)
0011
0012
0013
DO2J=1,IB
DE V=(XB AR-XU ,J ,K)) * * 2+DEV
2 CONTINUE
0014
0015
0016
0017
0018
0019
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0001
0002
0003
0003
0004
0005
0006
0007
0008
0009
C
C
C
C
STD=SQRT(DEV/(AI(I)-1.))
STE=STD/SQRT(AI(I))
KA=IB
CALL TDIST (KA.T)
CL$=T*STE
END
SUBROUTINE STDEV2 (STD,STE,CL$)
REAL *8X( 10,60,7)
DIMENSION TABLE(30), IA(2), AI(2).
COMMON X,TABLE,IA,IB,I,K,KD,ID
AI3=IA(1)+IA(2)
STE=STD/SQRT(AI3)
KA=IA(l)-i-IA(2)
CALL TDIST (KA,T)
CL$=T*STE
RETURN
END
SUBROUTINE NCOMP
DIMENSION TABLE(30),MONTH(12),ALOC(2,60,6),TOT(10,8),STD(23,8)
1, STE(23,8),CL95(23,8),IA(2), AI(2)
DIMENSION IYRVAL(5),ITOTDA(5)
REAL *4M(17,8)
REAL *8X(10,60,7),NAME(23)
INTEGER CST(2,60),CDATE(2,60,3)
INTEGER SUMDA(2)
COMMON X,TABLE,IA,IB,I,K,KD,ID
COMMON /BLK1/ NAME,TOT,M,STD,STE,CL95,ALOC,YR,CST)CDATE,MONTH,L1
IJ=0
90
-------
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
C
C
C
C
C
C
C
C
C
DO 131=1,2
K=0
SUMDA(I)=0
ITOTDA(I)=0
STORE INITIAL TIME
MO1=CDATE(I,1,1)
IDA1=CDATE(I,1,2)
IYR1=CDATE(I,1,3)
IYRVAL(1)=365
A=IYR1
IZ=A/4.
IF (IZ*4.EQ.IYR1.AND.MO1.GT.2) IYRVAL(1)=366
INT1=MONTH(MO1)+IDA1
FIND TIME INTERVAL OF FIRST DATE TO END OF FIRST YEAR
INT2=IYRVAL(1)-INT1
IB=IA(I)
DO4J=1,IB
MO2=CDATE(I,J,1)
IDA2=CDATE(IJ,2)
IYR2=CDATE(IJ,3)
COMPUTE YEAR VALUES-365 OR 366
IF (IYR1.EQ.IYR2) TO TO 3
K STORES NUMBER OF INTERVENING YEARS
K=IYR2-IYR1
DO 1 L=1,K
IYRVAL(L+1)=365
A-IYR1+L
IZ=A/4.
IF (IZ*4.EQ.IYR1+L) IYRVAL(L+1)=366
1 CONTINUE
COMPUTE INTERVAL OF LAST YEAR
LAST1=MONTH(MO2)+IDA2
CHECK FOR LEAPYR OF LAST YEAR
IF (IYRVAL(K).EQ.366.AND.MO2.GT.2) LAST=LAST+1
LAST2=IYRVAL(K)-LAST1
COMPUTE TOTAL DAYS OF DATA SET
INT=FIRST YEAR
K=K+1
K= NUMBER OF YEARS
91
-------
0040 DO 2 L=1,K
0041 2 ITOTDA(I)=ITOTDA(I)+IYRVAL(L)
C
C SUM ALL DAYS OF YEARS INVOLVED
0042 ITOTDA(I )=ITOTD A( I )-INT 1-LAST2
0043 SUMDA(I)=SUMDA(I)+ITOTDA(I)
0044 GO TO 4
0045 3 INT2=MONTH(MO2)+IDA2
0046 ITOTDA(I)=INT2-INT1
0047 SUMDA(I)-SUMDA(I)+ITOTDA(I)
0048 4 CONTINUE
C
C COMPUTE MEAN OF TIME
0049 MEANT=SUMDA(I)/IA(I)
C SUBTRACT FIRST YEAR
0050 IX=MEANT+INT1
0051 IF (K.EQ.O) GO TO 7
C
0052 DO 5 L=1,K
005 3 IF (IX.LT.IYRVAL(L)) GO TO 6
0054 IF (IX.EQ.IYRVAL(L)) TO TO 6
0055 IX=IX-IYRVAL(L)
0056 5 CONTINUE
C
C COMPUTE YEAR
0057 6 IYR=L-1+IYR1
0058 IF (IYRVAL(L).EQ.366.AND.IX.GT.59) IX=IX-1
0059 GO TO 8
0060 7 IF (IYRVAL(1).EQ.366.AND.IX.GT.59) IX=IX-1
0061 IYR=CDATE(I,1,3)
C
0062 8DO9N=1,12
C LOCATE MONTH
0063 IF (IX.LT.MONTH(N+1)) GO TO 10
0064 IF (IX.EQ.MONTH(N)) GO TO 10
0065 9 CONTINUE
C
0066 10 IMON=N
0067 IDAY=IX-MONTH(N)
0068 IF(I.EQ.1)IC=IA(2)
0069 IF(I.EQ.2)IC=IA(1)
C
0070 DO12J = 1,IC
0071 IJ=IJ+1
007 2 CALL LEAPYR (J ,IMON ,IDAY ,IYR)
92
-------
0073
0075
0076
0077
11
DO 11K=1,ID
X(5,IJ,K)=YR
CONTINUE
12 CONTINUE
13 CONTINUE
0078
0079
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
C
C
RETURN
END
SUBROUTINE LEAPYR (J,IMON,IDAY,IYR)
DIMENSION TABLE(30), MONTH(12), ALOC(2,60,6), TOT(10,8), STD(23,8)
1,STE(23,8), CL95(23,8), IA(2), AI(2)
REAL MMEAN, M(17,8)
REAL *8X(10,60,7),NAME(23)
INTEGER YR1,YR2,DA1(DA2,DAYS
INTEGER CST(2,60),CDATE(2,60,3)
COMMON /BLK1/ NAME,TOTJVl)STD,STE>CL95,ALOC,YR,CST,CDATE JVIONTH.L1.L
COMMON X,TABLE,IA,IB,I,K1KD,ID
DAYS=0
NT=0
IF (I.EQ.2) GO TO 1
MO1=IMON
DA1=IDAY
YR1=IYR
MO2=CDATE(2,J,1)
DA2=CDATE(2,J,2)
YR2=CDATE(2,J,3)
GO TO 2
1 MO2=IMON
DA2=IDAY
YR2=IYR
MO1=CDATE(1J,1)
DA1=CDATE(1,J,2)
YR1=CDATE(1J,3)
2 AMO=MO1
DO6IY=YR1,YR2
A=IY
LEAP=0
IZ=A/4.
Z=IZ
Z=Z*4.
IF (IY.EQ.YR1) GO TO 3
GO TO 4
93
-------
0034 3 DAYS=365-(MONTH(MO1)+DA1)
0035 IF (Z.EQ.A.AND.AMO.LT.3.) LEAP=1
0036 GO TO 5
0037 4 IF (Z.EQ.A) LEAP=1
0038 5 NT=DAYS+LEAP+NT
0039 6DAYS=365
C
0040 IF(LEAP.EQ.1)GOTO7
0041 GO TO 8
0042 7 IF (MO2.LT.3) NT=NT-1
0043 8 YR=NT-365+MONTH(MO2)+DA2
0044 YR=YR/365.
0045 RETURN
0046 END
94
-------
c
c
0001
0002
0003
0002
0003
0004
0005
0006
0007
C
c
c
c
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
SUBROUTINE FOR PROGRAM FOR ESTIMATING SYSTEM RATES
BASED ON SAMPLE VALUES PAIRED TO MEAN VALUES.
SUBROUTINE FACTOR
DIMENSION TABLE(30), MONTH(12), ALOC(2,60,6), TOT(10,8), STD(23,8)
1, STE(23,8), CL95(23,8), VAR1(7), VAR2(7), VAR3(7), SUM1(7), SUM2(
27), SUM3(7), SUM4(7), COV1(7), COV2(7), IA(2), AI(2)
REAL *4MEAN,M(17,8),MR(7)
REAL'8X(10,60,7),V2(60,7),NAME(23),V1(7)
INTEGER CST(2,60),CDATE(2>60,3)
COMMON X,TABLE,IA,IB,I,K,KD,ID
COMMON /BLK1/ NAME,TOT,M,STD,STE,CL95,ALOC,YR,CST,CDATE)MONTH,LI,L
COMMON /BLK2/ MR
CALCULATE CORRECTION FACTOR FOR STANDARD DEVIATION
IA3=IA(1)+IA(2)
AI3=IA3
DO 14 K=l, ID
IF (M(8,K» 1,4,5
1 JX=(AI3*M(7,K)-AI3*M(10,K))/M(8,K)
VJ=JX
V=( ((AI3*M(10,K) )/VJ)**2)*VJ
IF ((M(8,K)+( (AI3*M(10,K))/VJ))*VJ-(AI3*M(7,K))) 3,3,2
2 V=V+(A!3*M(7,K)-VJ*(M(8,K)+( (AI3*M(10,K) )/VJ))-M(8,K))**2
V=V+(((AI3*M(6,K))/AI3-VJ-1.0))-M(8,K))**2*(AI3-VJ-1.0)
GO TO 8
3 V=V+( ((AI3*M(6,K))/(AI3-VJ))-M(8,K))**2*(AI3-VJ)
GO TO 8
4 JX-AI3/2.0
VJ=JX
V=( (AI3*M(6,K)/VJ)**2)*VJ
V=V+( (AI3*M(7,K)/VJ)**2)*VJ
GO TO 8
5 JX=(AI3*M(6,K)-AI3*M(9,K))/M(8,K)
VJ=JX
V=(((AI3*M(9,K))AVJ)**2)*VJ
IF ((M(8,K)+((AI3'M(9,K))/VJ))*VJ-(AI3*M(6,K))) 6,7,7
6 V=V+(AI3*M(6,K)-VJ*(M(8,K)+((AI3*M(9,K))/VJ))-M(8,K))**2
V=V+(((AI3*M(7,K))/(AI3-VJ-1.0))-M(8,K))**2*(AI3-VJ-1.0)
GO TO 8
7 V=V+( ((AI3*M(7,K))/(AI3-VJ))-M(8,K))**2*(AI3-VJ)
8 V=V/(AI3-1.0)
W=0.0
C
C
95
-------
0035 DO9J=1,IA3
0036 9 W=W+(X(8,J,K)-M(8,K))**2
C
c
0037 W=W/(AI3-1.0)
0038 C=((W-V)/W)**2
C
C CALCULATE CORRECTED STD,6,7,AND 8
C STD(15,K) IS CORRECTED STD(6,K)
C
0039 STD(15,K)=SQRT(C*STD(6,K)**2)
0040 CALL STDEV2 (STD(15,K),STE(15,K),CL95(15,K))
C
C STD(16,K)IS CORRECTED STE(7,K)
C
0041 STD(16,K)=SQRT(C*STD(7,K)**2)
0042 CALL STDEV2 (STD(16,K),STE(16,K),CL95(16,K))
C
C STD(17,K) IS CORRECTED STD(8,K)
C
0043 STD(17,K)=SQRT(C*STD(8,K)**2)
0044 CALL STDEV2 (STD(17,K),STE(17,K),CL95(17,K))
C
C CALCULATE CORRECTED STD(9,K) AND STD(10,K)
C
0045 V=0.0
C
C
0046 DO11J=1,IA3
C
0047 IF (X(9,J,K)) 11,11,10
0048 10 V=V+X(9,J,K)**2
0049 11 CONTINUE
C
C
C
0050 V=V/(AI3-1.0)
0051 W=0.0
C
C
0052 DO13J=1,IA3
C
0053 IF (X(10J,K)) 12,13,13
0054 12 W=W+(X(10,J,K))**2
0055 13 CONTINUE
C
C
C
96
-------
0056 W=W/(AI3-1.0)
C
C STD(18,K) IS CORRECTED STD (9,K)
C
0057 STD(18,K)=SQRT((V/(V+W))**2*C*(STD(8,K)**2))
0058 CALL STDEV2 (STD(18,K),STE(18,K),CL95(18,K))
C
C
C STD(19,K) IS CORRECTED STD(10,K)
C
0059 STD(19,K)=SQRT((W/(V+W))**2*C*(STD(8,K)**2))
0060 CALL STDEV2 (STD(19,K),STE(19,K))CL95(19,K))
C
C CALCULATE CORRECTED STD(20,K) CORRECTED STD( 11 ,K)
C
0061 STD(20,K)=SQRT(STD(16,K)**2*(M(9,K)/M(6,K))**2)
C
C
0062 CALL STDEV2 (STD(20,K),STE(20,K),CL95(20,K))
C CALCULATE STD(21,K) CORRECTED STD(12,K)
0063 STD(21,K)=SQRT(STD(16,K)**2*(1.0-M(9,K)/M(6,K))**2)
0064 CALL STDEV2 (STD(21,K),STE(21,K),CL95(21,K))
C
C CALCULATE STD(22,K) CORRECTED STD(13,K)
C
0065 STD(22JC)=SQRT(STD(21,K)**2*(1.0/M(12,K)**2)**2)
0066 CALL STDEV2 (STD(22,K))STE(22,K))CL95(22,K))
C
C CALCULATE STD(23,K) CORRECTED STD(14,K)
C
0067 STD(23,K)=SQRT(STD(16,K)**2*(1.0/M(7,K)**2)**2)
0068 CALL STDEV2 (STD(23,K),STE(23,K),CL95(23,K))
C
C
0069 14 CONTINUE
C
C
0070 DO15K=1,ID
C
0071 WRITE (6,16) NAME(K)
0072 WRITE (6,17) NAME(K),M(8,K),STD(17,K),STE(17,K))CL95(17,K)
0073 WRITE (6,18) NAME(K),M(9,K),STD(18,K),STE(18)K),CL95(18,K)
0074 WRITE (6,19)
0075 WRITE (6,20) NAME(K),M(10,K))STD(19,K),STE(19,K),CL95(19,K)
0076 WRITE (6,21)
0077 WRITE (6,22) NAME(K),M(6,K),STD(15,K),STE(15,K),CL95(15,K)
0078 WRITE (6,23) NAME(K),M(7,K),STD(16,K),STE(16,K),CL95(16,K)
97
-------
0079 WRITE (6,24) NAME(K),M(11,K),STD(20,K),STE(20)K),CL95(20,K)
0080 WRITE (6,25) NAME(K),M(12,K),STD(21,K),STE(21)K),CL95(21,K)
0081 WRITE (6,26) NAME(K),M(13,K),STD(22,K),STE(22,K),CL95(22,K)
0082 WRITE (6,27)
0083 WRITE (6,26) NAME(K),M(14,K),STD(23,K),STE(23)K),CL95(23,K)
0084 WRITE (6,28)
0085 15 CONTINUE
C
C
0086 16 FORMAT (T,40X,'CORRECTED STANDARD DEVIATIONS'//,2X, 'RATES OF CHA
INGE FOR,' A8,30X,'S.D.',7X,'S.E.',4X,'95% LIMIT'//)
0087 17 FORMAT (2X,'MEAN OF K',13X,'=NET',3X,A8,'=',3X,4F11.4/)
0088 18 FORMAT (2X/MEAN OF + ( K - NET ) = T',5X,A8,<=',3X,4F11.4)
0089 19 FORMAT (27X.T/)
0090 20 FORMAT (2X/MEAN, OF - ( K - NET ) = T',5X,A8,'=',3X,4F11.4)
0091 21 FORMAT (27X/OV)
0092 22 FORMAT (2X/MEAN OF+ K',llX,'=r,5X,A8>',3X)4F11.4//)
0093 23 FORMAT (2X/MEAN OF - K',11X,'= 0+D',1X,A8,'=',3X,4F11.4//)
0094 24 FORMAT (26X)'O',5X,A8,'=',3X,4F11.4/)
0095 25 FORMAT (26X,'D',5X,A8,'=>,3X,4F11.4/)
0096 26 FORMAT (26X,'T')5X)A8,'=',3X,4F11.4)
0097 27 FORMAT (27X,'LV)
0098 28 FORMAT (27X,'R')
C
0099 RETURN
0100 END
98
-------
TECHNICAL REPORT DATA
I'lfusr rratl lnttfucttons on the rcivrse
i in rim i NCI.
EPA-660/3-75-013
I, I ifl. I. ANU I'UII II I LK
"An Analysis of the Dynamics of DDT in Marine
Sediments."
3. RECIPIENT'S ACCESSION-NO.
5. REPORT DATE
MAY 1975
6. PERFORMING ORGANIZATION CODE
/ AUTMOHI3)
Phillips, John H., Eugene E. Haderlie, and Welton L. L
8. PERFORMING ORGANIZATION REPORT NO.
e
y. I'l Rl OHMINO OHO -VNIZATION NAME AND ADDRESS
Hopkins Marine Station
Stanford University
Pacific Grove, CA 93950
A SI'l'NKOFMNG Af.LNCY NAME AND ADDRESS
Environmental Protection Agency
200 S. W. 35th St.
Corvallis, OR 97330
10. PROGRAM ELEMENT NO.
11. CONTRACT/GRANT NO.
13. TY.PE OF REPORT AND PERIOD COVERED
14. SPONSORING AGENCY CODE
16. SUPPl.C.MtNTARY NOTES
AT ACT The concentrations of DDT, ODD and DDE were measured in sediments at 57
stations in Monterey Bay on the Central California coast during 1970-1971. Mean
concentrations in parts per billion were DDT 3.1, ODD 2.3, and DDE 5.4. During
1973 nineteen of the original stations were sampled. Mean concentrations were
DDT 15.5, ODD 2.3, and DDE 5.4 PPB.
Two approaches to ithe estimation of annual system rates for input, I, output,
0, decay, D, and internal translocation, T, and TQ expressed as decimal fractions
of existing concentraticns were developed, and fraction programs that permit
rapid estimations were written. The mean annual rates in South Monterey Bay
obtained were for DDT, I +1.30, 0 -0.059, D -0.036 Tj and TQ +0.80 with a residence
time of 11 years and life time of 29 years. An I of 1.30 means the amount of input
is 130% of the existing concentration per year. Rates for ODD were, I +0.25, 0
-0.11, D -0.025, Tj and T0 +0.20 with residence time of 7 years and life time of 44
years. Rates for DDE were I +0.28, 0 -0.10, D -0.027, T0 and Tj_ +0.22 with residence
time of 8 years and life time of 39 years.
Laboratory assays were developed to determine the relative rate of decomposition
in sediment under conditions selective for various physiologically different kinds of
microorganisms. Decay under aerobic conditions was greater than under anaerobic con-
ditions. Nitrate increased the rate °f decomposition under anaarobic conditions.
The QlQ for decay was 2.5.
17.
it.
KEY WORDS AND DOCUMENT ANALYSIS
DESCRIPTORS
Degradation, DDT, ODD, DDE, Coastal,
Pollution
b.lDENTIFIERS/OPEN ENDED TERMS C. COSATI Field/Group
Microbial degradation,
chemical degradation,
rates of change,
pesticide residues,
chlorinated hydrocarbon
pesticides
organic pesticides
12 and 33
m. iiir,rHInurION STATEMENT
19. SECURITY CLASS (Thli Report}
21. NO. OF PAGES
20. SECURITY CLASS (This page I
22. PRICE
CPA form 2220-1 (9-73)
O U.S. GOVERNMENT PRINTING OFFICE: 1975-698-255/118 REGION 10
------- |