The kinetics of two-step photoionization of barium atoms through excited
P
1
1
,
3
levels was studied by numerical simulation method in the optically dense approximation. The first laser at wavelength
λ
1
=
553.5
nm (
λ
1
=
791
nm) pumps transition
S
0
1
↔
P
1
1
(
S
0
1
↔
P
1
3
). The second laser with
λ
2
=
266
nm (
λ
2
=
310
nm) ionizes atoms from the level
P
1
1
(
P
1
3
). The third probe laser with wavelength
λ
i
=
493.41
nm causes resonant fluorescence of barium ions. It is shown that the metastable
D
2
3
from the
D
1
,
2
,
3
3
triplet in barium is strongly populated through
P
1
1
,
3
levels. Therefore, the applying of the ionizing laser tuned to a continuum from this metastable
D
2
3
can yield a significant increase in the ionization efficiency. The laser-induced fluorescence method makes it possible to estimate the photoionization yield of barium ions with an error of up to 8%. The loading rate of the barium ion beam into the ion trap, obtained numerically, is consistent with the experimental data [Opt. Exp. 19(17), 16438 (2011)]. The deceleration of ultracold plasma recombination due to the heating of electrons in superelastic collisions for the characteristic times of a two-stage process is predicted.