We prove full Szegő-type large-box trace asymptotics for selfadjoint Z d -ergodic operators Ω ∋ ω → Hω acting on L 2 (R d ). More precisely, let g be a bounded, compactly supported and real-valued function such that the (averaged) operator kernel of g(Hω) decays sufficiently fast, and let h be a sufficiently smooth compactly supported function. We then prove a full asymptotic expansion of the averaged trace Tr h(g(Hω) [−L,L] d ) in terms of the length-scale L.