Deardorff model for estimating energy balance 325 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 VALIDATION OF THE DEARDORFF MODEL FOR ESTIMATING ENERGY BALANCE COMPONENTS FOR A SUGARCANE CROP Glauco de Souza Rolim1*; João Francisco Escobedo1; Amauri Pereira Oliveira2 1 UNESP/FCA - Depto. de Irrigação e Drenagem, C.P. 237 - 18610-307 - Botucatu, SP - Brasil. 2 USP/IAG - Depto. de Ciências Atmosféricas, C.P. 66318 - 05508-900 - São Paulo, SP - Brasil. *Corresponding author ABSTRACT: The quantification of the available energy in the environment is important because it determines photosynthesis, evapotranspiration and, therefore, the final yield of crops. Instruments for measuring the energy balance are costly and indirect estimation alternatives are desirable. This study assessed the Deardorff’s model performance during a cycle of a sugarcane crop in Piracicaba, State of São Paulo, Brazil, in comparison to the aerodynamic method. This mechanistic model simulates the energy fluxes (sensible, latent heat and net radiation) at three levels (atmosphere, canopy and soil) using only air temperature, relative humidity and wind speed measured at a reference level above the canopy, crop leaf area index, and some pre-calibrated parameters (canopy albedo, soil emissivity, atmospheric transmissivity and hydrological characteristics of the soil). The analysis was made for different time scales, insolation conditions and seasons (spring, summer and autumn). Analyzing all data of 15 minute intervals, the model presented good performance for net radiation simulation in different insolations and seasons. The latent heat flux in the atmosphere and the sensible heat flux in the atmosphere did not present differences in comparison to data from the aerodynamic method during the autumn. The sensible heat flux in the soil was poorly simulated by the model due to the poor performance of the soil water balance method. The Deardorff’s model improved in general the flux simulations in comparison to the aerodynamic method when more insolation was available in the environment. Key words: energy budget, net radiation, latent heat flux, sensible heat flux, microclimate VALIDAÇÃO DO MODELO DE DEARDORFF PARA CÁLCULO DO BALANÇO DE ENERGIA DURANTE UM CICLO DE CANA-DE-AÇÚCAR RESUMO: A quantificação da energia disponível no ambiente é importante porque ela afeta a fotossíntese, a evapotranspiração e conseqüentemente a produtividade final dos cultivos. Instrumentos para medidas de balanço de energia são caros e alternativas para estimações são desejáveis. O presente trabalho procura avaliar a performance do modelo de Deardorff (1978) ao longo do desenvolvimento de uma cultura de cana-de-açúcar em Piracicaba, SP, Brasil, em comparação ao método aerodinâmico. Este modelo mecanístico simula os fluxos energéticos (calor sensível, latente e saldo de radiação) em três níveis: a atmosfera, o dossel vegetativo e o solo, usando somente a temperatura do ar, umidade relativa e velocidade do vento medidos num nível de referência acima do dossel, o índice de área foliar e alguns parâmetros previamente calibrados (albedo do dossel, emissividade do solo e a transmissividade atmosférica e características hidrológicas do solo). As análises dos resultados foram feitas em diversas escalas de tempo, condições de insolação, nas diferentes estações do ano (primavera, verão, outono). Analisando todos os dados de 15 minutos, o modelo apresentou boa performance na simulação de radiação líquida em diferentes condições de insolação e estações do ano. O fluxo de calor latente e o fluxo de calor sensível na atmosfera não apresentaram diferenças em comparação ao método aerodinâmico no outono. O calor sensível no solo foi pobremente simulado pelo modelo devido à baixa capacidade de estimação do balanço hídrico do solo. Geralmente as estimações pelo modelo de Deardorff foram melhoradas quando mais insolação era disponível no ambiente. Palavras-chave: balanço energético, radiação líquida, fluxo de calor latente, fluxo de calor sensível, microclima Rolim et al.326 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 INTRODUCTION Studies on energy balance of agricultural ar- eas allow the understanding of energy change pro- cesses in relation to evapotranspiration and photosyn- thesis, which directly affect the biomass accumulation and agricultural yield. Several studies can be cited, as those related to evapotranspiration estimations (Villa Nova et al., 2007), to leaf wetness duration and its consequences on pathogen infection processes (Marta et al., 2007), and to the development of mechanistic crop growth models (Pauwels et al., 2007), among oth- ers. Several methods are available to estimate the energy balance of the environment, also when there is data limitation. The models presented by Shuttleworth & Wallace (1985), Zapata & Martinez- Cob (2002) and Hemakumara et al. (2003) are ex- amples of them. At the mesoclimatic scale, the Deardorff’s (1978) model has been used with suc- cess to simulate energy fluxes (Soares et al., 1996; Targino & Soares, 2002). Vogel et al. (1995) ob- served that Deardorff’s (1978) model presented a better performance simulating the sensible heat flux (in air and soil) than other models, when applied to vegetated surfaces. At the microclimatic scale, few agricultural studies can be found using the Deardorff (1978) model. Tattari et al. (1995) tested Deardorff’s (1978) model and compared it to the Bowen’s ratio method to estimate barley crop evapotranspiration, but the re- sults were not very good due to problems with mea- surements made under adverse climatic conditions, such as high humidity, gusts of wind and rain. How- ever, McCumber (1980) and Garret (1982) stated that among the methods found in the literature, the one pro- posed by Deardorff (1978) is the most efficient to esti- mate energy balances with quality and simplicity. Fi- nally, Oliveira et al. (1999) working with maize and grass crops at Candiota, RS, Brazil, obtained good re- sults in simulations of sensible and latent flux and net radiation. The advantage of using the model proposed by Deardorff (1978) is the significant reduction of the number of instruments required to estimate the energy balance of a given crop, since it only re- quires air temperature, relative humidity, wind speed, leaf area index and some hydrological parameters of the soil. This study analyses the performance of Deardorff’s (1978) model to estimate net radia- tion, latent heat flux, sensible heat flux in the atmo- sphere and in the soil, during a complete cycle of a sugarcane crop, in comparison to the aerodynamic method. MATERIAL AND METHODS Site descriptions and measurements The sugarcane cultivar IAC 87-3396 used in this study was sown on 12/08/2001 with 1 meter row spacing in the east-west orientation. The experimental field is characterized by a flat area of about 1 km2, located in Piracicaba, Sate of São Paulo, Brazil (22º40’ S, 47º38’ W and altitude 514 m). The meteorological data, were collected each second and integrated ev- ery 15 seconds, from 08/17/2001 to 05/31/2002 us- ing a Campbell 21X datalogger and sensors previously calibrated. The sensors were mounted on a tower at the center of the area and the measurements were made at six levels: 5 m from the surface (rain); 2 m from the top of the canopy (temperature, relative humidity, wind speed, net radiation, global and reflected solar radiation); ¾ of canopy height (temperature, relative humidity and wind speed); on the soil surface (heat flux); in the soil, 5 and 20 cm below the surface (tem- perature and humidity). The levels: 2 m from the top of the canopy and ¾ of canopy height were adjusted monthly to follow crop growth. The 2 m level men- tioned above, contains the necessary data for Deardorff’s (1978) model inputs (temperature, relative humidity and wind speed) and the extra observations were necessary for model comparisons. The sugarcane crop was chosen due to its great economic interest in State of São Paulo, besides its long cycle which makes possible performance analyses of the model during several seasons of the year. Determination of the energy balance components by the aerodynamic method The estimation of the energy balance compo- nents (Rn – Net radiation, LE – Latent heat flux, H – sensible heat flux in the atmosphere and G- Sensible heat flux in the soil) by the aerodynamic method us- ing data observed in the field was made to compare the simulated data obtained by the Deardorff's (1978) model. From data obtained every 15 seconds the Richardson number (Ri) was calculated in order to compare the thermal forces responsible for the free convection, with the mechanical forces responsible for the forced convection. The Richardson number was calculated by the following equation: 2)/.( /. dzdum dzdgRi Θ Θ = (1) where: g is the gravitational acceleration (9.8 m s-1), dΘ/dz is the potential temperature gradient (K m-1), du/ dz is the wind speed gradient (s-1). The potential tem- Deardorff model for estimating energy balance 327 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 perature is calculated by Q = air temperature (K).(1000/800)0.2857. In order to estimate the effect of Ri on the tur- bulent diffusion coefficient (Km), which makes pos- sible to quantify the transport of atmospheric activity (sensible and latent heat flux) the function φm was cal- culated for the following cases: unstable conditions (when Ri < 0): φm = (1–16.Ri)-0.25; stable conditions (when Ri > 0): φm = 1 + 5.Ri; and neutral conditions (when Ri = 0): φm = 1 The turbulent diffusion coefficient (Km) is given by: 2 22 . .).( mdz dudzkKm φ −= (2) where: k is the von Karman constant (= 0.41), z the reference sensor height and d is the zero place dis- placement (considered ¾ of the canopy height). The coefficient Km was used to estimate the sensible heat flux: dz dTKmCpH ...ρ−= (3) where: ρ is air density (=1,292 kg m-3), Cp the spe- cific heat of dry air at constant pressure (= 1,005 kJ kg-1 K-1) and dT/dz the temperature gradient between the two sensors. The values of the sensible heat flux in the soil (G) and net radiation (Rn) were obtained directly from field measurements and when H was possible to be calculated, the latent heat flux value (LE) was obtained from: LE = Rn – H – G (4) Deardorff’s (1978) model The Deardorff’s (1978) model comprises sev- eral stages, however some relevant aspects can be briefly described: the model considers that there is only one vegetation layer, whose thermal capacity is insig- nificant and characterized by a coefficient (σf) that is associated with the degree to which the canopy pre- vents the shortwave radiation to reach the soil, which is calculated using the leaf area index (LAI). σf = 0 refers to absent vegetation cover and σf = 1 to com- plete vegetation cover (theoretical). With the coefficient σf and considering the wind logarithmic profile, three heat transfer coeffi- cients can be calculated for the following levels: bare soil surface (CH0 = k2.ln(z/zo)-2, where k = 0.4; z = 2 m above the canopy and zo = roughness parameter); high immediately above the canopy top (CHh = (1/ k).ln[(z-d)/zoh]), where d = zero plane displacement = ¾. Canopy high; zoh = (canopy high – d)/3); soil surface (CHg = (1-σf).CH0 + σf.CHh ). Using the coefficient CHh at the reference level (ua), the mean speed inside the canopy (uaf) is calcu- lated by: uaf = 0.38s f cHh 1/2 ua + (1 – s f ) ua (m s-1) (5) where: ua is the wind speed (m s-1) measured at refer- ence level It was assumed that the inner part of the canopy presents average conditions between the atmo- sphere and the soil, thus the air temperature inside the vegetation (Taf) and the correspondent humidity (qaf) were calculated by the following equations: Taf = (1 – σf) Ta + σf (0.3Ta + 0.6Tf + 0.1Tg) (ºC) (6) qaf = (1 – sf) qa + σf (0.3qa + 0.6qf + 0.1qg) (cm3 cm-3) (7) where: Tf = average leaf temperature, Tg = soil sur- face temperature, Ta = temperature at the reference level (measured) and qaf , qg , qa and qf are analogous specific humidities. In a simple way, Tf, and qf have been estimated from the energy balance at the level of a single leaf inside the canopy. This is an iterative calculation pro- cedure that aims to minimize the differences in the fol- lowing equation: )( LgSgLgSgLhShLhSh RRRRRRRR =−−+−−−+ ↑↑↓↓↑↑↓↓ )( ghsgsh EEHH −+−= λ (8) where: RS, RL, Hs and E are short and long wave ra- diative fluxes, sensible heat flux and evapotranspira- tion, respectively, and l is the evaporation latent heat. The subscribed letters h and g represent values at the canopy top and on the soil surface, respectively. Fi- nally, the arrows indicate the radiative flux directions. When the Stefan-Boltzmann’s equation is used in each side of the equation 8 and other changes mentioned by Deardorff (1978) are introduced, the result is equa- tion 9, which must be solved to find the representa- tive canopy leaf temperature: g gfgf gf LhfShff TRR σ εεεε εε εασ − −+ ++− ↓↓ )( )1[( 4 ffff gfgf gfgf EHsT λσε εεεε εεεε += −+ −+ − ] )( )2( 4 (9) where: αf and εf are the albedo and the foliage emis- sivity, εg is the soil emissivity and s is the Stefan-- Boltzmann constant (5.67 × 10-8 W m-2 K-4). Another peculiarity of Deardorff’s (1978) model is the use of the F-R method (Force Restore) proposed by Bhumralkar (1975) to calculate Tg. It con- sists in a resolution of an equation which depends on Rolim et al.328 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 the energy flux sum in the atmosphere and contains a mechanism that causes an influence of the deep soil layers on the surface temperature. Thus, the soil is di- vided in two layers: a superficial one, for which the soil temperature is influenced by the daytime cycle, and a deep one, of infinite depth in which the temperature changes in an annual scale, as follows: 1 2 1 1 2 τρ TTgc dc Hc t T ss Ag − −−= ∂ ∂ (10) where: T2 is the average temperature of the soil in the deep layer; c1 and c2 are dimensionless constants; ρs is the soil density; cs the soil specific heat; d1 the soil depth under the influence of daytime cycle; and τ1 is the photoperiod; HA is the sum of the atmospheric fluxes at the surface as a result of the energy balance at the soil surface. Finally, the Deardorff’s (1978) model, besides temperature, relative humidity and wind speed data, requires some parameters that were previously cali- brated (Table 1). Basically, the calibration of these pa- rameters aimed to an approach between the measured and simulated data by the model, using the iterative method Simplex with initial values of εf = 0.95; εg = 0.95; αf = 0.2 according to results of Targino & Soares (2002) and τ = 0.7. A total of 12 days were used for this process (4 days of each season with different in- solation conditions), considering the available 249 days. These parameters are similar to those found in literature, as Monteith & Unsworth (1990), who re- lated an emissivity for a sugarcane surface of 97.7 ± 0.4%. Gates (1980), on the other hand, found an al- bedo of 0.17 during the complete development of a given sugarcane crop. Finally, Unsworth & Monteith (1972), who made atmospheric transmissivity measure- ments in England, obtained values of 0.95 for very clear days and 0.5 for days with too much pollution or for cloudy sky conditions. The Deardorff’s (1978) model also requires three soil hydrologic parameters: water contents at wiltting point (wwilt), at field capacity (Wk) and at satu- ration (Wmax). In relation to this a physico-hydric analysis of the soil was accomplished giving empha- sis on the determination of the soil water characteris- tic curve. The results were: wwilt = 0.15; Wk = 0.21; Wmax = 0.32 (cm3 cm-3). Biometrical measurements of the sugarcane crop The parameters that affect the logarithmic wind profile above canopy, leaf area (LA) and aver- age height of the canopy, were measured every 15 days from 08/17/2001 to 05/30/2002, considering five rep- lications. The leaf area index (LAI) was determined by collecting and measuring leaves from a randomized area of 1 m2, using a paper guide with accuracy of 1 mm. Evaluation of Deardorff’s (1978) model The data simulated by the model (SIM): Taf, uaf, canopy relative humidity (RHcanopy), Rn, LE, H, G, were compared to those calculated by the aerodynamic method (AERO) using the index of agreement id sug- gested by Willmott (1981): � 2 1 1 2)( 1 � � � � ��� � �� N i ii N i ii d OOOP OP i � (11) where: id is dimensionless and ranges from 0 (zero) to 1, considering that the value 1 indicates complete agreement between observed and estimated values; Pi are SIM values; Oi are AERO values; the mean of AERO values and N the number of observations. RESULTS AND DISCUSSION Meteorological elements inside the canopy were directly affected by crop development. The in- crease of the leaf area index and canopy height (Fig- ure 3) affected the wind speed inside the canopy (uaf), which was reduced in relation to the reference level (Figure 1C). Because of that, the difference uref - ua was smaller during the spring than in the autumn. Fig- ure 1D shows that soil temperature (5 and 20 cm) was reduced during the crop growth. The five-day moving average of Rn ranged from 500 W m-2 to 700 W m-2 from summer to spring and from 450 W m-2 to 600 W m-2 during autumn (Fig- ure 2). The first summer weeks presented maximum values of LE, up to 490 W m-2, whereas in middle of autumn they reached only 150 W m-2. The maximum LAI and height found for sug- arcane crop were 4.6 and 2.8 m, respectively (Figure 3). These results also correspond to those found in the literature, e.g. Leme et al. (1984), who observed a maximum LAI of 4.5 for the sugarcane cultivar CB47-355, and Bull & Glasziou (1975) and Machado et al. (1982) that also observed that the crop reached a maximum height of 2.6 m. α f ε g =( ε f ) τ gnirpS 71.0 79.0 95.0 remmuS 12.0 69.0 75.0 nmutuA 02.0 89.0 85.0 Table 1 - Calibrated dimensionless parameters used in the Deardorff’s (1978) model: αf - canopy albedo, εg - soil emissivity, εf - canopy emissivity, and t - atmospheric transmissivity. O Deardorff model for estimating energy balance 329 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 Figure 1 - Five-day moving averages of measured meteorological variables inside the canopy (3/4 of canopy height) and at the reference level (2 meters above canopy): A) Temperature, B) Relative humidity and accumulated rain, C) Wind speed, D) Soil temperature (5 and 20cm below the surface). 18 19 20 21 22 23 24 25 26 27 26 6 28 4 30 2 32 0 33 8 35 6 9 27 45 63 81 99 11 7 13 5 Julian Day A ve ra ge a ir te m pe ra tu re , º C Ta Taf Spring Summer Autumn A 40 50 60 70 80 90 100 26 6 28 3 30 0 31 7 33 4 35 1 3 20 37 54 71 88 10 5 12 2 13 9 Julian Day R el at iv e H um id ity , % 0 50 100 150 200 250 300 R ai n, m m Rain RHaf RHa Spring Summer Autumn B 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 26 6 28 3 30 0 31 7 33 4 35 1 3 20 37 54 71 88 10 5 12 2 13 9 Julian Day W in d sp ee d, m s -1 uaf ua Spring Summer Autumn C 30 32 34 36 38 40 26 6 28 4 30 2 32 0 33 8 35 6 9 27 45 63 81 99 11 7 13 5 Julian Day So il Te m pe ra tu re , º C Tsoil 20cm Tsoil 5cm Spring Summer Autumn D Figure 2 - Five-day moving average of maximum hourly values of energy balance components (Rn – Net radiation, LE – Latent heat flux, H – Sensible heat flux at the atmosphere and G – Sensible heat flux in the soil) estimated by the aerodynamic method. 0 100 200 300 400 500 600 700 800 26 6 27 9 29 2 30 5 31 8 33 1 34 4 35 7 5 18 31 44 57 70 83 96 10 9 12 2 13 5 14 8 Julian Day En er ge tic fl ux , W m -2 Rn LE H (-)G Spring Summer Autumn The best meteorological variable simulated by the Deardorff (1978) model was the air temperature inside the canopy (Taf). The smallest errors for all sea- sons were found next to midday. At other periods, the model had a tendency to underestimation, however, mean absolute differences of no more than 3oC were found. As expected, the minimum hourly tempera- ture inside the canopy always occurred between 6h00 and 7h00 a.m., with values around 16ºC during spring and autumn, and around 18ºC during summer, while the mean temperature inside the canopy remained around 27ºC in spring, 30ºC in summer and 28ºC in autumn. Deardorff’s (1978) model also presented a good performance for simulating temperature at the canopy level (Taf) in days with different insolation con- ditions. The id index was always higher than 0.83, de- spite of some variability noticed by R2 values around 0.49, or moreover, by the coefficient of variation (CVSIM), around 22% (Table 2). The wind speed inside the canopy (uaf), mainly in autumn, was well simulated by the model, showing that the parametrization of the wind aerodynamic profile Rolim et al.330 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 Figure 3 - A) Leaf area index variation (LAI); B) canopy height variation of the variety IAC 87-3396 of sugarcane between August 17, 2001 and May 30, 2002. 0 1 2 3 4 5 22 9 24 5 26 3 28 1 29 9 31 7 33 4 35 2 5 23 41 59 Julian day Le af A re a In de x (m 2 m -2 ) Adjusted Observed A 0 1 2 3 22 9 24 5 26 3 28 1 29 9 31 7 33 4 35 2 5 23 41 59 Julian day C an op y he ig ht (m ) B for vegetated surfaces was efficient to represent field conditions. The simulated uaf did not differ from the observed uaf, despite to the low values of R2 (Table 2), since the greatest daily mean differences found were of 0.57 m s-1 in summer, 0.51 m s -1 in autumn and a little higher in spring, reaching 0.78 m s-1. The small values of R2 and high values of id, indicated that the model did not follow instantaneous variations of the wind speed but the magnitude of adjustments in subsequent moments was good. The relationships between AERO and SIM data were less significant in spring, when the sugarcane crop was in the first stages of development. Still considering the canopy level, the last me- teorological variable analyzed was the canopy relative humidity (RHcanopy). The SIM values were not close to those found by AERO due the low values of R2, de- spite of the high values of id, and mainly, because the obtained mean differences were of 8.96%, 12.04% and 6.5% in the spring, summer and autumn, respectively. For all seasons there was a high variability of AERO for 15 minute averages of relative humidity, re- sulting in small R2 values when compared with SIM. The model presented a tendency of underestimation under high relative humidity, mostly before sunrise, and a tendency of overestimation under low relative hu- midity, close to 2h00 p.m. Finally, when SIM and AERO are compared for all meteorological variables when extra insolation was available in the environment, there was an increase of id and R2 (Table 2). In summary, high values of id indicate that SIM values are close to AERO, that is, close to the 1:1 line. On the other hand, the small values of R2 indicate that SIM and AERO data often have a high dispersion. The simulations of soil water content were worse than those at the canopy level, which is related to the low ca- pacity of the Deardorff (1978) model in simulating the soil water balance for successive days, as shown in Figure 4A. The model was not able to simulate the dry- ing of the soil (Figure 4), even though after long drought periods like at the end of spring and autumn, yielding different values from the AERO soil moisture. Thus, 57% of the days in spring and summer and 86% of the days in autumn, the simulated soil moisture in the root zone was above the observed values. Figure 4 - Daily mean soil moisture (A) and soil thermal diffusivity (ks) (B) simulated by Deardorff’s (1978) model (SIM) and calculated by the aerodynamic method (AERO) during spring, summer and autumn seasons. 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 26 6 28 2 29 8 31 4 33 0 34 6 36 2 14 30 46 62 78 94 11 0 12 6 14 2 Julian day SIM Soil moisture OBS Soil moisture Spring Summer AutumnA w at er c on te nt (m 3 m -3 ) 0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 26 6 28 2 29 8 31 4 33 0 34 6 36 2 14 30 46 62 78 94 11 0 12 6 14 2 Julian day So il th er m al d iff us iv ity - K s (c m 2 s -1 ) Ks SIM Ks OBS Spring Summer AutumnB Deardorff model for estimating energy balance 331 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 Ta bl e 2 - D et er m in at io n co ef fic ie nt (R 2 ), W ill m ot t c oe ffi ci en t ( i d), a ve ra ge , s ta nd ar t d ev ia tio n (S D ) a nd c oe ffi ci en t o f v ar ia tio n (C V ) f or 1 5 m in ut e da ta s im ul at ed b y D ea rd or ff ’s (1 97 8) m od el (S IM ) a nd ca lc ul at ed b y th e a er od yn am ic m et ho d (A ER O ) i ns id e t he ca no py : a ir te m pe ra tu re (T af , º C ), re la tiv e h um id ity (R H , % ) a nd w in d sp ee d (u af , m s-1 ) f or d iff er en t s ea so ns o f t he y ea r a nd d iff er en t i ns ol at io n le ve ls . nosae S noitalosnI )ruoh( ºn syad ) Cº( faT,erutarep met ria ypona C u ,deeps dni w ypona C fa s m( 1- ) H R ,ytidi mu H evitale R yponac ) %( i d R 2 .grev A MI S .grev A O R E A D S MI S D S O R E A V C MI S V C O R E A i d R 2 .grev A MI S .grev A O R E A D S MI S D S O R E A V C MI S V C O R E A i d R 2 .grev A MI S .grev A O R E A D S MI S D S O R E A V C MI S V C O R E A gnirp S ll A 98 88.0 86.0 86.12 42.32 42.5 94.4 42.0 91.0 94.0 31.0 45.1 55.0 39.0 54.0 06.0 38.0 86.0 82.0 27.08 67.17 09.41 26.91 81.0 72.0 3 < x = <0 91 68.0 56.0 47.12 56.32 57.4 81.4 22.0 81.0 94.0 80.0 94.1 26.0 59.0 64.0 46.0 47.0 85.0 61.0 02.58 96.37 44.21 90.71 51.0 32.0 7 < x = <3 81 48.0 85.0 56.12 01.32 91.4 99.3 91.0 71.0 94.0 61.0 36.1 45.0 00.1 34.0 16.0 08.0 55.0 61.0 81.58 27.57 50.11 79.61 31.0 22.0 01 < x = <7 62 88.0 76.0 53.12 57.22 94.5 78.4 62.0 12.0 84.0 11.0 15.1 84.0 28.0 63.0 55.0 57.0 46.0 81.0 23.77 74.07 18.41 24.91 91.0 82.0 x = <01 62 09.0 77.0 89.12 45.32 39.5 95.4 72.0 02.0 25.0 81.0 55.1 75.0 59.0 35.0 26.0 49.0 87.0 64.0 77.77 19.86 51.71 74.22 22.0 33.0 re m mu S ll A 98 58.0 85.0 66.22 47.32 16.4 88.3 02.0 61.0 15.0 70.0 40.1 25.0 57.0 83.0 27.0 37.0 04.0 30.0 69.29 29.08 60.9 42.51 01.0 91.0 3 < x = <0 91 78.0 16.0 34.22 15.32 33.4 89.3 91.0 71.0 94.0 70.0 51.1 35.0 28.0 73.0 27.0 07.0 15.0 70.0 59.09 70.08 53.11 79.51 21.0 02.0 7 < x = <3 81 38.0 94.0 39.22 17.32 63.4 17.3 91.0 61.0 94.0 30.0 60.1 65.0 77.0 24.0 37.0 47.0 53.0 10.0 72.39 64.18 69.7 69.31 90.0 71.0 01 < x = <7 62 48.0 75.0 39.22 48.32 50.5 18.3 22.0 61.0 05.0 70.0 40.1 64.0 17.0 13.0 86.0 76.0 43.0 10.0 93.49 84.18 52.7 01.51 80.0 91.0 x = <01 62 88.0 96.0 13.22 68.32 85.4 40.4 12.0 71.0 65.0 31.0 49.0 35.0 86.0 04.0 37.0 57.0 04.0 20.0 48.29 94.08 73.9 69.51 01.0 02.0 m mutu A ll A 98 98.0 07.0 95.12 58.22 09.4 52.4 32.0 91.0 56.0 52.0 32.1 57.0 31.1 27.0 29.0 69.0 67.0 63.0 12.38 17.67 35.71 75.02 12.0 72.0 3 < x = <0 91 68.0 17.0 14.12 22.32 15.4 95.3 12.0 51.0 46.0 12.0 52.1 28.0 80.1 67.0 78.0 29.0 47.0 93.0 68.88 46.08 09.31 06.81 61.0 32.0 7 < x = <3 81 19.0 57.0 23.22 34.32 43.4 20.4 91.0 71.0 66.0 82.0 70.1 66.0 21.1 86.0 50.1 20.1 77.0 14.0 33.68 66.77 90.61 93.02 91.0 62.0 01 < x = <7 62 09.0 86.0 21.12 80.22 30.5 07.4 42.0 12.0 36.0 32.0 34.1 18.0 02.1 67.0 48.0 49.0 17.0 32.0 80.87 18.47 56.61 03.02 12.0 72.0 x = <01 62 88.0 27.0 10.22 34.32 51.5 88.3 32.0 71.0 96.0 13.0 30.1 76.0 10.1 46.0 89.0 96.0 18.0 05.0 25.58 54.67 76.91 29.12 32.0 92.0 egarev A ll A 78.0 56.0 89.12 82.32 29.4 12.4 22.0 81.0 55.0 51.0 72.1 16.0 49.0 25.0 57.0 48.0 16.0 22.0 36.58 64.67 38.31 84.81 61.0 42.0 Rolim et al.332 Sci. Agric. (Piracicaba, Braz.), v.65, n.4, p.325-334, July/August 2008 Soil moisture tends to damp the soil against sud- den changes in temperature, since the heat moves from the soil to the water about 150 times easier than from the soil to the air (Brady, 1989). Therefore, if the simu- lated soil moisture by Deardorff’s (1978) model is higher than expected, the simulated thermal capacity of this soil is also higher than the observed one. Conse- quently, it makes higher temperature variations in the soil more difficult, indicating that the simulated thermal difusivity (Figure 4B) of this soil may be lower as well. This supposition was confirmed by Rolim (2003), with the similar data, because the amplitude of soil tempera- ture variations (on surface and at the 20 cm depth) were actually lower than the observed ones, resulting in low values of id and R2 for all seasons. Despite these differences found between the observed meteorological elements and those simulated by Deardorff’s model, one can see (Table 3) that simu- lations of the energy balance components – Net Ra- diation (Rn), Latent heat flux (LE) and Sensible heat flux at the atmosphere (H) (Table 4) – do not present distinction in relation to the AERO fluxes during all the seasons of the year and in different insolation levels for the sugarcane crop, except the sensible heat flux in the soil (G). Differences between Rn from AERO and the SIM data at every 15 minute interval were small, con- sidering the high values of id and R2 for all seasons (Table 3), with exception to summer (d = 0.79, R2 = 0.50) due to a greater nebulosity and occurrence of rain (Figure 1B), because there was as well a tendency to improve Rn estimates as the insolation availability increased (Table 3). This similarity between AERO and SIM data of Rn is also due, among other reasons, to the fact that Deardorff’s model made a good approach to the total albedo data of the surface (maximum val- ues: spring = 0.34; summer = 0.36 and autumn = 0.33), that directly affects the amount of radiation available in the system soil-vegetation-atmosphere (Figure 5). Relationships between the AERO and SIM data at every 15 minute intervals of latent heat flux in the atmosphere (LE) (Table 3) and sensible heat flux in the atmosphere (H) (Table 4) also were high, taking into account the values of id found for all seasons. The lower values of id during the summer were, on account of the same reasons already discussed in relation to Rn, due to the sensibility of the Deardorf’s (1978) model in relation to the available insolation. Table 3 - Determination Coefficient (R2), Willmott Coefficient (id), average, standart deviation (SD) and Coefficient of Variation (CV) of 15 minute data simulated by Deardorff’s (1978) model (SIM) and calculated by the aerodynamic method (AERO) of Net Radiation (Rn, W m-2) and Latent Heat Flux (LE, W m-2) in different seasons of the year and for different insolation levels. nosaeS noitalosnI )ruoh( ºn syad mW(nR,noitaidaRteN 2- ) mW(EL,xulftaehtnetaL 2- ) i d R 2 .grevA MIS .grevA OREA DS MIS DS OREA VC MIS VC OREA i d R 2 .grevA MIS .grevA OREA DS MIS DS OREA VC MIS VC OREA gnirpS llA 98 78.0 66.0 45.581 91.701 41.062 45.412 04.1 00.2 96.0 05.0 00.551 56.95 93.712 12.301 04.1 37.1 3