###### Abstract  
This study collected acceleration records and station reports from 379 strong-motion stations triggered by the M\({}_{\rm w}\)7.7 earthquake in Pazarcik, Turkiye, on 6 February 2023. A total of 451 horizontal and 194 vertical strong-motion acceleration records with complete waveforms and reasonable data were selected to establish appropriate equations for predicting seismic parameters. Two iterative algorithms, specifically the optimization algorithm and orthogonal distance regression, were employed to formulate prediction equations for PGA and EPA in the NS (north-south), EW (east-west), and vertical directions. Furthermore, a comparative analysis of the attenuation curves for PGA and EPA was conducted within the vertical and horizontal directions. The results indicated that the observed and predicted values of vertical PGA and EPA at the same fault distance were approximately 0.577-0.673 times those of horizontal PGA and EPA. The ratio range calculated in our study aligns closely with previous research results, slightly more than 2/3 only in the vicinity of fault zones. The analysis demonstrated that the horizontal PGA value of the soil site at the same fault distance was approximately 1.46 times that of the bedrock site, while the vertical PGA value was only 1.12 times. The comparison of attenuation relationships revealed that within the fault distance range of 20 km to 100 km, the horizontal PGA of the Wenchuan earthquake was approximately 1.26-2.16 times that of the Pazarcik earthquake, and the ratio increased with an increasing fault distance. Pazarcik earthquake in Turkiye; peak ground motion; prediction equations; attenuation relationship; site classification 2023 13 11436 stations and the improvement of strong-motion networks, Turkiye has gradually gathered a sufficient number of strong-motion observation records. Some researchers developed regional ground motion attenuation models for Turkiye and fitted ground motion prediction equations using various parameters [5; 6; 7; 8]. To establish prediction equations for the horizontal peak ground acceleration (PGA) and acceleration response spectrum (Sa) for different site categories, [PERSON] and Kalkan utilized strong-motion records of earthquakes with magnitudes greater than 5 that occurred in Turkiye between 1976 and 1999 [5]. They found that the PGA and Sa values were overestimated for distances of approximately 15-20 km compared to the research results of [PERSON] et al. [4]. [PERSON] et al. employed strong-motion records from 22 stations triggered by the Izmit earthquake to fit the attenuation relationships of PGA and Sa values at periods of 0.3, 1.0, and 2.0 s and established ground motion prediction equations for three different site categories: rock, stiff soil, and soft soil [9]. They analyzed the influence of site classification on ground motion attenuation relationships. [PERSON] and [PERSON] collected 100 vertical strong-motion records from earthquakes within Turkiye between 1976 and 2002 and established a set of vertical ground motion prediction equations to predict vertical PGA and vertical Sa values [10]. Through strong-motion data, [PERSON] and [PERSON] developed a new set of regional ground motion prediction models for Turkiye [8]. They fitted attenuation relationships for horizontal PGA, peak ground velocity (PGV), and Sa values and compared them with the prediction equations of European and global ground motion models, which led to the identification of the differences between the Turkish models and global prediction models [11]. [PERSON] and [PERSON] optimized the NGA-W1 model using strong-motion network data in Turkiye and calculated spectral ratios of vertical-to-horizontal acceleration attenuation curves to better estimate vertical response spectrum values [12]. [PERSON] et al.