The characterization of the covariance function of the solution process to a stochastic partial differential equation is considered in the parabolic case with multiplicative Lévy noise of affine type. For the second moment of the mild solution, a well-posed deterministic space-time variational problem posed on projective and injective tensor product spaces is derived, which subsequently leads to a deterministic equation for the covariance function.