A bivariate generating function for zeta values and related supercongruences

By using the Wilf-Zeilberger method, we prove a novel finite combinatorial identity related to a bivariate generating function for $\zeta (2+r+2s)$ (an extension of a Bailey-Borwein-Bradley Apery-like formula for even zeta values). Such identity is then applied to show several supercongruences.