The island-arc crustal growth rate (IACGR) is the island-arc magma production volume per 1 km width along the arc strike within one million years, and its variations are highly related to slab dehydration and mantle wedge melting. A novel method that includes Earth density modelling, gravity forward and inversion, and arc crustal growth thickness integration is designed to estimate the IACGR. This method can not only estimate the IACGR along the entire arc length but also assess the crustal growth of both remnant and active arcs. Therefore, the estimation result has high coverage and low uncertainty. Here, the Izu-Bonin-Mariana (IBM) subduction zone is taken as a case study region. The estimated time-averaged IACGR along the IBM arc changes between 16 and 59 km³ km^-1 Myr^-1, with a mean value of 40 km³ km^-1 Myr^-1, and this result matches the findings of previous studies well. The uncertainties due to crustal thickness inversion are relatively larger than those from flow line reconstruction. The rate results of the Mariana part have lower uncertainties than those of the Izu-Bonin parts since the arc boundaries can be delineated more accurately. The IACGR of the region where a plateau approaches the trench tends to be overestimated because the collision of the plateau with the island arc thickens the island arc crust and bias the uniform pre-existing crustal thickness assumption.