We investigate via Molecular Dynamics simulations the propagation of solitons in a twodimensional many-body system characterized by Yukawa interaction potential. The solitons are created in an equilibrated system by the application of electric field pulses. Such pulses generate pairs of solitons, which are characterized by a positive and negative density peak, respectively, and which propagate into opposite directions. At small perturbation, the features propagate with the longitudinal sound speed, form which an increasing deviation is found at higher density perturbations. An external magnetic field is found to block the propagation of the solitons, which can, however, be released upon the termination of the magnetic field and can propagate further into directions that depend on the time of trapping and the magnetic field strength.