In this paper, a bivariate power Lomax distribution based on Farlie-Gumbel-Morgenstern (FGM) copulas and univariate power Lomax distribution is proposed, which is referred to as BFGMPLx. It is a significant lifetime distribution for modeling bivariate lifetime data. The statistical properties of the proposed distribution, such as conditional distributions, conditional expectations, marginal distributions, moment-generating functions, product moments, positive quadrant dependence property, and Pearson's correlation, have been studied. The reliability measures, such as the survival function, hazard rate function, mean residual life function, and vitality function, have also been discussed. The parameters of the model can be estimated through maximum likelihood and Bayesian estimation. Additionally, asymptotic confidence intervals and credible intervals of Bayesian's highest posterior density are computed for the parameter model. Monte Carlo simulation analysis is used to estimate both the maximum likelihood and Bayesian estimators.