ABSTRACT: The dynamics of many viral infections, including rotaviral infections (RIs), are known to have a complex non-linear, non-stationary structure with strong seasonality indicative of virus and host sensitivity to environmental conditions. However, analytical tools suitable for the identification of seasonal peaks are limited. We introduced a two-step procedure to determine seasonal patterns in RI and examined the relationship between daily rates of rotaviral infection and ambient temperature in cold climates in three Russian cities: Chelyabinsk, Yekaterinburg, and Barnaul from 2005 to 2011. We described the structure of temporal variations using a new class of singular spectral analysis (SSA) models based on the "Caterpillar" algorithm. We then fitted Poisson polyharmonic regression (PPHR) models and examined the relationship between daily RI rates and ambient temperature. In SSA models, RI rates reached their seasonal peaks around 24 February, 5 March, and 12 March (i.e., the 55.17 ± 3.21, 64.17 ± 5.12, and 71.11 ± 7.48 day of the year) in Chelyabinsk, Yekaterinburg, and Barnaul, respectively. Yet, in all three cities, the minimum temperature was observed, on average, to be on 15 January, which translates to a lag between the peak in disease incidence and time of temperature minimum of 38-40 days for Chelyabinsk, 45-49 days in Yekaterinburg, and 56-59 days in Barnaul. The proposed approach takes advantage of an accurate description of the time series data offered by the SSA-model coupled with a straightforward interpretation of the PPHR model. By better tailoring analytical methodology to estimate seasonal features and understand the relationships between infection and environmental conditions, regional and global disease forecasting can be further improved.