Lomax distribution has been studied by many statisticians due to its important role in reliability modeling and lifetime testing. The bivariate Lomax distribution is an important lifetime distribution in survival analysis. In this article, a bivariate Lomax distribution is constructed based on Clayton copula. The joint probability density function and the joint cumulative distribution function are derived in closed forms. Some properties of this bivariate distribution are discussed. The maximum likelihood and Bayes estimators for the unknown parameters are derived. Also, the maximum likelihood and Bayesian two-sample prediction for the future observations are obtained. The performance of the proposed bivariate distribution is examined using a simulation study. Finally, one real data set under the proposed distribution is considered to illustrate its flexibility and applicability for real-life applications.