A benchmark for two-dimensional numerical subduction modeling using COMSOL Multiphysics (R)

Subduction has been the focal point of numerical studies for decades because it plays an important role in the Earth's mass and energy circulations and generates earthquakes and arc volcanoes which are closely related to the human lives. Among the studies on subduction, numerical modeling has been broadly applied to the quantitative studies on the subducting slab in the mantle which cannot be directly observed. In this study, we benchmark the numerical subduction modeling using a finite element package, COMSOL Multiphysics (R) and the results are consistent with the previously reported benchmark results.