Abstract: Seismic return period is one of the important concepts in the field of earthquake engineering, which has been widely used in the seismic design and assessment of structures, but there is a lack of research on the generalized conditional intensity measures and vector generalized conditional intensity measures based on multivariate return period. In this paper, three multivariate seismic return period calculation methods were summarized, the theories of generalized conditional intensity measures and vector generalized conditional intensity measures were summarized. For the example site, scalar probabilistic seismic hazard analysis and disaggregation were carried out to construct the conditional distribution of the generalized conditional intensity measures, and the results of vector disaggregation based on Kendall function were obtained based on scalar hazard analysis and disaggregation. Additionally, the conditional distribution of the vector generalized conditional intensity measures were constructed, and the joint multivariate seismic return periods of multiple intensity measures were calculated. It is found that the multivariate seismic return periods based on simultaneous exceedance and arbitrary exceedance are either too strict or too loose for GCIM and VGCIM, while the multivariate seismic return period based on the Kendall function is between the above extremes and more stable than the other two definitions; when the conditional IM disaggregation results are similar to the disaggregation of the other intensity measures calculated by GCIM and VGCIM, the multivariate seismic return period based on the Kendall function is close to the set multivariate seismic return period; the MRP of the ground motion intensity measures in the VGCIM is closer to the set multivariate seismic return period than the MRP of the ground motion intensity measures in the GCIM.