## 1 Introduction  
In deterministic seismic hazard analysis, the worst-case scenario or maximum credible earthquake is used to estimate the seismic ground-motion intensities [1, 2]. The current approach for obtaining deterministic ground-motion parameters is to employ only one set of modeling parameters for ground-motion simulations or to calculate average values using ground-motion prediction equations (GMPEs) [3, 4, 5]. However, because of the scarcity of near-field ground-motion data on moderate-to-large earthquakes, such GMPEs hardly consider the influence of near-fault effects [6, 7, 8, 9, 10, 11, 12]. The stochastic finite-fault method has been proven to enable reliable simulations of the near-field ground-motion parameters of large earthquakes, which can effectively synthesize Fourier amplitude spectra, response spectra, and the time history of acceleration [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. Nevertheless, obtaining accurate ground-motion parameters for simulations is difficult owing to inevitable aleatory and epistemic uncertainties. Therefore, uncertainties in the modeling parameters must be considered when calculating the peak ground acceleration (PGA) and spectral acceleration (SA) of large earthquakes at near field. In an uncertainty analysis, [PERSON] et al. [24] introduced a comprehensive evaluation system to obtain ground-motion parameters by applying the multi-scheme stochastic finite-fault simulation method (multi-SFFSM). They analyzedthe potential of near-field seismogenic tectonics at an engineering site by considering the impact of different source parameters, such as slip models, high-frequency attenuation, and quality factor, for ground-motion simulations. [PERSON] et al. [25] gave the value range of the key input parameters of stochastic finite-fault simulations for China's mainland, such as the fault size, rupture velocity, shear wave velocity, stress drop, high-frequency attenuation, and quality factor. However, the uncertainties of parameters when input into the multi-SFSFM have not been investigated. Since the multi-SFSM helps to statistically analyze the parameter uncertainty of ground-motion simulation through thousands or even more results, the multi-SFSM has recently been used in engineering practice for dams to evaluate the ground-motion intensity of the maximum credible earthquake. In addition to parameter uncertainty, the location and size of asperity in the source model is also challenging for earthquake scenarios in deterministic seismic hazard analysis. Either empirical relationships or a random slip model were usually adopted; in a random slip model, the slip distribution on the rupture plane is generated using a random seed. [PERSON] et al. [26] analyzed the impact of different source models on the simulation of large earthquakes at a near-field site. Their results showed that a random slip model was not suitable when conducting near-field ground-motion simulations because it cannot represent the maximum ground-motion scenario for some sites. Furthermore, for important engineering sites, the PGA, peak ground velocity (PGV), and SA are usually determined using the combined influence of multiple faults. [PERSON] et al. [27] found that considering earthquake nucleation on multiple receiver fault orientations can substantially change the predicted spatial patterns and total number of triggered events. Their result indicates that multiple faults may cause high-frequency or long-period variations in the SA or other ground-motion intensity measures. [PERSON] et al. [28] analyzed the seismogenic structures of the Xiluoudu dam site and determined three seismogenic structures and their maximum credible earthquake magnitude that can affect the dam site. These three faults influence the SA of the dam site differently; therefore, an envelope of the statistical values of the pseudo-spectral acceleration response spectra (PSA) is recommended.