In this work, we consider the (3+1)-dimensional Burgers equation with variable coefficients, which is frequently used to define the motion of solitary waves. Abundant lump waves are constructed by taking the ansatz as a rational function. Furthermore, mixed solutions utilizing lump waves, rogue waves, and kink solitons are obtained by combining the rational function with an exponential function, resulting in fission and fusion phenomena.