I think there is something odd here. For example the following code gives the same values for residuals and innovations:
fit <- us_change %>% model(ARIMA(