This paper presents a novel study on boxplot-valued data in climatological applications. Our methodologies are applied to the Berkeley Earth Surface Temperature Study. We validate our approaches through comprehensive simulations, comparing Bayesian and frequentist estimators for efficiency and accuracy. The results provide robust insights into climatic trends, particularly summer average temperatures across European countries.