In this work, a computationally efficient Bayesian framework for the reduction and characterization of parametric uncertainty in computationally demanding environmental 3-D numerical models has been developed. The framework is based on the combined application of the Stochastic Response Surface Method (SRSM, which generates accurate and computationally efficient statistically equivalent reduced models) and the Markov chain Monte Carlo method. The application selected to demonstrate this framework involves steady state groundwater flow at the U.S. Department of Energy Savannah River Site General Separations Area, modeled using the Subsurface Flow And Contaminant Transport (FACT) code. Input parameter uncertainty, based initially on expert opinion, was found to decrease in all variables of the posterior distribution. The joint posterior distribution obtained was then further used for the final uncertainty analysis of the stream base flows and well location hydraulic head values.