We investigate how well the SPHEREx all-sky survey can constrain local primordial non-Gaussianity beyond the parameter f NL using galaxy power spectra. We forecast joint constraints on the parameters f NL, g NL and τ NL obtained assuming a simple two-field curvaton model of inflation. The parameters f NL and g NL characterise the squeezed limits of the primordial bispectrum and trispectrum respectively, and lead to a characteristic scale-dependence of the galaxy bias that increases out to arbitrarily large scales. Values of the parameter τ NL > (6/5f NL)2 cause the galaxy power spectrum to have a stochastic component which also increases out to arbitrarily large scales. Our MCMC forecasts indicate that SPHEREx can provide joint constraints on any two of the three parameters f NL, g NL and τ NL. Due to strong degeneracies among these parameters, measurements of the galaxy power spectra alone may not be sufficient to jointly constrain all three. Constraints on f NL, g NL and τ NL obtained from galaxy power spectrum observations depend on the modelling of underlying nuisance parameters. We study the robustness of our forecast constraints to modelling choices and note that even with relatively conservative modelling assumptions, SPHEREx galaxy power spectra can provide strong evidence of local non-Gaussianity, even if the particular values of f NL and g NL cannot be measured precisely.