  This article reports global fits of short-baseline neutrino data to
oscillation models involving light sterile neutrinos. In the commonly-used 3+1
plane wave model, there is a well-known 4.9$\sigma$ tension between data sets
sensitive to appearance versus disappearance of neutrinos. We find that models
that damp the oscillation prediction for the reactor data sets, especially at
low energy, substantially improve the fits and reduce the tension. We consider
two such scenarios. The first scenario introduces the quantum mechanical
wavepacket effect that accounts for the source size in reactor experiments into
the 3+1 model. We find that inclusion of the wavepacket effect greatly improves
the overall fit compared to a 3$\nu$ model by $\Delta \chi^2/$DOF$=61.1/4$
($7.1\sigma$ improvement) with best-fit $\Delta m^2=1.4$ eV$^2$ and wavepacket
length of 67fm. The internal tension is reduced to 3.4$\sigma$. If reactor-data
only is fit, then the wavepacket preferred length is 91 fm ($>20$ fm at 99\%
CL). The second model introduces oscillations involving sterile flavor and
allows the decay of the heaviest, mostly sterile mass state, $\nu_4$. This
model introduces a damping term similar to the wavepacket effect, but across
all experiments. Compared to a three-neutrino fit, this has a $\Delta
\chi^2/$DOF$=60.6/4$ ($7\sigma$ improvement) with preferred $\Delta m^2=1.4$
eV$^2$ and decay $\Gamma = 0.35$ eV$^2$. The internal tension is reduced to
3.7$\sigma$.
  For many years, the reactor event rates have been observed to have structure
that deviates from prediction. Community discussion has focused on an excess
compared to prediction observed at 5 MeV; however, other deviations are
apparent. This structure has $L$ dependence that is well-fit by the damped
models. Before assuming this points to new physics, we urge closer examination
of systematic effects that could lead to this $L$ dependence.
