We investigate adaptivity issues for the approximation of Poisson equations via radial basis
function-based partition of unity collocation. The adaptive residual subsampling approach
is performed with quasi-uniform node sequences leading to a flexible tool which however
might suffer from numerical instability due to ill-conditioning of the collocation matrices.
We thus develop a hybrid method which makes use of the so-called variably scaled kernels.
The proposed algorithm numerically ensures the convergence of the adaptive procedure
