Soil heat flux calculation for sunlit and shaded surfaces under row crops: 1. Model development and sensitivity analysis

Soil heat flux at the surface (G₀) is strongly influenced by whether the soil is shaded or sunlit, and therefore can have large spatial variability for incomplete vegetation cover, such as across the interrows of row crops. Most practical soil-plant-atmosphere energy balance models calculate G₀ as a function of either total (R_N) or soil net radiation (R_N,S). Even though R_N,S includes sunlit and shaded conditions, this is seldom considered, even at spatial scales of a few m. In order to improve the utility of surface energy balance models designed for row crops at relatively small spatial scales, a method was developed to calculate G₀ as a function of shaded, partially sunlit, or fully sunlit R_N,S. Calculation of R_N,S was derived using a geometric approach, and G₀ was derived by the calorimetric method using measurements of soil temperature and volumetric soil water content under upland cotton (Gossypium hirsutum L.) over a wide range of canopy cover conditions. Calorimetric G₀ and calculated R_N,S were related by assuming their normalized values were equal at 24h time steps and in phase. The method required only a single empirical parameter, which related the 24h minimum G₀ (G_0,MIN) to the 24h maximum R_N,S (R_N,S,MAX) as G_0,MIN=a�R_N,S,MAX, and a=-0.31 was found by simple linear regression (p<0.01). Model sensitivity (S_M) of calculated G₀ was calculated for sparse, medium, and full canopy cover; nighttime, shaded, partially sunlit, and fully sunlit surface conditions; and north-south and east-west row orientations, where input values of agronomic, shortwave, and longwave input variables were varied �25% of their base values. The method was most sensitive (1.0<S_M<36) to canopy width, canopy height, leaf area index, row spacing, canopy and soil emittances, and canopy and soil temperatures for medium to full canopy cover. Also, there was generally greater sensitivity for shaded and partially sunlit surfaces compared with sunlit and nighttime surfaces, and north-south rows compared with east-west rows. However, the method had little sensitivity (S_M usually <0.50) to input variables used to calculate the shortwave components of R_N,S.