The importance of a good model for the γ-ray energy spectrum from the radiative thermal neutron capture on Gadolinium (Gd) is specially increased in the present era of Gd-enhancedν e -search detectors. Its an essential prerequisite for MC studies to evaluate the neutron tagging efficiency, in order to enhance signal sensitivity in the Gd-loadedν e -search detectors. The γ-ray spectra produced from the thermal neutron capture on enriched gadolinium targets ( 155 Gd, 157 Gd and Natural Gd) in the energy range 0.11 MeV to 8.0 MeV, were measured using the ANNRI Germanium Spectrometer at MLF, J-PARC. Based on the data acquired and a GEANT4 simulation of the ANNRI detector, we reported the energy spectrum of 157 Gd(n, γ) and developed a γ-ray emission model of 157 Gd(n, γ) in our previous publication. We now present the analysed data of 155 Gd(n, γ) and nat Gd(n, γ) reactions, the energy spectra of γ-rays and an improved model for 155 Gd(n, γ), 157 Gd(n, γ) and nat Gd(n, γ) reactions. The consistency of the results from the devised model is checked among all the 14 germanium crystals, at the level of 15% spectral shape deviation at 0.2 MeV binning.