Modeling Apple Surface Temperature Dynamics Based on Weather Data

The exposure of fruit surfaces to direct sunlight during the summer months can result in sunburn damage. Losses due to sunburn damage are a major economic problem when marketing fresh apples. The objective of this study was to develop and validate a model for simulating fruit surface temperature (FST) dynamics based on energy balance and measured weather data. A series of weather data (air temperature, humidity, solar radiation, and wind speed) was recorded for seven hours between 11:00–18:00 for two months at fifteen minute intervals. To validate the model, the FSTs of “Fuji” apples were monitored using an infrared camera in a natural orchard environment. The FST dynamics were measured using a series of thermal images. For the apples that were completely exposed to the sun, the RMSE of the model for estimating FST was less than 2.0 °C. A sensitivity analysis of the emissivity of the apple surface and the conductance of the fruit surface to water vapour showed that accurate estimations of the apple surface emissivity were important for the model. The validation results showed that the model was capable of accurately describing the thermal performances of apples under different solar radiation intensities. Thus, this model could be used to more accurately estimate the FST relative to estimates that only consider the air temperature. In addition, this model provides useful information for sunburn protection management.


Introduction
Sunburn on apple surfaces is a quality defect that results in culled apples. The exposure of apple surfaces to direct sunlight results in sunburn when the area becomes hot enough to cause cell tissue damage. Sunburn is a major economic problem in marketing fresh apples and can cause significant economic losses for some apple varieties, such as "Fuji", "Jonagold", and "Granny Smith". Many orchardists apply overhead irrigation water to cool the fruit in an attempt to avoid sunburn. Temperature influences many biological processes in fruit development and affects the fruit size, colour, sugar content, acid content, nutritional quality, sunburn injury, and pest development [1][2][3][4][5]. Hence, temperature is an important parameter that affects fruit quality.
Methods for directly measuring the fruit surface temperature (FST) include pushing the bulb of a thermometer beneath the surface of the fruit, inserting thermocouples into the fruit, and placing thermocouples on the fruit surface. In addition, infrared thermometry and thermographic imaging are used for non-contact measurements.
Pushing the sensory bulb of a thermometer beneath the fruit's skin has commonly been used to monitor apple surface temperature [6][7][8][9][10]. However, this process wounds the fruit surface, affects the measurement and does not provide a direct measure of surface temperature. Therefore, this method should not be considered a decisive method for measuring fruit surface temperature [11]. Although the insertion of thermocouples into the fruit could make temperature measurements easier, this process wounds the fruit [12], and the measurements cannot be repeated on the same fruit [13,14] at the same location. To solve these problems, researchers have placed thermocouples on the fruit surface using small fabric adhesive bandages [15,16]. However, the adhesive bandages change the fruit surface conditions, including the reflective radiation properties, and result in less accurate temperature measurements. Furthermore, the above measurements do not consider the spatial heterogeneity of the FST. The gradually changing position of the sun throughout the day causes the location of the highest temperature point on the fruit surface to change. Therefore, to measure the peak fruit surface temperature throughout the day, it is also important to change the measurement location on the fruit surface. Infrared thermometry and thermography can be used for non-destructive and remote determinations of FSTs because they detect the long-wave infrared radiation emitted from the fruit [17] and can be used to dynamically observe the highest temperature location on the fruit as it changes throughout the day. In addition, changes in the FST due to varying solar radiation can change the physiological and biochemical processes of the fruit due to active, natural, protective mechanisms (e.g., induction of HSPs genes). The associated changes in the form of heat balance on the fruit surface can be monitored instantly and remotely by using thermographic imaging [18]. Thus, infrared thermography shows great potential for detecting temperature differences on apple surfaces [2].
Although the temperature dynamics several fruits within a canopy could be monitored via infrared thermometry, thermography, or by an infrared camera, monitoring the temperature dynamics of all fruits throughout the entire growing season remains challenging. Modelling fruit surface temperatures promises to help overcome such difficulties. In addition, FST models could help guide orchardists when applying overhead cooling water. Thus, modelling fruit surface temperatures could save significant amounts of water and pumping energy.
Although modelling fruit surface temperature is meaningful and useful to orchardists, it is rarely performed. In some studies [19,20], air temperature has been considered rather than fruit temperature. However, theoretical and experimental evidence shows that the fruit temperature can be 10 °C higher than the air temperature under sunny conditions [9,10,21]. Evans [22] and Saudreau et al. [8] modelled the temperature dynamics of detached fruits based on climatic factors and used measured fruit temperature data from inserted thermocouples to verify their results. A linear regression of Evans's model with measured results for detached apples resulted in an R 2 value of 0.73. The root mean square error between Saudreau's model estimation and the thermocouple measurements at a 1 min time step was approximately 0.8 °C.
In this study, a model for estimating the highest FST using measured weather data was developed. This model was validated by using data collected from thermal images. Because the solar radiation and maximum FST were highly correlated (r = 0.65) between 11:00 h and 17:00 h [16], this model considered the effects of sun exposure and leaf shading.

Model Description
Energy budgets in a biosphere have been used for various purposes [23], such as modelling animal body temperature [23], analysing human comfort [24], and designing greenhouses [25]. In addition, the temperature of growing fruit is governed by these energy budgets: where ELh is the latent heat, ESh is the sensible heat, EHs is the heat storage, ENp is the energy produced by net plant photosynthesis and ENr is the net radiation. The heat energy balance of the fruit surface during growth is modelled by specifying the normal heat flux at any point on the fruit surface: where Einput is the input energy, Eoutput is the output energy, Rabs is the total incoming radiation (W/m 2 ), Re is the emitted thermal radiation by the fruit (W/m 2 ), E λ is the loss of energy by evaporation (latent heat loss), and H is the loss or gain of sensible energy by convection. Furthermore, f E is the total heat transfer within the fruit, including that from the plant to the fruit, the fruit temperature gradient, and the metabolism activity. The flux of water from the plant to the fruit may induce heat fluxes (Dufour effect) or may occur because of a temperature gradient (Soret effect). However, the Soret and Dufour effects act at large time scales or in response to high temperature gradients, such as those in drying wood.
Regarding the heat flux at the fruit surface, these two heat fluxes are usually small compared with conduction and radiation. Compared with the exchange between the fruit surface and the surrounding environment, the heat release from metabolic activity within the fruit and the energy exchange between the plant and the fruit are small. In addition, f E includes the heat fluxes (Hc) that are produced during heat transfer by conduction (Fourier effect) and can be calculated by using Equation (A.1) in Appendix A. The highest FST represents only one point temperature on the fruit surface at some time, and the apple thermal conductivity is approximately 0.5 Wm −1 K −1 . The heat fluxes (Hc) within the fruit are not considered relative to other fluxes. Thus, f E is assumed to be small enough to be neglected at any particular point in time, and the input energy is approximately equal to the output energy. Therefore, the heat energy balance equation of the fruit surface could be simplified as follows: The difference between Rabs and Re is the net absorbed radiation (Rn) (W/m 2 ). The FST is influenced by many environmental factors, such as solar radiation, air temperature, wind, and humidity. These factors are described below.

Radiation
The total incoming radiation Rabs includes the incoming net short-wave radiation (Rns) and the total incoming net long-wave radiation (Rnl), as defined in the following set of equations: where Rns is the incoming short-wave radiation (Wm −2 ), α is the reflectance of the apple surface (approximately 0.6 [22]), Rnl is the incoming long-wave radiation (Wm −2 ), σ is the Stephan-Boltzmann constant (5.67E −8 Wm −2 T −4 , with T in Kelvin), Ta is the surrounding air temperature (°C), Tgmean is the mean temperature (°C) of the surrounding ground and leaves, A is the estimated maximum projected sunlight on the fruit surface area in the incident sunlit direction, Ar is the projected area exposed to reflected radiation from the ground and canopy, and A' is the total surface area of the fruit (e.g., 2 4πr for a sphere). Furthermore, Ad is an estimate of the projected sunlit fruit surface area in the direction of the incident sunlight. Because this study is used for preventing sunburn, we were primarily interested in estimating the maximum temperature of the fruit surface (i.e., relatively small areas where sunburn occurs). In this study, the Ad/A values were estimated from RGB exposed images (as illustrated in Figure 1). The lighting scenario on the fruit surface could be the surface that is perpendicular to the light direction or is partially or completely shaded by an object [26]. To estimate the maximum temperature of the fruit surfaces under different illumination, the Ad/A ratios were set at 1, 0.8, 0.6, 0.5, 0.4, 0.2, and 0.0. A ratio of 0.0 represents a fruit that is completely hidden by leaves with no exposure to sunlight. The effects of different Ad/A ratios on FST were analysed. In addition, g ε is the emissivity of the ground. The Ar/A' ratio (i.e., exposed to reflected radiation from the ground and canopy) is an empirical parameter (maximum is 0.6). In addition, a ε is the emissivity of the atmosphere and is given by Equations (A.2) and (A.3) in the Appendix.
The emitted radiation ( e R ) is calculated as follows: where f ε is the emissivity of the fruit, 0.95 [8], and fs T is the temperature (°C) at the fruit surface.
where p c is the specific heat capacity of the air at 29.3 J mol −1°C−1 . In addition, a g is the boundary layer conductance of the turbulent heat flow transfer from the fruit surface to the air and can be estimated as where d is a characteristic length estimated by multiplying the fruit diameter by 0.84, which is an empirical parameter when the apple is treated as a sphere [22]. Furthermore, u is the wind speed. The constant 1.4 was empirically determined for application in outdoor environments that were influenced by air temperature and humidity [23].

Transpiration
The latent heat flux ( E λ ) due to transpiration is expressed as follows [23]:

Model Implementation
The general objective of the model is to simulate the dynamics of the highest fruit surface temperatures fs T based on the corresponding variations in the air characteristics and incoming radiation (solar radiation and long-wave radiation). The fs T can be estimated by substituting (Equations (4)- (8)) by (Equation (3)) and combining (Equations (A.2)-(A.9)) to obtain the following a quartic equation: The optical and thermal properties were provided in Section 2.4. The meteorological data for the model are detailed in Section 2.5. The apple characteristics that were required for the model include an estimation of the maximum projected sunlit fruit surface area in the incident sunlit direction and an estimation of the projected sunlit fruit surface area in the direction of the incident sunlight.

Fruit Material
Field data collection was conducted in August and September, 2013, at Washington State University at the Roza Research Farm near Prosser, WA (46.29 N, 119.73 W). Tree rows were arranged in a north-south direction. A set of 41 "Fuji" apples at heights of 1.5 m to 2.5 m on ten trees were randomly selected in the orchard. Thermal images of each apple (an example is shown in Figure 2) were captured using a FLIR A615 infrared camera (FLIR Systems Inc., Portland, OR, USA) to detect the temperature of the fruit surface. These images were calibrated using a BB701 Blackbody Calibrator (OMEGA Engineering, Inc., Stamford, CT, USA).
Apple sunburn normally occurs between 11:00 h and 17:00 h. To acquire information regarding sunburn development, fruit images were acquired from 11:00 h to 18:00 h (midday at 13:00 h) at 15 min intervals in August and from 11:00 h to 17:00 h (midday at 12:00 h) in September to reflect the earlier sunset. The distance of the camera's lens to the monitored apple surfaces was readjusted manually to approximately 0.7 m, and the camera was automatically refocused before each image acquisition. To obtain the highest FST value, the shooting angles were manually adjusted every hour following direct sunlight to maintain consistent angles for the images with the incident sunlight while avoiding casting a shadow over the viewing area. The highest fruit surface temperature measured by the thermal images was more precise than the highest fruit surface temperatures measured by the temperature sensors, including the thermistors and thermocouples. Thermal images could record a large number of temperature points at the same time, ideally allowing users to obtain the maximum temperature on the fruit surface. Figure 2. RGB image (left) and thermal image (right) in the field. In the RGB image, the red region is an estimation of the projected sunlit fruit surface area in the direction of the incident sunlight, (Ad) and the blue region is an estimation of the maximum projected sunlit fruit surface area in the incident sunlit direction (A). In the thermal image, the black region is the fruit region. The RGB images were acquired using a digital colour camera (GC1290, Allied Vision Technologies Inc., Newburyport, MA, USA) from the same location to record the illumination conditions of the monitored apples. The illumination levels were estimated by the ratio between the projected sunlit fruit surface area in the direction of the incident sunlight (Ad) and the maximum projected sunlit fruit surface area in the incident sunlit direction (A). The Ad/A ratio equalled 1, while the ratio was 0.9 to 1, and so on. Based on the exposure image of RGB, the selected exposure region was used to estimate Ad, and the contour region was used to estimate A. Our experiment was conducted at the Roza Research Farm at Washington State University near Prosser, WA (46.29 N, 119.73 W) under mainly sunny conditions. Because sunburn occurs on sunny days, we used the exposure image of RGB to estimate the Ad.
By using thermal image processing software (FLIR EXAMINIR; FLIR Systems Inc.), the fruit surface on the thermal image was manually obtained and the surface temperature (including the mean temperature, maximum temperature and minimum temperature) was obtained from those thermal images. Improved surface temperature measurement accuracy also requires inputs from the surrounding environment, including parameters like distance, atmospheric temperature, and relative humidity. For example, based on the thermal image in Figure 2 (right), measuring the highest FST step includes the following five steps: (1) open the thermal image in the thermal processing software (FLIR EXAMINIR); (2) input the surrounding environmental parameters, including the distance between the fruit and the thermal camera, atmospheric temperature, relative humidity and fruit emissivity; (3) select the region of each fruit, as shown in Figure 2 (right); (4) calculate the highest FST of each fruit region; (5) calculate the highest average FST of the fruits at the same level according to the illumination levels classified by the illumination area based on the image. For example, f1, f2, f3 and f4 are at the same level, and the highest FST at a given level and at a given time is the average value of the highest FST.

Optical and Thermal Properties
Two parameters were used in the model that were related to radiation at the fruit surfaces and must be known prior to simulation. These parameters include the emissivity ( f ε ) for long-wave IR radiation (Equation (4)) and the albedo or reflectance ( α ) for short-wave radiation (Equation (4)). For apples, f ε is generally accepted as 0.94-0.97 [27]. To investigate the influences of these values on the model, the minimum and maximum values of the fruit emissivity ( f ε ) (0.90 and 1.00, respectively) corresponded to relative variations of 5% ± when compared with an f ε of 0.95. The α value is approximately 0.6. However, this value varies with the incidence angle of sunlight, and according to the literature [8,10], its relative variation is approximately ±20% of 0.6.
Another important property for the model is the surface conductance to water vapour diffusion gw. It is difficult to measure gw, and measured gw values are often inaccurate. In this case, a gw value of 5 × 10 −5 ms −1 was used [28,29]. The surroundings have a large influence on gw, which varied from 0 (no transpiration process) to +200%. To investigate the influences of these parameters on the model predictions, a sensitivity analysis of these parameters was conducted in Section 3.2.2.

Measurements of the Climatic Variables
An automatic weather station [30] located approximately 400 m upwind from the testing orchard provided meteorological data, including the air temperature, relative humidity, wind speed, wind direction and soil temperature at 2 m above the ground. The short-wave radiation Rs, from which the direct radiation was estimated, was recorded using a CS300-L Pyranometer (CS300, Campbell Scientific, Inc. Logan, UT, USA) every 10 s, and the averaged values of the 15 min readings (90 data points) were used as the measured data.

Difference Comparisons of Air Temperature and FST
Evidence suggests that fruits can be heated to temperatures that are significantly higher than air temperatures in the field. Figure 3 presents the difference between the air temperature and the highest fruit surface temperature under different area ratios of illuminated surface at different times of the day. According to Figure 3, this difference varied from 0 to 12.0 °C during the 11:00 h to 18:00 h time window. The largest differences reached 12.0 °C around solar noon for the fruit with an Ad/A ratio near 1.0, 9.8 °C for the fruits with an Ad/A ratio of 0.5, and 3.2 °C for the fruits with an Ad/A ratio of 0. These results suggest that different proportions of illumination largely influence the FST. The maximum differences occurred at approximately 14:00 h (i.e., one hour after solar noon), mainly because the solar radiation was highest between midday and two hours after midday, which would cause the highest FST. However, the validation measurement indicated the FST reached its peak value if the wind speed was not very high. In almost all of our field measurements, we did not observed any cases where high wind affected the peak FST value. The FST peaked while the air temperature did not significantly increase. Next, the temperature differences gradually decreased. For shaded fruit (Ad/A = 0), the total difference between air temperature and fruit surface temperature was not notable, and the fruit surface temperature was estimated using the air temperature. Figure 3. The mean difference between the hottest FST and the air temperature from 11:00 h to 18:00 h for three illumination levers (Ad/A = 0, 0.5, 1).  Figure 4. The RMSE was less than 2.0 °C (Figures 4 and 5). Compared with previous analytical studies (R 2 of 0.73) [8,22], the presented study estimated peak FST more accurately. The use of thermal images more precisely measured the effects of a realistic environment (i.e., unsteady and non-homogeneous heat fluxes) on the peak temperature dynamics of the fruit surface.

Model Evaluation
To highlight the temperature dynamics of the fruit surface, the highest FSTs under Ad/A values of 1, 0.5 and 0 on DOY 232 are shown in Figure 5. The largest standard deviation was 1.9 °C within the allowable error of 2.0 °C, which validated the overall accuracy of the developed model for simulating  FST. Statistically, we observed that the simulated FST was more precise when the Ad/A value was 1 for an Ar/A' value of 0.6, when the Ad/A value was 0.5 for an Ar/A value of 0.3, and when the Ad/A value was 0 for an Ar/A value of 0.6.    Figure 6 presents the results of a sensitivity analysis of the model with selected fruit physical parameters. A negative value indicates a decrease in temperature for a corresponding decrease in the input parameter or vice versa. The reference fruit surface temperature (Tref) of the output value of the model is shown. While some parameters were changed, other parameters were held constant at reference values of α = 0.6, gw = 5 × 10 −5 ms −1 , and f ε = 0.95. The changes of α ( α : 20% ± ) resulted in larger changes in FST (i.e., approximately 2.0 °C). Changes in physical parameters, such as the emissivity of the fruit ( f ε : 5% ± ) and the surface conductance to water vapour diffusion (gw:0%-+200%), only resulted in minor variations in FST (less than 0.5 °C). When the emissivity ( f ε ) was increased by 5%, the fruit emitted more energy, which decreased the FST. When the fruit surface reflectance α was reduced by 20%, the amount of energy absorbed by the fruit increased, which caused the FST to increase. Additionally, when the surface conductance to water vapour diffusion gw increased by 200%, the evaporative cooling increased and caused a decrease in FST. This situation was reversed when gw was set to 0.

Figure 6.
Sensitivity analysis of the model to fruit physical parameters (the surface conductance to water vapour diffusion gw, the fruit surface reflectance α and the emissivity ( f ε ) of the fruit surface).
As expected, the energy received by the apples from direct sunlight is a major factor that affects FST [31]. As shown in Figure 7, different illumination areas largely influenced the peak FST. The reference temperature corresponded to the FST when the Ad/A ratio was equal to 1. When the Ad/A ratio was equal to 0.6 or 0.4, the temperature difference was 1.5 °C. The difference between all of the shaded areas (i.e., Ad/A = 0) and all of the sunny areas (i.e., Ad/A = 1) reached 7.4 °C.

Impact Analysis of Fruit Temperature and Thermal Parameters on FST
The model can be used to analyse the relevance of different thermal parameters (the surface conductance to water vapour diffusion gw, and the emissivity ( f ε ) of fruit surface) for predicting fruit temperature by using a sensitivity analysis. This point is critical because it shows which parameter is relevant when making decisions regarding conducting potential experiments. The sensitivity analysis ( Figure 6) showed that only the parameters involved in the heat fluxes at the fruit surface were required to correctly predict the temperature on the fruit surface. The values of gw obtained from previous studies [28] were not well characterised because the relative variations could reach +200%. The three following factors could explain this large variation: (1) the difficulty in measuring gw because its order of magnitude was 0.001 ms −1 ; (2) the deduction of this parameter from the VPD and other parameters according to (Equation (8)) (thus accurate measurements were not straightforward [29]); and (3) the gw value varied with the fruit growth conditions and with fruit development [28]. As expected, increasing gw (+200% of the original value) enhanced the transpiration process and cooled the entire fruit. This scenario decreased the FST by 0.2 °C. Although increasing gw by +200% resulted in the apparent surface conductance to water vapour diffusion, no substantial decrease in the FST was observed ( Figure 6). This finding suggested that the evaporation process could be neglected.
The sensitivity analysis revealed that the fruit emissivity f ε was also an important parameter in the FST. Increasing f ε (+5%) resulted in a decrease (−0.8 °C) in the FST. In contrast with the gw value, the f ε value was accurately estimated in previous studies [32]. The value of f ε increased from 0.9 to 1.0, which resulted in a 1.5 °C decrease in FST. Thus, fruit emissivity had a relatively substantial effect on the FST compared with the gw.
The sensitivity of the model to thermal parameters (gw and f ε ) was relatively low. This phenomena was useful for investigating the effects of microclimates on fruit responses in terms of physiological disorders. Moreover, the FST is not sensitive to thermal parameters (gw and f ε ). Thus, the amount of energy received or released, as well as the time and length scale ratios between the heat fluxes at the fruit surface, were very important to the FST.

Relationship between Microclimate and FST
To investigate the relationships between the environment and the FSTs, all of the heat fluxes at the fruit surface were calculated based on the simulated temperature of the model and by using (Equations (4), (6) and (8)). Figure 8 shows the amplitudes of the heat fluxes, including the total incoming radiation (Rabs), convection (H), transpiration ( E λ ), and emitted radiation of the apple surface (Re) on DOY 232. A positive value of heat flux corresponds to a gain of heat (i.e., the surface becomes warmer). A negative value of the heat flux corresponds to a loss of heat (i.e., the surface becomes cooler).
The loss of energy caused by evaporation (i.e., latent heat loss) was negligible and smaller than the loss caused by the total incoming radiation and convection processes. Between 11:00 h and 18:00 h, the changes in the latent heat loss and emitted radiation were relatively flat. A similar trend was observed regarding changes in the heat fluxes on the fruit surface for the other days observed. In addition to estimating the FST dynamics, the model simulated the heat flux dynamics on the fruit surface. Radiation and convection govern the dynamics of fruit surface temperature, and should be accurately implemented. The results validated the methods of estimating FST, as shown by (Equations (4), (8) and (14)).
The forced convection used in the model (Equation (8)) was caused by the difference between the FST and air temperature. Radiation and convection remained the driving processes that governed the FST and were substantially affected by the distribution of tree foliage. Therefore, the characterisation or modelling of the air within a tree canopy is critical for accurate FST modelling.
In addition, the maximum FST was only simulated on clear, sunny days because fruit sunburn occurs under these conditions. In future research, the model will need to be modified to estimate the FST under different weather conditions (such as cloudy, overcast, or rainy conditions). In addition, the influences

Conclusions
A physical model simulating the variations of the maximum fruit surface temperature was developed in this study. The model was based on an energy balance and field data. Model outputs were compared with fruit surface temperatures that were measured from thermal images. The linear regression of the simulated FST and the measured FST had a coefficient of determination (R 2 ) of 0.90 and a mean absolute error of less than 2.0 °C. The obtained results showed that the model accurately described the thermal performance of the apples and was more accurate than the estimates that are usually made based on air temperature. This finding suggested that the model could be used to study microclimatic effects on fruits, such as sunburn, and to provide useful information for sunburn protection management. Next, this model should be used to define temperature and time thresholds to automate fruit evaporative cooling systems and avoid sunburn. where a e is the actual vapour pressure (kPa). The actual vapour pressure is the saturation vapour pressure at the dew point temperature (Tdew) (°C) and is expressed as follows: The air density air ρ can be determined using the following equation [34]: where P is the atmospheric pressure at elevation Z (kPa), P0 is the atmospheric pressure at sea level (101.3 kPa), Z is the elevation (m), Z0 is the elevation at the reference level (m), g is the gravitational acceleration (9.807 m s −2 ), R is the specific gas constant (287 J kg −1 K −1 ), 1 α is the constant lapse rate of moist air (0.0065 K m −1 ), and Tk0 is the reference temperature (k) at elevation Z0. The elevation Z0 is generally considered as sea level for standardised calculations.
The American Society of Civil Engineers (ASCE) assumed the standard atmosphere to be P0 = 101.3 kPa at Z0 = 0 and Tk0 = 293 K for Tmean = 20. Therefore,(Equation (A.5)) simplifies to 5 where Δ is the increasing rate of the saturation vapour pressure at the dew point temperature (kPa°C −1 ).

Conflicts of Interest
The authors declare no conflict of interest.