With increasing availability and improving spatial and temporal resolution of geodetic data, more details of earthquake fault geometries and slip can be determined. Detailed estimates of fault model parameters are beneficial for better understanding of the earthquake mechanics at the different fault systems in the world. However, for the same earthquake, notably dissimilar coseismic fault-slip models have been produced by different authors depending on their diverse estimation methods and modeling assumptions, for example, regarding the fault geometry, elastic layering, smoothing parameters, etc. (Lay, 2018;Mai & Thingbaijam, 2014;Razafindrakoto et al., 2015). Bayesian inference of earthquake sources allows for constraining the posterior probability distributions of the different fault model parameters and can also take uncertain knowledge of the elastic layering, elastic parameters, etc. into account through model covariances or a priori constraints (