In previous work, the first and third authors introduced staircase tableaux, which they used to give combinatorial formulas for the stationary distribution of the asymmetric simple exclusion process (ASEP) and for the moments of the Askey-Wilson weight function. The fact that the ASEP and Askey-Wilson moments are related at all is quite surprising, and is due to Uchiyama-Sasamoto-Wadati. The ASEP is a model of particles hopping on a one-dimensional lattice of N sites with open boundaries, particles can enter and exit at both left and right borders. It was introduced around 1970 and is cited as a model for both traffic flow and translation in protein synthesis. Meanwhile, the Askey-Wilson polynomials are a family of orthogonal polynomials in one variable, they sit at the top of the hierarchy of classical orthogonal polynomials. So we have the relationship ASEP -- staircase tableaux -- Askey-Wilson moments It is well-known that Askey-Wilson polynomials can be viewed as the one-variable case of the multivariate Koornwinder polynomials, also known as the Macdonald polynomials for the type BC root system. It is natural then to ask whether one can generalize the relationships among the ASEP, Askey-Wilson moments, and staircase tableaux, in such a way that Koornwinder moments replace Askey-Wilson moments. In a recent work, we demonstrated a close connection between Koornwinder moments and the two-species ASEP (a particle model involving two species of particles with different "weights"). In this article we introduce rhombic staircase tableaux, and show that we have the relationship 2-species ASEP -- rhombic staircase tableaux -- Koornwinder moments In particular, we give formulas for the steady state distribution of the two-species ASEP and for Koornwinder moments, in terms of rhombic staircase tableaux.