Modeling atmospheric CO2 concentration profiles and fluxes above sloping terrain at a boreal site

CO2 fluxes and concentrations were simulated in the planetary boundary layer above subarctic hilly terrain using a three dimensional model. The model solves the transport equations in the local scale and includes a vegetation sub-model. A WMO/GAW background concentration measurement site and an ecosystem flux measurement site are located inside the modeled region at a hilltop and above a mixed boreal forest, respectively. According to model results, the concentration measurement at the hill site was representative for continental background. However, this was not the case for the whole model domain. Concentration at few meters above active vegetation represented mainly local variation. Local variation became inseparable from the regional signal at about 60-100 m above ground. Flow over hills changed profiles of environmental variables and height of inversion layer, however CO2 profiles were more affected by upwind land use than topography. The hill site was above boundary layer during night and inside boundary layer during daytime. The CO2 input from model lateral boundaries dominated in both cases. Daily variation in the CO2 assimilation rate was clearly seen in the CO2 profiles. Concentration difference between the hill site and the forest site was about 5 ppm during afternoon according to both model and measurements. The average modeled flux to the whole model region was about 40% of measured and modeled local flux at the forest site.