Hydrochemical tracers in the middle Rio Grande Basin,USA: 2. Calibration of a groundwater-flow model

The calibration of a groundwater model with the aid of hydrochemical data has demonstrated that low recharge rates in the Middle Rio Grande Basin may be responsible for a groundwater trough in the center of the basin and for a substantial amount of Rio Grande water in the regional flow system. Earlier models of the basin had difficulty reproducing these features without any hydrochemical data to constrain the rates and distribution of recharge. The objective of this study was to use the large quantity of available hydrochemical data to help calibrate the model parameters, including the recharge rates. The model was constructed using the US Geological Surveys software MODFLOW, MODPATH, and UCODE, and calibrated using ¹⁴C activities and the positions of certain flow zones defined by the hydrochemical data. Parameter estimation was performed using a combination of nonlinear regression techniques and a manual search for the minimum difference between field and simulated observations. The calibrated recharge values were substantially smaller than those used in previous models. Results from a 30,000-year transient simulation suggest that recharge was at a maximum about 20,000 years ago and at a minimum about 10,000 years ago.

---

Resumen La calibración de un modelo de aguas subterráneas con el apoyo de datos hidroquímicos ha demostrado que la recarga relativamente baja en la cuenca media del Río Grande es probablemente responsable de una depresión de aguas subterráneas en el centro de la cuenca y de la presencia de una cantidad considerable de agua del Río Grande en el acuífero del Grupo Santa Fe. Los modelos propuestos con anterioridad para la cuenca tenían dificultades para reproducir estas características ya que no tenían datos hidroquímicos que permitieran delimitar los ritmos y distribución de recarga. El objetivo del presente estudio consistió en utilizar una gran cantidad de datos hidroquímicos disponibles para ayudar a calibrar los parámetros del modelo, incluyendo los ritmos de recarga. El modelo se construyó utilizando los modelos MODFLOW, MODPATH, y UCODE del USGS, mientras que la calibración se realizó en base a concentraciones de ¹⁴C y a la posición de ciertas zonas definidas con los datos hidroquímicos. La estimación de parámetros se realizó en base a una combinación de técnicas de regresiones no lineares y a una búsqueda a viva fuerza del error mínimo entre los datos observados y los simulados. Los valores de recarga calibrados fueron significativamente más bajos que los estimados en los modelos anteriores. Los resultados de una simulación transitoria de 30,000 años sugieren que la recarga durante la última glacial máxima (LGM) fue diez veces el ritmo moderno, pero que la recarga que ocurrió inmediatamente después de la LGM fue más baja que el ritmo moderno.

---

Résumé Le calibrage dun modèle hydrogéologique avec laide de données hydrochimiques a démontré que la recharge relativement faible dans le Grand Bassin du Middle Rio est vraisemblablement responsable dune dépression des eaux souterraines dans le centre du bassin et de la présence dune quantité substantielle deau du Rio Grande dans laquifère du Groupe de Santa Fe. Les modèles antérieurs avaient des difficultés à reproduire ses conclusions sans laide de données hydrochimiques pour contraindre les taux et la distribution de la recharge. Lobjectif de cette étude était dutiliser une grande quantité de données hydrochimiques permettant de calibrer les paramètres du modèle, et notamment les taux de recharge. Le modèle a été construit avec les logiciels MODFLOW, MODPATH et UCODE, et calibré en utilisant les concentrations en ¹⁴C et la position de certaines zones définies par les données hydrochimiques. Lestimation de certains paramètres a été réalisée en utilisant une combinaison de techniques de régression non linéaire et une méthode de recherche exhaustive (Brute Force Search) de lerreur minimum entre les résultats des observations et les simulations. Les valeurs de la recharge calibrée sont substantiellement plus basses que celles estimées dans les modèles antérieurs. Les résultats dune simulation en régime transitoire sur 30.000 ans suggèrent que la recharge au maximum de la dernière glaciation (last glacial maximum, LGM) était 10 fois supérieure au taux actuel, mais que la recharge qui a suivit la LGM était plus bas que la recharge actuelle.