Underground formation characterization and uncertainty quantification using transitional Markov chain Monte Carlo sampling with an efficient parallel implementation