We use the linear programming algorithm introduced by Akulin et al. [V. M. Akulin, G. A. Kabatiansky, and A. Mandilara, Phys. Rev. A 92, 042322 (2015)] to perform best separable approximation on two-qutrit random density matrices. We combine the numerical results with theoretical methods in order to generate random representative families of positive partial transposed bound entangled (BE) states and analyze their properties. Our results are disclosing that for the two-qutrit system the BE states have negligible volume and that these form tiny `islands' sporadically distributed over the surface of the polytope of separable states. %We devise a method for estimating numerically the average thickness of these formations and their frequency of occurrence. The detected families of BE states are found to be located under a layer of pseudo one-copy undistillable negative partial transposed states with the latter covering the vast majority of the surface of the separable polytope.