In a planar smoothly bounded domain
$\Omega$, we consider the model for oncolytic virotherapy given by
$$\left\{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v) - uz, \\[1mm] v_t = - (u+w)v, \\[1mm] w_t = d_w \Delta w - w + uz, \\[1mm] z_t = d_z \Delta z - z - uz + \beta w, \end{array} \right.$$ with positive parameters
$ D_w $,
$ D_z $ and
$\beta$. It is firstly shown that whenever
$\beta \lt 1$, for any choice of
$M \gt 0$, one can find initial data such that the solution of an associated no-flux initial-boundary value problem, well known to exist globally actually for any choice of
$\beta \gt 0$, satisfies
$$u\ge M \qquad \mbox{in } \Omega\times (0,\infty).$$ If
$\beta \gt 1$, however, then for arbitrary initial data the corresponding is seen to have the property that
$$\liminf_{t\to\infty} \inf_{x\in\Omega} u(x,t)\le \frac{1}{\beta-1}.$$ This may be interpreted as indicating that
$\beta$ plays the role of a critical virus replication rate with regard to efficiency of the considered virotherapy, with corresponding threshold value given by
$\beta = 1$.