Abstract: The creep characteristics of frozen soil are the determinant factors for judging the stability of frozen soil structures. In this study, the clay in a mining area of Shandong Province was taken as the research object. Uniaxial compression tests and creep tests were conducted on the frozen clay, through the tests we obtained the creep law of frozen clay at different freezing temperatures and stress load levels. The Singh-Mitchell (S-M) model was modified with the introduction of a hyperbolic function and non-integer order calculus, creating a non-integer order modified S-M creep model that can reflect the influence of freezing temperature and stress load levels. Further a damage variable was introduced, we established a formula reflecting the relation between damage variable and the temperature, so we proposed a non-integer order modified S-M creep-damage model that can describe the non-stable creep stage of frozen clay more accurately. Research results showed that as the freezing temperature fell, the strength of the frozen clay significantly increased, and the creep deformation reduced. The change of stress load levels also dramatically affected the creep deformation of the frozen clay. Compared with the calculated values of the classical S-M model and the experimental values, the established model could accurately reflect the creep law of frozen clay under different stress load levels in the context of temperature effects, and demonstrated a higher degree of fit. The parameters of this model have both physical and mathematical meanings, and the relative less number of parameters makes it more suitable for practical engineering applications.