I would use Statsmodels SARIMAX with data having seasonality greater than one year. In my case, I have seasonality equal 4 years, therefore I have tried:
smt.SARIMAX(endog, order=(p, d, q), seasonal_order=(P,D,Q,4))
Is it correct set
seasonality = 4 when the seasonality is equal 4 years?