Performance assessment of the Yucca Mountain unsaturated zone (UZ) as the site for an underground repository of high-level radioactive waste relies on the crucial assumption that water percolation processes in the unsaturated zone can be approximated as a steadystate condition. Justification of such an assumption is based on temporal damping effects of several geological units within the unsaturated tuff formation. In particular, the nonwelded tuff of the Painbrush Group (PTn unit) at Yucca Mountain, because of its highly porous physical properties, has been conceptualized to have a significant capacity for temporally damping transient percolation fluxes. The objective of this study is to investigate these damping effects, using a three-dimensional (3-D) mountain-scale model as well as several one-dimensional (1-D) models. The 3-D model incorporates a wide variety of the updated field data for the highly heterogeneous unsaturated formation at Yucca Mountain. The model is first run to steady state and calibrated using fieldmeasured data and then transient pulse infiltrations are applied to the model top boundary. Subsequent changes in percolation fluxes at the bottom of and within the PTn unit are examined under episodic infiltration boundary conditions. The 1-D model is used to examine the long-term response of the flow system to higher infiltration pulses, while the damping effect is also investigated through modeling tracer transport in the UZ under episodic infiltration condition. Simulation results show the existence of damping effects within the PTn unit and also indicate that the assumption of steady-state flow conditions below the PTn unit is reasonable. However, the study also finds that some fast flow paths along faults exist, causing vertical-flux quick responses at the PTn bottom to the episodic infiltration at the top boundary.