Tissue covered by thin liquid layers treated by atmospheric pressure plasmas for biomedical applications ultimately requires a reproducible protocol for human healthcare. The desired outcomes of wet tissue treatment by dielectric barrier discharges (DBDs) depend on the plasma dose which determines the integral fluence of radicals, ions, electric fields and UV/VUV photons incident onto the tissue. These fluences are controlled by power, frequency and treatment time. To first order, these parameters determine the energy deposition (J cm −2) onto the tissue. However, energy deposition may not be the only parameter that determines the fluences of reactants to the underlying tissue. In this paper, we report on a computational investigation of multipulse DBDs interacting with wet tissue. The DBDs were simulated for 100 pulses at different repetition rates and liquid thicknesses followed by 10 s or more of afterglow. Two schemes were investigated-stationary and random. In the stationary scheme, the DBD plasma streamer continues to strike at the same location on the liquid layer, whereas in the random scheme the plasma streamer strikes at random locations on the liquid layer. These differences in streamer locations strongly affect the spatial distribution of solvated species such as OH aq and H 2 O 2aq (' aq ' represents an aqueous species), which have high rates of solvation. The spatial distribution of species such as NO aq, which have low rates of solvation, are less affected by the location of the streamer due to the remediating effects of diffusion in the air. The end result is that fluences to the tissue are sensitive to the spatial location of the streamer due to the ensuing reactions in the liquid between species that have low and high rates of solvation. These reactions can be controlled not only through location of the streamer, but also by repetition rate and thickness of the liquid layer.