Regression analysis and variable selection to determine the key subduction-zone parameters that determine the maximum earthquake magnitude

Abstract

Large variations in the maximum earthquake magnitude (Mmax) have been observed among the world’s subduction zones. There is still no universal relationship between Mmax and a given subduction-zone parameter, such as plate age, plate dip angle, or plate velocity, which suggests that multiple parameters control Mmax . Here, we conduct exhaustive variable selections that are based on three evaluation criteria; leave-one-out cross-validation errors (LOOCVE), Akaike information criterion (AIC), and Bayesian information criterion (BIC) to determine the combination of subduction-zone parameters that best explains Mmax . Multiple linear regression analyses are applied using 18 subduction-zone parameters as potential candidates for the explanatory variables of Mmax . The minimum BIC is obtained when five variables (trench sediment thickness, existence of an accretionary prism, upper-plate crustal thickness, bending radius of the subducting oceanic plate, and trench depth) are selected as explanatory variables; each variable contributes positively to Mmax . Minimum LOOCVE and AIC values are obtained when eight variables (the five parameters for BIC, plus the along-strike plate convergence rate, age of the subducting plate, and maximum depth of the subducting plate) are selected. Our selection of the trench sediment thickness and plate bending radius contributing to Mmax is consistent with previous studies. The results show that increasing upper-plate crustal thickness results in a large Mmax . In addition to smoothing the subducting-plate interface via subducted sediments, along-dip extension of the crustal area along the convergent plate boundary would be important for generating a large earthquake. Graphic Abstract: [Figure not available: see fulltext.].