The Zarshuran Carlin-like gold deposit is located at the Takab Metallogenic belt in the northern part of the Sanandaj-Sirjan zone, NW Iran. The high-grade ore bodies are mainly hosted by black shale and cream to gray massive limestone along the NNE-trending extensional fault/fracture zones. The aim of this investigation was to determine and separate the gold mineralized stages based on the surface litho-geochemical Au, Hg, and As data using the Concentration-Area (C-A) fractal model and stepwise factor analysis in the Zarshuran gold deposit. Three mineralized stages were determined by the C-A fractal modeling and factor analysis, which were correlated with the mineralized stages from geological studies. The main stage of Au mineralization was higher than 1.995 ppm, which was correlated with the main sulfidation stage, whereas the As and Hg highly intense anomalies (higher than 6409 and 19 ppm, respectively) were associated with the quartz-sulfide veins and veinlets. The results obtained by the C-A fractal model and stepwise factor analysis showed that the main gold mineralized stage occurred in the southern part of the Zarshuran deposit, which was correlated with the geological particulars.