Hub facility location problems involve the establishment of strategic hub facilities and allocation of demand nodes to them. Meanwhile, with the emergence of new transport and distribution networks with multi-level structures, the design of such networks has been evolved. The diversity in transportation system has added multimodality to these problems. Considering the strategic nature and long-term implications of decision-making in this field, the decisions shall be of high reliability. The present study is an attempt to consider reliability of communication paths in proportion to the volume of transaction through them while covering all of the above-mentioned issues in hub location problems at the same time. Incorporating hierarchical characteristics of the problem into the model, one can obtained a significantly enhanced model in terms of multimodality and reliability. A Lagrangian solution method was developed considering the strategic level of the problem and the importance of the solution accuracy. The model was then validated in terms of time and quality.