este trabajo se centró en el estudio del comportamiento del sistema de acuífero de manglaralto a través del desarrollo de un modelo 3d con densidad variable en régimen estacionario, utilizando como software feflow v. 5.3, diersch (2006), el mismo que se basa en el método de elementos finitos para resolver problemas de flujo subterráneo y transporte con densidad variable en 2d y 3d, considerando la advección y dispersión hidrodinámica dentro del proceso de transporte de solutos. previo a la modelización se desarrolló un modelo conceptual, en donde se utilizó herramientas de hidrodinámica, hidroquímica e isotopía, con el fin de definir el modelo matemático, conocer el tipo de facie, el comportamiento de las aguas en las épocas favorables y desfavorables, el tipo de reacción agua-roca, los procesos de salinización, la fuente y zona de recarga.