Modelos de flujo y transporte de solutos en las gravas del aluvial del río Gállego en el entorno del vertedero de Sardas (Sabiñánigo, Huesca)

Abstract

[Resumen:] Se ha elaborado un modelo 2D en planta de flujo subterráneo en el aluvial del río Gállego desde la presa de Sabiñánigo hasta la confluencia con el río Aurín 2 km aguas arriba. El dominio se ha discretizado con una malla 2D de elementos finitos triangulares. Se han calibrado las entradas de agua procedentes de la terraza fluvioglaciar en la margen derecha y los coeficientes de goteo de las descargas al embalse de Sabiñánigo y por debajo de la presa. Los resultados del modelo confirman que el caudal de agua procedente del vertedero fluye por las gravas del aluvial y descarga al embalse a través de los sedimentos de colmatación depositados en el fondo del embalse, especialmente en el antiguo cauce del río Gállego. Por otro lado, se ha realizado un modelo de flujo y transporte del hexaclorociclohexano (HCH) disuelto total. El modelo se basa en la hipótesis de que el producto disuelto está en equilibrio con el presente en la fase densa o sorbido en las fases sólidas del acuífero. La relación de equilibrio entre las concentraciones del HCH disuelto total y el HCH en la fase sólida se ha simulado mediante un coeficiente de distribución, Kd, constante. Los resultados de la simulación indican que la conductividad hidráulica vertical y el coeficiente de distribución del HCH en los sedimentos de colmatación son los parámetros que más influyen en la descarga de este producto desde el aluvial hacia el embalse a través de los sedimentos y los limos del aluvial.

[Abstract:] We present a 2D groundwater flow model of the Gállego alluvial aquifer, which extends from the Sabiñánigo dam to the mouth of the Aurín river located 2 km upstream. The domain has been discretized with a 2D triangular finite element mesh. The water inflows from the fluvioglacial terrace on the right bank and the leakage coefficients of the discharges to the Sabiñánigo reservoir and underneath the dam were calibrated. The results of the model confirm that the water coming the Sardas landfill discharges into the reservoir through the silting sediments deposited at its bottom, especially along the bed of the Gállego river prior to the construction of the dam. The solute transport model is based on the hypothesis that total dissolved hexaclorociclohexane (HCH) is in equilibrium with the HCH in the nonaqueous phase (DNAPL) and the contaminants sorbed in the solid materials of the aquifer. The equilibrium relationship between the concentrations of dissolved and sorbed concentrations has been simulated by using a constant distribution coefficient, Kd. The simulation results indicate that the vertical hydraulic conductivity and the distribution coefficient of HCH in the silting sediments are the parameters that most influence the discharge of HCH from the alluvial aquifer towards the reservoir through the silting sediments.