Adaptive Management Framework for Evaluating and Adjusting Microclimate Parameters in Tropical Greenhouse Crop Production Systems

High operational costs of greenhouse production in hot and humid climate condition due to the initial investments on structure, equipment, and energy necessitate practicing advanced techniques for more efficient use of available resources. This chapter describes design and concepts of an adaptive management framework for evaluating and adjusting optimality degrees and comfort ratios of microclimate parameters, as well as predicting the expected yield in greenhouse cultivation of tomato. A systematic approach is presented for automatic data collection and processing with the objective to produce knowledge-based information in achieving optimum microclimate for high-quality and high-yield tomato. Applications of relevant computer models are demonstrated through case-study examples for use in an iterative way to simulate and compare different scenarios. The presented framework can contribute to future studies for providing best management decisions such as site selection, optimum growing season, scheduling efficiencies, energy management with different climate control systems, and risk assessments associated with each task.


Introduction
The increasing market demand for high-quality food products have replaced open-field cultivations of Solanaceae and Cucurbits crops with modern plant production systems for more efficient use of available resources.Closed-field cultivations by means of commercial greenhouses have been changed over the last three decades from basic structures to advanced controlled environments for optimizing plant's productivity and producing high yields at low expenses.The higher costs of greenhouse operation due to the initial investments on structure, equipment, and energy necessitate practicing advanced techniques of automation for efficient control of the microenvironment.Research trends in this field are directed toward developing innovative solutions for shifting from energy-consuming to energy neutral greenhouses with the ultimate objective of increasing profits.This is, however, challenging due to the lack of accurate information about interactions between crops and environment at different growth stages, as well as the complexity of the dynamic system that is subjected to changes with internal and external factors.Plant-based engineering has helped researchers with proper management policies to embrace these uncertainties through modeling and integrated-learning approaches.Several uncertainties with greenhouse cultivation include climate variability, expected yield, optimum references of microclimate parameters, comfort ratios, insecurity of resources, complexity of the system states, lack of accurate information about interactions between plants and environment, and the relationships between biological and ecological system.
Greenhouse microclimate control has been a large field of study for many years.Much work has been done for moderate cold climate conditions as opposed to tropical lowlands.In contrast to cold arid climate, the main objective of a greenhouse in hot and humid regions such as lowlands of south-east Asia (Figure 1) is to protect crop against fluctuations of external conditions such as extreme winds, heavy seasonal rainfalls, typhoons, extreme solar radiation, occasional water shortage, high air temperature, high humidity, and invasion of pests and diseases.The major concern with greenhouses in these regions is the crop stress due to the adverse microclimate that reduces plant evapotranspiration rate and causes production failure.Evaporative cooling systems by means of misting, pad-and-fan, and swamp cooling are widely used in tropical greenhouses of south-east Asia for manipulating crop growth microclimate; however, these systems have not reached their optimum potential due to their conventional automation and control methods.If properly managed, tropical greenhouses can provide suitable growth condition for tomato cultivation by maintaining inside microclimate close to the outside, with an expected yield that varies between 30 and over 100 tons/ha (vs.open-field yield of 15-30 tons/ha) depending on soil culture or hydroponics medium.
Profitability and investment returns of commercial greenhouses are tightly linked to management decisions.One of the main factors to be considered in this context is the sustainability of Plant Engineering operations through proper management of available resources.Modern greenhouses are required to exhibit integration of automation, cultural practices, and environmental control using object-oriented analysis of the subsystems.The primary concepts and methods of automation-culture-environment system analysis (ACESYS) in controlled environments plant production (CEPP) have been introduced and expanded in the works of [1][2][3].Some of the earliest examples of object-oriented analysis and modeling applications including optimal control strategies and decision-support software in advanced CEPP systems can be found in the works of [2,[4][5][6].The purpose of object-oriented system analysis approach according to Ref. [3] is to develop a set of foundation classes that can be used to effectively describe the components of the automation system.This, however, requires a comprehensive understanding of the interaction between crop's growth response and environment characteristics.Some of the specific applications and benefits of system analysis in greenhouse production includes integrated energy-efficient strategies, extracting unique and new knowledge that provides valuable insight to local growers and beyond, understanding limitations of resources and balancing between input and output expectancies, improving technology and increasing returns, providing business attraction for local stakeholders, minimizing energy requirements and eliminating tedious operations, increasing production quality and quantity to satisfy market demands, and technology adaptation by balancing between fixed and flexible automation for various crop production.With this perspective, the convolution of several possible scenarios and combination of culture classes (i.e., climate control parameters) and objects (i.e., tomato crop at different growth stages) necessitates computer-based analysis program within the concepts of a systematic framework approach such as adaptive management.
Adaptive management was initially introduced at the University of Florida [7] as an iterative method for managing natural resource in the systems with wide range of responses to management choices and to help manager's difficulty in understanding the systems' dynamics [8] and plant's responses [9].It is defined as "a systematic process for continually improving management policies and practices by learning from the outcomes of operational programs" [10].Adaptive management was created based on the needs of environment and ecosystem managers with an iterative processing tool that acknowledged complexity and uncertainty, with a focus on learning and for continuous inputs [11].It has been widely used as a new design technique for large database that manages and assists the immense data collection, data analysis, and data storage of distributed sets of experiments associated with environmental, meteorological, biological, and medical research problems or other technical and experimental assessments that utilize large-scale data sources within multiple and separate engineering or laboratory facilities.Examples include the work of Refs.[12][13][14][15].The principles of adaptive management according to Ref. [16] suggest using the best available knowledge to design and implement management plans, while establishing an institutional structure that enables learning from outcomes to adjust and improve future decision making.This structured approach is an efficient method in developing decision-support tools for systems design, management, and operation by recognizing the importance of natural variability in contributing to ecological resilience and productivity.
This chapter provides a systematic process of incorporating new and existing knowledge that can be used in developing management decisions for achieving optimum microclimate.
It describes design and concepts of an adaptive management framework for evaluating optimality degrees and comfort ratios of air temperature (T), relative humidity (RH), and vapor pressure deficit (VPD), as well as prediction of the expected yield in greenhouse cultivation of tomato.The presented framework was designed to allow production managers to ask "what-if" type of questions for further quantitative inclusion to avoid possible detriment decisions.It also provides an in-depth rigorous analysis tool for decision making or decision procrastination when facing uncertainties.It can assist in enhancing scheduling efficiency and guiding investments through different simulated scenarios that are based on information analysis to support optimal restoration strategies.In the rest of this chapter, we refer to the term "microclimate parameters T, RH, or VPD" by M. We also use Opt M ðÞ ¼ α, and Cf t M,t,α ðÞ ¼ β, to refer to the terms "optimality degree" and "comfort ratio," respectively, defined in Sections 3.1 and 3.2.

Adaptive greenhouse model
The key to an adaptive greenhouse is the computer model that drives specific implementations of other components.An adaptive management framework for microclimate evaluation and control in greenhouse production systems is proposed in Figure 2. A diagram of the steps in the analysis process is shown by rectangles.The arrows are the direction of the process, and the central spiral highlights the goal of arriving at a compromising decision based on a shared set of objectives developed through the iterative process.The three essential elements in this structure are (i) data entry and retrieval, (ii) computer model (expressed by mathematical equations), and (iii) data analysis components.The data entry component may be implemented by direct interfacing with real-time data acquisition system or by using web-based and desktop application software.Computer model is application specific; it can be updated and is usually condensed and produced from previous extensive research works in crop physiology.Data analysis comprises implementing relevant techniques within the retrieval component (i.e., programmable spread sheets) or by integrating with third-party applications (i.e., Simulink blocks).The proposed framework can be adapted to new research projects for working with different culture classes and objects by which many specific scenarios may be modeled and analyzed.It carefully monitors the possible outcomes of the system for better understanding of the process in order to adjust control parameters through an iterative learning process.
The framework utilizes a custom-designed data acquisition, and control system [17] that was built using Arduino Uno prototype microcontroller board for monitoring and manipulating of the microclimate parameters.Three computer models were employed by the framework for evaluation and adjusting of optimality degrees Opt M ðÞ , comfort ratio Cf t M,t,α ðÞ , and prediction of the expected yield.The framework was implemented in MATLAB® (The MathWorks Inc, Natick, MA, USA) environment through Simulink blocks and coding of various main functions and sub-functions that were stored as "m-files."Different toolboxes were developed for the immense data-analyzing tasks as shown in Figure 3.The framework structure was designed in a way that end users can create (or update) entries in database, select report type (1-day or multidays report), and proceed with a specific analysis procedure.The database is a dynamic flat file type that can be created by entering collected data, either manually from previously stored sources such as excel sheets or directly from the hardware interface.The computer models presented in this chapter are focused on tomato (Lycopersicon esculentum); however, with slight modification, the framework can be reprogrammed to work with other greenhouse crops provided that their yield prediction and growth response models are available.Results of microclimate evaluation and set-point manipulation discussed in Sections 3 and 4 can contribute to dynamic greenhouse climate control strategies [18] such as the one in Ref. [19].An example is provided by comparing a model reference-adaptive greenhouse microclimate controller with conventional closed-loop feedback shown in  Plant Engineering parameter (i.e., air temperature), causing other microclimate parameter (i.e., humidity) to be actually controlled via set-point manipulation.

Optimality degrees of microclimate
Optimality degree of a microclimate parameter denoted by Opt M ðÞ ¼ α is a quantitative value between 0 and 1 that represents how close a microclimate measurement (T, RH, or VPD) is to its ideal value as required by the greenhouse crop at specific growth stage and climate condition.This value can be computed from experimental models that correlate different levels of microclimate parameters with yield and quality of the greenhouse crop.An example of such models is the one developed for air temperature and relative humidity by the Ohio Agricultural Research and Development Center [20,21].These models define optimality degrees of T, RH for greenhouse cultivation of tomato with independent trapezoid membership-function growth response plots that are specific for different growth stages and three light conditions (night, sun, and cloud).These plots were originated using utility theory with the goal of simultaneously achieving high-yield and high-quality fruit.The knowledge behind these plots was condensed from extensive scientific literature and peer-reviewed published research on greenhouse tomato production and physiology.Mathematical expressions and plots of membership functions for defining optimality degrees of T and RH are available in Ref. [22].The sets of membership functions for defining optimality degrees of VPD are presented in the work of Ref. [23].According to this model, a membership function for specific growth stage and light condition on the universe of discourse is defined as Opt M ðÞ GS, Light ðÞ : where M : T, RH, and VPD is the universe of discourse (input).In other words, each M reading in the greenhouse at time t m, n , is mapped to a value between 0 and 1 that quantifies its optimality for tomato production.The two indexes m and n refer to specific minute and date of a measurement.In this model, an optimality degree equal to 1 refers to a potential yield with marketable value high-quality fruit.For example, Opt T ðÞ ¼ 1 is associated with T ∈ 24, 27 ½ C at the vegetative to mature fruiting growth stage during sun hours.For the same growth stage and light condition, a wider reference border, that is, T ∈ 18:4, 32:2 ½ C, is associated with a lower range of optimality degrees, Opt T ðÞ ∈ 0:6, 1 ½ .In other words, a greenhouse air temperature equal to 32:2 C during sun hours is 60% optimal for tomato production in the vegetative to mature fruiting growth stage.The reference values corresponding to the optimal, marginal, and failure T and RH are summarized in Table 1.These values for VPD depend on the range of T and RH and are discussed in Ref. [23].The optimality-degree model was implemented in the framework as a toolbox and was successfully used in evaluating microclimate parameters.Results of an actual case study on a net-screen-covered greenhouse in tropical lowlands of Malaysia are provided in Figures 6 and 7 [22].

Comfort ratio of microclimate
Comfort ratio of a microclimate parameter, denoted by Cf t M, t, α s ðÞ GS ¼ β, is defined as the percentage of M data collected during time frame t that falls inside reference borders of M Vegetative to termination (GS2-5) Flowering to termination (GS3-5) Indices are: L: lower border, H: higher border, N: night, C: cloud, S: sun, α 0 : index of failure, α 0:5 : index of Opt¼0.5, α 1 : index of Opt¼1.
Table 1.Reference values of optimal and failure T and RH at different growth stages and light conditions.
Figure 6.Plots of daily averaged air temperature, RH, and associated optimality degrees from a tropical greenhouse experiment (Source: [22]).associated with α s at a specific growth stage.A 100% ideal microclimate growth condition is therefore defined as Cf t M,t,1 ðÞ ¼ 1.The notation α s refers to user-preferred optimality degree for adjusting the reference borders that is desired for microclimate evaluation or control.The reference borders for a given α s are calculated from available simulation models (i.e., from the membership function growth response models of [21,23]).For the purpose of this chapter, mathematical descriptions of Ref. [21] model for defining reference borders of air temperature and relative humidity are adapted and provided in Table 2.An example is demonstrated in Figure 8 for constructing reference borders of air temperature associated with α s ¼ 0:8 at the vegetative to mature fruiting growth stage.The procedure is similar for other microclimate parameters (RH and VPD) at other growth stages and for any selection of 0 ≤ α s ≤ 1.The framework algorithm automatically selects proper membership functions from database according to the light condition and growth stage and computes the reference borders for the given α s .The light condition in this demonstration belongs to a random day, date: December 15, 2013.The reference borders corresponding to α s ¼ 0, α s ¼ 0:8andα s ¼ 1 are shown in red, blue, and green colors, respectively.The framework plots data inside each reference border in different colors (black for α s ¼ 0, blue for a preferred α s , and green for α s ¼ 1).If a measurement lies outside marginal reference borders (α s ¼ 0), it will be plotted in red.
The main purpose of introducing comfort ratio and corresponding graphical demonstration is to address deviation of microclimate responses with respect to different reference borders and to compare it for different cultivation days or greenhouse designs.A practical example is provided in Figure 9 for air temperature collected from a naturally ventilated greenhouse in two random days, one at the early growth and the other at the mature fruiting stage.The reference borders associated with a preferred optimality degree (i.e., α s ¼ 0:7) are shown in blue color-dashed lines.Moreover, the reference borders corresponding to failure air temperature (α s ¼ 0) and optimum air temperature (α s ¼ 1) are, respectively, shown in redand green-dashed lines.In this example, the percentage of data that falls inside these three reference borders (α s ¼ 0, 0.7 and 1) are 100, 92 and 41% for the early growth stage, and 100, 73, and 3% for the mature fruiting stage.These values are expressed on the plots of Figure 9 as Cf t T, 24, 0 ðÞ GS1 ¼ 1, Cf t T, 24, 0:7 ðÞ GS1 ¼ 0:92, Cf t T, 24, 1 ðÞ GS1 ¼ 0:41, Cf t T, 24, 0 ðÞ GS5 ¼ 1, Cf t T, 24, 0:7 ðÞ GS5 ¼ 0:73, and Cf t T, 24, 1 ðÞ GS5 ¼ 0:03.In other words, Cf t T, 24, 0:7 ðÞ GS1 ¼ 0:92 and Cf t T, 24, 0:7 ðÞ GS5 ¼ 0:73 imply that for nearly 22 h (92% of the entire 24 h) of the random day at the early growth, and for 17.5 h (73% of the entire 24 h) of the random day at the mature fruiting stage, the climate controller (for this example, natural ventilation) provided the greenhouse with air temperature that was at least 70% optimal for tomato cultivation.Moreover, Cf t T, 24, 1 ðÞ GS1 ¼ 0:41 implies that at the early growth stage, the greenhouse was controlled

Reference function
Preferred optimality with 100% optimal air temperature for a total of 9.6 h (41% of the total 24 h, shown by green color between hours of 00:00-11:00 on the left plot of Figure 9).For the random day at the mature fruiting stage, it can be seen that only 3% of the air temperature response is inside α s ¼ 1 reference borders (around hour 8:00 to 8:30 am).
The discussion for comfort ratio is extended to compare VPD response in three different greenhouses for a random data collection day during the flowering growth stage (GS3).The greenhouses had different covering materials and climate control system (labeled by A, B, and CinFigure 10, respectively, covered with net-screen mesh, polyethylene film, and polycarbonate panels).The preferred reference border for this evaluation is α s ¼ 0:6 (blue-color borders).It can be observed that VPD response never crossed α ¼ 0 or the failure reference borders in greenhouses A and C.This can be expressed by saying that Cf t VPD, 24, 0 ðÞ GS3 was never less than 1 in greenhouses A and C. It should be mentioned that these two greenhouses were, respectively, operating on natural ventilation and evaporative cooling system during the experiment.According to the plots of the three greenhouses in Figure 10, no significant difference can be observed in their VPD responses between 0.1 and 1.2 kPa (corresponding to air temperature between 20 and 30 C, and RH between 80 and 100%); however, as air temperature starts rising above 30 C, differences in the environments start growing nonlinearly.The hourly averaged values of microclimate parameters for this experiment reveal that the major differences between these greenhouses occur between hours of 11:30 am to 4:00 pm.The mean VPD value for greenhouses B and C was equal to 2.9 and 1.19 kPa, respectively, which are less desirable for plant growth compared with the 0.97-kPa value observed from greenhouse A.  This observation indicates that as long as the outside temperature is less than 30 C, no major differences between the three greenhouses resulted.This example indicates that for this particular day of experiment, the net-screen-covered greenhouse operating on natural ventilation had a comfort ratio equal to 1 at α s ¼ 0:6, which is slightly higher than Cf t T, 24, 0:6 ðÞ GS5 ¼ 0:95 of the polycarbonate panel greenhouse with evaporative cooling system.It should be noted that greenhouse C was constructed with more expensive materials, including polycarbonate panels to reduce direct sun radiation, and was operating on evaporative cooling system with large fans that consume substantial amount of electricity.This example clearly shows the potential of natural ventilation in providing more desirable response for tomato cultivation under tropical climate conditions.

Simulation of expected yield
A peer-reviewed published state-variable tomato growth model, developed by Ref. [24] in Microsoft Excel spreadsheets, was studied and implemented in MATLAB Simulink (shown by Simulink blocks in Figures 11-13).The objective was to provide a standalone application in a way that end users unfamiliar with programming languages and/or crop modeling would have an easier access to yield prediction in different greenhouse environments.Data from spreadsheet version of the model were used for testing the Simulink blocks and validation of the results [25].The five state variables included in the tomato growth model of Ref. [24] were node number (N), leaf area index (LAI), total plant weight (W)orbiomass, total fruit weight (W F ), and mature fruit weight (W M ).Vegetative node development is calculated on an hourly time step using greenhouse temperature (T).The state-variable equation for the rate of node development (dN=dt)isexpressedbydN=dt ¼ N m :f N T ðÞ ,where N m is the maximum rate of node appearance per day and f N T ðÞ , is a function to reduce node development under nonoptimal temperatures on an hourly basis.Based on studies of tomato phenology, N m was established to be 0.02083 nodes:d À1 in the model, and the function, f N T ðÞ ,i sf N T ðÞ ¼ min 1, min 0:25 þ 0:025T, 2:5 À 0:05T ðÞ ðÞ ,w h e r eT is the hourly greenhouse temperature in C. Gross hourly photosynthesis (P h ) was calculated as a function of hourly temperature, incoming solar radiation, and LAI using Eq. ( 1) developed by Ref. [26].The Simulink blocks for hourly node development and hourly photosynthesis are shown in Figure 11.Here, D is a coefficient to convert P h from μmol CO 2 ðÞ m À2 : s À1 to gC H 2 O ðÞ m À2 : d À1 , K is the light extinction coefficient, m is the leaf light transmission coefficient, LF max is the maximum leaf photosynthetic rate, Q e T ðÞ is the leaf quantum efficiency and a function of temperature, PPFD is the photosynthetic photon flux density or the level of incoming solar radiation, and PGRED T ðÞ is a function to modify P h under suboptimal temperatures.Based on previous work with tomato growth models [24], D, K, m,andLF max were set to 0.108, 0.58, 0.1, and 26, respectively.The function for Q e T ðÞ can be expressed by Q e T ðÞ ¼ 0:084 : 1 À 0:143 exp 0:0295: T À 23 ðÞ ðÞ ÀÁ .
The function for PGRED T ðÞwas disregarded for this model because environmental conditions inside a greenhouse will not fluctuate significantly enough such that this function would have an effect on tomato growth simulations.Temperature and incoming solar radiation information necessary for computation of P h were obtained from hourly measured data in the greenhouses under study, and LAI was obtained using a feedback loop in the model.Gross daily photosynthesis (P g ) was found by integrating over the 24-hourly photosynthesis calculations during each day.Hourly maintenance respiration (R h ) was computed as ,w h e r er m and Q 10 are maintenance respiration coefficients for tomato with values of 0.019 and 1.4, respectively.Daily maintenance respiration (R m ) was computed by integrating over the 24-hourly respiration calculations during the day.Vegetative node development was the only state variable computed on an hourly time step.The remaining state variables were calculated on a daily time step.The state-variable equation for computing LAI was derived from the work of [27,28].This state-variable equation is expressed by Eq. ( 2), where ρ is the plant density, λ T d ðÞ is a function to reduce the rate of leaf area expansion for nonoptimal temperatures, and δ, β,a n dN b are coefficients in the expolinear growth equation developed by Ref. [27].For this work, the values for ρ, δ, β,andN b were 3.12 plants m À2 ,0.038m À2 node À1 ,0.169node À1 ,and 16 nodes, respectively.The function, λ T d ðÞ , was not necessary for this model because temperatures within a greenhouse will not fluctuate enough for this function to significantly affect leaf area expansion simulations.The value for N is the node count at the end of the previous day, and dN=dt is the change in node count during the current day.The model assumes that when LAI reaches LAI max , any additional leaf growth will be either pruned or senesced to maintain LAI at a constant value for the remainder of the growing period.For this work, the value of LAI max was set to 4 as recommended by Ref. [24].The state-variable equation for computing the accumulation of aboveground biomass (W) is based on the equation for daily plant growth (GR net ), that is, . Here, (W À W M ) is the difference between the total aboveground biomass and the total mature fruit, and this difference represents the growing and respiring plant mass.This difference is multiplied by the daily respiration rate (R m )togetthe amount of carbon necessary for plant maintenance.Subtracting this value from the total carbon assimilated during the day (P g ) gives the total carbon available for plant growth.The coefficient, E, represents the efficiency at converting photosynthate to crop biomass, and this value was set to 0.75 in this work.The function, f R N ðÞ , determines the proportion of carbon that is partitioned to roots as a function of the number of nodes, and it can be expressed as f R N ðÞ ¼ max 0:02, 0:18 À 0:0032:N ðÞ .The function allows a relatively large portion of carbon to be allocated to roots when the plant is young, and this portion tapers off to 0.02 as the plant matures.The state-variable equation for computing the accumulation of aboveground biomass (W)i sdW=dt ¼ GR net À p 1 :ρ:dN=dt, where p 1 is the dry matter weight of leaves removed per day due either to senescence or to pruning after LAI max is achieved.For this work, the value of p 1 was 0 g:node À1 before LAI max was reached and 2 g:node À1 after LAI max was reached.The state-variable equation to calculate the total fruit weight (W F ) is expressed by Eq. ( 3).Simulink blocks for daily biomass accumulation and senescence are shown in Figure 12.
Here, α F is the maximum partitioning of new growth to fruit, f F T d ðÞ is a function to modify partitioning to fruit according to the average daily temperature (T d ), ϑ is the transition coefficient between vegetative and full fruit growth, N FF is the nodes per plant when the first fruit appears, and gT daytime ÀÁ is a function to reduce fruit growth due to high daytime temperature.
For this work, α F, ϑ, and N FF were 0.95 d À1 , 0.2 node À1 , and 10 nodes, respectively.The function f F T d ðÞ is expressed as f F T d ðÞ ¼ max 0, min 1, 0:0625: T d À T min ðÞ ðÞ ðÞ , where T min is the minimum temperature below which no fruit growth occurs.The function gT daytime ÀÁ is expressed by Eq. ( 4) where T daytime is the average temperature during daylight hours and T crit is the temperature above which fruit abortion begins.For tomato, T and T crit are 8.5 and 24.4 C, respectively.The state-variable equation to calculate the total weight of mature fruit or the total tomato yield is expressed by Eq. ( 5) where D F T d ðÞ is a function for the rate of fruit development according to the average daily temperature, and κ F is the development time from first fruit to first ripe fruit.For this work, κ F was five nodes, and the function, D F T d ðÞ , is expressed as D F T d ðÞ ¼ 0:04 :max 0, min 1, 0:0714: T d À 9 ðÞ ðÞ ðÞ .Mature fruit is assumed to be harvested from the plants immediately upon ripening, as shown by the subtraction of W M during each time step from net crop growth explained by GR net equation.Simulink blocks for daily mature fruit weight and daily fruit growth are shown in Figure 13.This description completely explicates the reduced state-variable tomato model implemented in Simulink for this project, and the state-variable equations for LAI, total biomass accumulation (dW=dt), total fruit weight (dW F =dt), and mature fruit weight ((dW M =dt)) are highlighted.The implemented model was validated [25] using the Lake City experiment datasets of Ref. [24] to show that the Simulink version of the model is an exact replication of the original spreadsheet version.It was then used in yield prediction from the three greenhouses shown in Figure 10.Results of the prediction are summarized in Figure 14, showing that the net-screen greenhouse operating on natural ventilation (greenhouse labeled A) had the highest yield compared with the polycarbonate panel and polyethylene film greenhouses.This result is completely consistent with results of the optimality degrees and comfort ratios obtained in the previous sections.

Set-point manipulation for optimum climate control 4.1. Critical reference borders
The comfort ratio curve, denoted by Cf t-curve, refers to the plot of Cf tðM,t,α s ) values calculated for all α s ¼ 0 : d α : 1.It shows how much close a microclimate parameter can be controlled to different preferred reference borders.The horizontal blue-dashed line at Cf t M,t,α s ðÞ ¼ 1 represents 100% satisfied control objective; that is, parameter M is always inside reference borders of α s .The Cf t-curve can be used as a tool to demonstrate the behavior of Cf tðM,t,α s ) in different greenhouses or at different cultivation days for decision making in set-point manipulation for the climate controller.For example, it can be used in finding the largest α s for which Cf t M, t, α s ðÞ ¼ 1 (in other words, finding α max corresponding to the narrowest achievable reference border by the climate controller).An example is provided in Figure 15 by plotting air temperature response for 2 consecutive days of an experiment inside a tropical greenhouse.It can be observed that the narrowest reference borders of air temperature that was completely satisfied by the climate controller in these two days are, respectively, equal to α s ¼ 0:55 and α s ¼ 0:67.After these points, comfort ratio starts decreasing until it arrives at its lowest value of 0.42 for both days at α s ¼ 1.
Another application of the Cf t-curve includes finding critical reference borders, denoted by is maximum (reference borders that cause significant loss in comfort ratio).To further explain, comfort ratios of air temperature for two distinct cases are plotted in Figure 16.In the first case, increasing α s from 0.3 to 0.65 has not caused significant loss in the resulting comfort ratio.The values of Cf t T, t, 0:3 ðÞ and Cf t T, t, 0:75 ðÞ for this case are nearly the same and equal to 0.8 and 0.77.In other words, by increasing α s from 0.3 to 0.75 to provide air temperature response that is more favored by tomato plants, performance of the controller in achieving the extra accuracy was not decreased.In a greenhouse with natural ventilation, this means that the extra 0.35 increase in α s comes at no cost (no significant loss of response).In the case of an energy-consuming climate controller (i.e., pad-and-fan-evaporative system or swamp cooler), it means that the cooler can be set to maintain air temperature inside a narrower reference border (by selecting α s ¼ 0:75 rather than 0.3) without imposing additional energy cost.On the right plot of Figure 16, this situation is, however, different.Significant loss in Cf t T, t, α s ðÞ can be observed for a slight increase from α s ¼ 0:7t oα s þ E ¼ 0:75.Here, increasing α s for as little as 0.05 has led to a sudden drop in the comfort ratio by 50% (from 1 to 0.5).The α s at which the largest loss appear is referred to α Crit and can be calculated by differentiating Cf t -curve with respect to α as α crit ¼ d=d α Cf t M,t,α ðÞ ðÞ .

Performance of climate controller
Plots of measured optimality degrees of a response parameter, denoted by Opt M ðÞ ¼ α y , corresponding to the preferred α s reference borders can provide a useful graphical tool to monitor performance of the climate control system.For the sake of demonstration, Cf t-curves and performance curves of the climate controller for T, RH, and VPD are shown in Figure 17.For a perfectly control task with a preferred α s , the control system must achieve microclimate parameter M that has optimality degree of at least α s .For example, if reference borders of air temperature control are set at α s ¼ 0:8, it is expected that the optimality degree of air temperature response inside the greenhouse is at least α y ¼ 0:8 at any measured time.As mentioned earlier, in a 100% perfectly controlled greenhouses, the measured optimality degrees are at least equal to the preferred optimality degrees of the reference border (α y ¼ α s ).This is shown by the perfect control line (line of α y ¼ α x ) on the response plot of Figure 17.It should be noted that α y can also be calculated by integrating Cf t M,t,α ðÞ curve over α ¼ 0t oα ¼ α s (Eq.( 6)), indicating that α y is equal to α s only when Cf t M,t,α s ðÞ ¼ 1.In other words, performance of a climate control system in achieving preferred reference borders of M is considered 100% perfect only when 100% of M-response falls inside the α s preferred optimal reference borders.
In controlled greenhouses, both Cft curve and performance curve provide a graphical assessment tool for comparing different control strategies and scenarios (i.e., microclimate responses due to different greenhouse designs, cooling systems, and covering materials at different Figure 16.Comparison between comfort versus α s 2 days of experiment in a greenhouse with evaporative cooling system for demonstration of α Crit .Left: significant increase in α s from 0.3 to 0.75 resulting in significant loss in Cft, right: slight increase in α s from 0.7 to 0.75 causing significant loss in Cft. growth stages).The performance curve in fact reveals how much a greenhouse microclimate parameter deviates from a perfectly controlled response.Deviation of the greenhouse from this ideal line at any α s can be used as an index factor of the perfect climate control task.The lesser deviation means the more perfect control task.Adaptability factor of the controller for microclimate parameter M at a preferred α s , denoted by ADPðM, α s Þ, is then defined as the ability of the controller to adapt itself with different preferred references and is calculated using Eq. ( 7).

Optimum reference borders
The optimum preferred reference border for parameter M, denoted by α Opt , is defined as the largest possible α s value for which the largest Cf t M,t,α ðÞ can be achieved.In other words, it is the value of an unknown α i for which Cf t M,t,α i ðÞ ¼ β i has the minimum distance to Cf t M,t,1 ¼ 1.In that sense, the cost function for this optimization problem is defined as , which is the Euclidean distant between the unknown point (α i and β i ) on the Cft curve and the point of ideal microclimate (α ¼ 1 and β ¼ 1).The objective is therefore to minimize this cost function by finding 0 ≤ α i ≤ 1 value that leads to the shortest Euclidean distant (minimum Dist i ) to the Cf t M,t,1 ðÞ ¼ 1.An example is demonstrated in Figure 18 for VPD response in a random day of experiment with α Opt ¼ 0:77.The plot on the right side of Figure 18 shows the values of D i versus 0 ≤ α i ≤ 1, and the of α Opt is shown as the global minimum point.

Conclusion
An adaptive management was designed, developed, and introduced in this research to respond to the needs for an iterative processing tool that acknowledge complexity and uncertainty in microclimate control and management.A systematic approach was presented for automatic data collection and processing with the objective to produce knowledge-based information in achieving optimum microclimate for producing high-quality and high-yield tomato.Applications of computer models were demonstrated through case-study examples for measuring and adjusting optimality degrees, comfort ratios, and prediction of the expected yield.Several applications of the framework toolboxes were demonstrated through case-study examples for evaluating and comparing microclimate parameters as well as yield prediction in different greenhouse environments.Specific applications of the optimization toolbox of the framework were discussed for evaluating and adjusting greenhouse climate controller through manipulated set points.It was shown that using adaptive greenhouse model for tropical climate condition, efficient use of natural ventilation, or shading will cause up to 70% savings on other energy-consuming cooling systems without sacrificing fruit quality or yield.The presented approach can be used in cost-benefit analysis for providing best management decisions such as site selection, optimum growing season, scheduling efficiencies, energy management with different climate control systems, and risk assessments associated with each task.Results of microclimate evaluation and yield prediction that are generated by this framework can be used in other crop models that estimate plant responses to the environment, or contribute to task-planning algorithms for hierarchical decomposition of climate management, and in economic models of tomato for energy conservation and energy efficient greenhouse crop productions.The framework can also be used as a research tool in future studies such evaluating effects of different greenhouse designs and shapes on comfort ratios of microclimate parameters, or finding optimum combination of ventilation and evaporative cooling systems for best fruit quality and yield.

Technical data
The custom-designed data acquisition and control system [17] for monitoring and manipulating of the microclimate parameters was built using Arduino Uno prototype board utilizing ATmega328P (Atmel®, San Jose, CA) microcontroller on the open source Arduino Uno prototyping platform programmable in Arduino sketch environment software with C (C Compiler, Brookfield, WI), a liquid crystal display, power supply, and serial port RS-232 communication cable (bidirectional with a maximum baud speed up to 115,200 bites per seconds) for transferring and storing collected data in PC.All vital components (i.e., clock generator, 2 KB of RAM, 32 KB of flash memory for storing programs and 1 KB of EEPROM for storing parameters, a 16-MHz crystal oscillator, digital input/output pins, USB connection, power regulator, power jack, and a reset button) for operating the microcontroller, as well as direct programming and access to input/output pins, available on the prototype board.Four arrays of HSM-20G-combined sensors modules (Shenzhen Mingjiada Electronics LTD, Futian Shenzhen, China), external micro-secure digital (SD) cardboard for storing larger amount of sensor data, output connection, sensor input, and relay circuit board for on/off control purposes were used.The data acquisition interface was tested for accuracy and reliability with available commercial models, and with a control sample data collected from airport weather station at Sultan Abdul Aziz Shah-Subang in Malaysia.

Figure 1 .
Figure 1.Outside and inside view of tropical greenhouses in the lowlands of Malaysia.

Figure 2 .
Figure2.Diagram of the adaptive management framework for monitoring, data processing, evaluating, and adjusting greenhouse microclimate with an iterative analysis approach for scenario analysis with greenhouse crop production.

Figure 3 .
Figure 3. Arrangement of the process in the framework toolboxes.

Figure 4 .Figure 5 .
Figure 5. Adaptive control of greenhouse air temperature based on manipulated set point as discussed in Ref.[19].

Figure 4 .
Figure 4. Demonstration of conventional greenhouse controller (left) versus model reference adaptive controller (right).

Figure 7 .
Figure 7. Demonstration of real-time measured air temperature and RH (left) and corresponding optimality degrees (right) for a random cultivation day at the flowering to mature fruiting growth stage (date: March 11, 2015) in a tropical greenhouse.Each color represents a light condition, back: night, red: sun, blue: cloud (Source: [22]).

Figure 8 .
Figure 8. Demonstration of adjusting reference borders with light condition and a preferred optimality degree of α ¼ 0:8 for air temperature control and evaluation in a random day at the flowering to mature fruiting growth stage.

Figure 9 .
Figure 9. Demonstration of air temperature response and corresponding comfort ratios for two random days of experiment at the early growth (left) and mature fruiting stage (right) in a tropical greenhouse.

Figure 11 .
Figure 11.Simulink blocks for hourly node development and hourly photosynthesis.

Figure 12 .
Figure 12.Simulink blocks for daily biomass accumulation and senescence.

Figure 13 .
Figure 13.Simulink blocks for daily mature fruit weight and daily fruit growth.

Figure 14 .
Figure 14.Simulated results with TOMGRO model for three experimental greenhouses.

Figure 15 .
Figure 15.Comparison between air temperature responses from a tropical greenhouse in 2 days of experiment showing raw data (left), and comfort ratios (right).The controller did not satisfy 100% optimal references.

Figure 17 .Figure
Figure 17.Comfort ratio of microclimate parameters (left) and response of the climate controller (right) at 0 ≤ α s ≤ 1.

Table 2 .
Membership function model for adjusting reference borders of air temperature and RH.