Generalized linear models (GLMs) were applied to assess the variables that predicted the relationships between the environmental parameters and δ13CPOM and δ15NPOM as in Kodama et al. (2021). We assumed that the error distributions of δ13CPOM and δ15NPOM were normally distributed with a linear link function in the GLMs. The full GLM models are as follows: