Este estudo teve como objetivo elaborar uma metodologia para identificação das ondas de frio (OF) mais intensas ocorridas no centro-norte do Estado de São Paulo entre 1961 e 2017. Para tanto, foram usados dados de temperatura mínima e máxima medidos em 10 estações meteorológicas convencionais. O método de identificação das OF consistiu em compilar as informações das 10 séries temporais em apenas uma e na aplicação de um percentil muito restritivo, de 0,2%, para seleção dos casos mais intensos de OF. Foram identificadas 12 OF intensas, com duração média de 5 dias. Todos os casos provocaram geada na região de estudo, que muitas vezes causaram perdas na agricultura. Para descrição do padrão sinótico da atmosfera associado com os casos de OF, foi realizada uma análise de composição (média de uma dada variável atmosférica em datas específicas). Em superfície, foi observado que a incursão de ar polar para latitudes mais baixas ocorre devido às características de circulação do Anticiclone Subtropical do Pacífico Sul (ASPS) e de uma baixa pressão no oceano Atlântico, que interagem favorecendo a advecção de ar frio. Em médios e altos níveis da atmosfera, observa-se um padrão de onda amplificado, com crista no Pacífico e cavado no Atlântico e América do Sul, que propicia o avanço do ar frio para o norte e dá suporte dinâmico ao centro de alta pressão migratório.