Laser welding processes are widely used for fabrications in many engineering applications such as aerospace and automotives. In this paper, a moving distributed heat source model based on Goldak's method [1] has been implemented into finite volume thermal simulations in order to predict temperature distributions during the welding process of a magnesium alloy and to study the effects of variations in thermal properties, absorption coefficient and gas shielding on the computed temperature distributions and weld pool dimensions. The main conclusion is the significant effects of varying the thermal conductivity and absorption coefficient of magnesium. Also, it has been seen that the shielding gas, besides its main role of protection against oxidation, has a significant effect on the width of the weld pool. Finally, the obtained results have been compared to the experimental ones and a satisfactory correlation has been observed, indicating the reliability of the model developed in this study.