Multiple linear regression models were developed for 1–3-day lead forecasts of maximum and minimum temperature for two locations in the city of Lima, on the central coast of Peru (128S), and contrasted with the operational forecasts issued by the National Meteorological and Hydrological Service-SENAMHI and the output of a regional numerical atmospheric model. We developed empirical models, fitted to data from the 2000–13 period, and verified their skill for the 2014–19 period. Since El Niño produces a strong low-frequency signal, the models focus on the high-frequency weather and subseasonal variability (60-day cutoff). The empirical models outperformed the operational forecasts and the numerical model. For instance, the high-frequency annual correlation coefficient and root-mean-square error (RMSE) for the 1-day lead forecasts were 0.378–0.538 and 0.748–1.768C for the empirical model, respectively, but from around-0.058 to 0.248 and 0.888–4.218C in the operational case. Only three predictors were considered for the models, including persistence and large-scale atmospheric indices. Contrary to our belief, the model skill was lowest for the austral winter (June–August), when the extratropical influence is largest, suggesting an enhanced role of local effects. Including local specific humidity as a predictor for minimum temperature at the inland location substantially increased the skill and reduced its seasonality. There were cases in which both the operational and empirical forecast had similar strong errors and we suggest mesoscale circulations, such as the low-level cyclonic vortex over the ocean, as the culprit. Incorporating such information could be valuable for improving the forecasts.