This work is concerned with zeta functions of two-dimensional shifts of finite type. A two-dimensional zeta function ζ0(s) which generalizes the Artin-Mazur zeta function was given by Lind for Z2-action φ. The n-th order zeta function ζn of φ on Zn×∞, n ≥ 1, is studied first. The trace operator Tn which is the transition matrix for x-periodic patterns of period n with height 2 is rotationally symmetric. The rotational symmetry of Tn induces the reduced trace operator τn and ζn = (det (I − snτn))−1. The zeta function
ζ =∞ Qn=1 (det (I − snτn))−1 in the x-direction is now a reciprocal of an infinite product of polynomials. The zeta function can be presented in the y-direction and in the coordinates of any unimodular transformation in GL2(Z). Therefore, there exists a family of zeta functions that are meromorphic extensions of the same analytic function ζ0(s). The Taylor series at the origin for these
zeta functions are equal with integer coefficients, yielding a family of identities which are of interest in number theory. The method applies to thermodynamic zeta functions for the Ising model with finite range interactions.