Traumatic brain injury such as that developed as a consequence of blast is a complex injury with a broad range of symptoms and disabilities. Computational models of brain biomechanics hold promise for illuminating the mechanics of traumatic brain injury (TBI) and for developing preventive devices. However, reliable material parameters are needed for models to be predictive. Unfortunately, the properties of human brain tissue are difficult to measure, and the bulk modulus of brain tissue in particular is not well-characterized. Thus, a wide range of bulk modulus values are used in computational models of brain biomechanics, spanning up to three orders of magnitude in the differences between values. However, the sensitivity of these variations on computational predictions is not known. In this work, we study the sensitivity of a 3D computational human head model to various bulk modulus values. A subject-specific human head model was constructed from T1-weighted MRI images at 2 mm voxel resolution. Diffusion tensor imaging provided data on spatial distribution and orientation of axonal fiber-bundles for modeling white-matter anisotropy. Non-injurious, full-field brain deformations in a human volunteer were used to assess the simulated predictions. The comparison suggests that a bulk modulus value on the order of GPa gives the best agreement with experimentally measured deformation in the human brain. Further, simulations of injurious loading suggest that bulk modulus values on the order of GPa provide the closest match with the clinical findings in terms of predicated injured regions and extent of injury.