Open access peer-reviewed chapter

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

By Redmond R. Shamshiri, Muhammad Razif Mahadi, Kelly R. Thorp, Wan Ishak Wan Ismail, Desa Ahmad and Hasfalina Che Man

Submitted: December 6th 2016Reviewed: May 31st 2017Published: November 17th 2017

DOI: 10.5772/intechopen.69972

Downloaded: 1782


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.


  • greenhouse
  • climate control
  • microclimate evaluation
  • tomato
  • ventilation
  • evaporative cooling

1. 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.

Figure 1.

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

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 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 [13]. 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, 46]. 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. [1215]. 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 Cft(M,t,α)=β, to refer to the terms “optimality degree” and “comfort ratio,” respectively, defined in Sections 3.1 and 3.2.


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.

Figure 2.

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.

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 Cft(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 Figure 4. In this scheme, the control law is adapted with the new greenhouse states based on the optimized set points as shown in the diagram of Figure 5 [19] for a specific microclimate parameter (i.e., air temperature), causing other microclimate parameter (i.e., humidity) to be actually controlled via set‐point manipulation.

Figure 3.

Arrangement of the process in the framework toolboxes.

Figure 4.

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

Figure 5.

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


3. Microclimate evaluation with manipulated set points

3.1. 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, RHfor 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):M[0,1], where M: T,RH, and VPDis the universe of discourse (input). In other words, each Mreading in the greenhouse at time tm,n, is mapped to a value between 0 and 1 that quantifies its optimality for tomato production. The two indexes mand nrefer 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)=1is associated with T[24,27]°Cat 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°Cduring 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].

Figure 6.

Plots of daily averaged air temperature, RH, and associated optimality degrees from a tropical greenhouse experiment (Source: [22]).

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]).

TemperatureRelative humidity
Growth stageReference borderValue (°C)Growth stageReference borderValue (%)
Early growth (GS1)T1α0L9Early growth (GS1)RH1α0L60
Vegetative to termination(GS2‐5)T2α0L10Vegetative (GS2)RH2α0L40
T2α1H,N20Flowering to termination (GS3‐5)

Table 1.

Reference values of optimal and failure T and RH at different growth stages and light conditions.

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.

3.2. Comfort ratio of microclimate

Comfort ratio of a microclimate parameter, denoted by Cft(M,t,αs)GS=β, is defined as the percentage of Mdata collected during time frame tthat falls inside reference borders of Massociated with αsat a specific growth stage. A 100% ideal microclimate growth condition is therefore defined as Cft(M,t,1)=1. The notation αsrefers to user‐preferred optimality degree for adjusting the reference borders that is desired for microclimate evaluation or control. The reference borders for a given αsare 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.8at 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αs1. 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=1are 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.

Figure 8.

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

Reference functionPreferred optimality
T(α)G1A={T1α0L  T1α0Hα=0
α(T1α1LT1α0L)+T1α0L  α(T1α1HT1α0H)+T1α0H0<α<1
T(α)G2S={T2α0L  T2α0Hα=0
α(T2α1L,ST2α0L)+T2α0L  α(T2α1H,ST2α0H)+T2α0H0<α<1
T(α)G2C={T2α0L  T2α0H α=0
α(T2α1L,CT2α0L)+T2α0L  α(T2α1H,CT2α0H)+T2α0H0<α<1
T(α)G2N={T2α0L  T2α0H α=0
RH(α)G2A={RH2α0L  RH2α0Hα=0
α(RH2α1LRH2α0L)+RH2α0L  α(RH2α1HRH2α0H)+RH2α0H0<α<1
RH(α)G3A={RH3α0L  RH3α0Hα=0
α(RH3α1LRH3α0L)+RH3α0L  α(RH3α1HRH3G0,max)+RH3α0H0<α<1

Table 2.

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

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 red‐ and 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 Cft(T,24,0)GS1=1, Cft(T,24,0.7)GS1=0.92, Cft(T,24,1)GS1=0.41, Cft(T,24,0)GS5=1, Cft(T,24,0.7)GS5=0.73, and Cft(T,24,1)GS5=0.03. In other words, Cft(T,24,0.7)GS1=0.92and Cft(T,24,0.7)GS5=0.73imply 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, Cft(T,24,1)GS1=0.41implies that at the early growth stage, the greenhouse was controlled 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=1reference borders (around hour 8:00 to 8:30 am).

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.

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 C in Figure 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 α=0or the failure reference borders in greenhouses A and C. This can be expressed by saying that Cft(VPD,24,0)GS3was 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 Cft(T,24,0.6)GS5=0.95of 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.

Figure 10.

A comparison between comfort ratio of VPD at reference borders ofα=0,α=0.6,andα=1in three different greenhouses. Date of data collection March 18, 2013.

3.3. 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 1113). 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) or biomass, total fruit weight (WF), and mature fruit weight (WM). 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) is expressed by dN/dt=Nm.fN(T), where Nmis the maximum rate of node appearance per day and fN(T), is a function to reduce node development under nonoptimal temperatures on an hourly basis. Based on studies of tomato phenology, Nmwas established to be 0.02083 nodes.d1in the model, and the function, fN(T), is fN(T)=min(1,min(0.25+0.025T,2.50.05T)), where Tis the hourly greenhouse temperature in °C. Gross hourly photosynthesis (Ph) 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 Phfrom μmol(CO2)m2. s1to g(CH2O)m2. d1, Kis the light extinction coefficient, mis the leaf light transmission coefficient, LFmaxis the maximum leaf photosynthetic rate, Qe(T)is the leaf quantum efficiency and a function of temperature, PPFDis the photosynthetic photon flux density or the level of incoming solar radiation, and PGRED(T)is a function to modify Phunder suboptimal temperatures. Based on previous work with tomato growth models [24], D, K, m, and LFmaxwere set to 0.108, 0.58, 0.1, and 26, respectively. The function for Qe(T)can be expressed by Qe(T)=0.084 . (10.143exp(0.0295.(T23))).

Figure 11.

Simulink blocks for hourly node development and hourly photosynthesis.

Figure 12.

Simulink blocks for daily biomass accumulation and senescence.

Figure 13.

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

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 Phwere obtained from hourly measured data in the greenhouses under study, and LAI was obtained using a feedback loop in the model. Gross daily photosynthesis (Pg) was found by integrating over the 24‐hourly photosynthesis calculations during each day. Hourly maintenance respiration (Rh) was computed as Rh=rm.Q10(T20)/10, where rmand Q10are maintenance respiration coefficients for tomato with values of 0.019 and 1.4, respectively. Daily maintenance respiration (Rm) 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, λ(Td)is a function to reduce the rate of leaf area expansion for nonoptimal temperatures, and δ, β, and Nbare coefficients in the expolinear growth equation developed by Ref. [27]. For this work, the values for ρ, δ, β, and Nbwere 3.12 plants m2, 0.038 m2node1, 0.169 node1, and 16 nodes, respectively. The function, λ(Td), 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 Nis the node count at the end of the previous day, and dN/dtis the change in node count during the current day. The model assumes that when LAI reaches LAImax, 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 LAImaxwas 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 (GRnet), that is, GRnet=E.[PgRm(WWM)]. [1fR(N)]. Here, (WWM) 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 (Rm) to get the amount of carbon necessary for plant maintenance. Subtracting this value from the total carbon assimilated during the day (Pg) 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, fR(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 fR(N)=max(0.02,0.180.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) is dW/dt=GRnetp1.ρ.dN/dt, where p1is the dry matter weight of leaves removed per day due either to senescence or to pruning after LAImax is achieved. For this work, the value of p1was 0 g.node1before LAImaxwas reached and 2 g.node1after LAImaxwas reached. The state‐variable equation to calculate the total fruit weight (WF) is expressed by Eq. (3). Simulink blocks for daily biomass accumulation and senescence are shown in Figure 12.

Here, αFis the maximum partitioning of new growth to fruit, fF(Td)is a function to modify partitioning to fruit according to the average daily temperature (Td), ϑis the transition coefficient between vegetative and full fruit growth, NFFis the nodes per plant when the first fruit appears, and g(Tdaytime)is a function to reduce fruit growth due to high daytime temperature. For this work, αF, ϑ, and NFFwere 0.95 d1, 0.2 node1, and 10 nodes, respectively. The function fF(Td)is expressed as fF(Td)=max(0,min(1,0.0625.(TdTmin))), where Tminis the minimum temperature below which no fruit growth occurs. The function g(Tdaytime)is expressed by Eq. (4) where Tdaytimeis the average temperature during daylight hours and Tcritis the temperature above which fruit abortion begins. For tomato, Tminand Tcritare 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 DF(Td)is a function for the rate of fruit development according to the average daily temperature, and κFis the development time from first fruit to first ripe fruit. For this work, κFwas five nodes, and the function, DF(Td), is expressed as DF(Td)=0.04 .max(0,min(1, 0.0714.(Td9))). Mature fruit is assumed to be harvested from the plants immediately upon ripening, as shown by the subtraction of WMduring each time step from net crop growth explained by GRnetequation. 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 (dWF/dt), and mature fruit weight ((dWM/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.

Figure 14.

Simulated results with TOMGRO model for three experimental greenhouses.


4. Set‐point manipulation for optimum climate control

4.1. Critical reference borders

The comfort ratio curve, denoted by Cft‐ curve, refers to the plot of Cft(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 Cft(M,t,αs)=1represents 100% satisfied control objective; that is, parameter Mis always inside reference borders of αs. The Cft‐ curve can be used as a tool to demonstrate the behavior of Cft(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 αsfor which Cft(M,t,αs)=1(in other words, finding αmaxcorresponding 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.55and α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.

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.

Another application of the Cft‐ curve includes finding critical reference borders, denoted by αCritat which Δ=Cft(M,t,αs)Cft(M,t,αs+ϵ)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 αsfrom 0.3 to 0.65 has not caused significant loss in the resulting comfort ratio. The values of Cft(T,t,0.3)and Cft(T,t,0.75)for this case are nearly the same and equal to 0.8 and 0.77. In other words, by increasing αsfrom 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 αscomes 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.75rather than 0.3) without imposing additional energy cost. On the right plot of Figure 16, this situation is, however, different. Significant loss in Cft(T,t,αs)can be observed for a slight increase from αs=0.7to αs+ϵ=0.75. Here, increasing αsfor as little as 0.05 has led to a sudden drop in the comfort ratio by 50% (from 1 to 0.5). The αsat which the largest loss appear is referred to αCritand can be calculated by differentiating Cft ‐curve with respect to αas αcrit=d/dα(Cft(M,t,α)).

Figure 16.

Comparison between comfort ratios versusαsin 2 days of experiment in a greenhouse with evaporative cooling system for demonstration ofαCrit. Left: significant increase inαsfrom 0.3 to 0.75 resulting in significant loss in Cft, right: slight increase inαsfrom 0.7 to 0.75 causing significant loss in Cft.

4.2. Performance of climate controller

Plots of measured optimality degrees of a response parameter, denoted by Opt(M)=αy, corresponding to the preferred αsreference borders can provide a useful graphical tool to monitor performance of the climate control system. For the sake of demonstration, Cft‐ 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 Mthat 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.8at 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 αycan also be calculated by integrating Cft(M,t,α)curve over α=0to α=αs(Eq. (6)), indicating that αyis equal to αsonly when Cft(M,t,αs)=1. In other words, performance of a climate control system in achieving preferred reference borders of Mis considered 100% perfect only when 100% of M‐response falls inside the αspreferred optimal reference borders.

Figure 17.

Comfort ratio of microclimate parameters (left) and response of the climate controller (right) at0αs1.

In controlled greenhouses, both Cft curveand 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 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 αscan 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 Mat 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).

4.3. Optimum reference borders

The optimum preferred reference border for parameter M, denoted by αOpt,is defined as the largest possible αsvalue for which the largest Cft(M,t,α)can be achieved. In other words, it is the value of an unknown αifor which Cft(M,t,αi)=βihas the minimum distance to Cft(M,t,1)=1. In that sense, the cost function for this optimization problem is defined as Di=(αi1)2+(βi1)2, which is the Euclidean distant between the unknown point (αiand βi) on the Cft curveand the point of ideal microclimate (α=1and β=1). The objective is therefore to minimize this cost function by finding 0αi1value that leads to the shortest Euclidean distant (minimum Disti) to the Cft(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 Diversus 0 ≤ αi≤ 1, and the position of αOptis shown as the global minimum point.

Figure 18.

Demonstration of the algorithm for finding optimum preferred reference border for adjusting the climate controller. Data belongs to VPD response from a random data collection day in a tropical greenhouse experiment.


5. Conclusion

An adaptive management framework 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 as 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.


6. 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, were 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.

Ph=D . LFmax . PGRED(T)K.ln[ (1m).LFmax+Qe(T) . K. PPFD(1m) . LFmax+Qe(T). K. PPFD . exp(k . LAI)]E1
{d(LAI)dt=ρ.δ.λ(Td).exp[β.(NNb)] 1+exp[β.(NNb)].dNdt: LAILAImaxd(LAI)dt=0: LAILAImaxE2
dWFdt=GRnet.αF.fF(Td).[1exp(ϑ(NNFF))].g(Tdaytime) if N>NFFE3
dWMdt=DF(Td).(WFWM), if N>NFF+κFE5
ADP(M,αs) =12(α=0α=αsα.dαα=0α=αsOpt(M).dα)E7


The financial support provided by the University of Putra Malaysia, Grant Number GP‐IPB/2013/9415600, and the scientific comments and suggestions from Professor Warren Dixon, Professor Jim Jones, and Professor Ray Bucklin at the University of Florida, and Professor Jan Bontsema at the Wageningen University and Research Center are duly acknowledged.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Redmond R. Shamshiri, Muhammad Razif Mahadi, Kelly R. Thorp, Wan Ishak Wan Ismail, Desa Ahmad and Hasfalina Che Man (November 17th 2017). Adaptive Management Framework for Evaluating and Adjusting Microclimate Parameters in Tropical Greenhouse Crop Production Systems, Plant Engineering, Snježana Jurić, IntechOpen, DOI: 10.5772/intechopen.69972. Available from:

chapter statistics

1782total chapter downloads

7Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Introductory Chapter: Plant Engineering for the Future Ahead

By Snježana Jurić

Related Book

First chapter

Introductory Chapter: Technological Innovation as Tool for Products Qualification

By Innocenzo Muzzalupo

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us