The primary purpose of this study is to understand quantitative characteristics of mobile, residual, and dissolved CO 2 trapping mechanisms within ranges of systematic variations in different geologic and hydrologic parameters. For this purpose, we conducted an extensive suite of numerical simulations to evaluate the sensitivities included in these parameters. We generated two-dimensional numerical models representing subsurface porous media with various permutations of vertical and horizontal permeability (k v and k h ), porosity (φ), maximum residual CO 2 saturation (S max gr ), and brine density (ρ br ). Simulation results indicate that residual CO 2 trapping increases proportionally to k v , k h , S max gr and ρ br but is inversely proportional to φ. In addition, the amount of dissolution-trapped CO 2 increases with k v and k h , but does not vary with φ, and decreases with S max gr and ρ br . Additionally, the distance of buoyancy-driven CO 2 migration increases proportionally to k v and ρ br only and is inversely proportional to k h , φ, and S max gr . These complex behaviors occur because the chosen sensitivity parameters perturb the distances of vertical and horizontal CO 2 plume migration, pore volume size, and fraction of trapped CO 2 in both pores and formation fluids. Finally, in an effort to characterize complex relationships among residual CO 2 trapping and buoyancy-driven CO 2 migration, we quantified three characteristic zones. Zone I, expressing 123 808 W. S. Han et al.the variations of S max gr and k h , represents the optimized conditions for geologic CO 2 sequestration. Zone II, showing the variation of φ, would be preferred for secure CO 2 sequestration since CO 2 has less potential to escape from the target formation. In zone III, both residual CO 2 trapping and buoyancy-driven migration distance increase with k v and ρ br .