In this study, a Volterra system for the variations of metacentric height (GM) in waves is employed to simulate the parametric roll phenomena of ships in head sea condition. Using the present Volterra system, the transfer function of each harmonic component in the GM variation is computed for different ship models, including mathematical models and a real containership, and the results are validated through the comparison with the values obtained using the direct calculations based on a weakly nonlinear time-domain method. Then, a semi-analytic approach employing a 1-degree of freedom equation for roll motion is developed to simulate the parametric roll motions in irregular waves. In the derived approach, the nonlinear and time-varying restoring forces in the waves are approximated using the Volterra system. Through simulations of the parametric roll for different sea states, the effects of the 1st and 2nd-order harmonic components of the variations in the occurrence and amplitude of the parametric roll motions are investigated. Because of the strong nonlinearities in the phenomena, a stochastic analysis is conducted to examine the statistical properties of the roll motions in consideration of the sensitivities and uncertainties in the computations.