Summary: The system \[ \begin{cases} u_t = \Delta u - \chi \nabla \cdot(\frac{u}{v} \nabla v) - u v + B_1(x, t), \\ v_t = \Delta v + u v - v + B_2(x, t), \end{cases} \tag{\(\star\)} \] is considered in a disk \(\Omega \subset \mathbb{R}^2\), with a positive parameter \(\chi\) and given nonnegative and suitably regular functions \(B_1\) and \(B_2\) defined on \(\Omega \times(0, \infty)\). In the particular version obtained when \(\chi = 2\), (\(\star\)) was proposed in [M. B. Short et al., Math. Models Methods Appl. Sci. 18, 1249–1267 (2008; Zbl 1180.35530)] as a model for crime propagation in urban regions. Within a suitable generalized framework, it is shown that under mild assumptions on the parameter functions and the initial data the no-flux initial-boundary value problem for (\(\star\)) possesses at least one global solution in the case when all model ingredients are radially symmetric with respect to the center of \(\Omega\). Moreover, under an additional hypothesis on stabilization of the given external source terms in both equations, these solutions are shown to approach the solution of an elliptic boundary value problem in an appropriate sense. The analysis is based on deriving a priori estimates for a family of approximate problems, in a first step achieving some spatially global but weak initial regularity information which in a series of spatially localized arguments is thereafter successively improved. To the best of our knowledge, this is the first result on global existence of solutions to the two-dimensional version of the full original system (\(\star\)) for arbitrarily large values of \(\chi\).