At the international congress of mathematicians in 1900, Hilbert claimed that the Riemann zeta function \(\zeta (s)\) is not the solution of any algebraic ordinary differential equations on its region of analyticity. Let T be an infinite order linear differential operator introduced by Van Gorder in 2015. Recently, Prado and Klinger-Logan [9] showed that the Hurwitz zeta function \(\zeta (s,a)\) formally satisfies the following linear differential equation \(\begin{aligned} T\left[ \zeta (s,a) - \frac{1}{a^s}\right] = \frac{1}{(s-1)a^{s-1}}. \end{aligned}\) Then in [6], by defining \(T_{p}^{a}\) , a p-adic analogue of Van Gorder’s operator T, we constructed the following convergent infinite order linear differential equation satisfied by the p-adic Hurwitz-type Euler zeta function \(\zeta _{p,E}(s,a)\) \(\begin{aligned} T_{p}^{a}\left[ \zeta _{p,E}(s,a)-\langle a\rangle ^{1-s}\right] =\frac{1}{s-1}\left( \langle a-1 \rangle ^{1-s}-\langle a\rangle ^{1-s}\right) . \end{aligned}\) In this paper, we consider this problem in the positive characteristic case. That is, by introducing \(\zeta _{\infty }(s_{0},s,a,n)\) , a Hurwitz type refinement of Goss zeta function, and an infinite order linear difference operator L, we establish the following difference equation \(\begin{aligned} L\left[ \zeta _{\infty }\left( \frac{1}{T},s,a,0\right) \right] =\sum _{\gamma \in \mathbb {F}_{q}} \frac{1}{\langle a+\gamma \rangle ^{s}}. \end{aligned}\)
