Future climate projections provide an opportunity to evaluate cultivar climate classification and preferred styles of wine production for a wine grape growing region. However, ensemble selection must account for downscaled archive model skills and interdependence rather than be arbitrary and subjective. Relatedly, methods for generalizing climate model choice remain uncertain, particularly for identifying optimal ensemble subsets. In this study we consider the complete archive of the thirty-two Coupled Model Intercomparison Project Phase 5 (CMIP5) daily Localized Constructed Analogs (LOCA) downscaled historic datasets and their observational data that were used for downscaling and bias corrections. We apply four model averaging methods to determine optimal ensembles for the computation of six common climate classification indices for the Willamette Valley (WV) American Viticultural Area (AVA). Among the four methods evaluated, elastic-net regularization consistently performed best with identifying optimal ensemble subsets. Variation exists among the optimal ensembles computed for each of the six bioclimatic indices. However, a subset of approximately seven to ten climate models were consistently excluded across all six indices’ ensembles. While specific to the archive and wine region, optimal ensemble sizes were noticeably larger than ensemble sizes commonly employed in published studies. Results are reported such that they can be used by researchers to independently perform analyses involving any one of the six bioclimatic indices throughout the WV AVA while using historic and future LOCA CMIP5 climate projections. The data and methods employed herein are applicable for other wine regions.