This dataset contains shapefiles and rasters that summarize the results of a stochastic analysis of temperatures at depth in the Appalachian Basin states of New York, Pennsylvania, and West Virginia. This analysis provides an update to the temperature-at-depth maps provided in the Geothermal Play Fairway Analysis of the Appalachian Basin (GPFA-AB) Thermal Quality Analysis (GDR repository 879: https://gdr.openei.org/submissions/879). This dataset improves upon the GPFA-AB dataset by considering several additional uncertainties in the temperature-at-depth calculations, including geologic properties and thermal properties. A Monte Carlo analysis of these uncertain properties and the GPFA-AB estimated surface heat flow was used to predict temperatures at depth using a 1-D heat conduction model. In this data submission, temperatures are provided for depths from 1-5 km in 0.5 km increments. The mean, standard deviation, and selected quantiles of temperatures at these depths are provided as shapefiles with attribute tables that contain the data. Rasters are provided for the mean and standard deviation data. Figures and maps that summarize the data are also provided. For the pixel corresponding to Cornell University, Ithaca, NY, a .csv file containing the 10,000 temperature-depth profiles estimated from the Monte Carlo analysis is provided. These data are summarized in a figure containing violin plots that illustrate the probability of obtaining certain temperatures at depths below Cornell.