This article studies the non-homogeneous quadratic Bessel zeta function ζRB(s, v, a), defined as the sum of the squares of the positive zeros of the Bessel function Jv(z) plus a positive constant. In particular, explicit formulas for the main associated zeta invariants, namely, poles and residua ζRB(0, v, a) and ζRB(0, v, a), are given.