We present ALMA dust polarization and molecular line observations toward 4 clumps (I(N), I, IV, and V) in the massive star-forming region NGC 6334. In conjunction with large-scale dust polarization and molecular line data from JCMT, Planck, and NANTEN2, we make a synergistic analysis of relative orientations between magnetic fields (θ B ), column density gradients (θ NG ), local gravity (θ LG ), and velocity gradients (θ VG ) to investigate the multi-scale (from ∼30 pc to 0.003 pc) physical properties in NGC 6334. We find that the relative orientation between θ B and θ NG changes from statistically more perpendicular to parallel as column density (N H2 ) increases, which is a signature of trans-tosub-Alfvénic turbulence at complex/cloud scales as revealed by previous numerical studies. Because θ NG and θ LG are preferentially aligned within the NGC 6334 cloud, we suggest that the more parallel alignment between θ B and θ NG at higher N H2 is because the magnetic field line is dragged by gravity. At even higher N H2 , the angle between θ B and θ NG or θ LG transits back to having no preferred orientation or statistically slightly more perpendicular, suggesting that the magnetic field structure is impacted by star formation activities. A statistically more perpendicular alignment is found between θ B and θ VG throughout our studied N H2 range, which indicates a trans-to-sub-Alfvénic state at small scales as well and signifies an important role of magnetic field in the star formation process in NGC 6334. The normalised mass-to-flux ratio derived from the polarization-intensity gradient (KTH) method increases with N H2 , but the KTH method may fail at high N H2 due to the impact of star formation feedback.