  We present lattice non-relativistic QCD calculations of bottomonium
correlation functions at temperatures $T \simeq 150-350$ MeV. The correlation
functions were computed using extended bottomonium operators, and on background
gauge-field configurations for 2+1-flavor QCD having physical kaon and
nearly-physical pion masses. We analyzed these correlation functions based on
simple theoretically-motivated parameterizations of the corresponding spectral
functions. The results of our analyses are compatible with significant
in-medium thermal broadening of the ground state S- and P-wave bottomonia.
