The 22 March 2020 Zagreb earthquake of local magnitude ML = 5.5 has caused great damage in the city of Zagreb and in the local area. The most affected were old buildings (built before 1963) and cultural heritage buildings. Therefore, the thesis presents seismic site response analysis for the Zagreb ML5.5 earthquake scenario that was calculated on realistic soil profiles, which are based on data from seismic and geological microzonation of the part of the city of Zagreb area and from additional measurements with Multichannel analysis of surface waves on three soil profiles. 1D equivalent linear site response analysis was done for two cases, 1.) case in which all soil profiles are in epicenter area (1 km from epicenter) and 2.) case in which epicentral distance was taken into account (input rock motions are attenuated with distance). These calculated values were used to construct an interpolated 2D profile showing the variation of peak ground acceleration with depth, and thus the amplification of seismic motion due to spreading through different soil layers. The profile extends from the epicenter of the Zagreb ML5.5 earthquake to the accelerometer stations QUHS and QARH. Furthermore, obtained results were compared with empirical values and with results from similar research. The amplification factor is estimated to be about 2 for the scenarios of the Zagreb ML5.5 earthquake, and the estimated values of peak ground acceleration are 0.35 g in the epicentral area and 0.20 g at a distance of 12 km from the epicenter. The paper contributes to a better understanding of the consequences of the Zagreb ML5.5 earthquake, it shows the importance of local soil conditions on the damage extent and the importance of microzonation and seismic soil response analysis for building earthquake-resistant structures.