Presentamos un marco computacional que permite investigar la respuesta Termo-Hidro-Mecánica de una parte representativa de una cuenca sedimentaria durante un ciclo de glaciación. Abordamos la complejidad del problema, que surge de la interacción mutua entre varios fenómenos, mediante un modelo multifísico y multiescalar en relación tanto al espacio como al tiempo. Nuestra contribución aborda tanto la generación de la malla computacional como el algoritmo para la solución numérica del problema. En particular, presentamos un enfoque multiescalar que tiene en cuenta el campo de deformación global de la litosfera acoplado con la retroalimentación Termo-Hidro-Mecánica de la carga de hielo en una parte representativa del dominio a una escala más fina. En el modelo de escala fina también incluimos la erosión posiblemente causada por el derretimiento del hielo. Esta metodología permite investigar la evolución de la cuenca sedimentaria como respuesta al ciclo de glaciación a una escala fina, teniendo también en cuenta el movimiento de gran escala espacial de la litosfera debido a la isostasia. Los experimentos numéricos se basan en el análisis de un escenario simple y muestran la aparición de efectos debido a la naturaleza multifísica del problema que apenas son capturados por enfoques más simples.