Eu iria para aditivo também. Como seu sinal aparente parece de baixa frequência, você pode ir um pouco além, pelo menos empiricamente. Você pode verificar, por exemplo, a homocedasticidade das diferenças finitas dos dados (primeira ou segunda ordem). Isso funcionaria como um filtro passa-alto muito bruto, onde você poderia esperar que o ruído fosse dominante.
Se o seu sinal for muito mais longo, mover janelas e transformadas de Fourier pode ajudar.
No entanto, quanto à previsão, você pode executar os dois modelos em paralelo e decidir qual deles aplicar com base, por exemplo, no melhor desempenho de um deles com base nas estatísticas anteriores. Este é um método heurístico que recentemente usei na previsão de resultados para a co-simulação de sistemas híbridos, onde nenhum modelo é conhecido: execute diferentes extrapolações em paralelo, muito rápido e decida. Não é muito teórico, mas funciona bem em nossos dados.
Se estiver interessado, eu poderia me desenvolver. A referência é chamada: CHOPtrey: extrapolação polinomial contextual on-line para co-simulação aprimorada de múltiplos núcleos de sistemas complexos
Como os dados são bastante curtos e não tenho certeza se temos um período sazonal completo, tentei realizar algumas análises de Fourier sobre os dados, seu gradiente e Laplaciano. A flutuação parece ser bastante periódica; portanto, no gráfico inferior, tentei projetar uma média móvel de "filtragem". O resíduo não varia muito em amplitude. Realmente não parece ser aleatório.