A computer-implemented method of estimating reservoir attributes as a function of position in a subsurface volume of interest, the method being implemented in a computer system that includes one or more physical computer processors, non-transient electronic storage, and a graphical user interface, comprising: a. obtaining , from the non-transient electronic storage, a geophysical dataset and a geological dataset representative of a subsurface volume of interest; b. obtaining, from the non-transient electronic storage, a parameter model, the parameter model having been trained using training data, wherein the geological dataset includes well data of one or more wells in the subsurface volume of interest and the training data includes the well data; c. picking , with the one or more physical computer processors, at least one surface in the geophysical dataset; d. assigning stratigraphic meaning to the at least one surface based on the geological dataset; e. identifying at least one region of interest on the at least one surface based on the stratigraphic meaning; f. generating , with the one or more physical computer processors, statistical seismic attributes for the at least one region; g. identifying , with the one or more physical computer processors, statistically significant critical attributes among the statistical seismic attributes by applying the parameter model to generate response variable maps for the at least one region; h. generating , with the one or more physical computer processors, uncertainty maps of the at least one region for each of the statistically significant critical attributes and uncertainty maps for the response variables; and i. displaying, via the graphical user interface, the response variable maps and the uncertainty maps, and further comprising using the response variable maps and the uncertainty maps to manage hydrocarbon production from the subsurface volume of interest. a. obtaining , from the non-transient electronic storage, a geophysical dataset and a geological dataset representative of a subsurface volume of interest; b. obtaining, from the non-transient electronic storage, a parameter model, the parameter model having been trained using training data, wherein the geological dataset includes well data of one or more wells in the subsurface volume of interest and the training data includes the well data; c. picking , with the one or more physical computer processors, at least one surface in the geophysical dataset; d. assigning stratigraphic meaning to the at least one surface based on the geological dataset; e. identifying at least one region of interest on the at least one surface based on the stratigraphic meaning; f. generating , with the one or more physical computer processors, statistical seismic attributes for the at least one region; g. identifying , with the one or more physical computer processors, statistically significant critical attributes among the statistical seismic attributes by applying the parameter model to generate response variable maps for the at least one region; h. generating , with the one or more physical computer processors, uncertainty maps of the at least one region for each of the statistically significant critical attributes and uncertainty maps for the response variables; and i. displaying, via the graphical user interface, the response variable maps and the uncertainty maps, and further comprising using the response variable maps and the uncertainty maps to manage hydrocarbon production from the subsurface volume of interest. <EOS>