We study variational inequalities which are governed by a strongly monotone and Lipschitz continuous operator F over a closed and convex set S. We assume that is the nonempty solution set of a (multiple-set) split convex feasibility problem, where C and Q are both closed and convex subsets of two real Hilbert spaces and respectively, and the operator A acting between them is linear. We consider a modification of the gradient projection method the main idea of which is to replace at each step the metric projection onto S by another metric projection onto a half-space which contains S. We propose three variants of a method for constructing the above-mentioned half-spaces by employing the multiple-set and the split structure of the set S. For the split part we make use of the Landweber transform.