In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium’s operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER’s 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.