A seismic risk analysis of harp type cable stayed bridges is presented using the concept of damage probability matrix. The cable stayed bridge is modeled as a two dimensional system with the deck idealized as a continuous beam subjected to bending action and axial compression. The response of the bridge is obtained by the frequency domain spectral analysis. A double filtered power spectral density function with seismic intensity parameter taken as the magnitude of earthquake and a correlation function between the support excitations are considered as a seismic inputs. For a given magnitude of earthquake, the damage probability matrix is determined by defining three damage states namely, major, moderate and minor. The seismic risk index is determined by combining the dam- age probability matrix with the probability of occurrence of different magnitudes of earthquake. As an illustrative example, a three span cable stayed bridge is analyzed for an extensive parametric study. The parameters include degree of correlation, angle of incidence of earthquake, ratio of the components of ground motion and soil condition. Some of the important conclusions of the study indicate that i) Longitudinal component of ground motion has considerable effect on the probability of failure of the bridge deck; ii) Probability of failure is significantly more for the soft soil condition; iii) Fully correlated ground motion between support excitation provides less value of the probability of failure as compared to the uncorrelated ground motion; and iv) Probability of failure is not very sensitive to the variation of angle of incidence of earthquake.