We describe a procedure that allows us to solve the three-dimensional time-dependent Schrödinger equation for an atom interacting with a quantized one-mode electromagnetic field. Atom-field interaction is treated in an ab initio way prescribed by quantum electrodynamics. We use the procedure to calculate probability distributions of absorbed photons in the regime of tunneling ionization. We analyze evolution of the reduced photon density matrix describing the state of the field. We show that non-diagonal density matrix elements decay quickly, as a result of the decoherence process. A stochastic model, viewing ionization as a Markovian birth-death process, reproduces the main features of the calculated photon distributions.