We study the flow in the lower atmosphere in response to a sudden surface cooling. Our numerical approach aims to resolve all relevant scales of the flow but only along a vertical column. Complex turbulent motions are modeled by simple random mappings. We show that the numerical model accurately captures some relevant features of near-surface turbulent winds and temperature fluctuations. The model offers new opportunities for atmospheric chemistry and polar boundary layer application cases.