La modélisation inverse est devenue une approche fréquemment utilisée pour l'estimation des paramètres en hydrogéologie. Fondamentalement cette technique est basée sur les méthodes de contrôle optimal qui nécessitent des observations et un modèle pour le calcul des dérivées du premier ordre. Le modèle adjoint du modèle de Richards est construit pour obtenir le gradient de la fonction coût par rapport aux paramètres de contrôle. Les paramètres hydrodynamiques sont pris comme paramètres de contrôle; leurs valeurs optimales sont trouvées en minimisant la fonction coût ceci en utilisant un algorithme de minimisation de type descente quasi-Newton. Cette approche est utilisée pour l'identification des paramètres hydrodynamiques sur un modèle d'écou-lement souterrain en zone non saturée, ainsi que les études de sensibilité du modèle.