Cortical bone is a heterogeneous, composite medium with a porosity from 5-10%. The characterization of cortical bone using ultrasonic techniques is a complicated procedure especially in numerical studies as several assumptions must be made to describe the concentration and size of pores. This study presents numerical simulations of ultrasound propagation in two-dimensional numerical models of cortical bone to investigate the effect of porosity on: a) the propagation of the first arriving signal (FAS) velocity using the axial transmission method, and b) the displacement and scattering amplitude in the backward direction. The excitation frequency 1 MHz was used and different receiving positions were examined to provide a variation profile of the examined parameters along cortical bone. Cortical porosity was simulated using ellipsoid scatterers and the concentrations of 0-10% were examined. The results indicate that the backscattering method is more appropriate for the evaluation of cortical porosity in comparison to the axial transmission method.