Optimal Aggregation of Prediction Intervals under Unsupervised Domain Shift