We applied a stochastic method for the finite-fault modeling of strong ground motions to the 2016 Meinong, Taiwan earthquake. Newly developed attenuation models in Southern Taiwan with the frequency-dependent Q = 86.4f 0.73 and the high-frequency decay factor κ 0 were used in the synthetic model. The horizontal-to-vertical spectral ratios (HVSR) were calculated from weak motions and the Meinong mainshock and used for the site amplification correction of the synthetic waveforms produced by the stochastic ground motion simulation. Simulations incorporating the attenuation models and site correction improved the prediction of the S-wave envelope, duration, and peak ground acceleration (PGA). The nonlinear site response during the Meinong mainshock was identified by the degree of nonlinear site response (DNL), which is a summation of HVSR differences between weak motions and the Meinong mainshock as recorded by the Taiwan Strong Motion Instrument Program. The DNL showed a positive correlation with ground motion intensity. The surface site conditions influenced DNL strength. The percentage of PGA reduction calculated in this study can be an indicator of the spatial distribution of the degree of nonlinear soil effects on the Meinong earthquake in the time domain. Areas that had high levels of PGA reduction overlap with areas that had high liquefaction potential. Based on the residual analysis, forward directivity was identified in a 105° range in the northwestward direction. The amplification of forward rupture directivity was three times greater than the backward rupture directivity. which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.