Two-dimensional simulation of airflow and carbon dioxide transport over a forested mountain: Part I: Interactions between thermally-forced circulations

Local-scale wind regimes over an idealized two-dimensional (2D) mountain having similar horizontal and vertical scales as Vancouver Island were investigated using a high-resolution mesoscale model. The model-generated flow outputs were then used as ersatz data to assess the impact of limiting assumptions in the eddy-covariance (EC) method as well as the effectiveness of tower flux data time-filtering for the main tower site of the Fluxnet-Canada British Columbia flux station. In this paper, we describe the simulated mesoscale and local-scale flow regimes, and in Part II we describe their use in assessing tower-data analyses of CO2 fluxes. The numerical model was enhanced to include parameterizations for: tree drag, radiation effects of forest canopies on the surface energy budget, and soil heat conduction. Simulations were performed both over an idealized bare hill and over a forested hill. The simulated flow involved interactions between the land/sea breeze, convective thermals, and mountain circulations under fair-weather conditions. The resulting simulated winds over the forested slope were much weaker than those over a bare slope. The nocturnal drainage flow over the forested slope was separated into sub-canopy and above-canopy regimes due to the temperature inversion at treetop. The strongest downslope flows occurred above the canopy with a minimum velocity occurring in the upper canopy where canopy drag is the strongest. The overturning of large convective eddies during daytime resulted in intermittent downslope flows under the canopy during the day. Wind regimes had a rapid shift from upslope to downslope flow above the canopy just after sunset. Further testing will be made using 3D simulations with higher resolution.