Three temporary Automatic Weather Stations measured summertime surface layer climate over the Darwin–Hatherton Glacial system. These data were used to test a Polar optimized Weather Research and Forecasting model (Polar-WRF) simulation for December as a case study. Observations show differences in hourly averaged solar and net all-wave radiation between white ice and blue ice areas (BIAs). Although the down-welling solar radiation is higher over the white ice region, the net all-wave energy is higher over the BIA. Derived albedo for each surface type confirms that the blue ice areas have lower albedo. Also, the hourly averaged temperatures are higher at lower elevation stations, creating a gradient towards the Ross Ice Shelf. Analysis shows that there is a diurnal oscillation in strength and intensity of the katabatic wind. The two lower stations register a distinct reversal of wind direction in the early afternoon due to intrusion of an anabatic circulation. Anabatic winds are not prevalent further up the Darwin Glacier. A high-resolution Polar-WRF simulation as a case study shows good agreement with observations. The December 2008 case study is characterized by a strong south-westerly katabatic wind over Hatherton, whereas the flow over Lower Darwin was diurnally reversing. Polar-WRF shows that the katabatic front advanced and retreated periodically between Hatherton and Lower Darwin.

Direct numerical simulation (DNS) is applied to investigate properties of katabatic and anabatic flows along thermally perturbed (in terms of surface buoyancy flux) sloping surfaces in the absence of rotation. Numerical experiments are conducted for homogeneous surface forcings over infinite planar slopes. The simulated flows are the turbulent analogs of the Prandtl (1942) one-dimensional laminar slope flow. The simulated flows achieve quasi-steady periodic regimes at large times, with turbulent fluctuations being modified by persistent low-frequency oscillatory motions with frequency equal to the product of the ambient buoyancy frequency and the sine of the slope angle. These oscillatory wave-type motions result from interactions between turbulence and ambient stable stratification despite the temporal constancy of the surface buoyant forcing. The structure of the mean-flow fields and turbulence statistics in simulated slope flows is analyzed. An integral dynamic similarity constraint for steady slope/wall flows forced by surface buoyancy flux is proposed and quantitatively verified against the DNS data.