The purpose of this paper is to present a method for the ultrasonic characterization of air-saturated porous media, by solving the inverse problem using only the reflected waves from the first interface to infer the porosity, the tortuosity, and the viscous and thermal characteristic lengths. The solution of the inverse problem relies on the use of different reflected pressure signals obtained under multiple obliquely incident waves, in the time domain. In this paper, the authors propose to solve the inverse problem numerically with a first level Bayesian inference method, summarizing the authors' knowledge on the inferred parameters in the form of posterior probability densities, exploring these densities using a Markov-Chain Monte-Carlo approach. Despite their low sensitivity to the reflection coefficient, it is still possible to extract the knowledge of the viscous and thermal characteristic lengths, allowing the simultaneous determination of all the physical parameters involved in the expression of the reflection operator. To further constrain the problem and guide the inference, the knowledge of a particular incident angle is used at one's advantage in order to more precisely define the thermal length, by effectively yielding a statistical relationship between tortuosity and characteristic length ratio.