One possibility is to obtain forecasts as the weighted average of the forecasts based on the parameter estimates of each regime. The weights are the smoothed probabilities of each regime, as obtained for example via the Kim's smoothing algorithm.
In the case of the ARMA model that you give (and given $ns=2$ states):
$$
y_{T+1} = \sum_{i=1}^{ns} Pr(S_{T+1}=j \,|\, \Psi_T) \left(\mu_s + \phi_{s,1} y_{T} + \theta_{s,1} \epsilon_T\right) \,,
$$
where the index $s$ denotes the regime, for example:
$\mu_1=0.00839$, $\mu_2=-0.05615$, $\phi_{1,1}=0.26694$, $\phi_{2,1}=0.70866$; $Pr(S_{T+1}=j \,|\,\Psi_T)$ is the smoothed probability that observation at time $T+1$ is in state $j$, given the entire set of information $\Psi_T=y_1,...,y_T$.
After obtaining each one-step-ahead forecast, the smoothed probabilities
$Pr(S_t=j \,|\, \Psi_T) = \sum_{i=1}^{ns} Pr(S_t=j, S_{t-1}=i \,|\, y_1,y_2,...,y_T)\,$ should be updated adding the forecasts to the set of information $\Psi_T$, $\Psi_{T+1}$,...
I don't know if it would be straightforward to apply Kim's algorithm in this case with an MA term; with an AR model the algorithm can be used as described in the reference paper.
As an alternative to the smoothed probabilities, Boot and Pick propose the usage of other weights (a link to the document is given below).
References:
Kim, CJ (1994). "Dynamic linear models with Markov-switching".
Journal of Econometrics, 60(1), pp. 1-22.
Boot, T. and Pick, A. (2014).
"Optimal forecasts from Markov switching models".
DNB Working Paper. No. 452.
URL:
http://www.dnb.nl/binaries/Working%20Paper%20452_tcm46-316416.pdf