The primary objective of this research was to analyze the spatial pattern and interactions of dendritic rill channels and surface/subsurface channels (pipes)/collapsed pipes. To achieve this, photogrammetry drones were employed on over 70 ha in the protected region of Takhtsoltan in Razavi Khorasan Province. This approach pinpointed areas impacted by erosions from subsurface channels (pipes)/collapsed pipes and dendritic rill channels and surfaces. Subsequent field visits documented 76 instances of subsurface channels (pipes)/collapsed pipes and 58 dendritic rill channels and surfaces, using orthophotos derived from the aerial drone images. In the analytical phase, topographical, hydrological, soil, and biological factors were examined as independent variables, whereas the types of erosion from subsurface channels (pipes)/collapsed pipes and dendritic rill channels were viewed as dependent variables. Spatial patterns of the subsurface channels (pipes)/collapsed pipes and dendritic rill channels were discerned using univariate functions. Further, bivariate functions were deployed to probe the interrelations, revealing that the distribution pattern of the study area's subsurface channels (pipes)/collapsed pipes is predominantly clustered, in contrast to the dispersed spatial pattern of the dendritic rill channels. This analysis also confirmed significant positive correlations between dendritic rill channel erosions and subsurface channels (pipes)/collapsed pipes. Ultimately, by identifying the erosional landforms of subsurface channels (pipes)/collapsed pipes and dendritic rill channels and understanding the spatial and processual relationships between them, a deeper insight into the natural processes inherent in their spatial structure was achieved, paving the way for devising suitable strategies for their timely management.