Abstarct: The Four Gonbad Copper Deposit is located 65 kilometers northeast of Sirjan and 110 kilometers southwest of Kerman. Geologically, it is situated in the southern part of the Central Iranian Zone and within the Urumieh-Dokhtar Magmatic Zone. The Urumieh-Dokhtar Magmatic Zone is a long strip extending northwest-southeast from Sahand to Bazman, primarily composed of Cenozoic volcanic and pyroclastic rocks. The mineralization resulted from the intrusion of granitic, granodioritic, and quartz dioritic bodies into Eocene and Oligo-Miocene volcanic-sedimentary rocks, such as basalt, andesite, tuff, sandstone, and conglomerate, which have undergone propylitic, phyllic, and argillic alterations. The mineralization occurs in vein form, with ore present as veins, veinlets, and disseminated within the host rock. The Four Gonbad deposit exhibits both hypogene and supergene mineralization, containing chalcopyrite, pyrite, galena, hematite, goethite, and limonite, accompanied by comb quartz and calcite. Modeling, reserve estimation, and final boundary determination are among the most critical stages in the design of any open-pit mine. One of the most significant challenges faced by geologists and mining engineers has been defining the boundaries of the ore deposit and estimating the quantity and quality of the reserves. This study presents the results of geostatistical modeling and reserve estimation for the Four Gonbad Copper Mine. The Leapfrog software was used to carry out these stages. Initially, using available data such as exploration drill holes, geological cross-sections, and other relevant information, a model of the deposit and the mineralization zone of the Four Gonbad Copper Mine was prepared. Subsequently, statistical and geostatistical analyses were performed on the exploration data, and the grade for the block model was estimated using the kriging method. The reserve of the deposit was then calculated. baxsed on the 3D models created in Leapfrog and combined with the grade model, a suitable grade classification in three-dimensional space was prepared, allowing for the separation of high-grade zones from waste material and effective management of the feed to the processing plant.