The contribution of this paper is twofold. First, we exploit copula methodology, with two threshold GARCH models as marginals, to construct a bivariate copula threshold GARCH model, simultaneously capturing asymmetric nonlinear behaviour in univariate stock returns of spot and futures markets...