The neutron scattering properties of water ice are of interest to the nuclear criticality safety community for the transport and storage of nuclear materials in cold environments. The common hexagonal phase ice Ih has locally ordered, but globally disordered, H2O molecular orientations. A 96-molecule supercell is modeled using the VASP ab initio density functional theory code and PHONON lattice dynamics code to calculate the phonon vibrational spectra of H and O in ice Ih. These spectra are supplied to the LEAPR module of the NJOY2012 nuclear data processing code to generate thermal neutron scattering laws for H and O in ice Ih in the incoherent approximation. The predicted vibrational spectra are optimized to be representative of the globally averaged ice Ih structure by comparing theoretically calculated and experimentally measured total cross sections and inelastic neutron scattering spectra. The resulting scattering kernel is then supplied to the MC21 Monte Carlo transport code to calculate time eigenvalues for the fundamental mode decay in ice cylinders at various temperatures. Results are compared to experimental flux decay measurements for a pulsed-neutron die-away diffusion benchmark.