Abstract. Climate indicators are useful tools to synthesise climate information from multiple station time-series into a single national indicator. The method applied should be spatially representative and robust over time. We introduce a new method, based on Empirical Orthogonal Functions (EOF) during the calibration period 1961–2018, in order to reconstruct the climate indicator for temperature in Sweden for the full 1860–2020 period of available observations. The new method delivers results in good overall agreement with the reference method (i.e. arithmetic mean from selected stations in the reference network). Discrepancies are found prior to 1900, primarily related to the reduced number of active stations: the robustness of the indicator estimation is assessed by an ensemble computation with added random noise, which confirms that the ensemble spread increases significantly prior to 1880. The present study establishes that the 10-year running averaged temperature indicator rose from −1.03 °C in 1903 to +1.19 °C in 2010 (with respect to a mean value of 4.64 °C over the 1961–2018 calibration period), i.e. an increase by +2.22 °C in a century. The temperature difference between 1860 and 2020 was largest for winter (DJF) averages (+3 °C) and minimal for summer (+2 °C). The leading EOF patterns illustrate the spatial modes of variability for climate variability, with a predominantly homogeneous, mono-modal distribution for temperature. For precipitation, the first EOF pattern displays more pronounced regional features (maximum over the West coast), which is completed by a north-south seesaw pattern for the second EOF. We illustrate that EOF patterns calculated from observation data reproduce the major features of EOF calculated from GridClim, a gridded data-set over Sweden, for annual and seasonal averages. The leading EOF patterns vary significantly for seasonal averages (DJF, MAM, JJA, SON) for temperature. Finally, future developments of the EOF-method are discussed for calculating regional aggregated climate indicators, their relationship to synoptic circulation patterns and the benefits of homogenisation of observation series. The EOF-based method to compute a spatially aggregated indicator for precipitation is presented in a companion article (Sturm, 2024a). The code and data for this study is available on Zenodo (Sturm, 2024b).