Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A framework of artificial light management for optimal plant development for smart greenhouse application

  • João Pereira,

    Roles Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation School of Mechanical, Aerospace and Automotive Engineering, Coventry University, Coventry, United Kingdom

  • Abdul Mounem Mouazen,

    Roles Formal analysis, Writing – original draft, Writing – review & editing

    Affiliation Department of Environment, Ghent University, Ghent, Belgium

  • Mathias Foo,

    Roles Conceptualization, Formal analysis, Funding acquisition, Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliation School of Engineering, University of Warwick, Coventry, United Kingdom

  • Hafiz Ahmed

    Roles Conceptualization, Formal analysis, Writing – original draft, Writing – review & editing

    hafiz.ahmed@bangor.ac.uk

    Affiliation Nuclear Futures Institute, Bangor University, Bangor, United Kingdom

Abstract

Smart greenhouse farming has emerged as one of the solutions to global food security, where farming productivity can be managed and improved in an automated manner. While it is known that plant development is highly dependent on the quantity and quality of light exposure, the specific impact of the different light properties is yet to be fully understood. In this study, using the model plant Arabidopsis, we systematically investigate how six different light properties (i.e., photoperiod, light offset, intensity, phase of dawn, duration of twilight and period) would affect plant development i.e., flowering time and hypocotyl (seedling stem) elongation using an established mathematical model of the plant circadian system relating light input to flowering time and hypocotyl elongation outputs for smart greenhouse application. We vary each of the light properties individually and then collectively to understand their effect on plant development. Our analyses show in comparison to the nominal value, the photoperiod of 18 hours, period of 24 hours, no light offset, phase of dawn of 0 hour, duration of twilight of 0.05 hour and a reduced light intensity of 1% are able to improve by at least 30% in days to flower (from 32.52 days to 20.61 days) and hypocotyl length (from 1.90 mm to 1.19mm) with the added benefit of reducing energy consumption by at least 15% (from 4.27 MWh/year to 3.62 MWh/year). These findings could provide beneficial solutions to the smart greenhouse farming industries in terms of achieving enhanced productivity while consuming less energy.

Introduction

The increase in the global population together with the surge of demand in the food industry has created concerns about food security, prompting a need to explore sustainable agriculture practices to solve this problem [1]. One of the emerging solutions to address this problem is the concept of precision and smart greenhouse farming (see e.g. [24]), where the productivity of farming is managed and improved by using technologies involving sensors, actuators and control systems (see e.g. [57]). The application of these technologies in smart greenhouse farming has enabled variables such as light, temperature, soil moisture, etc to be controlled and managed to ensure maximum farming productivity can be achieved (see [8] and references therein for more details).

Among all the variables mentioned above, it is known that plant development and growth are greatly influenced by light [9]. Hence, this makes light one of the most important variables that need to be well-managed to ensure good yield and productivity in smart greenhouse farming environment [10, 11]. The evidence of this can be seen with the extensive studies on how to determine the efficient lighting properties for optimal plant growth that is applicable to smart greenhouse farming. In a review paper by Bian et al. [12] the photoperiod, light quality (colours) and intensity have been identified to be the key light properties in regulating plant growth. The authors also discussed the progress and advantages of LED technology to facilitate plant growth in a controlled environment. Nevertheless, due to the different scope of study, the discussion on the quantitative light properties for plant growth is not included. The effect of that same three light properties on the lettuce growth was investigated by Zhang et al., [13]. The authors considered two photoperiods and four different light intensities in their experimental studies and suggested that a 16-hour light/8-hour dark photoperiod with a light intensity of 250 μmol/m2s provides the ideal condition for lettuce growth with the lowest energy consumption. In a similar study, Hiroki et al. [14] used the same 16-hour light photoperiod but varied the light period from 18 to 24 hours and showed that a shorter period has the best growth in lettuce albeit at the expense of high energy consumption. On the specific effect of photoperiod on lettuce growth, [15] and [16] both reported that the preferred photoperiod range is between 16- to 20-hour light. The speed breeding method [17] introduced by Watson et al., marks the culmination of the importance of photoperiod in improving crop productivity, where the productivity of crops such as wheat and barleys that are subjected to extended light exposure (up to 22-hour light) have been experimentally shown to accelerate by at least twofold [17]. More importantly, this speed breeding protocol can be readily extended to growth chambers and greenhouses [18], thus making it an attractive avenue for integrating optimal lighting management in smart greenhouse farming for increasing crop productivity.

While significant efforts have been made to identify lighting configuration for maximum plant productivity [19] with lowest energy consumption [20], a systematic framework to provide quantitative optimal lighting recommendations is still lacking. Most studies have been primarily focusing on the effect of light properties like photoperiod, period, and light intensity despite there being other light properties such light offset, phase of dawn, etc. that could also influence plant development. Moreover, as mentioned above, often there is a wide range of these light properties (e.g., photoperiod of 16 to 20 hours) being reported that could improve productivity across different plant species.

The lack of a systematic framework is what forms the motivation of this study. We anticipate that the reported wide range of light properties can be further narrowed down and made more precise through a systematic approach of exploring different aspects of light properties using computational models. The use of computational models can also serve as a complement to existing experimental approaches, which is often time and labour consuming. By providing a more precise quantitative range of these light properties, we could also obtain the added benefit of further energy consumption reduction associated with these lighting operations. This thereby could reduce the operational cost and carbon emission making the food production system more environment friendly (see e.g. [10, 21, 22]).

In this study, using the model plant Arabidopsis, we present a systematic framework of artificial light management across six light properties (i.e., photoperiod, light offset, intensity, phase of dawn, duration of twilight and period) for optimal plant development with the added benefit of improved energy efficiency. This is done by varying these six light properties individually and then collectively using an established mathematical model of the Arabidopsis circadian system relating light input to flowering time and hypocotyl elongation outputs. The circadian system is considered here following the nascent research focus on circadian agriculture, where chronobiology is utilised in agriculture for improved food productivity, security and sustainability (see [23] and references therein). The results from this analysis could facilitate better decision making for smart greenhouse farming practitioners with a more specific quantitative range of light properties, which at the same time is able to reduce energy consumption.

The main contributions of this study are as follows: providing for the first time a systematic analysis on the effect of six different light properties, (i.e., photoperiod, light offset, intensity, phase of dawn, duration of twilight and period) on plant development and recommending the best combination of light properties that not only ensures optimum plant development but with minimum energy consumption.

Materials and methods

Input light function

Most of the modelling of plant circadian literature uses a simple binary representation, i.e. ‘1’ for ON and ‘0’ for OFF to represent a light function. This binary representation is an adequate representation of the light function used in smart greenhouse farming as light usually turns on and off almost instantaneously. Nevertheless, the binary representation has limited light properties to be analysed. Therefore, in this study, we consider a more comprehensive input light function, which is used in Seaton et al [24] that incorporates six light properties as shown in Eq (1), (1) where Pr is the period, Ph is the photoperiod, Dw is the phase of dawn, Tw is the duration of twilight, A is the light intensity (amplitude), Δl is the light offset, tanh is the hyperbolic tangent function and ⌊.⌋ represents the floor function.

Output phenotype calculation

The two output phenotypes are represented by gene expression of ATHB2 and FT for hypocotyl elongation and flowering time regulatory pathways, respectively. To convert these two gene expressions into measurable outputs, the following equations presented in [24] are used.

The hypocotyl length measured in mm is calculated using (2) where z(t) is the gene expression of ATHB2. If ATHB2 < a3, z(t) = ATHB2 and if ATHB2a3, we have z(t) = a3. The variable a3 represents a saturation term to limit the effect of ATHB2 to further downstream genes within the gene network.

The days to flower measured in days is calculated using (3) where d0 is a parameter that fits a sigmoid function to experimental data used in [25]. The parameters a0 to a5 are estimated from data obtained across different photoperiods [26, 27]. Following [24], we use d0 = 16.55; a1 = 0.9; a2 = 0.03103; a3 = 0.8; a4 = -2308.141; a5 = 0.02. The expression FTarea is the area under the curve for the gene expression of FT over a period. Note that in [25], FTarea is approximated by a quadratic equation, whereas in our study, FTarea is calculated using the above equation.

Simulation model

The simulation model used in this study is from Seaton et al [24]. The Arabidopsis circadian mathematical model incorporating the flowering time and hypocotyl elongation pathways consist of 48 ODEs, where 30 ODEs are associated with the Arabidopsis circadian clock [28] and 18 ODEs are associated with the two phenotype pathways. For the complete ODEs, see Supplementary Information of [24].

The MATLAB scripts used in our simulation are obtained from S1 File of [24], with minor modifications on the data processing and figure plotting to cater to our analysis. The two main MATLAB scripts from S1 File of [24] that are required for the analysis are light_conditions.m and simulate_model.m, which simulate the input light function and the output gene expressions of the phenotypes, respectively.

In the original MATLAB script, only three light properties are varied. Here, we modify the MATLAB script such that we can vary the six light properties to obtain the output gene expression of ATHB2 and FT phenotypes simulated using simulate_model.m. From the two gene expressions of ATHB2 and FT, the hypocotyl length and days to flower are calculated using Eqs (2) and (3) given above using the file FT_ATHB2_simscript.m. All the simulation files can be downloaded from https://github.com/mathiasfoo/lightmanagement.

Effective light duration for flowering calculation

The “Effective Light Duration for Flowering” refers to the duration over the period when the light is turned ON across the total number of days to flower. This can be calculated as follows: (4) where Ph is the photoperiod and Pr is the period.

Energy consumption calculation

LED lights used for smart greenhouse farming are usually powered by the LED driver, which converts AC grid voltage to DC voltage. The luminous flux of the LED is controlled by the current [29] and often the light intensity is assumed to be proportional to the current. Here, we consider the LED lights used for smart greenhouse farming are powered by the 500W Mean Well Model RSP-500-48 LED driver. The LED driver provides a constant voltage of 48 V, a current that varies between 0 to 10.5 A and the AC-DC conversion efficiency is 90.5% [30]. With that, following standard energy consumption derivation (see e.g., [31]), the energy consumption for the LED driver can be calculated as follow: (5) where Vout is the rated output voltage of the LED driver, Iout is the maximum rated current of the LED driver, A is the intensity of the LED, ηLED is the energy conversion efficiency of the LED driver, Ph is the photoperiod, Pr is the period, N is the number of days and the multiplication by 24 is to convert days to hours.

Results and discussion

Mathematical model relating light to plant development

All living beings are embedded with a biological clock called the circadian clock that can synchronise the organism’s biological functions with the 24-hour day cycles. In plants, the circadian system, which governs the optimal coordination of biological timing has been identified to be responsible for most of the plant development [3234] such as flowering time, hypocotyl elongation, petal opening, roots growth etc. (see e.g. [35, 36]). The nascent circadian agriculture [23] has led plant biologists to look at crop productivity and sustainability from the circadian system perspective, thereby justifying the use of the plant circadian mathematical model in our analysis.

One of the key influential inputs driving the circadian system that affects plant development is light [37, 38] and to complement experimental studies, many plant circadian mathematical models (see e.g. [24, 28, 3942]) have been developed to provide better insights into the light-plant development mechanisms. In order for us to carry out a systematic analysis of the light management on plant development, we employ a well-established mathematical model of the plant circadian system developed by Seaton et al., [24], as this model is the only known model to date that comprehensively relates plant circadian system to two phenotypes namely, flowering time and hypocotyl elongation, which will act as a proxy for plant development. The dynamical behaviour of all the regulatory genes involved in flowering time and hypocotyl elongation is modelled using Ordinary Differential Equations (ODE) as discussed in Materials and methods section.

Fig 1 shows the summary of the regulatory genetic pathways relating the input light, circadian system and the two aforementioned phenotypes. The genetic pathways that characterise the plant physiological properties are mathematically modelled taking the input light condition and providing output gene concentrations of FT and ATHB2 that can be converted to days to flower and hypocotyl length respectively (see Materials and methods section). For more details on this model, see [24].

thumbnail
Fig 1. Regulatory pathways relating light input, circadian clock and output phenotypes.

The six light properties are varied to study their effect on the output phenotypes viz hypocotyl length and days to flower. The effect of temperature is not investigated in this study, and it is fixed at 22-degree Celsius following [24].

https://doi.org/10.1371/journal.pone.0261281.g001

We note that most of the plant circadian mathematical models are developed based on the plant model Arabidopsis thaliana grown in a controlled laboratory condition. Although not all knowledge from Arabidopsis can be transferable to crop-based plants grown in the smart greenhouse, the findings from using Arabidopsis can still provide the relevant fundamental knowledge about general plant behaviour [43], thus warranting the validity of using Arabidopsis circadian mathematical model in this study.

The input light function given by Eq (1) used in Seaton et al., is taken from [44], where the light function is governed by six properties, namely photoperiod (i.e., length of light and dark in a period) Ph, phase of dawn (i.e., start time of light) Dw, light offset Δl, light intensity (amplitude) A, duration of twilight (i.e., light/dark transitions regime) Tw and period (i.e., length of one full light cycle) Pr. The combination of these six properties results in a distribution of light as a function of time and by individually and collectively manipulating the different properties of the light condition, the role and importance of each property can be well understood.

Characterisation of nominal phenotypes

To have a reference for comparison in terms of the two phenotypes (days to flower and hypocotyl length), a nominal value for each of the two phenotypes associated with the light function that is commonly used in a typical smart greenhouse environment needs to be first calculated. With that we set the light function (Eq (1)) with the following parameters (Table 1); Ph = 12 h, Pr = 24 h, Tw = 0.05 h, Dw = 0 h, A = 1 and Δl = 0. The output gene concentrations of FT and ATHB2 are then converted into days to flower and length in mm using Eqs (2) and (3), respectively following [24]. With these light properties, the nominal value for days to flower is 32.52 days and the hypocotyl length is 1.90 mm (Table 1). By comparing these two nominal values against the one subjected to the variation of different light properties, we can then evaluate which light properties are most influential to plant development. More importantly, we want to determine which of these light properties can be manipulated to improve the two phenotypes viz reduction in days to flower and relatively short hypocotyl length.

thumbnail
Table 1. Days to flower and hypocotyl length associated with nominal and recommended light properties.

https://doi.org/10.1371/journal.pone.0261281.t001

Flowering time is greatly related to plant development that results in flowers, fruits, and seeds production [45, 46]. Therefore, an increase in flowering time is not preferred as it impedes productivity due to longer generation times for flowers, fruits, and seeds production [47]. The elongation of the hypocotyl is usually advantageous when plants need to compete for more sunlight when growing in a dense environment. Nevertheless, under the absence of competition, a long hypocotyl is disadvantageous as this induces a higher risk of having mechanical damages to the stems of the plants due to the increase of the centre of gravity [4851]. In the following section, we will vary each light property individually to determine which property can better enable the improvement to the days to flower and hypocotyl length.

Effect of varying photoperiod

A photoperiod Ph = 12 h means the input light has 12 hours of light and 12 hours of dark (see Fig 2A). The variation of photoperiod modifies the total hours of light and dark within a 24-hour period. A Ph < 12 h usually represents a short day (e.g., winter season) while a Ph > 12 h usually represents a long day (e.g., summer season). The value of Ph is varied between 0 and 24 hours with increments of 1 hour and the results are shown in Fig 2B and S1 Table of S1 Appendix.

thumbnail
Fig 2. Effect of varying photoperiod Ph.

(A) Variation of Ph on the light function given by Eq (1). (B) Top row: Simulated gene expressions of FT and ATHB2. Bottom row: Days to flower and hypocotyl length associated with the simulated gene expressions of FT and ATHB2 calculated using Eqs (2) and (3). The dotted lines represent the nominal value of the days to flower and hypocotyl length.

https://doi.org/10.1371/journal.pone.0261281.g002

The top row of Fig 2B shows the simulated gene expression levels for FT and ATHB2 while their corresponding hypocotyl length and days to flower calculated using Eqs (2) and (3) are given at the bottom row of Fig 2B. Note that these results are the identical results obtained and shown in [24].

Here, we can clearly see the dependencies of these two phenotypes with their corresponding gene concentrations, where the days to flower are inversely proportional to FT expression levels and the hypocotyl length is proportional to ATHB2 expression levels. For Ph = 0 h and 1 h, we obtain negative values for the days to flower. These negative values occur because the relative expression levels of FT are so abnormally low that the second term of Eq (3) becomes more negative than the first term, (i.e., FTarea < a5). Note that a4 and a5 in Eq (3) are empirical parameters estimated under the experimental conditions, where the photoperiod ranges from 8 to 24 hours [26]. As the considered photoperiod does not account for a low Ph value, this leads to the negative values of days to flower. Nevertheless, for the purpose of this study, this result is insignificant since the small Ph values correspond to longer days to flower (>100 days), which is undesirable for plant development.

As Ph values increase, we observe the hypocotyl length and days to flower decrease exponentially. The days to flower significantly reduced from 113 days (Ph = 8 h) to 32 days (Ph = 12 h), which corresponds to a drop of 72% with a 50% increase in the number of photoperiods, while the hypocotyl length also significantly drops by almost two-fold from 3.58 mm (Ph = 8 h) to 1.90 mm (Ph = 12 h). By exposing the plant to a longer light (Ph > 16 h) the days to flower and hypocotyl length can be further reduced by 10 days and 0.7 mm respectively. Our analysis shows that Ph = 21 h yields the shortest days to flower with 20.1 days and the hypocotyl length of 1.18 mm. While this Ph = 21 h is close to the one suggested in the speed breeding approach [17], our analysis also reveals that Ph = 18 h would yield very similar phenotype behaviours (i.e. days to flower = 20.6 days and hypocotyl length = 1.19 mm), whilst using up to three hours lesser light per 24 hours period. Utilising three hours less light but with a very similar number of days to flower is favourable as this leads to a reduction of energy consumption (S1 Table of S1 Appendix).

Effect of varying light offset

The light offset Δl shifts the light function up or down. For Δl > 0, this represents the presence of background light with minimum light intensity (Fig 3A). We vary Δl from 0 to 5% with increments of 1%. Fig 3B and S2 Table of S1 Appendix show the days to flower and hypocotyl length against the change in Δl, with the black dashed line in Fig 3B representing the nominal values. We observe that light offset brings improvement to both phenotypes. The days to flower decreases exponentially and a change in Δl from 0 to 1% can reduce the days to flower by ~6 days (i.e., from the nominal value of 32.5 days to 26.7 days). Further increase in Δl does further reduce the days to flower albeit not by much. Our variation of Δl ends at 5% because for Δl > 5%, the days to flower remain close to 23 days. For hypocotyl length, only at the onset of Δl do we see changes in the hypocotyl length from the nominal value of 1.90 mm to 1.18 mm at Δl = 1%. Further increase in Δl does not affect the hypocotyl length due to plants being consistently exposed to minimal background light thereby negating the need for the plant to elongate its hypocotyl to seek the light. While Δl = 0.05 yields the best phenotypic values, we recommend Δl = 0.04 instead as the difference in the phenotypes is minimal, at the same time, could potentially reduce energy consumption due to a lower background light intensity.

thumbnail
Fig 3. Effect of varying light offset Δl and light intensity A.

(A, C) Variation of Δl and A on the light function given by Eq (1). (B, D) Days to flower and hypocotyl length associated with the variation of Δl and A. The dotted lines represent the nominal value of the days to flower and hypocotyl length.

https://doi.org/10.1371/journal.pone.0261281.g003

Effect of varying light intensity (amplitude)

The light intensity (amplitude) A (Fig 3C) corresponds to the amount of photosynthetic photon flux density (PPFD) being received by the plant [52]. We vary A within ±10% from the nominal value and the results are shown in Fig 3D and S3 Table of S1 Appendix. Like before, the dashed line represents the nominal values for the two phenotypes. Similarly, to the effect of varying light offset, light intensity is shown to be influential to the days to flower in an exponential manner. For lower light intensities (A < 1), the days to flower are significantly larger (i.e., > 50 days) compared to the case of higher light intensities (A > 1), where the days to flower ranges between 23–25 days. For A < 1, the days to flower increase approximately 20 days for every 2% decrease in the light intensity from the nominal value. For A > 1, the days to flower initially reduced by approximately 8 days from the nominal value but then remain unchanged with the increase of light intensity, which is in agreement with some of the reported findings that there is a limit on how much light intensity can affect the circadian clock and its downstream phenotypes (see [53] and the references therein). On the other hand, the hypocotyl length does not seem to be greatly affected by light intensity as we observe the hypocotyl length ranges between 2.0–2.4 mm with A varied within ±10%, which is consistent with the findings in [54]. Based on these analyses, we recommend an increase of light intensity by 6% to provide benefits to plant development albeit requiring higher energy consumption due to the increase in light intensity (see e.g., [55]).

Effect of varying phase of dawn

The phase of dawn, Dw is the time when the light first shines. As shown in Fig 4A, varying Dw affects the start position of the light function with Dw > 0 h and Dw < 0 h, respectively, representing a later and earlier start position of the light function. Notably, varying Dw does not alter the shape of the light function but only the start position and thus should not have any influence on the two phenotypes. In Fig 4B, we plot the days to flower and hypocotyl length for different values of Dw, with the dashed line representing the nominal values with Dw = 0 h. As expected, there is no significant change in the days to flower and hypocotyl length across different values of Dw (S4 Table of S1 Appendix). The shape of FT and ATHB2 expression levels are practically identical but have their peak expressions at different times of the day (S1 Fig of S1 Appendix). While these expression levels translate into no change in days to flower and hypocotyl length, having different peak expressions could potentially provide other advantages to plants such as temporal coordination of other circadian clock events [42] and reducing cost for protein synthesis associated with waveform generation [56, 57]. Nevertheless, from the perspective of plant development, varying Dw bears no significant impact.

thumbnail
Fig 4. Effect of varying phase of dawn Dw and duration of twilight Tw.

(A, C) Variation of Dw and Tw on the light function given by Eq (1). (B, D) Days to flower and hypocotyl length associated with the variation of Dw and Tw. The dotted lines represent the nominal value of the days to flower and hypocotyl length.

https://doi.org/10.1371/journal.pone.0261281.g004

Effect of varying duration of twilight

The duration of twilight, Tw is related to the decadence of light between extremes of intensity. The increase in the value of Tw results in a slower transition from maximum light to minimum dark (Fig 4C). The variation of Tw also influences the shape of the light function, thereby influencing the total amount of light intensity input to the plant. We vary the value of Tw from 0.01 h to 3 h with a logarithmic increment for Tw < 1 h and linear increment of 0.5 h for Tw > 1 h. Fig 4D shows the two phenotypes plotted against different Tw with the dashed line representing the nominal value. Despite having a significant influence on the days to flower (S5 Table of S1 Appendix), variation of Tw does not reduce the number of days to flower. In fact, the increase in Tw increases the number of days to flower, which is not desirable from the plant development point of view. This is because, with increasing value of Tw, reduces the total average amount of light intensity required for the plant to achieve earlier flowering. As for hypocotyl length, we observe a decreasing trend with an increase in Tw prior plateauing to a constant value of 1.18 mm (S5 Table of S1 Appendix). This is because an increase in Tw also introduces background light as this causes the light presence to linger towards the dark cycle that resulting in a behaviour similar to introducing light offset. These results indicate that Tw should not be altered from the perspective of plant development.

Effect of varying period

The period Pr manipulation corresponds to the modification of the diurnal period, whereby the light cycle will repeat after a particular number of hours instead of the usual 24 h. It has been reported in the literature that plants grow best in an environment that matches the natural 24 h period with an equal light-dark cycle (12L12D), when compared to a plant growing in Pr = 20 h and Pr = 28 h with their respective equal light-dark cycle i.e., 10L10D and 14L14D respectively [55]. Here, we are interested in the effect of varying the period on plant development by considering an extended range of Pr < 20 h and Pr > 28 h. We vary Pr from 16 to 32 hours with 2 hours increment and, in each period, we set Ph = Pr/2 h, i.e., there is always an equal amount of light-dark cycle within the considered period. The illustration of the light function is shown in Fig 5A.

thumbnail
Fig 5. Effect of varying period Pr.

(A) Variation of Pr on the light function given by Eq (1). (B) Days to flower and hypocotyl length associated with the variation of Pr. The dotted lines represent the nominal value of the days to flower and hypocotyl length.

https://doi.org/10.1371/journal.pone.0261281.g005

While the increase in Pr value to 32 hours reduces the days to flower from nominal value by ~10 days (Fig 5B (left) and S6 Table of S1 Appendix), the hypocotyl length increases by ~10 mm (Fig 5B (right), and S6 Table of S1 Appendix), which is not desirable. For smaller values of (Pr < 24 h), when compared to the nominal values, while no significant change in the hypocotyl length is observed, we notice an increase in the days to flower, which is also not desirable. Taken together, this result indicates that the phenotypes are sensitive to changes in Pr and the recommended Pr is 24 hours, which is in agreement with the finding of [58].

Effect of collective manipulation of light properties

Having analysed the effect of varying individual light properties, we are now in a good position to investigate the effects of these light properties collectively. The goal here is to find the combination of light input properties that could further improve plant development (i.e., reduced number of days to flower and relatively short hypocotyl length compared to nominal values). From the six light properties analysed previously, the phase of dawn is not sensitive to plant development, while the duration of twilight deteriorates the plant development and thus these two light properties will be retained at their respective nominal values. Next, the collective manipulation will be carried out between the remaining four light properties, and they are done in a 2-stage manner:

  • Stage 1: Finding the best combination of Ph and Pr
  • Stage 2: Applying variations to A and Δl with the optimal combination from Stage 1.

The search for the collective manipulation of Ph and Pr is carried out first because the analyses given above show that these two light properties are most influential to the phenotypes. S2 Fig of S1 Appendix shows the days to flower and hypocotyl length for several practical combinations of Ph and Pr, where we focus only on the combinations that produce the days to flower and hypocotyl length that are not exceeding the nominal values of 32.52 days and 1.90 mm respectively.

Defining the Ph/Pr ratio as the metric to account for the amount of light exposure, where the higher the Ph/Pr ratio, the more the plant is exposed to the light per Pr cycle, we see that in general, the trends of the days to flower and hypocotyl length in Stage 1 are as expected, where days to flower and hypocotyl length reduce with the increase in the Ph/Pr ratio. The collective manipulation of Pr and Ph enables us to investigate the effect of extended combinations of Ph and Pr on the two phenotypes that are not previously seen when they are individually manipulated. When Ph is manipulated individually, we find that the best phenotypes are obtained with Ph = 18 h with Pr = 24 h. When Pr is manipulated individually, we find that the best value of days to flower is obtained with Pr = 32 h with Ph = 16 h, albeit this induces undesirable hypocotyl length. Interestingly, when considered collectively, due to the wider range and combination of Pr and Ph being considered, we see that Pr = 32 h and Ph = 28 h, produces the best phenotypes in terms of the shortest days to flower of 19.4 days and hypocotyl length of 1.60 mm (Fig 6A and S2 Fig of S1 Appendix), which is respectively, a reduction of ~13 days and 0.3 mm from the nominal values.

thumbnail
Fig 6. Effect of collective manipulation of light properties.

(A) Nominal, best and alternate combinations of Pr and Ph on days to flower and hypocotyl length following S2 Fig of S1 Appendix. (B) Optimal Pr = 24 h and Ph = 18 h from Stage 1, best and alternate combination of A and Δl on days to flower and hypocotyl length following S3 Fig of S1 Appendix. The colours used are associated with the same colours used in S2 and S3 Figs of S1 Appendix.

https://doi.org/10.1371/journal.pone.0261281.g006

Despite Pr = 32 h with Ph = 28 h bringing the best improvement to the two phenotypes, having Ph = 28 h means a longer duration of light is required, which may not be beneficial from an energy savings point of view. A detailed look at S2 Fig of S1 Appendix reveals that significant reduction in days to flower can also be achieved with alternate combinations with Pr between 22 to 24 hours and Ph between 18 to 21 hours, where the days to flower is ~20 days and the hypocotyl length is ~1.13 mm. Comparing these alternate combinations with Pr = 32 h with Ph = 28 h, (Fig 6A and S2 Fig of S1 Appendix) we see an almost negligible difference in terms of the phenotype values but a substantial reduction in the photoperiod, from 28 to 18 hours, which could account for substantial energy savings. We will further discuss this in the next section where we introduce the concept of Effective Light Duration for Flowering. Additionally, having light with Pr > 24 h is not a natural light cycle compared to Pr = 24 h and this may affect other phenotype developments (see e.g. [59]). In view of this coupled with a detailed inspection of S2 Fig of S1 Appendix, the recommended optimal combination from Stage 1 is Pr = 24 h and Ph = 18 h.

Following Stage 1, we set Pr = 24 h and Ph = 18 h and proceed to Stage 2 where we collectively manipulate A and Δl. S3 Fig of S1 Appendix shows the days to flower and hypocotyl length for several practical combinations of light intensities and light offsets and we note several interesting observations. In our previous analysis when varying A and Δl individually, the improvement in the phenotypes can be achieved with A > 1 and Δl > 0. While this trend holds when collective manipulation is considered, introducing light offset seems to be able to improve the phenotypes even when A < 1. Achieving improvement with reduced light intensity is desirable as this could help with reducing energy consumption [55]. As shown in Fig 6B and S3 Fig of S1 Appendix, when A = 0.99 and Δl = 0.04, we obtain further improvement to the phenotype outputs, with a reduction to the days to flower by ~1 day, while virtually no change in the hypocotyl length.

Having a light offset means the implementation requires a constant background light, which incurs more energy consumption. A detailed look at S3 Fig of S1 Appendix suggests that if we consider A = 0.99 and Δl = 0, we increase the days to flower by 1 day but with no change in the hypocotyl length. Whether the energy consumption incurred by having a light offset would be higher than having A > 1 depends on various factors that are beyond the scope of this study. Nevertheless, given that the increase in the days to flower is merely by 1 day when A = 0.99 and Δl = 0, we are more in favour of these light properties, primarily due to the non-requirement of constant background light due to the light offset.

Based on the findings from both stages of analysis, the recommended light properties (see also Table 1) corresponding to the optimal plant development is given by Pr = 24 h, Ph = 18 h, A = 0.99 and Δl = 0 with Dw = 0 h and Tw = 0.05 h. Compared to the nominal value, these light properties improve the days to flower by ~11 days (~34% improvement) and produce shorter hypocotyl length by 0.71 mm (~37% improvement).

Effective light duration for flowering

Here, we introduce the concept of Effective Light Duration for Flowering to aid our choice of light property for optimal output phenotypes with energy savings consideration. As the light function turns ON and OFF across the total number of days to flower, the effective light duration for flowering is defined as duration when only the light is turned ON across the total number of days to flower. The effective light duration for flowering is of particular importance especially when correlating the days to flower with the collective manipulation of photoperiod and period in particular the light exposure, i.e., the Ph/Pr ratio.

As an illustration, let us consider the following two configurations of photoperiod and period. In the first configuration, we have days to flower of 32.52 days obtained with Pr = 24 h and Ph = 12 h and for the second configuration, we have days to flower of 24.79 days obtained with Pr = 18h and Ph = 12 h (see also S2 Fig of S1 Appendix). At a first glance, we may think the second configuration is preferred due to the fewer days to flower compared to the first configuration. However, when we compute the Effective Light Duration for Flowering using Eq (4), the duration of light ON are 390.24 hours and 396.64 hours, for the first and second configurations, respectively. Here, we can see that the second configuration in fact uses a longer light duration than the first configuration. This is because despite the two configurations having the same Ph = 12 h, their Ph/Pr ratio differs, and this affects the Effective Light Duration for Flowering. By comparing the Effective Light Duration for Flowering for the best and alternate combinations of the Ph/Pr ratios given in S2 Fig of S1 Appendix, we, therefore, recommend Pr = 24 h and Ph = 18 h as the preferred light property for optimal plant development with the added benefit of reducing energy consumption (S7 Table of S1 Appendix).

Recommended light properties and energy consumption

Here, we explore the added benefit from the energy consumption perspective based on the recommended light properties that we suggest in the previous section. The calculation of the energy consumption is given in Eq (5), which incorporate three light properties: photoperiod, period, and intensity. We will compare the energy consumption of the recommended light properties against two other lighting recommendations. The first recommendation follows the speed breeding approach [17] and the second recommendation is based on the combination of light properties that produce the best output phenotypes as given in S2 Fig of S1 Appendix. Table 2 shows the energy consumption for these three recommendations across one year. Our recommendation utilises the least energy consumption with 3.62 MWh/year compared to the speed breeding and best output phenotype where they consume 4.47 MWh/year and 4.27 MWh/year, respectively, which is a respective 19% and 15.1%, reduction in energy consumption. While the calculation of energy consumption given in Eq (5) is not comprehensive, it suffices to demonstrate the added benefit of lower energy consumption based on our recommended light properties.

thumbnail
Table 2. Energy consumption for different recommended light properties calculated using Eq (5).

https://doi.org/10.1371/journal.pone.0261281.t002

Conclusions

The impact of light on plant development has been subjected to extensive studies as evident by the copious literature on this topic. In view of the huge laborious and timely effort, many of those aforementioned studies often could only consider a fraction of light property before making recommendations of the light input. In this study, using a well-established Arabidopsis circadian mathematical model [24], we have systematically investigated the effect of different light properties to develop a framework for artificial light management for optimal plant development for smart greenhouse applications. In particular, we want to determine whether any further improvement on plant development could be achieved by covering a larger range of light properties in our analysis as compared to previous studies and to suggest a quantitative range of light properties with the added benefit of energy savings.

The input light function considered in this study encompasses six different properties, i.e., photoperiod (Ph), period (Pr), phase of dawn (Dw), light intensity (amplitude) (A), light offset (Δl) and duration of twilight (Tw). We first vary each of the light properties individually while retaining the others and compare the effect with the nominal values of days to flower and hypocotyl length. From the individual light property variation, as expected the most influential light property is the photoperiod as it is associated with the duration of the input light to the plant. On the other hand, the least influential light property is the phase of dawn, where we observe no change in the phenotypes when this property is varied, indicating that the start of light does not affect the phenotypes (see Table 3 for the summary).

thumbnail
Table 3. Summary of the influence of each light property on days to flower and hypocotyl length when varied individually.

https://doi.org/10.1371/journal.pone.0261281.t003

The collective manipulation reveals a more complete influence of light properties on the phenotypes. The collective manipulation suggests to us that the light properties that produce the best value of the phenotype are not necessarily the suitable ones for implementation due to potential high energy usage. Our analysis shows that our recommended alternate light properties could also achieve comparable phenotypes values (i.e., days to flower of 21.62 days and hypocotyl length of 1.19 mm) to the best one (i.e., days to flower of 20.61 days and hypocotyl length of 1.19 mm) but with lower energy consumption with savings up to 15%. The outcome of this study would benefit practitioners of smart greenhouse farming by providing them with an efficient way of artificial light management for improved plant development with the added benefit of improved energy consumption.

So far, our analysis involves only white light due to the absence of comprehensive Arabidopsis circadian mathematical model that incorporates different light qualities (colours). There have been many studies reporting the enhancement of plant development when different light qualities are used instead (see e.g. [5962]). With more smart greenhouse farming begin adopting the use of light qualities in their operation, it would be of great interest to extend our analysis to include the effect of light qualities, which is currently part of our future works.

Supporting information

S1 Appendix. Supporting material: SI tables, SI figures.

https://doi.org/10.1371/journal.pone.0261281.s001

(DOCX)

Acknowledgments

MF would like to thank Laura Roden (Coventry University) for the useful discussion regarding the physiology of plant growth and development. HA would like to thank Rishad Ahmad (University of Nottingham) for the useful discussion regarding the energy consumption of LED lights.

References

  1. 1. Food and Agriculture Organization of the United Nations. The future of food and agriculture: Trends and challenges. Food & Agriculture Org.; 2018.
  2. 2. De Pinto A, Cenacchi N, Kwon H-Y, Koo J, Dunston S. Climate smart agriculture and global food-crop production. PLoS One. 2020;15: e0231764. pmid:32348336
  3. 3. Panayi E, Peters GW, Kyriakides G. Statistical modelling for precision agriculture: A case study in optimal environmental schedules for Agaricus Bisporus production via variable domain functional regression. PLoS One. 2017;12: e0181921. pmid:28961254
  4. 4. Al-Gaadi KA, Hassaballa AA, Tola E, Kayad AG, Madugundu R, Alblewi B, et al. Prediction of potato crop yield using precision agriculture techniques. PLoS One. 2016;11: e0162219. pmid:27611577
  5. 5. Writer G. IoT applications in agriculture. 30 Jun 2020 [cited 2 Jun 2021]. https://www.iotforall.com/iot-applications-in-agriculture/
  6. 6. Gómez-Chabla R, Real-Avilés K, Morán C, Grijalva P, Recalde T. IoT applications in agriculture: a systematic literature review. ICT for Agriculture and Environment. Cham: Springer International Publishing; 2019. pp. 68–76.
  7. 7. Lee J, Kang H, Bang H, Kang S. Dynamic greenhouse supplemental light source control with wireless sensor network. Proceedings of IEEE International Conference on ICT Convergence; Jeju, South Korea, 2012.
  8. 8. Rayhana R, Xiao G, Liu Z. Internet of things empowered smart greenhouse farming. IEEE J Radio Freq Identif. 2020;4: 195–211.
  9. 9. Kami C, Lorrain S, Hornitschek P, Fankhauser C. Light-regulated plant growth and development. Curr. Top. Dev. Biol. 2010;91: 29–66. pmid:20705178
  10. 10. Mahdavian M, Wattanapongsakorn N. Optimizing greenhouse lighting for advanced agriculture based on real time electricity market price. Math Probl Eng. 2017; 6862038: 1–11.
  11. 11. Tewolde FT, Shiina K, Maruo T, Takagaki M, Kozai T, Yamori W. Supplemental LED inter-lighting compensates for a shortage of light for plant growth and yield under the lack of sunshine. PLoS One. 2018;13: e0206592. pmid:30383825
  12. 12. Bian ZH, Yang QC, Liu WK. Effects of light quality on the accumulation of phytochemicals in vegetables produced in controlled environments: a review. J Sci Food Agric. 2015;95: 869–877. pmid:24930957
  13. 13. Zhang X, He D, Niu G, Yan Z, Song J. Effects of environment lighting on the growth, photosynthesis, and quality of hydroponic lettuce in a plant factory. Int J Agr Biol Eng. 2018;11: 33–40.
  14. 14. Hiroki R, Shimizu H, Ito A, Nakashima H, Miyasaka J, Ohdoi K. Identifying the optimum light cycle for lettuce growth in plant factory. International Symposium on New Technologies for Environment Control, Energy-Saving and Crop Production in Greenhouses and Plant; Jeju, South Korea, 2013.
  15. 15. Dorais M. The use of supplemental lighting for vegetable crop production: light intensity, crop response, nutrition, crop management, cultural practices. Proceeding of the Canadian Greenhouse Conference; Ontario, Canada, 2003.
  16. 16. Lee HI, Kim YH. Utilization efficiencies of electric energy and photosynthetically active radiation of lettuce grown under red LED, blue LED and fluorescent lamps with different photoperiods. J. Biosyst. Eng. 2013;38: 279–286.
  17. 17. Watson A, Ghosh S, Williams MJ, Cuddy WS, Simmonds J, Rey M-D, et al. Speed breeding is a powerful tool to accelerate crop research and breeding. Nat Plants. 2018;4: 23–29. pmid:29292376
  18. 18. Ghosh S, Watson A, Gonzalez-Navarro OE, Ramirez-Gonzalez RH, Yanes L, Mendoza-Suárez M, et al. Speed breeding in growth chambers and glasshouses for crop breeding and model plant research. Nat Protoc. 2018;13: 2944–2963. pmid:30446746
  19. 19. Paradiso R, Proietti S. Light-quality manipulation to control plant growth and photomorphogenesis in greenhouse horticulture: The state of the art and the opportunities of modern LED systems. J Plant Growth Regul. 2021.
  20. 20. Minh QT, Phan TN, Takahashi A, Thanh TT, Duy SN, Thanh MN, et al. A cost-effective smart farming system with knowledge base. Proceedings of the Eighth International Symposium on Information and Communication Technology; New York, USA, 2017.
  21. 21. Watson RT, Boudreau M-C, van Iersel MW. Simulation of greenhouse energy use: an application of energy informatics. Energy Informatics. 2018;1: 1–14.
  22. 22. Lork C, Cubillas M, Ng BKK, Yuen C, Tan M. Minimizing electricity cost through smart lighting control for indoor plant factories. Proceedings of the 46th Annual Conference of the IEEE Industrial Electronics Society; Singapore, 2020.
  23. 23. Hotta C, Article from crops to shops: how agriculture can use circadian clocks. J Exp Bot. 2021; erab371.
  24. 24. Seaton DD, Smith RW, Song YH, MacGregor DR, Stewart K, Steel G, et al. Linked circadian outputs control elongation growth and flowering in response to photoperiod and temperature. Mol Syst Biol. 2015;11: 776. pmid:25600997
  25. 25. Salazar JD, Saithong T, Brown PE, Foreman J, Locke JCW, Halliday KJ, et al. Prediction of photoperiodic regulators from quantitative gene circuit models. Cell. 2009;139: 1170–1179. pmid:20005809
  26. 26. Corbesier L, Gadisseur I, Silvestre G, Jacqmard A, Bernier G. Design in Arabidopsis thaliana of a synchronous system of floral induction by one long day. Plant J. 1996;9: 947–952. pmid:8696370
  27. 27. Kunihiro A, Yamashino T, Nakamichi N, Niwa Y, Nakanishi H, Mizuno T. Phytochrome-interacting factor 4 and 5 (PIF4 and PIF5) activate the homeobox ATHB2 and auxin-inducible IAA29 genes in the coincidence mechanism underlying photoperiodic control of plant growth of Arabidopsis thaliana. Plant Cell Physiol. 2011;52: 1315–1329. pmid:21666227
  28. 28. Pokhilko A, Fernández AP, Edwards KD, Southern MM, Halliday KJ, Millar AJ. The clock gene circuit in Arabidopsis includes a repressilator with additional feedback loops. Mol Syst Biol. 2012;8: 574. pmid:22395476
  29. 29. Lai R, Fang J, Jiang Q, Wei Y, Wang T, Liu Y, et al. Design and analysis of the average current-detection method for wide input voltage range constant-current lighting LED driver. Proceedings of 32nd International Symposium on Power Semiconductor Devices and ICs; Vienna, Austria, 2020.
  30. 30. Mean Well Web. [cited 10 Jul 2021]. https://www.meanwell-web.com/en-gb/
  31. 31. Hughes E. Electrical and electronic technology. Pearson; 2016.
  32. 32. Henriques R, Papdi C, Ahmad Z, Bögre L. Circadian regulation of plant growth. Annual Plant Reviews online. Chichester, UK: John Wiley & Sons, Ltd; 2018. pp. 1–29.
  33. 33. Martín G, Rovira A, Veciana N, Soy J, Toledo-Ortiz G, Gommers CMM, et al. Circadian waves of transcriptional repression shape PIF-regulated photoperiod-responsive growth in Arabidopsis. Curr Biol. 2018;28: 311–318.e5. pmid:29337078
  34. 34. Farré EM. The regulation of plant growth by the circadian clock. Plant Biol. 2012;14: 401–410. pmid:22284304
  35. 35. Hoff M. Time to grow: circadian clock controls plant hormone signaling and response. PLoS Biology. 2007;5: e227. pmid:20076684
  36. 36. Inoue K, Araki T, Endo M. Circadian clock during plant development. J Plant Res. 2018;131: 59–66. pmid:29134443
  37. 37. Wenden B, Kozma-Bognár L, Edwards KD, Hall AJW, Locke JCW, Millar AJ. Light inputs shape the Arabidopsis circadian system. Plant J. 2011;66: 480–491. pmid:21255161
  38. 38. Fankhauser C, Chory J. Light control of plant development. Annu Rev Cell Dev Biol. 1997;13: 203–229. pmid:9442873
  39. 39. Fogelmark K, Troein C. Rethinking transcriptional activation in the Arabidopsis circadian clock. PLoS Comput Biol. 2014;10: e1003705. pmid:25033214
  40. 40. De Caluwé J, Xiao Q, Hermans C, Verbruggen N, Leloup J-C, Gonze D. A compact model for the complex plant circadian clock. Front Plant Sci. 2016;7: 74. pmid:26904049
  41. 41. Foo M, Bates DG, Akman OE. A simplified modelling framework facilitates more complex representations of plant circadian clocks. PLoS Comput Biol. 2020;16: e1007671. pmid:32176683
  42. 42. Foo M, Somers DE, Kim P-J. Kernel architecture of the genetic circuitry of the Arabidopsis circadian system. PLoS Comput Biol. 2016;12; e1004748. pmid:26828650
  43. 43. Borrill P. Blurring the boundaries between cereal crops and model plants. New Phytol. 2020;228: 1721–1727. pmid:31571228
  44. 44. Pokhilko A, Hodge SK, Stratford K, Knox K, Edwards KD, Thomson AW, et al. Data assimilation constrains new connections and components in a complex, eukaryotic circadian clock model. Mol Syst Biol. 2010;6: 416. pmid:20865009
  45. 45. Cockram J, Jones H, Leigh FJ, O’Sullivan D, Powell W, Laurie DA, et al. Control of flowering time in temperate cereals: genes, domestication, and sustainable productivity. J Exp Bot. 2007;58: 1231–1244. pmid:17420173
  46. 46. Mallik M, Division of Vegetable Crops, ICAR- IIHR, Bengaluru-560089, India. Flowering control mechanisms in plants and its importance in crop production and breeding. Int J Pure Appl Biosci. 2018;6: 1033–1038.
  47. 47. Jung C, Pillen K, Staiger D, Coupland G, von Korff M. Recent advances in flowering time control. Front Plant Sci. 2016;7: 2011. pmid:28105041
  48. 48. Schäfer E, Nagy F. Photomorphogenesis in plants and bacteria: function and signal transduction mechanisms. Springer Science & Business Media; 2006.
  49. 49. Casal JJ, Smith H. The function, action and adaptive significance of phytochrome in light-grown plants. Plant, Cell Environ. 1989;12: 855–862.
  50. 50. Casal J. Anatomy, Growth and survival of a long-hypocotyl mutant of Cucumis sativus deficient in Phytochrome B. Ann Bot. 1994;73: 569–575.
  51. 51. Schmitt J. Is photomorphogenic shade avoidance adaptive? Perspectives from population biology. Plant, Cell Environ. 1997;20: 826–830.
  52. 52. Baligar VC, Elson MK, He Z, Li Y, de Q, Paiva A, Almeida AAF, et al. Light intensity effects on the growth, physiological and nutritional parameters of tropical perennial legume cover crops. Agronomy. 2020;10: 1515.
  53. 53. Oakenfull RJ, Davis SJ. Shining a light on the Arabidopsis circadian clock. Plant Cell Environ. 2017;40: 2571–2585. pmid:28732105
  54. 54. Chan AH, Chawla R, Johnston AS, Lee KY, Parmar MS. Effect of light intensity on the hypocotyl length of Arabidopsis thaliana during germination. Expedition. 2013;3. Available: http://ojs.library.ubc.ca/index.php/expedition/article/view/184807
  55. 55. Velasco MH, Mattsson A. Light quality and intensity of Light-Emitting Diodes (LEDs) during pre-cultivation of Picea abies (L.) Karst. and Pinus sylvestris L. seedlings—impact on growth performance, seedling quality and energy consumption. Scand J For Res. 2019;34: 1–49.
  56. 56. Jo H-H, Kim YJ, Kim JK, Foo M, Somers DE, Kim P-J. Waveforms of molecular oscillations reveal circadian timekeeping mechanisms. Commun Biol. 2018;1: 207. pmid:30511021
  57. 57. Lim R, Chae J, Somers DE, Ghim C-M, Kim P-J. Cost-effective circadian mechanism: rhythmic degradation of circadian proteins spontaneously emerges without rhythmic post-translational regulation. iScience. 2021:102726 pmid:34355141
  58. 58. Dodd AN, Salathia N, Hall A, Kévei E, Tóth R, Nagy F, et al. Plant circadian clocks increase photosynthesis, growth, survival, and competitive advantage. Science. 2005;309: 630–633. pmid:16040710
  59. 59. Naznin MT, Lefsrud M, Gravel V, Azad MOK. Blue light added with Red LEDs enhance growth characteristics, pigments concent, and antioxidant capacity in lettuce, spinach, kale, basil, and sweet pepper in a controlled environment. Plants. 2019;8: 93.
  60. 60. Dou H, Niu G, Gu M, Masabni J. Effects of light quality on growth and phytonutrient accumulation of herbs under controlled environments. Horticulturae. 2017;3: 36.
  61. 61. Pennis G, Blasioli S, Cellini A, Maia L, Crepaldi A, Braschi I, et al. Unraveling the role of red:blue LED lights on resource use efficiency and nutritional properties of indoor grown sweet basil. Front. Plant Sci. 2019; 10: 305. pmid:30918510
  62. 62. Jeong HW, Lee HR, Kim HM, Kim HM, Hwang HS, Hwang SJ. Using light quality for growth control of cucumber seedlings in closed-type plant production system. Plants. 2020; 9: 5. pmid:32429540