107-Texto Del Artículo-141-2-10-20160722

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 8

MODELACIÓN HIDROLÓGICA CON ESCENARIOS DE

CONTAMINACIÓN DIFUSA, PARA LA PLANIFICACIÓN


AMBIENTAL EN LA CUENCA DEL RÍO DAVID
Viccelda María Domínguez, Ph. D.
Universidad Tecnológica de Panamá
[email protected]

RESUMEN three different months (October, December de 2006 and


El objetivo principal de esta investigación es calibrar un May 2007), and the calibration of the hydrological model,
modelo hidrológico que permita simular escenarios de that represents the principal equation using by SWAT.
contaminación difusa por plaguicidas utilizados en una Compounds found in river water, were Propanil,
cuenca agrícola, con el propósito de que sea aplicado en Pendimethalin, Endosulfan, Dieldrin y Lindano; these
la Planificación de las mejores prácticas de manejo en la pesticide were analyzed using solid-phase “electrical
cuenca del Río David, una importante fuente de micro conductivity” extraction technique with detection by
abastecimiento de agua potable. Los cultivos de arroz y gas chromatography with mass spectrometry (GC-MS).
café han sido identificados, entre las principales Different scenarios of critical runoff were projected. The
actividades antropogénicas, siendo necesario evaluar los results of the hydrological model calibration have a R2 of
impactos de los plaguicidas utilizados en el cultivo de 0.80 and in the simulation of pesticide the principal
arroz y café, en el sistema acuático. El propósito de este impacts were generated by Pendimetaline and 2,4-D.
artículo es evaluar el fenómeno de transporte de Further monitoring studies are recommended, especially
plaguicidas utilizando el modelo “Herramienta de from different insecticides and their metabolites.
Evaluación para Suelo y Agua” (SWAT). Basados en los
resultados del muestreo de agua, en tres meses diferentes Keywords: Hydrological Model, Pesticides, Metabolites
(octubre, diciembre de 2006 y mayo de 2007), y la
calibración del modelo hidrológico, que representa la 1. INTRODUCCIÓN
principal ecuación utilizada por SWAT. Los compuestos
encontrados en aguas del río, fueron Propanil, La modelación de la contaminación no puntual,
Pedimentalina, Endosulfán, Dieldrin y Lindano; estos se compone de distintas ecuaciones matemáticas,
plaguicidas fueron analizados utilizando tecnología de de modo de representar los procesos de
extracción “microconductividad eléctrica” con detección precipitación-escorrentía, erosión transporte
por cromatografía de gases con espectrofotómetro de masa disuelto y adsorbido de los compuestos
(GC-MS). Diferentes escenarios críticos de máxima
escorrentía fueron proyectados. Se obtuvo una R2 de 0.80 (Maidment, 1993; Novotny & Olem, 1994). De
en la calibración del modelo hidrológico y en la acuerdo a la Organización de las Naciones
simulación de plaguicidas se obtuvo la mayor afectación Unidas para la Agricultura.(FAO, 2006) entre los
por Pendimetalina y 2,4-D. Se recomiendan realizar modelos de simulación de fuentes de
estudios a futuros, especialmente de distintos insecticidas y contaminación no puntuales o difusas más
sus metabolitos que se pueden generar en dichas
combinaciones.
reconocidos, por su importancia en la aplicación
de procesos hidrológicos, sedimentos,
Palabras claves: Modelo Hidrológicos, Plaguicidas, plaguicidas y nutrientes, además de su nivel de
Metabolitos precisión temporal y espacial; se encuentra el
ABSTRACT SWAT (Soil and Water Assessment Tool),
The overall objective of this work is to calibrate a Herramienta de Evaluación de Suelo y Agua,
hydrological model to allow the simulation of different desarrollado por el Departamento de Agricultura
scenarios of no point source pollution by pesticide in an
de los Estados Unidos (USDA) (Arnold et al.,
agricultural watershed to create a methodology for a
sustainable management of an important drinking water 2002). El SWAT integra un número considerable
source in Panama, David River. Rice and coffee culture de submodelos; sin embargo, su fundamento es el
has been identified as the main anthropogenic activity, balance hídrico para determinar la entrada, salida
been necessary to evaluate the impacts of pesticides used y almacenamiento de agua en la cuenca, este
in rice and coffee culture on the aquatic systems. The balance se calcula con la ecuación (1):
purpose of this paper is to evaluate pesticides transport
phenomenon using the Soil and Water Assessment Tool
(SWAT) software. Based on the results of water samples, in

RIDTECRIDTEC
– VOL –6VOL.
No. 16 yNº1
2 y Nº2 17
27
SWt = sw + ∑ (Ri - Qi - ETi - Pi - QRi) en la cuenca del río David, la cual es una
(1) importante fuente de agua potable en la
República de Panamá.
Donde: SWt es el contenido de agua en el suelo
en el día t, sw es el agua aprovechable por las 3. METODOLOGÍA:
plantas o el contenido de agua en le suelo a 0.1- Implementación de Arc View 3.2 en el
bar menos el contenido de agua a 15-bar, t es el Diagnóstico Físico de la Cuenca: De acuerdo al
tiempo en días, R es la precipitación diaria, Qi la Instituto de Investigación del Sistema Ambiental
cantidad de escorrentía diaria, ETi la por sus siglas en inglés ESRI (Environmental
evapotranspiración diaria, P la percolación diaria System Research Intitute, 1997), el
y QRi el flujo de retorno o flujo base; todas las procedimiento para implementar los Sistemas de
unidades en mm. Información Geográfica (SIG); fase preliminar a
la modelación hidrológica; se divide en cuatro
La presente investigación se basa en Maidment etapas, las cuales se definen como sigue:
(1993) y Legates and McCabe (1999); dividiendo • Se referencia y digitaliza la red hídrica de
la cuenca del Río David en Unidades de la cuenca
Respuesta Hidrológica (UHRs) y modela cada • Se referencia geográficamente los mapas
una de las UHRs, permitiendo conocer su o cartas en Arc View 3.2 para así
respuesta al final de cada UHRs, al igual que la digitalizar, y obtener uso de suelo
respuesta al final de la cuenca total (Domínguez (agrícola, ganadería de pastoreo y otros),
et al, 2005). Utilizando las UHRs se puede y tipo de suelo basado en la metodología
realizar escenarios de plaguicidas aplicados en abalada y descrita del Departamento de
distintas áreas de cultivos y su comportamiento Agricultura de los Estados Unidos
en suelo y su posterior llegada al río. Los (USDA) mediante el método de
compuestos encontrados en aguas del río y Bouyoucos.
analizados por un Cromatógrafo de Gases
• Se referencia y digitaliza las curvas de
Agilent Modelo 6890 con receptor de
nivel del área de estudio
microconductividad eléctrica, fueron los
• A partir de las curvas de nivel se genera
plaguicidas Propanil, Endosulfan, Pendimetalina,
el modelo de elevación digital (MDE), el
Dieldrin, Lindano y 2,4 D (Olmos 2007).
cual se comprueba con los resultados de
Diferentes escenarios críticos de máxima la red hídrica, observando la similitud
escorrentía fueron proyectados. Entre los entre las líneas de la red hídrica
resultados de las simulaciones en los distintos digitalizada y las depresiones observadas
escenarios, hay presencia de compuestos en el MDE en tres dimensiones.
carcinogénicos, siendo estos el caso del Dieldrin
y el 2,4-D. Se recomienda implementar el Modelación Hidrológica: Se realiza a través de
Sistema de Gestión de Cuencas en el río David, la recopilación de datos generados a partir del
en el cual en la próxima fase se calibren los SIG y de la base de datos de la estaciones de
plaguicidas encontrados en la cuenca, y además ETESA número 108-03-02 y Rovira, ubicadas en
se relacionen con datos epidemiológicos del la parte cuenca del Río David en Unidades de
Ministerio de Salud. Respuesta Hidrológica (UHRs) modelando cada
una de las UHRs, lo que permite su respuesta al
2. OBJETIVO final de cada una, al igual que la respuesta al
El objetivo principal de esta investigación es final de la cuenca total (Domínguez et al., 2005).
calibrar un modelo hidrológico que permita La información se obtuvo y se introdujo de
simular escenarios de contaminación difusa por acuerdo al Diagrama de Flujo 1.
plaguicidas utilizados en la cuenca alta, media y
baja, con el propósito de que sea aplicado en la
Planificación de las mejores prácticas de manejo
RIDTEC
28 – VOL. 6 Nº1 y Nº2 RIDTEC – VOL. 6 No. 118
y2
Requerimientos del Modelo SWAT-Arc View y Escenarios de

Modelo de Red Punto de Parámetros de Mapas de Precipitación,


Elevación Hídrica Descarga Almacenamiento Tipo y Uso Temperatura
Digital de suelos y
Generadores

Plano de Proyección del Conversión y


Referencia Plano análisis de
base de datos

Simulación y Calibración del Modelo SWAT

Escenarios de Plaguicidas

Diagrama de flujo 1. Requerimientos del Modelo SWAT - Arc View y Escenarios de Plaguicidas.

4. RESULTADOS Y DISCUSIÓN alejan la respuesta a la realidad y posterior


A partir del modelo hidrológico SWAT, se correlación obtenida de la simulación del modelo
obtiene una definición del área de la cuenca del hidrológico.
Río David a partir de la dirección y el flujo
representado según las pendientes del MDE. Como muestra la Figura 1, la correlación para los
Como resultado se obtienen 46 UHRs que años 1994 hasta 1997, es de R2 = 0.73, valor de
conforman la cuenca del río David. correlación aceptable, considerando los datos
faltantes y rellenos realizados en su mayoría para
los años 1994, 1995, 1996 y 1997 según lo
350.00 observado por Oñate (2003). A partir de estos
Caudal Simulado (m3/s)

300.00 resultados se trabaja en base al año 1996, el cual


250.00
y = 0.7159x + 9.2513 posee la mayor cantidad de datos reales. De los
200.00 2
R = 0.7268
150.00
cambios realizados a los parámetros de curva
Caudal Real vs. Caudal
100.00 Simulado número, surlag, tiempo de concentración, flujo
50.00 Lineal (Caudal Real vs.
Caudal Simulado)
base o flujo de recarga subterránea (GWQMN),
0.00 coeficiente de evapotranspiración (ESCO);
0.00 100.00 200.00 300.00
Caudal Real (m3/s)
erivados del análisis de sensibilidad combinado
2k; donde K es el número de parámetros de
Figura 1. Correlación para los años 1994 – 1997. mayor influencia, siendo en total 4 los
fundamentales en cuanto a la ecuación de
De la calibración que se realiza específicamente escorrentía; en el año 1996 se obtiene los
para el año 1996 (288 simulaciones) se obtiene resultados mostrados en la figura 1.2 para los
la Figura 1, la cual muestra la correlación para meses de mayor precipitación.
distintos eventos entre los años 1994 – 1997. Los
eventos utilizados son aquellos que han sido La Figura 2 muestra los cuatro escenarios más
poco intervenidos por rellenados de datos por el representativos de la calibración que permiten ir
método de Pearson, debido a que datos ajustando los parámetros de entrada al modelo
rellenados generan un caudal diferente al real y SWAT para la cuenca del río David, la cual
RIDTEC
RIDTEC– VOL. 6 Nº1
– VOL y Nº21 y 2
6 No. 1929
600.00
0.7981 ≈ 0.8, con un valor de flujo pico de
mayor similitud con el valor real de caudal. Por
500.00
otra parte, segregando de la Figura 2 el último
Qreal
Caudal (m3/s)

400.00
Qsimulado 15
evento de precipitación observado (el cual
300.00 Qsimulado 91 cubre 5 días de evento de lluvia entre el 1 al 5
200.00
Qsimulado 105 de octubre de 1996) se obtiene una correlación de
Qsimulado Final
0.8001 con respecto a lo medido en la estación
100.00
limnimétrica de ETESA (Figura 3).
0.00

En la Figura 3 se observa que el flujo pico ocurre


/2 9 6
/2 9 6
/2 9 6
/2 9 6
/2 9 6
/2 9 6
/3 9 6
/0 9 6
/0 9 6
/0 9 6
/0 9 6
/0 9 6
96
0 9 /1 9
09 19
0 9 /1 9
0 9 /1 9
0 9 /1 9
0 9 /1 9
09 19
1 0 /1 9
10 19
1 0 /1 9
10 19
1 0 /1 9
19
4/

9/

1/

3/

5/
3

5
6
7
8

4
para el día 4 de octubre de 1996; y que los
/2
09

Tiempo (mes/día/año) Qsimulado91 y Qsimuladofinal se asemejan al


Figura 2. Escenarios para la calibración. Qreal. Observándose una tendencia de cercanía
de valores simulados a valores reales. En la
muestra una similitud entre Qreal (caudal real) y Figura 4 se observan las diferentes líneas de
las diferentes simulaciones (Qsimulado) hasta tendencia para las cinco simulaciones realizadas
alcanzar una correlación R2 entre 0.7981 a en la Figura 3.
0.8005, los cuales se mantiene dentro de los Figura 3. Evento del 1 al 5 de octubre de 1996
sugeridos por Torres et al. (2005), Aitken,
(1973), Karvonen et al., (1999). En la figura 2 se
300.00
observa para el evento de precipitación-
escorrentía superficial, que los primeros días del 250.00

evento se mantienen cercano a los valores reales; Qreal


Caudal (m 3/s)

200.00
Qsimulado 15
sin embargo, el flujo de los valores pico 150.00 Qsimulado 91
presentes, no puede ser disminuido de mejor Qsimulado 105
100.00
manera, sin pérdida en la correlación. La Qsimulado Final

interpretación del modelo en cuanto a la 50.00

diferencia del flujo pico, puede ser causada por la 0.00


10/01/1996 10/02/1996 10/03/1996 10/04/1996 10/05/1996
implementación de usos de suelo actuales (año
Tiempo (m3/s)
2005 – 2006), mientras se trabaja en base a
comparaciones de caudales del año 1996. El último evento de calibración (QsimuadoFinal)
en la Figura 3 se encuentra representado en la
Lo planteado sugiere una diferencia de 10 años correlación mostrada en la Figura 4 d).
de desarrollo económico y social, lo cual se Observándose una correlación mayor a 0.8. Se
refleja como mayor área pavimentada, observa en la Figura 4 que para a) se tiene valor
crecimiento urbanístico, usos pecuarios en la de correlación de 0.9210, sin embargo existe un
región entre otros; que inciden en una menor flujo pico superior al flujo real u observado. En
permeabilidad del suelo, aumento en la cambio la d) de la misma Figura tiene una
escorrentía y subsiguiente aumento de caudal en correlación de 0.8005, con un valor pico de flujo
la región. Lo que coincide con autores como de mayor ajuste al valor real de caudal. Los
Mays (1996). Aunando la inconsistencia de los valores de caudal para la gráfica d)
valores de precipitación obtenidos de los datos de QsimuladoFinal son menores a los caudales
ETESA, lo cual disminuye el factor de reales presentes en la Figura 3.
correlación del modelo, que coincide con lo
mencionado por Ndomba & Mtalo (2005), El modelo SWAT muestra una correlación para
Arnold (2002) y Oñate (2003), el ultimo sobre ambos eventos escogidos para la calibración de
calibración del Modelo SWAT en la cuenca R2 ≈ 0.8, esto sugiere que el modelo se ajusta a
binacional del Río Catamayo entre Ecuador y las condiciones de estudio y su calibración para
Perú. Se observa sin embargo una correlación de el año 1996 se da por concluido este proceso.
RIDTEC
30 – VOL. 6 Nº1 y Nº2 RIDTEC – VOL. 6 No. 120 y2
Figura 4. Línea de Tendencia y Correlación por Simulación para el evento del 1 al 5 de octubre de 1996. a) R2
= 0.9210 para Qsimulado15 para el segundo evento; b) R2 = 0.8573 para Qsimulado91 para el ultimo evento; c)
R2 = 0.8558 para Qsimulado105 para el ultimo evento; d) R2 = 0.8005 para QsimuladoFinal para el último
evento.

A partir de los resultados de la calibración de los Comportamien Comportamie Cantidad que


Plaguicida to en Suelo nto en Agua llega al Río (mg)
eventos de lluvia, se procede a realizar los Se adhiere al Poco soluble
escenarios de plaguicidas, con el valor obtenido Endosulfan suelo en agua 0.7e-4 a 8.48 mg
en QsimuladoFinal. El objetivo de ejecutar estas Poca adherencia Soluble en
Propanil al suelo agua 3.76 a 27.85 mg
corridas es encontrar valores que serian menores
Pendimetali Se adhiere al Poco soluble
a los que pueden ser generados en la situación na suelo en agua 1511 a 13360 mg
real, y nos sugirieran la llegada o no al río de Se adhiere al Poco soluble
plaguicidas en una situación similar. Dieldrin suelo en agua 5.2e-4 a 246 mg
Se adhiere al Soluble en
Lindano suelo agua 0.12 a 241.5 mg
Resultados preliminares de muestras tomadas en Poca adherencia Soluble en
distintos sitios de la cuenca del Río David (medio 2,4-D al suelo agua 137 a 619 mg
agua en 6 puntos, y envase de mezcla de Tabla 1. Resultados de Escenarios de
plaguicidas en cultivos), y analizados por un Plaguicidas y su Comportamiento en los Medios
Cromatógrafo de Gases Agilent Modelo 6890 de Transporte.
con receptor de microconductividad eléctrica,
revelan la presencia de los plaguicidas Propanil, Según los resultados obtenidos para simulaciones
Endosulfan, Pendimetalina, Dieldrin, Lindano y de escenarios de plaguicidas bajo condiciones
2,4 D Olmos, (2007). Se realiza la ejecución y establecidas de cantidad del mismo en follaje y
análisis de los escenarios en base a estos enriquecimiento en suelo se obtiene los
plaguicidas de acuerdo a sus características siguientes resultados descrita en la Tabla 1.
físico-químicas.
De acuerdo a lo observado en la tabla 1 se
observa que los plaguicidas de menor llegada son
RIDTEC
RIDTEC– VOL.
– VOL6 Nº1 y Nº2
6 No. 1y2 21
31
Endosulfan y Propanil, ambos tienen Franco Arcillosa, aunque su llegada al río se da a
comportamiento parecidos en cuanto a llegada al lo largo y ancho de la cuenca la cual es de textura
río, teniendo similitud en las UHRs en la cual se Franco Arenosa como se observa en la Figura 4.
obtiene llegada al río. El Endosulfan y el
Propanil tienen comportamiento de transporte 5. CONCLUSIONES Y
distinto, el primero depende del sedimento es un RECOMENDACIONES
plaguicida que se adhiere fuertemente a suelos Se debe comprender claramente, que estos son
de partícula pequeña (textura de franco arcillosa). escenarios posibles, no son escenarios calibrados
Christodoulatos y Mohiuddin, (1996), y el último en campo con datos reales. La intención es
una moderada a alta solubilidad en agua conocer como se proyectaría el transporte, y si
dependiendo de la escorrentía superficial. Un son capaces de llegar los plaguicidas objeto de
comportamiento parecido se observa para los este estudio, bajo el escenario hidrológico
plaguicidas Dieldrin y Lindano, ambos se calibrado, características fisicoquímicas del suelo
adhieren al suelo, sin embargo el Lindano tiene y las condiciones de uso de suelo actuales.
mayor solubilidad en agua que el Dieldrin,
observándose una parecida llegada al río en Entre los plaguicidas simulados en los distintos
cantidades cercanas a 250 mg y un similar escenarios, hay presencia de compuestos
movimiento en UHRs debido principalmente a su carcinogénicos, siendo estos el caso del Dieldrin
moderada a fuerte adherencia en suelos. Son la y el 2,4-D. La presencia de estos plaguicidas en
Pendimetalina y el 2,4-D los plaguicidas de la cuenca ha sido comprobada por Olmos (2007),
mayor presencia y llegada a las aguas del río y su transporte es más que probable por medio de
David basándonos en las simulaciones, la escorrentía simulada por este estudio. Como
observándose una cantidad cercana a las 13360 otra fase de la investigación, se realizó una
para la Pendimetalina y 619 para el 2,4-D. El encuesta socioeconómica en la cuenca, que
comportamiento y transporte de ambos incluía, recabar información de tipo
plaguicidas es completamente distinto, la epidemiológica y de uso de plaguicidas. Los
Pendimetalina se adhiere al suelo y por ende resultados mostrados en Mojica y Rodríguez,
depende del sedimento para llegar al río, por su (2007), arrojaron una incidencia del 19 % de
lado el 2,4-D es soluble en agua y por ende su casos de cáncer en la cuenca del río David. Sin
transporte en la cuenca se debe a la escorrentía embargo, las evidencias no son suficientes para
superficial. relacionar directamente los casos.

Se recomienda implementar el Sistema de


Gestión de Cuencas en el río David, en el cual en
la próxima fase se calibren los plaguicidas
encontrados en la cuenca, y además se relacionen
con datos epidemiológicos del Ministerio de
Salud.

Con este proyecto se comprueba la posible


implementación de diversos modelos
hidrológicos en cuencas con características
volcánicas montañosas con condiciones
meteorológicas e hidrológicas similares o muy
Figura 5. Escenarios de transporte de parecidas a las observadas en la cuenca del río
Pendimetalina en la Cuenca del Río David. David.
Es la Pendimetalina el de mayor llegada al río,
debido en gran parte a la parte baja de la cuenca
del río David, la cual es de textura de suelo
RIDTEC
32 – VOL. 6 Nº1 y Nº2 RIDTEC – VOL. 6 No. 122
y2
Agradecimientos: A la Secretaria Nacional de • Mojica, V., Rodríguez, X. (2007).
Ciencia, Tecnología e Innovación (SENACYT) Modelación hidrológica con escenarios de
por financiar esta investigación a través del contaminación difusa, para la
Proyecto de Fomento a la Investigación y planificación ambiental en la cuenca del
Desarrollo ID 147-2005. río David. Tesis para optar al grado de
Licenciado en Ingeniero Ambiental.
6. BIBLIOGRAFÍA Universidad Tecnológica de Panamá.
• Aitken, A. (1973). Assessing Systematic • Ndomba, P., Mtalo, F. (2005). The
Errors In Rainfall-Runoff Models. Suitability of SWAT in Sediment Yield
Journal of Hydrology 20:121-136. Modeling for Ungauged Catchments: A
• Arnold, J. et al., (2002). Soil and Water Case of Simiyu River Subcatchetment,
Assessment Tool User’s Manual. Texas Tanzania. 3rd International SWAT
Water Resources Institute College Conference.
Station. Texas. • Novotny, V., Ollem, V. (1994).
• Christodoulatos, C., Mohiuddin, M. Prevention, Identification, and
(1996). Generalized models for prediction Management of Diffuse Pollution.Van
of pentachlorophenol adsorption by Nostrand Reinhold. New York.
natural soils. Water Environmental • Olmos, J. (2007). Avances de la Tesina,
Research 68:370-378. Determinación Analítica por
• Domínguez, V.M., Franco, N., Mojica, Cromatografía de Gases de Compuestos
V., Caballero, M., Santamaría, E., Orgánicos clorados de Media a Larga
Rodríguez, X. (2005). Informes de vida media en Agua, Tomando como
Avance (1,2) proyecto Implementación Caso de Estudio la Cuenca del Río David.
de un Sistema de Gestión de Cuencas Trabajo de Graduación para optar al
tomando como base de respuesta, la Título de Magíster Profesional en
unidad ecosistémica del Río David. Ingeniería Ambiental. Panamá, Panamá.
Proyecto presentado a la Secretaría • Oñate, F., Aguilar, G., (2003). Aplicación
Nacional de Ciencia y Tecnología del Modelo SWAT para la Estimación de
(SENACYT). Panamá. Caudales y Sedimentos en la Cuenca Alta
• Environmental System Research Institute. del Río Catamayo. Perú.
(ESRI). (1997). Understanding GIS. The • Organización de las Naciones Unidas
ARC/INFO Method. John Wiley & Sons, para la Agricultura. (FAO). (2006).
Inc., New York, U.S.A. Sistema de Información Geográfica
• Karvonen, T., Koivusalo, H., Jauhiainen, conexión con los modelos de simulación
Palko, J., Weppling, K. (1999). A y aplicación a las materias relacionadas
hydrological model for predicting runoff con la erosión Pág. 1-6
from different land use areas. Journal of • Torres, E., Oropeza, J., Fernandez, D.
Hydrology 217: 253 – 265. (2005). Aplicación del Modelo SWAT en
• Legates, D., Mcabe, G., (1999). la Cuenca el Tejocote Altacamulco.
Evaluating the use of “goodness-of-fit” México.
measures in hidrologic and hydroclimatic
model validation. Water Resourcers
Research 35:233-241.
• Maidment, D. (1993). Handbook of
hydrology. The McGraw-Hill
Companies, Inc., New York, U.S.A.
• Mays, L.W. (1996). Water Resources
Handbook. McGraw-Hill. New York,
U.S.A.
RIDTEC
RIDTEC– VOL.
– VOL6 Nº1 y Nº2
6 No. 1y2 23
33
RIDTEC
34 – VOL. 6 Nº1 y Nº2 RIDTEC – VOL. 6 No. 1 24
y2

También podría gustarte