ABSTRACT: A fully coupled model is used to simulate the results of large-scale laboratory tests that investigated the non-isothermal behaviour of geosynthetic clay liners in landfill basal liner applications. Results of numerical simulations of the laboratory tests in terms of water content, capillary pressure, temperature and stress distributions are presented, and encouraging agreement between the numerical and experimental results is achieved. Under the conditions examined, a 25°C/m temperature gradient led to the development of tensile stresses within the GCL, increasing the resultant likelihood of cracking. A Poisson's ratio of 0.25 resulted in predicted horizontal tensile stresses supporting qualitative observation of desiccation cracking in the GCLs. One important implication of the work reported herein is that elevated temperatures need not occur for extended periods to create the risk of desiccation.