1.1. Zona de estudio

El área de estudio se encuentra al sur de la provincia de Córdoba, entre los 33° 34' y 33° 40' latitud sur y 64° 34' y 64° 33' longitud oeste (Figura 1) y se eleva en promedio unos 400 metros sobre el nivel del mar (msnm). El área de Coronel Moldes y sus alrededores pertenecen a la provincia geológica argentina denominada Llanura Chacopampeana. Dado que la tendencia general de esta llanura fue marcadamente negativa, su relleno sedimentario está integrado por rocas de edad cámbrico-ordovícica a sedimentos holocenos (^{Chebli et al., 1999}). El sitio de disposición de residuos sólidos urbanos corresponde a una zona de médanos a 4 km al NE de Coronel Moldes.

Figura 1 Ubicación regional y local (de las perforaciones realizadas en el predio del basural) del área de estudio. 

1.2. Características hidrodinámicas e hidroquímicas del acuífero

El acuífero libre, conformado por sedimentos arenosos muy finos-limosos de edad cuaternaria y origen eólico (loess ), constituye la fuente de abastecimiento para la mayor parte de la población rural. La importancia de su estudio y protección se basa en que Coronel Moldes posee una economía basada en el uso exclusivo del agua subterránea. El clima de la zona es Mesotermal Subhúmedo-húmedo y presenta pendientes topográficas promedio del 0.55 % (^{Gomez et al., 2009a}).

Los materiales aflorantes y que componen el acuífero libre en la zona de estudio corresponden a materiales cuaternarios, dominantemente arenas muy finas y limos de origen eólicos (loess ) con diferente grado de cementación carbonática. La homogeneidad de este paquete sedimentario se interrumpe en profundidad con delgados niveles arenosos que no superan los 7-10 m de potencia. Hacia el norte, en la localidad de Bulnes, una perforación alcanzó basamento a los 131 metros, mientras que otra realizada en Coronel Moldes llegó a los 165 m sin alcanzarlo. De los materiales que conforman el acuífero freático de la zona de Coronel Moldes y sus alrededores, se observa que hasta los 36 m de profundidad, paquete que comprometió todas las perforaciones muestreadas, dominan sedimentos areno finos-limosos y limo-arenosos finos, con diferente grado de cementación carbonática, destacándose algunas capas delgadas de tosca a los 1.5 m o una más importante entre los 16 y los 21 m. Si bien se dispone de pocos datos, se estima que los materiales ubicados entre los 49 y 60 m de profundidad corresponden a la base del acuífero y que los limos yacentes por debajo tendrían cierto grado de semiconfinamiento (^{Gomez et al, 2009a}; ^{Gomez et al., 2009b}).

La profundidad del nivel freático se ubica entre los 7.80 y 5.60 m. El menor espesor de la zona no saturada se encontró justo debajo del basural, posiblemente vinculado a la presencia de la laguna artificial donde se vuelcan los efluentes líquidos (Figura 1). El mapa potenciométrico realizado a escala local, en los alrededores del basural municipal (Figura 1) muestra una circulación del agua NO-SE, con líneas equipotenciales que muestran una morfología suavemente ondulada de la capa freática en una zona de recarga local asociada a los médanos. Se midieron gradientes hidráulicos del orden de 0.6 %, y gradientes locales en la zona del basural del orden del 1 % y se calcularon velocidades reales de escurrimiento en esta zona de 0.1 m/día, asumiendo una K de 1 m/día y porosidad efectiva del 10 % (^{Gomez et al., 2009b}; ^{Gomez y Martínez, 2010}).

La composición química del agua subterránea es dominantemente bicarbonatada sódica y bicarbonatada-sulfatada sódica, resultando en general aguas dulces con salinidades menores a 2000 mg/L. La laguna de residuos líquidos urbanos analizada presentó condiciones reductoras, presencia de amonio, alta CE y pH, elevada concentración de Cl^-, Na⁺, HCO₃ ^- y concentración más alta de iones metálicos que las medidas en el agua subterránea. Los valores de DBO₅, DQO, y recuento de bacterias superan los límites impuestos para el vertido superficial de efluentes líquidos (Tabla 2). Actualmente los líquidos urbanos son pre-tratados, pero aún existen efluentes industriales (lácteas y frigoríficos) que son volcados directamente sobre el terreno, sin tratamiento previo (^{Gomez et al., 2009b}). Los modelos hidrogeoquímicos reflejan la influencia de estos lixiviados sobre el agua subterránea debajo del basural y su disminución 30 m aguas abajo. Los modelos obtenidos permiten asociar que el incremento del contenido de CE, Cl^-, Na⁺ y SO₄ ^-2debajo del basural, derivan principalmente del aporte de los lixiviados desde la laguna de efluentes, resultando en condiciones geoquímicas locales propicia para algunas reacciones (intercambio, precipitación, oxido-reducción y disolución).

Tabla 2 Datos de entrada para la simulación del transporte reactivo con elementos mayoritarios y metales en una columna de 220 m y 30 m. 

3.1. Principales procesos involucrados en el movimiento de solutos en el agua subterránea

El agua subterránea moviliza los contaminantes una vez que éstos han atravesado la zona no saturada y han alcanzado el nivel freático. A partir de aquí, el transporte de solutos en el acuífero se ve afectado por procesos hidrodinámicos, químicos y bioquímicos que modifican su concentración y características. En general, las sustancias disueltas no son sustancias inertes y durante el curso de su migración pueden verse involucrados en procesos fisicoquímicos que afecten su movilidad (^{Fetter, 1993}).

En general, se produce una transferencia de masa y no una transformación de los contaminantes durante el transporte. Los principales procesos de transporte son la advección, difusión, dispersión y adsorción. Estos mecanismos afectan y condicionan el transporte de solutos en un acuífero (^{Appelo y Postma, 1996}).

La advección es el proceso por el cual los solutos son transportados a la misma velocidad media del agua subterránea, es decir, por el movimiento de la masa de agua subterránea que fluye. La velocidad real del agua en los poros del acuífero será:

Donde v _D es la velocidad de Darcy, m _e es la porosidad eficaz, que representaría el porcentaje de porosidad interconectada (^{Domenico y Schwartz, 1990}), K la permeabilidad del material y dh /dl el gradiente hidráulico.

El proceso de dispersión hidrodinámica actúa diluyendo los contaminantes y disminuyendo su concentración (^{Zaporozec, 2002}). Es el mecanismo por el cual se producen fenómenos de dispersión o extensión de mezcla o dilución del soluto al moverse en el medio poroso. Consta a su vez de dos mecanismos que ocurren simultáneamente y que son claramente diferenciables, la dispersión mecánica y la difusión molecular.

La difusión molecular en una solución es el proceso por el cual los constituyentes iónicos o moleculares se mueven bajo la influencia de su actividad cinética en la dirección de su gradiente de concentración. La difusión se da incluso en ausencia de movimiento de la masa de agua y sólo cesa cuando se anulan los gradientes de concentración. La difusión cobra mayor relevancia dentro de la dispersión hidrodinámica para pequeñas velocidades de flujo (^{Cerón et al., 2005}).

Este proceso es descrito por la primera ley de Fick, en estado estable (p.ej. la concentración no varía con el tiempo):

La masa de sustancia que se difunde y pasa a través de un área por unidad de tiempo es proporcional al gradiente de concentración, donde F es el flujo difusivo, (m²/s) y el signo negativo indica que el transporte es en la dirección de decrecimiento de la concentración, D _m es coeficiente de difusión efectiva (m²/s) para el medio poroso, C es la concentración de soluto, X es la distancia desde la fuente contaminante y ∂C /∂X es el gradiente de concentración en la dirección x (^{Domenico y Schwartz, 1990}). En general los valores de coeficientes de difusión de una partícula en el agua D _w están disponibles en la literatura y depende de la presión, de la temperatura y de la composición del sistema, los cuales deben ser convertidos a valores de D _m, coeficiente de difusión efectiva, para medios porosos, dado que la difusión de una partícula depende tanto del líquido en cuestión como de la porosidad del sólido.

Donde n es la porosidad (o fracción del volumen de agua en la zona saturada) y τ es la tortuosidad del medio.

La dispersión mecánica, que ocurre como consecuencia de las variaciones locales de la velocidad alrededor de la velocidad media del flujo, es un proceso advectivo y no un proceso químico (^{Domenico y Schwartz, 1990}). La dispersión se refiere a los procesos de mezcla debido a la estructura porosa del medio, que origina variaciones locales del módulo y dirección del vector de velocidad con respecto a su valor medio. Esto a su vez es causa de mezcla y expansión del contaminante (^{Carrera, 1984}).

La principal causa de las variaciones en la dirección y rangos de transportes es el medio poroso no ideal. La variable más importante es la conductividad hidráulica y la dispersión mecánica es generada como consecuencia de la falta de idealidad en todas las escalas (^{Domenico y Schwartz, 1990}).

La dispersión transversal reducirá las concentraciones en todas las direcciones detrás del frente advectivo, mientras que la dispersión longitudinal sólo lo hará en la porción frontal de la pluma (^{Domenico y Schwartz, 1990}). Generalmente, la dispersión longitudinal es mayor a la transversal, resultando en una distribución de la concentración en forma de elipsoide (^{Gelhar et al., 1985}).

La naturaleza del fenómeno de la dispersión y su formulación matemática son temas todavía en discusión. La dispersión está afectada por lo que se llama un Factor de Escala, es decir, los valores de la dispersión varían desde unos centímetros en los experimentos de laboratorio, a algunos metros, para estudios de trazadores, y decenas o centenares de metro para problemas de contaminación regional. La explicación comúnmente aceptada es que el coeficiente de dispersión hidrodinámica longitudinal es una medida de la escala de heterogeneidad que no está incluida en el análisis (y a mayor área, mayor heterogeneidad) (^{Cerón et al., 2005}). La tendencia general del gráfico de ^{Gelhar et al. (1985}) sugiere que la dispersividad longitudinal es alrededor del 10 % de la distancia recorrida:

Parece lógico que grandes dispersividades en campo están dadas por la variación en los tiempos de llegada al pozo de observación, lo que significa que la distribución de la permeabilidad en el acuífero es la principal propiedad que afecta el proceso (^{Appelo y Postma, 1996}).

El coeficiente de dispersión hidrodinámica D _L, en la dirección del flujo, es representado por el producto entre la velocidad del agua en el medio poroso y la dispersividad α, despreciando a menudo la contribución del coeficiente de difusión molecular (^{Domenico y Schwartz, 1990}):

Cuando v tiende a 0, p.ej. en sedimentos arcillosos o de baja permeabilidad, la dispersión hidrodinámica es igual al coeficiente de difusión molecular.

Si se combinan estos procesos en una sola ecuación que describa los cambios en la concentración debido al transporte, dispersión y adsorción, en medios porosos saturados, el transporte químico está gobernado por la ecuación:

Donde C es la concentración del soluto (mol/L), v representa la velocidad de agua en los poros (m/s), D _x el coeficiente de dispersión hidrodinámica (m²/s) en la dirección x, r tiene en cuenta la masa de soluto degradada, adsorbida o aportada por el medio sólido por unidad de tiempo y m es la porosidad.

3.2. Simulación de transporte reactivo unidimensional

La implementación de la modelación geoquímica en este punto pretende servir como otra herramienta más a fin de ajustar el modelo conceptual del acuífero poroso clástico, típico de la Llanura Chacopampeana, poder conocer mejor su funcionamiento, ajustar parámetros de transporte y estimar los mecanismos que gobiernan su respuesta frente a la contaminación.

El estudio de la evolución hidrogeoquímica usando los modelos de transporte unidimensional de PHREEQC (^{Parkhurst y Appelo, 1999}) ha sido implementado por muchos autores. Este programa provee un método numérico eficiente para la simulación del movimiento de soluciones a través de una columna o del flujo unidimensional con o sin el efecto de la dispersión. De este modo pueden modelarse los cambios en la composición de un agua subterránea debidos a la acción conjunta de la advección y dispersión. Los cambios pueden ser expresados en la siguiente ecuación:

Donde C es la concentración de la sustancia química en mol/L, t es el tiempo en segundos s, x es la distancia, en m, D _L es el coeficiente de dispersión hidrodinámica longitudinal, en m²/s; v es la velocidad del flujo, en m/s, y q es la concentración del sólido, en mol/L de agua pura.

Los objetivos de la aplicación de la modelación geoquímica utilizando el programa PHREEQC 2.0 (^{Parkhurst y Appelo, 1999}) se centraron en ajustar el parámetro de dispersividad. Se utilizó al ion Cl^-como trazador, además de elementos mayoritarios presentes en el lixiviado, en el agua del acuífero sin contaminar P1 (Figura 1) y en agua contaminada de P2, simulando en un primer momento el transporte unidimensional entre los pozos P2 y P3. Otro objetivo se basó en la simulación del flujo advectivo-dispersivo utilizando al efluente de la laguna y analizando el comportamiento de los elementos traza metálicos, considerando el valor de dispersividad asumido en el paso anterior.

4.1. Simulación con iones mayoritarios para una columna de 220 m

Se aplicó la modelación de transporte reactivo unidimensional de una solución de agua subterránea contaminada, y posteriormente de lixiviado usando en este caso la composición de la laguna de efluentes líquidos, atravesando una columna de sedimentos saturada en agua del acuífero. De este modo, el pozo P2 y el efluente de la laguna, corresponden a las soluciones que se infiltran en el acuífero cuyas características químicas naturales, no contaminadas, corresponden al agua de P1 (Figura 2). En todos los casos se usaron pares de muestras tomadas el mismo día (Tabla 1). Las fases consideradas en los modelos fueron: calcita, mirabilita, hornblenda, illita, fluorapatita, feldespato potásico (K-Spar), hematita, materia orgánica, CO₂Gas, intercambio Ca/Na y Mg/Na.

Tabla 1 Análisis químico y microbiológico a muestras pertenecientes a la laguna de efluentes localizada dentro del basural. Valores expresados en mg/L y conductividad eléctrica (CE) expresada en μS/cm. 

Figura 2 a) Simulación hidrogeoquímica con iones mayoritarios para una columna de 220 m. b) Simulación con iones metálicos para la misma columna. v = velocidad real del flujo, CIC = Capacidad de intercambio catiónico, Pe = Número de Péclet y Co = Coeficiente de Courant (Co). 

Inicialmente se simuló el flujo advectivo-dispersivo a lo largo de una columna de 220 m. Esta distancia incluye desde P2, justo debajo de la contaminación por la laguna de efluentes y el basural, hasta P3 localizado 30 m aguas abajo del basural. El flujo en el tubo fue dividido en 22 celdas de 10 m de largo cada una. Se consideró una velocidad real media del flujo de 0.1 m/día calculada para la zona del basural y un valor de CIC de 12.4 mEq/100 g correspondiente al valor medio entre los valores de CIC de la zona no saturada y saturada. El valor de dispersividad de 10 m utilizado fue tomado del diagrama de ^{Gelhar et al. (1985)}.

En un primer momento de la modelación el tubo es llenado con la solución correspondiente al agua no contaminada, aguas arriba del basural (en P1), a la que sigue luego la entrada del agua del acuífero contaminada, la correspondiente a P2. La concentración de los distintos iones utilizados en la modelación corresponde al análisis obtenido por ^{Gomez et al. (2009b)}. Para el análisis se consideró lo que ocurre en la última celda (celda 22) correspondiente al punto P3, localizado a 220 m de P2.

En este caso se satisface el criterio de número de Pe alto, de 1 y el coeficiente de Co igual a 1. Según estos valores y las características hidrogeológicas del acuífero, es correcto entonces considerar que para este acuífero la componente advectiva es dominante respecto a la dispersiva.

Los datos de entrada pueden verse en la Tabla 1 y los resultados de la modelación para los iones mayoritarios se presentan en las Figuras 3.

Figura 3 Concentración en mol/L de iones mayoritarios en la celda 22 a medida que P2 fluye en cada paso. En 22 pasos, se desplaza una vez el volumen de poro de la columna considerada. Co /Cf = 1 indica que el fluido original ha sido totalmente desplazado. 

En la Figura 3 puede observarse como Cl^- y SO₄ ^-2 muestran un comportamiento conservativo, alcanzando el equilibrio, esto es la concentración de P2 en P3, luego de 2 y 1/2 volúmenes de poro, lo que hace un tiempo de 15 años. En el caso de los cationes el avance es más lento y no logra alcanzarse en ninguno de los casos la concentración correspondiente a P2, incluso una vez alcanzados los 16 volúmenes de poro. En la Figura 3 puede verse la necesidad de mayor tiempo para que se alcance el equilibrio con Na⁺, Ca⁺², K⁺ y Mg⁺². En la Figura 4 se muestran las variaciones de estos cationes sobre las superficies de intercambio, que bien pueden ser las arcillas. Puede observarse cómo en los primeros pasos de la modelación, antes de la llegada de los contaminantes, el intercambio se produce con la liberación de Na⁺ y entrada de Ca⁺² y Mg⁺² (Figura 4), en coincidencia con los resultados de la modelación inversa realizada con NETPATH (^{Gomez et al., 2009b)}. A su vez puede verse que la variación del ion K⁺sobre las superficies de intercambio es prácticamente nula, involucrando entre 0.002 y 0.004 moles a lo largo de toda la columna.

Figura 4 Concentración de los cationes adsorbidos en mol/L sobre las superficies de intercambio, en la celda 22 a medida que P2 fluye en cada paso. En 22 pasos, se desplaza una vez el volumen de poro de la columna considerada. Aquí la simulación se realizó para un número mayor de pasos. 

Luego de un volumen de poro, una vez que el agua subterránea contaminada ha atravesado toda la columna, la concentración de los cationes aumenta en respuesta a procesos de disolución que ocurren en el acuífero en respuesta a cambios en las condiciones químicas y del aporte de las fuentes contaminantes. Este aumento lleva a que el sistema reaccione con un proceso de intercambio reverso, a partir del cual la gran concentración de Na⁺ que entra al sistema, principalmente derivado del lixiviado de la laguna del basural, comienza a ser consumida de la solución hacia las superficies de intercambio, con la consecuente liberación de Ca⁺² y Mg⁺²a partir de estas superficies y un aumento de ambos en la solución.

Así los procesos de intercambio retardan el avance de estos cationes en el agua contaminada. En el caso del Ca⁺², el aumento en los primeros pasos se vincula a fuertes procesos de disolución de calcita. Luego de los 12 volúmenes de poro la solución se sobresatura en este mineral y la disminución de Ca⁺²pasa a estar controlada por la precipitación de calcita (Figura 5a).

Figura 5 A) Variación del índice de saturación (IS ) de calcita en la celda 22, a medida que P2 fluye en cada paso. En 22 pasos, se desplaza una vez el volumen de poro de la columna considerada. B) Concentración en mol/L de HCO₃ ^- en la celda 22 a medida que P2 fluye en cada paso. 

Para el caso de los bicarbonatos (Figura 5b) se observa un gran aumento de su concentración luego de que la solución contaminada del acuífero atravesado una vez la columna completa. Luego de 12 volúmenes de poro, cuando la solución se vuelve sobresaturada en calcita, la concentración de bicarbonatos también comienza a disminuir hasta aproximarse al equilibrio.

El gran aumento inicial ocurrido para los HCO₃ ^- puede vincularse a los procesos de disolución de calcita, indicado por las condiciones de subsaturación de calcita que ocurren en los primeros pasos y en respuesta a la alta concentración de bicarbonatos aportados por la laguna de efluentes y el basural a partir de la oxidación de la materia orgánica como ha podido comprobarse en la modelación con NETPATH (^{Gomez et al., 2009b}).

Los procesos antes mencionados responden a los cambios no naturales que ocurren en el acuífero, debidos a la actividad antrópica como es la laguna de efluentes y el propio basural. De los modelos geoquímicos inversos realizados con NETPATH se sabe que ocurre un aumento de la fuerza iónica de las soluciones correspondientes a P2 y P3, y cambios en los índices de saturación de calcita, apareciendo con IS más elevados en P2.

Según esta primera aproximación los modelos indican que un 95 % de la contaminación por Cl^- y SO₄ ^-2son alcanzados luego de 11/2 volúmenes de poro y el total para estos iones luego de 21/2 volúmenes de poro.

En el caso de los bicarbonatos, el 100 % se alcanza luego de 14 volúmenes de poro, mostrando la influencia que tienen sobre él el aporte y oxidación de la materia orgánica junto a los procesos de disolución mineral que conforman el acuífero.

Luego de 16 volúmenes de poro se alcanza un 91 % de la concentración de Na⁺, y al igual que los demás cationes, la necesidad de un mayor tiempo para alcanzar el equilibrio, evidencian el retardo que los procesos de intercambio tienen en el avance de estos iones en el acuífero afectado por la contaminación.

4.2. Simulación con iones metálicos para una columna de 220 m

A fin de simular el flujo advectivo-dispersivo con metales, se consideró al lixiviado derivado de la laguna de efluentes, fluyendo sobre una columna con las mismas características a la considerada anteriormente. Se utilizaron los valores de dispersividad asumidos, y la composición del efluente que compone la laguna analizando los mismos iones mayoritarios e incorporando algunos iones metálicos (Figura 2). Para ello se utilizaron los análisis efectuados en febrero, que incluyen los resultados de metales pesados. Cabe señalar aquí que, si bien para la modelación con iones mayoritarios se utilizaron las muestras recogidas en el muestreo de diciembre, esto se hizo con el fin de utilizar más resultados y observar en este punto si existe alguna variación en el comportamiento de los mayoritarios.

Se consideró una columna de 22 celdas y los datos de entrada pueden verse en la Tabla 1. De este modo, el tubo es llenado con la solución correspondiente al agua no contaminada y llenando luego con la solución correspondiente a la laguna de efluentes. El análisis se efectuó considerando lo que ocurre en la última celda, luego de 220 m, en P3. La modelación con los iones mayoritarios resultantes pueden verse en la Figura 6.

Figura 6 Concentración en mol/L de Cl^-, SO₄ ^-2, Ca⁺², Na⁺, Mg⁺², K⁺, HCO₃ ^-, Pb y Zn a medida que P2 fluye en cada paso. En 22 pasos, se desplaza una vez el volumen de poro de la columna considerada. Puede observarse que tanto para el ion Cl^-, conservativo, y para el que exhibe el mismo comportamiento conservativo, los SO₄ ^-2. Los cationes evidencian el retardo de estos iones en el frente contaminante. Los bicarbonatos alcanzan el equilibrio luego de 14 volúmenes de poro, mostrando la influencia del proceso de disolución de calcita. Para la concentración de Pb y Zn en la celda 22 se desplaza una vez el volumen de poro de la columna considerada (ver en el texto el análisis detallado del comportamiento de ambos metales). 

Para los iones mayoritarios se observa el mismo comportamiento que en la simulación anterior. Puede observarse que tanto para el ion Cl^- conservativo, y para el que exhibe el mismo comportamiento conservativo, los SO₄ ^-2, ambos alcanzan el 96 % de la concentración del efluente luego de 11/2 volúmenes de poro y el 100 % luego de 21/2 volúmenes de poro (Figura 6), indicando la alta disipación del fluido contaminante una vez que se infiltró en el acuífero. Los cationes y las concentraciones sobre las superficies de intercambio muestran el mismo comportamiento que la modelación anterior, evidenciando el retardo de estos iones en el frente contaminante. Del mismo modo los bicarbonatos alcanzan el equilibrio luego de 14 volúmenes de poro, mostrando la influencia que tiene el proceso de disolución de calcita tienen sobre éstos.

Respecto al análisis de los metales puede verse que, especialmente con el Fe (Figura 7) el avance de éste en el frente contaminante está fuertemente controlado por el intercambio sobre las superficies que pueden ser tanto de compuestos orgánicos como inorgánicos. La concentración de Fe luego de 16 veces el volumen de poro alcanza apenas el 10 % de la concentración del efluente, resultando un elemento fuertemente controlado por los procesos de intercambio y adsorción (Figura 8a). El índice de saturación (IS ) de los hidróxidos de Fe (Figura 8b) muestra condiciones próximas al equilibrio y condiciones de sobresaturación respecto a esta fase luego de 13 volúmenes de poro. Al comienzo de la modelación, la solución, que representa el acuífero inmediatamente debajo del basural, se encuentra sobresaturada en respuesta al gran aporte del efluente, los cuales comienzan a ser rápidamente adsorbidos a las superficies de arcilla, materia orgánica y otros hidróxidos. El modelo conceptual indica que ocurren condiciones reductoras en la laguna de efluentes e inmediatamente debajo del basural, condiciones a partir de las cuales el Fe es movilizado desde los hidróxidos aumentando su concentración en la solución. Bajo las condiciones oxidantes del acuífero los rangos de disolución disminuyen y los procesos de precipitación y adsorción de Fe son dominantes.

Figura 7 Concentración en mol/L de Fe y Cu en la celda 22, a medida que el efluente de la laguna fluye en cada paso. En 22 pasos, se desplaza una vez el volumen de poro de la columna considerada. Co /Cf = 1 indica que el fluido original ha sido totalmente desplazado. 

Figura 8 a) Concentración de los iones metálicos adsorbidos, en mol/L sobre las superficies de intercambio (arcillas, materia orgánica e hidróxidos) en la celda 22, a medida que el efluente de la laguna fluye en cada paso. En 22 pasos, se desplaza una vez el volumen de poro de la columna considerada. b) Variación del índice de saturación (IS ) de algunos hidróxidos en la celda 22, a medida que el efluente de la laguna fluye en cada paso. 

Al inicio de la simulación el Cu (Figura 7) alcanza el 27 % de la concentración y luego de 1/2 volumen de poro ésta baja al 1 % de la concentración presente en el efluente, sin observarse cambios incluso luego de 16 volúmenes de poro. Esto muestra que el Cu en las condiciones químicas existentes en el acuífero, es un elemento que está fuertemente controlado por los procesos de adsorción y por lo tanto retrasando su avance en el frente de contaminación.

Para el Pb (Figura 8a) se alcanza el equilibrio luego de 13 volúmenes de poro. Los procesos de adsorción y de intercambio señalan que hasta los 13 volúmenes de poro el Pb es incluido en estos sitios y controlando hasta entonces la concentración de Pb en el agua del acuífero. Luego de los 13 volúmenes de poro comienza el aumento de Pb en la solución, teniendo en cuenta que el aumento que no supera los 0.0003 mg/L. Este pequeño aumento del Pb en la solución va acompañado de una disminución sobre las superficies de intercambio, aunque involucrando muy pocos moles (6×10^-9 a 2×10^-9 moles/L). Con esto, sólo los procesos de desorción no alcanzan para explicar el pequeño aumento de Pb observado en este modelo. En la Figura 8b puede observarse que los IS de los hidróxidos se encuentran subsaturados, indicando que parte de la concentración de Pb en el agua subterránea puede derivar de los procesos de disolución de los minerales fuente o del aporte mismo de los lixiviados.

En los primeros pasos del efluente sobre la columna con agua no contaminada del acuífero resulta en una concentración elevada de Zn, vinculado a que la concentración de este metal en el agua subterránea no contaminada resultó superior (Figura 6). Luego de 14 volúmenes de poro la concentración se equilibra con la concentración presente en el efluente. Luego de este tiempo la concentración de Zn adsorbido muestra muy poca variación suponiendo que son los procesos de disolución de los hidróxidos los que pueden estar controlando mayormente a este elemento en la solución, más que los procesos de intercambio y/o adsorción (Figura 8a).

Al analizar la concentración en P3 luego de 25 años (Paso, el tiempo de progreso de la reacción, 91, correspondiente a 4 veces el volumen de poro), que es el tiempo de funcionamiento del basural y la laguna, se observa que para el caso del Fe y el Zn la concentración predicha por los modelos (Fe = 0.03 mg/L y Zn= 0.02 mg/L) resultó algo mayor a la medida en P3 (Fe = 0.015 y Zn = 0). Mientras que en el caso del Cu la predicción del modelo resultó en una concentración menor (Cu = 0.0003) a la realmente medida en P3 (Cu = 0.0055), resultando una concentración similar en el caso del Pb (0.000005 mg/L y la medida en P3 de 0.000006 mg/L). Para el caso de los iones mayoritarios los modelos indican concentraciones similares para los cationes Na⁺ y K⁺ para 5 volúmenes de poro, es decir, sería necesario, según el modelo, un poco más de tiempo para alcanzar la concentración medida encontrada en P3, mientras que para Ca⁺², HCO₃ ^- y Mg⁺² las concentraciones resultaron notoriamente más altas a las medidas efectivamente en P3.

4.3. Simulación con iones mayoritarios y metálicos para una columna de 30 m

A fin de observar el efecto de la dispersividad a una distancia más corta se realizó una nueva modelación considerando una escala menor y con un nuevo valor de dispersividad.

Se simuló el transporte unidimensional sobre una columna de 30 m de largo, que es la distancia que separa al basural de P3 (Figura 9). La columna se dividió en 30 celdas de 1 m cada una, con una dispersividad de 1 m y considerando los mismo valores de velocidad del flujo, CIC y porosidad que la modelación anterior. Del mismo modo, el valor de dispersividad según la nueva escala utilizada fue tomado del diagrama de ^{Gelhar et al. (1985)}. Se consideró la composición de P2 como la solución localizada al borde del basural y por ende localizada a 30 m de P3.

Figura 9: Simulación hidrogeoquímica con iones mayoritarios para una columna de 30 m. v = velocidad real del flujo, CIC = Capacidad de intercambio catiónico, Pe = Número de Péclet y Co = Coeficiente de Courant. 

De este modo el agua subterránea contaminada atraviesa una columna de sedimentos de 30 m, saturada en agua del acuífero no contaminada de P1. Se observó el comportamiento de los iones mayoritarios en la última celda (celda Nº 30). Los datos de entrada y las figuras resultantes de la modelación para los iones mayoritarios pueden verse en la Tabla 1 y Figura 9.

Como se observa en la Figura 10, para los iones con comportamiento conservativo como el Cl^- y SO₄ ^-2, se compromete el mismo volumen de poro (21/2) hasta alcanzar el equilibrio, involucrando apenas 2 años, tiempo que está en relación a la distancia simulada. El mismo comportamiento se observa para los HCO₃ ^-y los cationes y sus concentraciones sobre las superficies de intercambio (Figura 10).

Figura 10 Concentración en mol/L de los iones mayoritarios, Cl^-, SO₄ ^-2, Ca⁺², Na⁺, Mg⁺², K⁺, en la celda 30, a medida que el agua de P2 fluye en cada paso. En 30 pasos, se desplaza una vez el volumen de poro de la columna considerada. Para los iones con comportamiento conservativo como el Cl^- y SO₄ ^-2, se involucra el mismo volumen de poro (21/2) hasta alcanzar el equilibrio, involucrando apenas 2 años, tiempo que está en relación a la distancia simulada. El mismo comportamiento se observa los cationes Ca⁺², Na⁺, Mg⁺², K⁺. 

Para la modelación del flujo advectivo-dispersivo con metales, se consideró al lixiviado derivado de la laguna de efluentes fluyendo sobre una columna de 30 m y dispersividad de 1 m. Los datos de entrada y las figuras resultantes de la modelación pueden verse en la Tabla 1 y Figura 11.

Figura 11 Concentración en mol/L de Fe, Cu, Pb y Zn, en la celda 30, a medida que el efluente de la laguna fluye en cada paso. En 30 pasos, se desplaza una vez el volumen de poro de la columna considerada. El análisis de los metales resulta igual que aquellos modelados en una columna de 220 m y dispersividad 10 m. Puede verse que la concentración de Fe luego de 16 veces el volumen de poro alcanza apenas del 10 % de la concentración del efluente, resultando un elemento fuertemente controlado por los procesos de intercambio y adsorción. Al inicio de la simulación el Cu alcanza el 27 % de la concentración y luego de medio volumen de poro ésta baja al 1 % de la concentración presente en el efluente, sin observarse cambios incluso luego de 16 volúmenes de poro. Para el Pb se alcanza el equilibrio luego de 13 volúmenes de poro señalan que el Pb es incluido en estos sitios y controlando hasta entonces su concentración en el agua del acuífero. Luego de los 13 volúmenes de poro comienza el aumento de Pb acompañado de una disminución sobre las superficies de intercambio. 

El análisis de los metales resulta igual que aquellos modelados en una columna de 220 m y dispersividad 10 m. Los procesos de adsorción (Figura 12a) controlan la presencia de los metales en el acuífero y su avance en el frente de contaminación. A partir de este análisis, además se pudo corroborar que la dispersividad resultó en función de la escala de análisis, como ya muchos autores han indicado. En la Figura 12b puede observarse que los IS de los hidróxidos se encuentran subsaturados, indicando que parte de la concentración los metales en el agua subterránea puede derivar de los procesos de disolución de los minerales fuente o del aporte mismo de los lixiviados.

Figura 12 a) Concentración de los iones metálicos adsorbidos, en mol/L, sobre las superficies de intercambio (arcillas, materia orgánica e hidróxidos), en la celda 30, a medida que el efluente de la laguna fluye en cada paso. En 30 pasos, se desplaza una vez el volumen de poro de la columna considerada. La concentración de Fe luego de 16 veces el volumen de poro alcanza apenas del 10 % de la concentración del efluente, resultando un elemento fuertemente controlado por los procesos de intercambio y adsorción. Luego de 14 volúmenes de poro, la concentración de Zn adsorbido muestra muy poca variación suponiendo que son los procesos de disolución de los hidróxidos los que pueden estar controlando mayormente a este elemento en la solución, más que los procesos de intercambio y/o adsorción (para los demás iones metálicos ver análisis detallado en el texto). b) Variación del índice de saturación (IS ) de algunos hidróxidos en la celda 30, a medida que el efluente de la laguna fluye en cada paso. Los IS de los hidróxidos de Fe muestran condiciones próximas al equilibrio y condiciones de sobresaturación respecto a esta fase luego de 13 volúmenes de poro (para los demás IS ver análisis detallado en el texto). 

Las concentraciones de los metales luego de 25 años para 900 pasos, resultan notoriamente más altas si se las compara con aquellas encontradas para el mismo tiempo en la simulación de una columna de 220 m. Del mismo modo la concentración de los iones mayoritarios resultó superior luego de simular 25 años al modelar una distancia de 30 m. Esto puede estar también relacionado con una carga contaminante más elevada respecto a la que realmente es aportada a una distancia 30 m aguas arriba de P3. Esto es, la laguna de efluentes se localiza en realidad a más de 200 m de P3.