Dar Zarouk Parameter - Nigeria
Dar Zarouk Parameter - Nigeria
Dar Zarouk Parameter - Nigeria
Research Article
ABSTRACT
This study was carried out to assess the direction of groundwater flow and the extent of interaction of the groundwater
within the premises of Micheal Okpara University of Agriculture, Umudike( MOUAU) and the Kwa Ibo River, Umudike -
Nigeria. The pattern of groundwater flow within the study locations was obtained through the use of hydraulic heads
which were determined using Global Positioning System (GPS). Fourteen (14) VES were taken in four traverses within
the premises and immediate environment of MOUAU in order to determine the aquifer parameters such as
conductivity, Dar Zarrouk parameters, transmissivity and the K values of the study area. The VES data were
interpreted using Resist software that works on the least-square inverse modeling technique which employs the
backward difference resistivity transforms to obtain the final model for each VES. The electrical conductivity of the
study area ranges from 0 .000014Ω-1 to 0.00014 Ω-1 , hydraulic conductivity K=8.00m/day, the transmissivity ranges
from 20.80m 2/day to 384.80m2/day. Water samples were taken from boreholes within the premises of MOUAU and
along the banks of Kwa Ibo River bordering the University. The water samples were analysed for physical and chemical
properties. Suspended particulate matter load at the College of Natural and Applied Sciences (CNAS) and Zennith bank
locations were 1200mg/L and 700mg/L respectively which is above World Health Organisation ( WHO ) standard. These
high values suggest contamination. Based on the pattern of flow and the determined values of particulate matter, the
interaction of the Kwa Iboe River and the boreholes were also found to occur at CNAS and the Zennith bank within the
study area. Generally, the flow is from north to south of the study area. The pattern and distribution of each parameter
was contoured using Surfer 8 softwares. The flow pattern and distribution actually show the direction of flow and the
actual position of groundwater interaction between the water in the boreholes and the water from the Kwa Ibo River.
Specifically, the analysed data constrained with borehole within the area suggests that the area has top laterite sands
underlain by sandy geomaterials.
Key words: Transmissivity, Electrical Conductivity, The Direction of Flow, Suspended Particulate Matter , Boreholes in
MOUAU, Kwa Ibo River , Umudike-Nigeria
Introduction
The knowledge of aquifer parameters is essential for the management of groundwater resource. The parameters,
which describe the dynamics of aquifers, include; geometry of the pore space, geometry of the rock particles
themselves, secondary geologic processes such as faulting and folding and secondary deposition. All these
parameters jointly affect the rate and pattern of groundwater flow.
SPM in earth sciences refers to suspended Particulate Matter. Suspended Particulate Matter (SPM) is
“…sediment carried in suspension by the turbulent components of the fluid or Brownian movement (Wilber 1983). It
is the residue in a well-mixed sample of water that will not pass through a standard, glass fibre (0.45µm) filter.
Suspended sediments can also be called, suspended solids or suspended particles. Generally, the suspended
loads in flowing water consist of grains less than 0.5mm in diameter. There are various types of suspended
sediment. Eroded soils produce the most important type of suspended solids on a large scale. Samples of eroded
soils include sand, silt and clay that are relocated by rainfall and overland flow and carried into rivers, forest, and
urban areas. Organic suspended particulates compose a significant part of suspended solid in most natural waters.
By convention, particulate matter (PM) in suspension is defined as the material that is retained on a 0.4 – 0.5 µm
pore size filter. Smaller material is considered to be dissolved. Muddy sediments consist of clays and some variable
www.gjournals.org 70
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
silt content and particle diameter (Pugh, 1987) is < 62.2 µm. Suspended particulate matter (SPM) regulates
the transport of all types of water pollutant in dissolved and particulate phases in Lakes, river and coastal areas,
regulates water clarity and sedimentation. Groundwater flow modeling is generally defined as the quantity of
groundwater available or direction of dissolved contaminant migration. It is also used to define the limit of a capture
zone, for a contamination recovery well (or well field), or for delineating a water well protection area (or recharge
area) for a water supply.( Igboekwe and Udoinyang, 2011)
The objective of this paper consists of presenting a suspended matter flow model to simulate the transport
of sediments by the Kwa Ibo River- Umudike into the borehole wells in MOUAU. This is an interesting problem
since a complex intrusion process occurs when a freshwater input runs into marine waters. To achieve this
objective, characterization of the geological formations within the area through interpretation of geophysical data,
hydrogeological analysis for aquifer characteristic and groundwater flow modeling of the area have been carried
out.
Kwa Ibo River, Umudike-Nigeria lies within the Ikwano and Umuahia South Local Government Areas and was
selected for detailed aquifer characterization and hydrogeological investigation. The location map of the study area
0 I 0 1 0 1
is in fig.1.1. The Kwa Ibo River watershed is located between latitudes 5 19 and 5 30 N and longitude 7 30 and
0 1 2
7 37 E . It covers a rectangular area of about 180Km . Also Michael Okpara University of Agriculture Umudike can
be located at the latitude of 5028. 6581N - 5029 .1761N and longitude of 7032 .2561E- 70 32 .8031E. The location
map of the University can be obtained in fig1.1
www.gjournals.org 71
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
Modelling of the suspended particulate matter by the kwa Ibo River was carried out using the interpretation of the
geophysical data within the area and some parameters from hydrogeological analysis of aquifer characteristics of
the wells located within the area. This enabled us to estimate the water quality of those wells and also the direction
of solute contaminant migration.
Groundwater quality reflects substances that are dissolved or suspended in the water. Suspended materials are not
transported far in most subsurface materials, but it is usually filtered out. In general, groundwater flow is very slow
and depends on the permeability (water transmitting ability) of the subsurface materials, as well as the hydraulic
gradient (slope of the water-table or pressure gradient for artesian conditions). Fig. 2 shows a typical modeling
layout.
The process of groundwater flow is generally assumed to be governed by the relations expressed in Darcy’s Laws
and the conservation of mass. Darcy’s Law is a phenomologically derived constitutive equation that describes the
flow of a fluid through a porous medium based on the results of experiments on the flow of water through beds of
sand. It states that the rate of flow is directly proportional to the drop in vertical elevation between two places in the
medium and directly proportional to the distance between them.
---------------------------1
The development of mathematical equations that describes the groundwater flow and transport process may be
developed from the fundamental principle of conservation of mass of fluid or of solute. Given a representative
volume of porous medium, a general equation for conservation of mass for volume mass is expressed as:
Rate of mass inflow – (rate of mass outflow + rate of mass production/consumption) = Rate of mass accumulation -
------------- 2
This statement of conservation of mass (or continuity equation) may be combined with a mathematical expression
of the relevant process to obtain a different equation describing flow of transport (Bear 1997, Domenico and
Schwartz, 1998; Freeze and Cherry, 1979).
The rate of flow of water through a porous media is related to the properties of the water, the properties of the
porous media and the gradient of the hydraulic head, as represented by Darcy’s Law which can be written as;
, ---------------- 3
A general form of the equation describing the transient flow of a compressible fluid in a non-homogeneous
anisotropic aquifer may be derived by combining Darcy’s Law with the continuity equation. A general groundwater
flow equation may be written in Cartesian coordinate as;
, !
----------------------------4
Where: Ss is the specific storage, LT-1 , t is the time, xi are Cartesian coordinates, W* is the volumetric flux per unit
volume.
If the aquifer is relatively thin compared to its lateral extent, it may be appropriate to assume that groundwater flow
is really two-dimensional. This allows the three-dimensional flow equation to be reduced to the case of two
dimensional areal flow for which several additional simplifications are possible. An expression similar to equation 4
may be derived for the two-dimensional areal flow of a homogeneous fluid in a confined aquifer and written as:
", !
---------------------- 5
W= ! % is the volumetric flux per unit area. When the equation 5 is applied to an unconfined (water table) aquifer
system, it must be assumed that the flow is horizontal and equipotential lines are vertical, that the horizontal
hydraulic gradient equals the slope of the water table and storage coefficient is equal to specific yield (Sy)
(Anderson and Woessner, 1992).
An equation describing the transport and dispersion of a dissolved chemical in flowing groundwater may be derived
from the principle of conservation of mass, equation (2), just as the general flow equation was derived (Bear, 1979;
Domenico and Schartz, 1998, Konikow and Groove, 1977, Igboekwe and Udoinyang ;2011, Reddel and Sunada,
1970).
A generalized form of the solute transport equation is presented by Groove (1976), in which terms are
incorporated to represent chemical reaction and solute concentration both in the pore fluid and on the solid
surfaces as
www.gjournals.org 73
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
Transmissivity is a measure of the amount of water that can be transmitted horizontally by a full-saturated thickness
of the aquifer under hydraulic gradient. The relationship between transmissivity and hydraulic conductivity is given
by the relation;
T = bk (m/day)m=m2day.------------------ 7
Where b = saturated thickness of the aquifer in meters in
K= hydraulic conductivity in m/day
The study area has aquifer thickness of about 88.0m, hydraulic conductivity of about 8.00m/day and a
transmissivity of about 704.00m2/day (Igboekwe, Okwueze and Okereke, 2006).
HYDRAULIC CONDUCTIVITY
Hydraulic conductivity is a measure of the ease with which sediments can transmit water,
i.e. a medium has a unit hydraulic conductivity, if it will transmit in unit from a unit volume of ground water at the
prevailing kinematics velocity.
To obtain a layer parameter, a unit square cross sectional area is cut out of a group of n- layers of infinite
lateral extent. The total transverse resistance R is given by
P P
5 Where N and 5are respectively the thickness and resistivity of the ith layer in the section. The total
longitudinal conductance S is
S =∑O< --------------------------- 9
Q
www.gjournals.org 74
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
Si=R*N Where R* is the layer conductivity. Conductivity in this case is analogous to the layer transmissivity T, given
by
T= Kihi ---------------------------------------------------- 10
Ki is the hydraulic conductivity of the ith layer of the thickness hi. R and S of the above equation are called the Dar
Zarrouk parameters, which have been shown to be powerful interpretational aid in groundwater surveys (Zohdy et
al, 1974)
According to the fundamental Darcy’s Law, the fluid discharge, , is given by
= KIA ------------------------------------ 11
Where K is the hydraulic conductivity, I is the hydraulic gradient, A is the Cross-sectional area per perpendicular
to the direction of flow. The differential form of Ohms Law gives
j= R1 ---------------------------------------- 12
Where j is the
Current density; and R is the electrical conductivity, which is the reciprocal of the resistivity, 5. For aquifer material
having unit cross-sectional area A and thickness h, the two fundamental laws can combine to give
S
" #R9 T
#N ------------------------------- 13
Where T is the transmissivity, R is the transverse resistance of the aquifer, K is the hydraulic conductivity and S is
the longitudinal Conductance. 46.97m.
Table 4.1 shows the values of aquifer transmissivity and the values of transverse resistance and longitudinal
conductance generally known as Dar Zarrouk parameters. Aquifer transmissivities were determined using the
analytical relationship established by (Niswass and Singhal 1981, Mbonu et al., 1981) between aquifer
transmissivity (Tr) and transverse resistance (T) on one hand and aquifer transmissivity and longitudinal
conductance (S) on the other hand. For an aquifer of electrical conductivity (R and hydraulic conductivity (K), the
relation given in equation 1 was used
S
"U KR" KSR --------------------------------- 14
T
Based on the equation above, the quantities (Kσ) and (k/σ) are assumed to be fairly constant within the watershed.
Thus knowing the values of K from existing boreholes and the conductivity, the inverse of resistivity from the
sounding interpretation around the borehole, one can estimate the transmissivity and its variation from place to
place using the Dar zarrouk parameter determined for each aquifer. On the whole, σ or 1/5 must be equivalent to
any of the sets kσT or KS or KS/5 to be used. For aquifer located within this study area, the mean of the hydraulic
conductivity K=8.00m/day, while the mean transmissivity Tr= 156.53m2/day the maximum Tr=384.80m2 /day while
the minimum Tr= 20.80m2/day. The values obtained by Uma (1989) for the Ogwashi – Asaba Formation are 1.00 to
2
14.00m/day for K and 54.98 to 702.60m day for T. For the aquifers within the Benin Formation, Kmean =
2 2 2
23.07m/day, Tmean =704.72m /day; Kmax = /day; Tmax=2036.21m /day; Kmin = 7.02m/day; Tmin = 87.05m /day. Uma,
(1989) obtained K values for Benin Formation that range from a minimum of 4.90m/day to a maximum of
43.99/day with a mean of 19.71m/day. This appears to be in range constant with the values obtained in this study.
The highest value of T was obtained about 200m from MOUAU gate towards the right, near the mechanic
workshop. This indicates that this location has the greatest potential for productive aquifers.
www.gjournals.org 75
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
www.gjournals.org 76
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
7.53 7.53 7.53 7.53 7.53 7.54 7.54 7.54 7.54 7.54 7.55
5.49 Transmissivity
5.49 360.00
340.00
L a t itu d e (d e g r e e )
320.00
5.48 300.00
280.00
260.00
5.48 240.00
220.00
200.00
5.48 180.00
160.00
5.48 140.00
120.00
100.00
5.48 80.00
60.00
40.00
7.53 7.53 7.53 7.53 7.53 7.54 7.54 7.54 7.54 7.54 7.55 20.00
Longitude (degree)
Fig 2: Contour Map of Transmissivity in the Study Area
Figure 3 show the contour map of transverse resistance which clearly reflects a high value of 194592.32Ωm2
recorded for the Bishop’s house off Lihman Hall at the back of the MOUAU stadium.
www.gjournals.org 77
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
5.49
L a t itu d e ( d e g r e e )
5.48 220000.00
200000.00
5.48 180000.00
160000.00
140000.00
5.48
120000.00
100000.00
5.48 80000.00
60000.00
5.48 40000.00
20000.00
7.53 7.53 7.53 7.53 7.53 7.54 7.54 7.54 7.54 7.54 7.55 0.00
Longitude (degree)
Fig 3: Contour Map of Transverse Resistance
Table 2: Result of groundwater quality analysis of samples collected from selected portions of Kwa Ibo
river.
www.gjournals.org 78
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
Hydraulic head or piezometric head is a specific measurement of water pressure above a geodetic datum. It is
usually measured as a water surface elevation, expressed in units of length, at the entrance (or bottom) of a
piezometer. In an aquifer, it can be calculated from the depth to water in a piezometric well (a specialized water
well), and given information of the piezometric elevation and screen depth. Water entering an unconfined or
confined well will stand at a particular level. This level is the hydraulic head and is actually the sum of three
components, the pressure head, elevation head and velocity head. Therefore the distribution of hydraulic head
through an aquifer system determines where groundwater will flow. Where the hydraulic head is constant, there is
no flow. However, if there is a difference in the hydraulic head from the top to the bottom due to draining from the
bottom, the water will flow downward, due to difference in heads, also called the hydraulic gradient.
Table.3: Result of groundwater quality analysis of samples collected from selected boreholes in the
study area
To determine the direction of flow of an aquifer in the study area, groundwater level measurement; relative
geographic position of the wells and elevation were collected and contoured in map perspectives, (using surfer 6
and 8 software), thereby generating groundwater surface maps. The direction of groundwater movement can be
understood in the fact that groundwater always flows in the direction of decreasing head. From the contour map
(showing direction) from samples collected from for Kwa Ibo River fig 2, it can be inferred that the water flows
downwards from the river(North) with a higher hydraulic head (elevation) towards the well (south) located within the
premises of MOUAU.
www.gjournals.org 79
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
Elevation-Kwa river
5.49
Elevation -borehole
5.49 165.00
160.00
5.48 155.00
Latitude (degree)
150.00
145.00
5.48 140.00
135.00
130.00
5.48 125.00
120.00
115.00
5.48 110.00
105.00
5.48 100.00
95.00
90.00
5.47
5.47
7.54 7.54 7.54 7.55 7.55
Longitude (degree)
www.gjournals.org 80
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
Similarly, the elevation map/direction of flow from modeling of water samples collected from borehole wells located
within the MOUAU campus and its environs indicates that water flows towards the south west direction interacting
with the groundwater from the Kwa Iboe River fig4.14
The rate of movement on the other hand is dependent on the hydraulic gradient, which is the change in
head per unit distance.
The modeling studies when compared with the contour maps from Kwa Ibo river water samples and that of
the borehole water samples reveals that there is a river water- aquifer water interaction at a particular geographic
location within the study area.
The water- table elevation in the water shed along the northern boundary was defined by the 97m (a. m. s.
l.) equipotential line at Government College in the north, whereas the lowest water -table elevation was 34m at
Usaga Elegu in the south (Igboekwe 2008).The top layer mostly consists of 5-120 m topsoil/laterite ,underlain by a
-1 -1
sandy zone 68-325m thick. The aquifer permeability varied from 3.15mday to 14.4m day in the water shed.
Higher permeability prevailed along the Kwa Ibo River course.
Table 3 shows the magnitude of loaded suspended particles in the borehole water samples analyzed in the
laboratory. The suspended particulates in the analyzed borehole water samples ranged from 100 to 1200mg/L with
an average value of 414mg/L. In comparison with the WHO standard, this pollutant seems to fall within the
acceptable WHO standard. The groundwater within the indicated depths and locations in the table does not seem
to be polluted/contaminated. The suspended solids particles from river water samples were also investigated and
the result is shown in table 2
The value of the suspended solids analyzed from Kwa Ibo River water samples ranged from 100 to 900
mg/L and the average value was found to be 518mg/L. Based on the flow interaction between the borehole water
samples and the Kwa Ibo River samples, the borehole water seems to be polluted with the highly loaded
suspended particles present in the Kwa Ibo River water samples.
By applying surfer 8 software, contour maps showing the distribution of the suspended solids were
obtained as shown in figs 8 and 9 , for borehole water samples and fig. 10 and 11 for Kwa Ibo River water
samples.
Kwa Ibo River according to the elevation contours in fig. 4 flows from north to south and in event acts as
the main recharge of the selected boreholes in Michael Okpara University. The fairly high distribution of the
suspended particles in Kwa Ibo River water samples is responsible for the high value of suspended particles in
borehole water samples.
5.48 900.00
Latitude (degree)
800.00
5.48 700.00
600.00
5.48
500.00
5.48 400.00
300.00
5.48 200.00
100.00
5.47
5.47
7.54 7.54 7.54 7.55 7.55
Longitude (degree)
www.gjournals.org 82
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
1000.00
900.00
800.00
700.00
600.00
500.00
400.00
300.00
200.00
Fig 8 and 9: 2-D and 3-D contour map distribution for suspended solids in Borehole water
850.00
5.48 800.00
750.00
700.00
Latitude (degree)
5.48 650.00
600.00
550.00
5.48 500.00
450.00
400.00
5.48 350.00
300.00
250.00
5.48 200.00
150.00
100.00
5.48
Fig 10 and 11: 2-D and 3-D contour maps distribution of suspended solids for Kwa Ibo River
www.gjournals.org 83
Greener Journal of Physical Sciences ISSN: 2276-7851 Vol. 2 (3), pp. 070-084, September 2012.
The summary of the geophysical data and the complementary laboratory analysis for boreholes and the Kwa Ibo
River are respectively shown in table 1, 2, and 3. The estimation of elevation from geophysical survey and within
the Kwa Ibo River and boreholes has enabled us to determine the pattern and direction of flow of the groundwater.
This was achievable by employing the idea that groundwater flows from topographically high elevation to a
topographically low elevation. In general water flows from the northern zone (Kwa Ibo River) with higher hydraulic
head to the southern zone (MOUAU) with somewhat low hydraulic head. The interpretation is that water from the
river recharges the boreholes thereby causing interaction due to the contaminant loads or plumes. This interaction
can be inferred from the elevation contour maps shown in Fig. 4 and 5 and the direction contour maps in 6 and 7.
In general, the combination of the geophysics and geology of the area indicated that it has lateritic topsoil,
sandy soil, medium to coarse grained sand, fine to medium grained sand, medium grained sand, clayey sand,
shaley sand and river sand with various thickness at various depths. These primary parameters were used to
determine the secondary parameters such as Dar Zarrouk parameters (Transverse resistance and longitudinal
conductance) and electrical conductivity. With the help of the unique hydraulic conductivity determined by
Igboekwe etal; 2008. The contour maps of some of these parameters were drawn to show the extent of
distributions in the mapped area. The higher values of resistivities of the aquifers indicate that the aquifer itself in
not contaminated. Moreover, the lower values of Kσ indicate that the water is not brackish. However, the observed
contamination is due to the steady flow and interaction between the boreholes and the Kwa Iboe River that has the
loaded contamination plumes.
The analysis of the borehole water at the MOUAU gave values for suspended solids, total dissolved solids,
Turbidity, Conductivity, PH, and total solids. Temperature and elevation at different locations were also determined
(see table 3). According to W.H.O standard the values of the suspended solids in mg/L indicates that the borehole
at college of Natural and Applied Science (CNAS) and Zenith bank are somehow contaminated since the value is
beyond the safe standard value. Moreover, the dissolved values of solids at CNAS, Zenith bank and VC’s lodge
indicates that the borehole is also contaminated. However, the reduced value of total suspended solids at the VC’s
lodge should not be a surprise as the VC’s lodge is about 2km away from the Kwa Ibo River.
REFERENCES
Anderson, M. P. and Woessrev, (1992) Applied Groundwater modeling. Academic Press, Inc., San Diego, CA., P
381.
Bear, J., (1979) Hydraulics of Groundwater. McGraw-Hill, New York; pp 569.
Domenico, P.A., Schwartz, F. W., (1998). Physical and chemical Hydrogeology. John Wiley & Sons, New York [2nd
Ed. J.pp 506].
Freeze, R. A. & Cherry, J. A. (1979). Groundwater Eagle Wood Cliffs, New Jersey, Prentice-Hall, 335-336.
Gelhar, L.W. (1993). Stochastic subsurface hydrology. Prentice Hall, Englewood Cliffs. Pp. 390.
Griffiths, D. H., (1976). Application of electrical resistivity measurements for the determination of porosity &
permeability in sandstones. Geo-exploration, 14, 207 – 213.
Igboekwe, M. U., Udoinyang , I. E, (2011). Modelling of suspended particulate matter by kwa Ibo River, Umudike,
South-eastern Nigeria. Implications for Pollutant dispersion. International Journal of Water Resources and
Environmental Engineering Vol. 3(6 ), pp.109-116
Igboekwe, M.U., Okwueze, E.E.and Okereke, C.S.(2006)" Delineation of Potential Aquifer Zones from Geoelectric
Soundings In Kwa Ibo River Watershed, South-easten, Nigeria", Journal Of Engineering and Applied Sciences,
Vol. 1,No.4,pp 410-421.
Konikow, L. F., Grove, D. B., (1977). Derivation of Equations Describing Solute Transport in Groundwater. Us Geol.
Survey Water- Res. Inv. 77-79; pp 30.
Mbonu, D. D. C., Ebeniro, J. O. Ofoegbu, C. O.& Ekine A. S. (1991). Geo-electrical sounding for the determination
of aquifer characteristics in parts of the Umuahia area of Nigeria of Nigeria. Geophysics 56 (5), 284 – 291.
Reddell, D.L., Sunada, D. K., (1970). Numerical Simulation of Dispersion in Groundwater Aquifers. Colorado State
University, Ft. Collins, Hydrology Paper 41. Pp 79.
Reilly, T.E., Franke, O.L., Buxton, H.T., Bennet, G.D. (1987). A conceptual framework for ground-water solute –
transport studies with emphasis on physical mechanics of solute movement. U.S. Geol. Survey Water – Res.
Inv. Rept. 87-4191. Pp. 44.
Wilber, C. G., (1983). Turbidity in an aquatic environment an environmental factor in fresh and oceanic waters.
Charles C. Thomas, Spring field IL.
www.gjournals.org 84