In 1927, Artin conjectured that any integer a which is not −1 or a perfect square is a primitive root for a positive density of primes p. While this conjecture still remains open, there has been a lot of progress in last six decades. In 2000, Moree and Stevenhagen proposed what is known as the two variable Artin's conjecture and proved that for any multiplicatively independent rational numbers a and b, the set{p⩽x:p prime, mamodp∈〈b〉modp} has positive density under the Generalised Riemann Hypothesis for certain Dedekind zeta functions. While the infinitude of such primes is known, the only unconditional lower bound for the size of the above set is due to Ram Murty, Séguin and Stewart who in 2019 showed that for infinitely many pairs (a,b)#{p⩽x:p prime, mamodp∈〈b〉modp}≫xlog2x. In this paper we improve this lower bound. In particular we show that given any three multiplicatively independent integers S={m1,m2,m3} such thatm1,m2,m3,−3m1m2,−3m2m3,−3m1m3,m1m2m3 are not squares, there exists a pair of elements a,b∈S such that#{p⩽x:p prime, mamodp∈〈b〉modp}≫xloglogxlog2x. Further, under the assumption of a level of distribution greater than x23 in a theorem of Bombieri, Friedlander and Iwaniec (as modified by Heath-Brown), we prove the following conditional result. Given any two multiplicatively independent integers S={m1,m2} such thatm1,m2,−3m1m2 are not squares, there exists a pair of elements a,b∈{m1,m2,−3m1m2} such that#{p⩽x:p prime, mamodp∈〈b〉modp}≫xloglogxlog2x.