Journal of Environmental Protection, 2011, 2, 903-914
doi:10.4236/jep.2011.27103 Published Online September2011 (http://www.SciRP.org/journal/jep)
Copyright © 2011 SciRes. JEP
903
Multivariate Chemometric Analysis of a Polluted
River of a Megalopolis
Alejandro Gabriel García-Reiriz1*, Jorge Federico Magallanes2, Marjan Vracko3, Jure Zupan3,
Silvia Reich4, Daniel Salvador Cicerone2,5
1Department of Analytical Chemistry, Faculty of Biochemical and Pharmaceutical Sciences, National University of Rosario, Rosario
Institute of Chemistry (IQUIR-CONICET), Rosario, Argentina; 2Gerencia Químic a, Comisión Nacional de Energía Atómica (CNEA),
Buenos Aires, Argentina; 3National Institute of Chemistry, Ljubljana, Slovenia; 4Universidad Nacional de General San Martín. San
Martín, Argentina; 5Instituto de Investigación en Ingeniería Ambiental – 3iA, Universidad Nacional de General San Martín,. San
Martín, Argentin.
Email: garciareiriz@iquir-conicet.gov.ar
Received April 28th, 2011; revised June 27th, 2011; accepted August 8th, 2011.
ABSTRACT
A chemometrical study regarding a 10-years water quality monitoring plan a t 15 sampling points along a section of the
Reconquista River and its stream channels, which embraces 21 campaigns, is presented. The original data were
pre-treated in or der to eliminate missing data and outliers , obtaining a final data matrix of 270 samples containing 26
physical-chemistry variables each. Multivariate statistical methods like multi curve resolution, canonical correlation
analysis and factor analysis methods, as well as current univariate statistics were applied. The interpretation was sim-
plified when varia bles were separated in groups containing environmentally and chemically related variables instead of
analyzing them all together. These methods have shown that the presence of metals likely come from at least 3 different
type of sources. Although the stream channels arriving to the main river course are highly polluted, their flow rates are
so low that do not significantly decrease its water quality. They mainly contribute to the high levels of biochemi-
cal-oxygen demand and chemical-oxygen demand as well as nitrogen-content species. Furthermore, regarding metals,
the pollutants coming from the upstream of the river is higher than those introduced by all channels.
Keywords: Surface Water, Water Pollution, Chemometric, Multivariate Statistic
1. Introduction
It is clear today that fresh-water flows are an essential
human resource for the mankind. About 20% of normal
flow of world’s rivers is used by humankind’s necessities
like agriculture, hydropower and domestic use. As a
consequence, several activities water-dependent re-
sources, like agricultu ral and industrial as well as heavily
dense populated cities with inadequate sanitary infra-
structure, have seriously compromised the quality of
surface waters because agricultural runoff, the introduc-
tion of heavy metals by industrial untreated waste, or-
ganic persistent compounds and pathogens. Furthermore,
the shanty towns waste the dirty water by channelling it
into the sewer to the surface water flows. The discharge
of untreated wastewater is a growing environmental
concern, with many rivers being turned into open sewers.
Only about 10 percent of wastewater in developing
countries is collected and only about 10 percent of exist-
ing wastewater treatment plants operate reliably and effi-
ciently [1]. As a consequence, many downstream popula-
tions receive water of low quality, unsuitable for domes-
tic use and protection of aquatic life.
When river water approaches to the lowlands their
transported contaminants impact the low basin area, par-
ticularly in the coastal zones, where more than half of the
world’s population lives. As a matter of fact, 12 out of 17
world’s mega cities, whose populations are grater than 10
millions are located at coastal zones, confirmed Envi-
ronment Matters (2003).
In Argentina, Buenos Aires City and its suburbs is one
of the coastal megalopolis having more than 11 millions
of inhabitants. The Matanza-Riachuelo and Reconquista
rivers are both the two most polluted rivers of the district.
Both of them flow into Rio de la Plata Estuary. A local
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis
904
study [2] reported that this estuary is 30.212 km2; with a
river flow rate of 18.000 m3/s; being the source of drink-
able water for Buenos Aires’s inhabitants.
The middle basin of the Reconquista River is being
studied by interest of the San Martín Municipality (see
Figures 1-3). Three affluent channels are located in this
sector, which collect non treated effluents and solid
waste coming from industries and populated irregular
settlements. These are low flow channels in comparison
with the river, but they are highly polluted.
The Arroyo Moron stream comes from highly industrial-
ized areas; it is larger than the channels and tributes to
the Reconquista Rive r next to the initial p art of the sector
under study. Another remarkable fact is the presence of
two landfills on the shore of this sector, one of them is
presently being filled and the other was closed in January
Figure 1. Location of Reconquista River in Bue nos Ai r es Province. It flows to Río de la Plata River.
Figure 2. Close sight of the sector under study in the Reconquista River.
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis905
Figure 3. Detailed sight of the sector under study show ing sample points and main geographical features.
31, 2004, but it is still biologically and chemically active.
Environmental systems are rather complex to interpret
because the amount of different variables required to be
taken into account and the complicate relationships
among them. These characteristics determine that the
data have to be processed by chemometrical techniques
and particularly by multivariate statistical methods. The
description and interpretation of environmental systems
requires the development of mathe matical models of two
different kinds. One of them is the mechanistic model [3],
which try to describe the system applying a combination
of known basic scientific laws, like mass transport,
chemical equilibrium, thermal and mechanical convec-
tion, etc. The other one, which is of our specific interest
here, applies statistical and computational methods based
on different strategies, but with the same objective. This
methodology of analysis includes several multivariate
statistical techniques like clusters, principal components
and factor analysis, partial least squares, artificial neural
networks, etc. All these methods are widely described,
for instance, in the following References: [4-7].
This work presents a first water-pollution analysis of
the Reconquista River and its present state. The concept
of studying a partial section of the river instead of the
complete basin meets the requirement of intensive sam-
pling with a low budget.
2. Sampling and Data Base
Samples were collected during a period of 10 years,
1995-2005. Sample points are shown in Figure 3. Two
or three sampling campaigns were carried out every year.
Samples were taken at all sites during the same day.
After elimination of variables with few measurements
or occasional sample points with lack of information, we
obtained a complete data base matrix M270×26 of 270
samples with 26 measured variables each, which are
shown in Table 1.
Analytical chemistry methods were carried out under
appropriate normalized procedures for every variable.
Most variables have a lognormal distribution and a
few of them are normally distributed. Those variables
values recorded as zero were changed to the average be-
tween the detection’s limit of th e method and zero. When
comparisons among them altogether were necessary, they
were scaled to 0 - 1 range using a minimax algorithm:
min
max min
s
xx
xxx
(1)
where xs is the scaled value, x is the original value, and
xmin and xmax are the extremes of its range. For instance,
Figure 4 shows spreads of all 26 variables. It would have
been impossible to show them altogether with their
original units because of their different scales. This fig-
ure also shows the extreme values of many variables,
which strongly affect the results of calculus, thus needing
them to be analyzed as outliers or not. Different criteria
have been used to remove (or not) extreme values shown
in Figure 4: robust Mahalanobis detection [8] and pro-
fessional judgmental. Further analyses have been made
on this data base, either scaled or with its original unities.
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis
906
Table 1. Physical and chemical measured variables and
used acronyms.
Variable Acronym
Alcalinity Al
Arsenic As
Cloride Cl
Conductivity Cl
Biochemical oxygen demand BOD
Chemical oxygen demand COD
Methylene blue active substances MBAS
Hardness
Phosphorous- as ortophosphate
Nitrogen- as amonium N_AMON
Nitrites NO2
Nitrates NO3
Organic nitrogen N_ORG
Total Kjeldahl nitrogen NTK
Disolved oxygen DO
pH
Total residues after evaporatión TR
Total disolved solids
Total suspended solids
Sedimentable solids at 10' SS10H
Sedimentable solids at 2 hs SS2H
Sulfate
Sulfide S
Phenols
Turbidity
Cadmium Cd
Zinc Zn
Copper Cu
Chromium Cr
Iron Fe
Manganese Mn
Nickel Ni
Lead Pb
3. Chemometric Methods of Analysis
Multivariate techniques are methods of analysis gener-
ally recognized as very well known tools to study envi-
ronmental problems [9,10]. From this kind of methods,
Multivariate Curve Resolution (MCR) [11,12] has been
selected as one of the most advantageous to have a close
look of our system. To apply the MCR method, the
minimax-scaled-matrix data base M270×26 was reorgan-
ized by assembling every sampling campaign as a
sub-matrix, which includes all its sampling sites as row
vectors. These row vectors embrace the 26 physi-
cal-chemistry parameters measured. Not necessarily all
these sub-matrixes have the same number of rows, but
the same number of columns. Twenty one of these
sub-matrixes were stacked to form the matrix D270×26,
which contain the data spatially and temporally organ-
ized. Then, the matrix D allows MCR to analyze the d ata
in space and time dimensions.
MCR is a bilinear based method of inverse calibration;
it can be basically described as a matrix decomposition:
T
DSL E (2)
where D is th e data matrix, S is the scores matrix related
with the objects and L is the load ings matrix related with
variables. Every vector of S is associated with a vector of
L through a product that represents a component. It is
supposed that each component represents a kind of
source (or combination of similar sources) which con-
tributes to spoil the system. A sketch of the calculus is
shown in Figure 5. The bilinear model in Equation (2)
assumes that the major sources of the experimental data
variance can be explained by a small number of compo-
nents defining the two reduced-size factor matrices
(scores and loadings). The model described by this equa-
tion assumes that the measured concentration of a con-
taminant (variable) in a particular sample is the sum of a
reduced number of contributions of this contaminant
coming from different sources. It is therefore a mixture
analysis problem with unknown sources which have to
be estimated from the analysis. Since the solution of
Equation (2) is ambiguous, the matrix decomposition in
this equation has to be performed under some constraints.
The decomposition of Equation (2) is similar to Principal
Components Analysis (PCA), but PCA decomposition is
performed under orthogonal constraints, loadings nor-
malization and maximum explained variance for the
successive extracted components. Under these con-
straints, PCA provides unique solutions. However, these
solutions are an abstract linear combination of the true
experimental variance sources and, although they are
very useful for data exploration and summary, in many
cases they can be too complicated in terms of environ-
mental interpretation. Although there are many good
textbooks about PCA, we address to the interested reader
to Jolliffe [13]. Unlike PCA, the matrix bilinear decom-
position performed by MCR alternating-least-squares
(ALS) uses softer natural constraints and as a result, the
interpretation of loading an d score p rofiles are more easy
and reasonable from an environmental point of view
[9,14]. Constraints used in this work during the MCR-
ALS bilinear matrix decomposition were non-negativity
and normalization of loadings to equal length as those
used in previous works [14].
In this study another applied method was canonical
correlation analysis (CCA). It is a way of measuring the
linear relationship between two subsets of multidimen-
sional variables. It begin when the original data set of n
variables and p objects are grouped into two data sets,
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis 907
Figure 4. Spreads of all 26 variables scaled according to minimax algorithm. See equation 1 in the text.
Figure 5. Sketch of the MCR calculus (see te xt).
X(p × nX) and Y(p × nY), where n = nX + nY, The method
searches for a linear combination of the variables in X
and another linear combination of the variables in Y to
optimize the correlation between both sets. Then, the
canonical correlations measure the strength of associa-
tion between the two sets of variables. Hotelling devel-
oped CCA [10], becoming then a standard tool in statis-
tical analysis, applied to economics, medical, and eco-
logical studies [15,16].
CCA proposes new va riables U(p × nX) and V(p × nY)
presented in Equations (3) and (4).
UXA (3)
VYB (4)
The matrix A (nX × nX) and B (nY × nY), defining the
transformation are chosen in orther that the correlation C
(U,V) is maximum. Thus, it is possible to find new vari-
ables as combinations of the original ones, which reveal
existing correlations between two different sets. Despite
the fact that canonical variables are arbitrary, they can be
then interpreted from the previous knowledge of the sub-
ject matter. Other methods, like principal component
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis
908
analysis (PCA) and univariate statistic, have been applied
to check particular results.
4. Results and Discussion
Most of the chemical variables considered in this study
overpass the concentration limits established as protec-
tion of aquatic life and consequently the water quality of
the river is in general of poor quality due to pollutants
coming from untreated domestic and industrial effluents
poured directly on the river. The distribution of metal-
cations concentrations along the river is shown in Fig-
ures 6 and 7. Figure 6 cannot shows the extreme values
(mg·L–1) of Zn (2.94), Cr (1.91), Ni (0.20) and Pb (0 .11),
but those of Cd and Cu (0.05 and 0.34 mg·L–1 respec-
tively). Fe and Mn have not been considered here be-
cause they have high concentrations due to their geo-
logical origin.
In these figures can be seen an increase of cation con-
centrations from RVI to RIII sampling points and then a
progressive fall. The rational explanation for the incre-
ment of concentrations is the entrance of the Moron
stream into the Reconquista River, with significant high
values of dissolved salts. High levels of conductivity
induce a compression of the double layer of water sus-
pended particles,, thus decreasing the repulsion among
them and consequentelyfavouring their coagulation and
further precipitation in the main course [17] with the
concomitant lowering of cation concentrations. The ex-
istence of a deposit of metals in the sediment between
points RII and RI has been checked. This behaviour has
been explained [18] through a water quality model that
contemplates the adsorption of metals onto suspended
particles and the precipitation into the sediments before
reaching RI. Moreover, the average of sum of cations
(without Fe and Mn) for the sampling points RV (en-
trance point) and RI (exit point) have similar values
without statistically significant differences (0.535 and
0.429 respectively); this means that the cations coming
from channels do not produce a significant increment of
metals content in the water quality of the river at RI.
Figure 7 shows the role of Fe and Mn. They are of
geological origin and follow the same path than cations
of different sources reinforcing the precipitation mecha-
nism explained in the previous paragr aph.
The interpretation of MCR’s results is simplified if in-
stead of analyzing the variables all together, they are
separated in groups containing environmentally and
chemically rela ted variables. W ith this purpose we formed
three groups of variables described in Table 2. Results for
the first group are shown in Figure 8(a) (singular value
decomposition (SVD) of the data matrix) and Figure 8(b)
(level of significance for each variable according to their
loading values). By selecting the first three factors of the
SVD allow to reach 58.57% of the cumulative variance.
The first row of Figure 8(b) shows a strong relationship
among metals, S and SS2H, this is a usual expected result
because metals that easily combine with sulphide form
very insoluble compounds, ending up associated to solid
suspended particles. Cappari [19] and Nader [18] have
arrived to similar conclusions in t he same site.
Figure 6. The distribution of metal-cation concentrations along the river in the sector under study. Fe and Mn are not shown.
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis909
Figure 7. MCR result for metal-related variables.
Figure 8. MCR result for nitrogen-related variables.
Table 2. Groups of variables for MCR analysis.
Group Variables
A
Metals related SS2H, Phenols, S, Cd, Zn, Cu, Cr, Mn, Ni, Pb
B
Nitrogen related MBAS, NO2, NO3, N_ORG, NTK, Phenols
C
Oxigen related Alkalinity, BOD, COD, DO, S, Phenols
(d)
(c)
(b)
(a)
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis
910
The second row of Figure 8(b) shows that Cd and Pb
are also associated, indicating a possible different source
(or sources) for these contaminants, phenols are included
in this group. The third row of Figure 8(b) shows the
contribution of the sed iment as a source of pollution. We
did not consider Fe because it is of geogenic origin and
its high concentrations are permanent around all the area.
Its levels correlate with Mn concentrations, demonstrat-
ing that the latter one has the same origin, although with
lower concentratio ns than Fe.
Figure 9 shows the MCR profiles of the nitrogen re-
lated variables: in this case three factors of the SVD
represents 80.49% of accumulative variance. The sig-
nificant contribution of MBAS, NO2, N_ORG, NTK and
in low proportion phenols, are shown in the first row of
Figure 9. These variables reveal the presence of reduced
forms of nitrogen species as well as other species like
phenols and surfactants. By relating the object scores
with the sampling sites, it is possible to show that this
component is mainly present in the lowest stream of the
channels. The same procedure helped us to determine
that the oxidized chemical species of nitrogen, NO3,
appears in the upper stream of the channels. Moreover,
the most reduced forms like NH3 and NTK are present in
the river (see Figure 9).
The third group of variables, those related with oxygen,
reinforces the previous result. It is shown in Figure 10
that dissolved oxygen is significantly present in the
channels (upper stream) whereas a reductant media exist
in the main river course, according to the weights of
BOD, COD, S and phenols.
Canonical correlation has been used to explore other
class of relationships, those that would be established
between cations and variables related to organic material.
Furthermore, we want to know if this possible relation-
ship is the same at the river and at the channels. Within
this frame, the X set includes variables of the B and C
groups of Table 2 (DBO, DQO, MBAS, N_ORG and
NTK); and the Y set includes metal concentrations (Cd,
Zn, Cu, Cr, Mg, Ni and Pb).
Figure 11(a) shows a high correlation between the X
and Y groups of variables with reference to the river,
whereas for the channels, a low correlation is observed
between the same groups of variables (Figure 11(b)).
(a)
(b)
(c)
(d)
Figure 9. MCR result for oxygen-related variables.
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis911
(a)
(
b
)
Figure 10. Correlation between groups of variables. Probably-complexating-agents and metallic (see text, Equations (5) and
(6)). At the channels (a) and at the river (b).
(a)
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis
Copyright © 2011 SciRes. JEP
912
(b)
Figure 11. PCA loadings show correlations at the river (a) and at the channels (b).
The structures of the U and V values for the river (UR
and VR) are :
UR = 1.23% DBO +1.91% DQO + 3.49% MBAS
+ 66.53% N_ORG + 26.83%NTK (5)
VR = 3.54% Cd + 1.85% Zn + 1.30% Cu + 0.01% Cr +
0.07% Mg + 90.14% Ni +3.07% Pb (6)
Since CAA also maximizes the correlation between
each group of variables, it can also be highlighted that
N_ORG and NTK are the most significant nitrogen re-
lated chemicals and, among metals, Ni contributes with
the main part of the total variability.
A different behaviour is observed for the channels,
showing a significant proportion of uncorrelated pres-
ence of chemicals associated with nitrogen and metals
respectively. When the structure of the U and V variables
for the channels (UC and VC) is analyzed for the first
canonical components we obtain:
UC = 6.37% DBO +0.77% DQO +0.01% MBAS +
40.22% N_ORG + 52.62%NTK (7)
VC = 23.99% Cd + 6.32% Zn + 4.53% Cu + 20.04% Cr +
8.16% Mg + 8.31% Ni +28.72% Pb (8)
The structure of the nitrogen related variables is very
similar for either channels or river sites but, the composi-
tion of the metal variables is completely different. Thus,
a mixture of metals, lightly correlated with nitrogen
variables is found in channels monitored sites while these
nitrogen variables show a strong correlation with Ni
concentrations for the river waters.
We looked for more evidence supporting these results
through principal components analysis (PCA). By de-
veloping two PCA with the same variab les used for CCA,
one for the river and the other one for channels, the ob-
tained results show agreements with those of CCA (see
Figures 12(a) and (b)). The PCA loadings at the river
(Figure 12(a)) show a good correlation among these
variables with exception of NTK, MBAS and Cadmium.
It is remarkable the correlation of all metal cations with
other variables related to organic materials and organic
nitrogen, possibly suggesting the co mplexation of cations
with organic matter. This is consistent with previous re-
sults that found high levels of dissolved metals in the
main course of the river [14]. The channels (Figure
12(b)) do not show the same correlation with organic
nitrogen, reinforcing the results found with CCA. The
difference between channels and river could be explained
because organic nitrogen have much higher median in
the river (11.80 g·L–1) than in the channels (6.85 g·L–1)
while the median of the sum of cations (without consid-
ering iron and manganese) keep almost the same: 0.10
g·L–1 for channels and 0.15 g·L–1 for the river.
5. Conclusions
Univariate and multivariate statistical techniques are both
of importance to analyze a multivariate problem. Com-
plex multivariate systems like environmental ones re-
quire a sequential application of both methods. It is pref-
erable to start with multivariate methods to arrive at gen-
eral and clear conclusions, and then, to check those pos-
sible doubtful points or conclusions wi t h u nivari a te to ols.
Because environmental systems involve many vari-
ables of different origins, like chemical, physical, mete-
orological, geological, etc., it is difficult to tray to relate
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis913
DBO
DQO
SAAM
N_ORG.
NTK
Cd
Zn Cu
Ni Pb
-1.0 -0.50.00.51.0
Factor 1 : 32.16%
-1. 0
-0. 5
0.0
0.5
1.0
Factor 2 : 15.74%
(a)
DBO
DQO SAAM
N_ORG .
NTK
Cd
Zn
Cu
Ni
Pb
-1.0 -0.50.00.51.0
Fac tor 1 : 2 4.98%
-1.0
-0.5
0.0
0.5
1.0
Fa c tor 2 : 1 7 .49%
(b)
Figure 12. Distribution of Fe and Mn along the river.
all variables at the same time. To facilitate the analysis,
the variables should be divided into groups according to
their environmental compatibility, as we exp lained in the
results section.
The river is already highly polluted at the point in
which it reaches the studied sector. Its contamination by
most of the measured pollutants overpasses the limits of
protection of aquatic life. Along the studied sector, the
presence of metals seems to come from at least three
different kinds of sources as the MCR result shows.
Regarding metal-cation concentrations, there are no
significant differences between the ones measured at the
entrance from those determined at the exit of the studied
sector. Even though the Arroyo Morón, one of the con-
sidered streams, significantly increases the concentration
of dissolved salts (and in this way the cation contents) in
the main river course from RVI to RII; a deposit of met-
als has been identified between RII and RI, probably due
to precipitation with sulphides and adsorption to sus-
pended particulate material, all of them falling to the
sediments. This situation requires the management of the
sediments in this area as hazardous waste.
Finally, at their sources, th e channels show compatible
levels of DO, BOD and COD and oxidized species of
nitrogen for aquatic life pro tection. The contamination of
channels increases along their courses, the content of
dissolved oxygen diminishes toward the river because the
presence of shanty towns in the area without sewer ser-
vices and the direct input of untreated industrial effluents.
When channels arrive to the main stream, they end up
polluting the river with organic matter coming from
waste, sewage, etc., although the river is in worse condi-
tion than channels. These chan nels contribute to increase
the pollution of the river, either of its water or its sedi-
ment. The lack of sewers is one of the main causes of this
type of pollution.
6. Acknowledgments
This study is part of an International Atomic Energy
Commission (IAEA) project (RLA /1/010) ‘Improvement
of the management of contamination of surface water
bodies contaminated with heavy metals’. As a comple-
mentary project, a bilateral Argentine-Slovenian agree-
ment (SLO/08/12) was established in order to analyze the
data provided by the water quality sampling plan. Sup-
port has also been received from ANPCyT project PIC07
01216 “Movilización de Contaminantes en Sistemas
Naturales”.
REFERENCES
[1] World Bank, “Environment Matters,” Annual Review
2003, Water Supply and Sanitation, p. 23.
[2] Aguas Argentinas, Obras Sanitarias de Buenos Aires,
Instituto Ringuelet y Servicio de Hidrografía Naval,
“Calidad de las Aguas de la Franja Costera Sur del Río de
la Plata, (San Fernado-Magdalena),” Consejo Permanente
para el monitoreo de la Calidad de Aguas de la franja
Costera Sur del Río de la Plata (ed.), Secretaría de Obras
Públicas, Dirección Nacional de Gestión de los Recursos
Hídricos, anexos I y II, Buenos Aires, 1997.
[3] G. Kiely, “Ingeniería Ambiental,” McGraw Hill, Boston,
1998.
[4] D. L. Massart, B. G. M. Vandeginste, S. N. Deming, Y.
Michotte and L. Kaufman, “Chemometrics: A Textbook,”
Elsevier, Amsterdam, 1988.
[5] D. L. Massart, B. G. M. Vandeginst, L. M. C. Buydens, S.
De Jong, P. J. Lewi and J. Smeyers-Verbeke, “Hanbook of
Chemometrics and Qualimetrics,” Elsevier, Amsterdam,
1997.
[6] J. W. Einax, “Chemometrics in Environmental Che-
mistry-Applications,” Springer-Verlag, Berlin, 1995.
[7] J. W. Einax, H. W. Zwanziger and S. Geiss, “Che-
mometrics in Environmental Analysis,” VCH-Weinheim,
Copyright © 2011 SciRes. JEP
Multivariate Chemometric Analysis of a Polluted River of a Megalopolis
Copyright © 2011 SciRes. JEP
914
Germany, 1997.
[8] S. Verboven and M. Hubert, “Libra: A Matlab Library for
Robust Analysis,” Chemometrics and Intelligent Labora-
tory Systems, Vol. 75, No. 2, 2005, pp. 127-136.
doi:10.1016/j.chemolab.2004.06.003
[9] M. Felipe-Sotelo, L. Gustems, I. Hernández, M. Terrado
and R. Tauler, “Investigation of Geographical and Tem-
poral Distribution of Tropospheric Ozone in Catalonia
(North-East Spain) during the period 2000-2004 Using
Multivariate Data Analysis Methods,” Atmospheric En-
vironment, Vol. 40, No. 38, 2006, pp. 7421-7436.
doi:10.1016/j.atmosenv.2006.07.013
[10] H. Hotelling, “Relations between Two Sets of Variates,”
Biometrika, Vol. 28, No. 3-4, 1936, pp. 321-377.
[11] R. Tauler, “Multivariate Curve Resolution Applied to
Second Order Data,” Chemometrics and Intelligent
Laboratory Systems, Vol. 30, No. 1, 1995, pp.133-146.
doi:10.1016/0169-7439(95)00047-X
[12] R. Tauler, A. Smilde and B. Kowalski, “Selectivity, Local
Rank, 3-Way Data Analysis and Ambiguity in Multivari-
ate Curve Resolution,” Journal of Chemometrics, Vol. 9,
No. 1, 1995, pp. 31-58. doi:10.1002/cem.1180090105
[13] T. Jolliffe, “Principal Component Analysis,” Springer-
Verlag, New York, 2002, p. 26.
[14] M. Terrado, D. Barceló and R. Tauler, “Quality Assess-
ment of the Multivariate Curve Resolution Alternating
Least Squares (MCR-ALS) Method for the Investigation
of Environmental Pollution Patterns,” Environmental
Science & Technology, Vol. 43, No. 14, 2009, pp.
5321-5326. doi:10.1021/es803333s
[15] C. J. F. Ter Braak, “Canonical Correspondence Analysis:
A New Eigenvector Technique for Multivariate Direct
Gradient Analysis,” Ecology, Vol. 67, No. 5, 1986, pp.
1167-1170.
[16] S. Becker, “Mutual Information Maximization: Models of
Cortical Self-Organization Network,” Computation in
Neural Systems, Vol. 7, No. 1, 1996, pp. 7-31.
[17] W. Stum and J. Morgan, “Aquatic Chemistry. Chemical
Equilibria and Rates in Natural Waters,” 3rd Edition,
John Wiley & Sons, Inc., New York, 1995.
[18] G. Nader, “Modelización del Transporte de Metales en el
río Reconquista (entre ex ruta 8 y Panamericana): Etapas
de Conceptualización, Formulación y Calibración,” Ph.D.
Dissertation, Universidad Nacional De San Martín, San
Martín, 2009.
[19] G. Cappari, “Evaluación de la Calidad de Agua
Superficial del Río Reconquista a la Altura del Partido de
San Martín Mediante Estudios Ecotoxicológicos en dos
Oportunidades de Muestreo (2007-2008) y Caracte-
rización de los Sedimentos Obtenidos en el Mismo Tramo
de la Cuenca,” Ph.D. Dissertation, Universidad Nacional
De San Martín, San Martín, 2008.