Temporal oscillation of the negative hydrogen ion (H−) beam in the Radio Frequency (RF) ion source is investigated by a simple 3D3V Particle-In-Cell (PIC) model. The model solves electron, proton, and H− transport processes in the extraction region. The calculation domain is from the vicinity of the beam aperture to the ground electrode. To understand the relation between the plasma density oscillation and the extracted H− beam characteristics, the input electron and proton fluxes from the driver region are varied parametrically with the fundamental J-PARC RF frequency (2 MHz). The numerical results give an idea of the main physical processes between the oscillations of the plasma parameters and the extracted H− ion trajectories depending on the plasma meniscus shape in the different RF phases.