Abstarct: in this study the probabilistic seismic hazard analysis has been conducted in the city of Babol. To achive the objectives, while studying the seismotectonic of the region, an area extending from 150 Km radius from Babol city centre was chosen for the study. By studying seismotectonic of the region, epicenter distribution and some information such as geological information of 11 earthquake source have been modeled. Linear sources or faults were used for deterministic estimation of the maximum value of earth motion parameters. There are several relationships for estimation of maximum magnitude of every linear source. Which five of them were chosen in order to do the study. The earthquakes greater than 6 Richter were collected from through country and the relation between earthquake magnitude value and the reported faulting value were controlled with each of five relationship. Finally Nowroozi and Wells & Cooper Smith relationships were introduced as proper relationships for this zone. The average results of this two relationship were used for determination of maximum magnitude resources. In order to chose proper attenuation relationship, among all attenuation relationships, four relationship, Joyner & Boore (1981), Ambraseys & Douglas (2003), Campbell and Bozorgnia (2003) and Khademi (2002) that were more important and applicable than others were chosen. By correlating these relationships with seismogram data of Baladeh earthquake (1383), two relations of Joyner & boore (1981) and Campbell and Bozorgnia (2003) , were chosen as proper relations for this region. . Two proper slowness relationships for this zone (tow relationship, Joyner & Boore and Campell & Bozorgnia) were used for deterministic estimation of maximum horizontal acceleration. On the basis of results, the maximum values are related to blocks 1 and 2 in Khazar faults and also block 1 in North Alborz fault which contained 0.36g, 0.30g and 0.21g respectively. Maximum vertical acceleration were also calculated for faulted zone blocks which maximum value of results are belong to blocks 1 and 2 in khazar fault and also blocks 1 and 2 in north alborz fault which maximum vertical acceleration values are 0.36g, 0.29g, 0.18 and 0.18g, respectively. Logic tree approach was used in order to consider uncertainty in danger analysis by deterministic method. For doing this, 11 regional sources were used as local seismic source. According to results, maximum horizontal acceleration are belongs to sources 1, 2 and 3 with values of 0.38g, 0.33g and 0.21g, respectively. In the part probabilistic seismic hazard analysis According to (Shapira et al 2007) studies and assuming that  is constant,  value for each earthquake source was normalized. In order to chose proper attenuation relationship, among all attenuation relationships, four relationship, Joyner & Boore (1981), Ambraseys & Douglas (2003), Campbell and Bozorgnia (2003) and Khademi (2002) that were more important and applicable than others were chosen. By correlating these relationships with seismogram data of Baladeh earthquake (1383), two relations of Joyner & boore (1981) and Campbell and Bozorgnia (2003) , were chosen as proper relations for this region. At the next steep, by using SEIS RISK I software, isoacceleration map of the design surfaces, equivalent to excess probability of 0.10 and 0.20 during 50 years (equivalent to 475 and 2500 years returning period, respectively) was prepared by using 7*9 network. Maximum horizontal acceleration for returning period of 475 and 2500 years was estimated at 0.27 and 0.4, respectively. Vertical acceleration values in different return periods were calculated by Campell and Bozorgnia (2003) relationship for concerned field. Maximum vertical acceleration values, in two return periods of 475 and 2475 years are 0.17g and 0.28g, respectively. Finally, considering the importance of response spectrum of baxsement, the values of this parameter for two returning periods of 200 and 2500 years were calculated by taking into account Joyner & Boore and Campbell and Bozorgnia data. According to obtained results in this return period, maximum response spectrum can be observed in period 0.5 secend with value of 0.28g, to other periods.