AbstractIn this manuscript we reply to the issues raised by G. Bagheri and C. Bonadonna in their comment 2019JB017697. We reiterate our definition of particle Reynolds number, which is appropriate for our data set of experimental measurements, and we show that the poor performance of Bagheri and Bonadonna (2016, https://doi.org/10.1016/j.powtec.2016.06.015) model discussed in Dioguardi et al. (2018, https://doi.org/10.1002/2017JB014926) was mainly due to the typos in the equations presented in their original manuscript. We believe that the iterative methodology for calculating the drag coefficient is not strictly necessary for our data set. In this reply, however, we also include results of the intercomparison study among different drag models using the iterative methodology. Results show that the performance of the different drag models considered in the intercomparison study is not significantly affected by the employed methodology (direct vs. iterative) and that, regardless the employed methodology and unlike what has been stated in the comment 2019JB017697, our model has the best performance in reproducing the experimentally measured terminal velocities.