We study the NERA model that describes the dynamic evolution of illicit drug usage in a population. The model consists of nonusers (N) and three categories of drug users: the experimental (E) category, the recreational (R) category and the addict (A) category. Two epidemic threshold terms known as the reproduction numbers, R0 and µ are defined and derived. Sensitivity analysis of R0 on the parameters are performed in order to determine their relative importance to illicit drug prevalence. The local and global stability of the equilibrium states are also analysed. We also prove that a transcritical bifurcation occurs at R0 = 1. It is shown that an effective campaign of prevention can help to fight against the prevalence of illicit drug consumption. We demonstrate persistence when R0 > 1 and conditions for the extinction of drug consumption are also established. Numerical simulations are performed to verify our model. Our results show that the NERA model can assist policy makers in targeting prevention for maximum effectiveness and can be used to adopt evidence-based policies to better monitor and quantify drug use trends.