Simulations of Drainage Flows with Topographic Shading and Surface Physics Inform Analytical Models

We perform large-eddy simulations (LESs) with realistic radiation, including topographic shading, and an advanced land surface model to investigate drainage flow dynamics in an idealized compound-slope mountain geometry. This allows an analysis not only of fully developed profiles in steady state—the subject of existing analytical solutions—but also of transient two- and three-dimensional dynamics. The evening onset of downslope flow is related to the duration of shadow front propagation along the eastern slopes, for which an analytic form is derived. We demonstrate that the flow response to this radiation pattern is mediated by the thermal inertia of the land through sensitivity to soil moisture. Onset timing differences on opposite sides of the peak are explained by convective structures that persist after sunset over the western slopes when topographic shading is considered. Although these preceding convective systems, as well as the presence of neighboring terrain, inhibit the initial development of drainage flows, the LES develops an approximately steady-state, fully developed flow over the finite slopes and finite nocturnal period. This allows a comparison to analytical models restricted to such cases. New analytical solutions based on surface heat flux boundary conditions, which can be estimated by the coupled land surface model, suggest the need for improved representation of the eddy diffusivity for analytical models of drainage flows.