Estimating Heat Sources And Fluxes In Thermally Stratified Canopy Flows Using Higher-Order Closure Models