Fractal geometry is used as a suitable tool for analysis of the sub-basin and the river networks geomorphology characteristics and simulation of complex natural phenomenon. The main objective of this article is the analysis and evaluation of Walnut Gulch Watershed based on multi-fractal analysis of the river network images and geomorphology characteristics. The basin is located at the south eastern Arizona State in the united states with total area of over 150 km2 that surrounds the historical western town of Tombstone (31o42' N, 110o03' W). To achieve this goal we inserted the required data in GIS template in the Arc Map software including Dem 10 topographic maps of the drainage network and sub-basins classification. Accordingly geo-morphologic characteristics including stream length gradient index (SL), river sinuosity (S), transverse topographic symmetry factor (T), drainage density (μ ), elongation ratio (Re) and form factor (Bs) are calculated. We also developed MATLAB codes for multi-fractal indices including fractal generalized dimension D (q), singularity spectrum f (α ) and the scaling exponent T(q). Besides, the fractal dimension is calculated using Box Counting method. Results depicted that the non-linear nature of the river network is mostly based on the multi-fractal characteristics rather than the fractal dimension. Results also showed that significant relations are recorded between the geomorphologic indices such as: μ , SL, S, T, BS and Re and the multi-fractal characteristics D (q), T (q) and f (α ). Results reported that sub-basins with most S, SL and μ values including sub-basins 3, 6, 9 and 10, also hold the greatest T (q) index and singularity angle, α . Meanwhile, the absolute difference between the fractal dimension and D (q), which is called Δ , was minimum and generalized dimension graph showed the most monotonic behavior. On the other hand for sub-basins 5, 7, 15, 2, and 1 for which the S, SL and μ indices were the smallest, T (q), D (q) and α also delivered their minimum values, Δ was maximum and generalized dimension diagram showed sharp variations. Results also showed that T has significant relation with the symmetry of the singularity spectrum. Sub-basins 7, 9, 5 and 11 with the smallest T have non-symmetric singularity spectrums and right asymmetric while sub-basins 3, 6, 8, 1 and 15 with the highest T values are more symmetric and sub-basin 10 with the maximum T is fully symmetric.