Hierarchies can be modeled by a set of exponential functions, from which we can derive a set of power laws indicative of scaling. The solution to a scaling relation equation is always a power law. The scaling laws are followed by many natural and social phenomena such as cities, earthquakes, and rivers. This paper reveals the power law behaviors in systems of natural cities by reconstructing the urban hierarchy with cascade structure. Cities of the U.S.A., Britain, France, and Germany are taken as examples to perform empirical analyses. The hierarchical scaling relations can be well fitted to the data points within the scaling ranges of the number, size and area of the natural cities. The size-number and area-number scaling exponents are close to 1, and the size-area allometric scaling exponent is slightly less than 1. The results show that natural cities follow hierarchical scaling laws very well. The principle of entropy maximization of urban evolution is then employed to explain the hierarchical scaling laws, and differences entropy maximizing processes are used to interpret the scaling exponents. This study is helpful for scientists to understand the power law behavior in the development of cities and systems of cities. © 2018 by the authors.