Based on the Hartree-Fock-Bogoliubov solutions in large deformed coordinate spaces, the finite amplitude method for quasiparticle random phase approximation (FAM-QRPA) has been implemented, providing a suitable approach to probe collective excitations of weakly-bound nuclei embedded in the continuum. The monopole excitation modes in Magnesium isotopes up to the neutron drip line have been studied with the FAM-QRPA framework on both the coordinate-space and harmonic oscillator basis methods. Enhanced soft monopole strengths and collectivity as a result of weak-binding effects have been unambiguously demonstrated.