The hot working behavior of Mg-Gd-Y-Nb-Zr alloy was investigated using constitutive model and hot processing maps in this work. Isothermal compression tests were conducted with temperature and strain rate range of 703–773 K and 0.01–5 s~(-1), respectively. Improved Arrhenius-type equation incorporated with strain compensations was used to predict flow behavior of the alloy, and the predictability was evaluated using correlation coefficient, root mean square error and absolute relative error. Processing maps were constructed at different strains for Mg-Gd-Y-Nb-Zr alloy based on dynamic materials model.The processing maps are divided into three domains and the corresponding microstructure evolutions are referred to the forming of straight grain boundaries, twinning, dynamic recrystallization and grain growth. Instability occurred mainly at the strain rate range of 0.3s~(-1)–0.5s~(-1). The optimum processing domain is mainly at the temperature range of 703–765 K with the strain rate range of 0.01–0.1 s~(-1). The hot working behavior of Mg-Gd-Y-Nb-Zr alloy was investigated using constitutive model and hot processing maps in this work. Isothermal compression tests were conducted with temperature and strain rate range of 703-773 K and 0.01-5 s ~ (-1), respectively. Improved Arrhenius-type equation incorporated with strain compensations was used to predict flow behavior of the alloy, and the predictability was evaluated using correlation coefficient, root mean square error and absolute relative error. Processing maps were constructed at different Dyna for Mg-Gd-Y-Nb-Zr alloy based on dynamic materials model. The processing maps are divided into three domains and the corresponding microstructure evolutions are referred to the forming of straight grain boundaries, twinning, dynamic recrystallization and grain growth. Instability Occurring mainly at the strain rate range of 0.3s ~ (-1) -0.5s ~ (-1). The optimum processing domain is at the temperature range of 703-765 K with the strain ra te range of 0.01-0.1 s ~ (-1).