Gloxinia – An Open-source Sensing Platform to Monitor the Dynamic Responses of Plants

Olivier Pieters Tom De Swaef Peter Lootens Michiel Stock Isabel Roldán-Ruiz Francis wyffels
Pre-print published on 25 March 2020
Keywords: plant monitoring real-time data acquisition sensor platform phenotyping
This is a pre-print version

Abstract

The study of the dynamic responses of plants to short-term environmental changes is becoming increasingly important in basic plant science, phenotyping, breeding, crop management, and modelling. These short-term variations are crucial in plant adaptation to new environments and, consequently, in plant fitness and productivity. Scalable, versatile, accurate, and low-cost data-logging solutions are necessary to advance these fields and complement existing sensing platforms such as high-throughput phenotyping. However, current data logging and sensing platforms do not meet the requirements to monitor these responses. Therefore, a new modular data logging platform was designed, named Gloxinia. Different sensor boards are interconnected depending upon the needs, with the potential to scale to hundreds of sensors in a distributed sensor system. To demonstrate the architecture, two sensor boards were designed: one for single-ended measurements and one for lock-in amplifier based measurements, named Sylvatica and Planalta, respectively. To evaluate the performance of the system in small setups, a small-scale trial was conducted in a growth chamber. Expected plant dynamics were successfully captured, indicating proper operation of the system. Though a large scale trial was not performed, we expect the system to scale very well to larger setups. Additionally, the platform is open-source, enabling other users to easily build upon our work and perform application-specific optimisations.

Introduction

Plants that grow in natural or agricultural environments are exposed to substantial short-term variations in environmental conditions. For instance, the occurrence of clouds or waving leaves can modify the light environment within seconds; relative humidity and temperature can change within minutes due to precipitation. Because fitness productivity of plants often lies in their ability to swiftly respond to these highly variable conditions, studying these dynamic responses is crucial [@lawsonStomatal2014]. As a result, the need for monitoring the response time of plant behaviour increases in different research fields of plant science. First of all, in basic plant science where the understanding of the short-term responses to environmental variation is central [@barillotWhat2010,windtMRI2006,behrensRapid1982,hubeauPlantPET2015]. Additionally, they are relevant in more applied phenotyping or breeding research, where the more dynamic behaviour of certain genotypes might be key to their success [@lawsonStomatal2014,rascherNoninvasive2011,caldeiraHydraulic2014]. Furthermore, in crop management where for instance irrigation scheduling can be done using continuous measurements of stem diameter [@deswaefStem2015,deswaefDetermining2009,fernandezIrrigation2010]. Finally, as crop and plant models become more dynamic and mechanistic, the need for dynamic data to calibrate and validate these parameter-rich models increases [@deswaefUse2019].

In high-throughput phenotyping, plant traits are extracted through image analysis and sometimes complemented with scoring from breeders. These techniques allow monitoring of a large number of plants both in controlled settings such as growth chambers and greenhouses, and in the field [@walterPlant2015; @liReview2014]. Depending on the set-up, the camera system can be mounted on drones [@arausField2014], tractors [@busemeyerBreedVision2013], ground-based robotic systems [@shafiekhaniVinobot2017; @andrade-sanchezDevelopment2014] or conveyor belt-based systems [@granierPHENOPSIS2006; @hartmannHTPheno2011]. These allow to extract very detailed data and specific plant traits, but their temporal resolution is generally low, ranging between once per day (e.g., a conveyor belt system) and once per week (e.g., drone flights). Consequently, they are not well suited for studying plant responses at a time scale of seconds or minutes. Other sensor types, such as porometers, cannot be permanently mounted on the plant and thus need manual intervention. Measurement devices for gas exchange (e.g., Li6800, LI-COR Biosciences, USA) can be attached to individual plants to record short-term changes of important features like stomatal conductance, photosynthesis, and transpiration automatically, but cannot be deployed on a larger number of plants due to the cost of these devices. Other contact sensors, such as sap flow sensors or Linear Variable Displacement Transducers (LVDT), are interesting solutions for short-term continuous monitoring [@deswaefStem2015; @clearwaterPotential1999; @steppeSap2015], although they require considerable time to install and might need some maintenance once installed.

Concurrent with the requirement for dynamic plant data, monitoring of a plant’s micro-environment is crucial for interpretation, considering that this is the driver for the plant’s response. Environmental data of a field trial is typically captured at a single measurement location. However, assuming the same environmental conditions for the whole field or greenhouse is not necessarily an accurate representation of reality [@xuEnvirotyping2016; @tardieuPlant2017]. Some extreme examples are the differences in microclimate between the top and bottom of a sloped field [@bennieSlope2008], shading caused by trees next to the field [@artruImpact2017], but also differences in temperature within greenhouses, or even growth chambers [@cabrera-bosquetHighthroughput2016; @granierPHENOPSIS2006].

Monitoring on the timescale of seconds or minutes, necessary to capture the dynamics of plant traits such as leaf temperature, stomatal conductance, photosynthesis, and transpiration, along with their environmental drivers, is not a trivial task [@inoueNormalized2008; @wallachDevelopment2010; @bakerChlorophyll2008; @maesEstimating2012; @jonesUse1999; @costaThermography2013]. First, interfacing and synchronising sensory readout is difficult when different sensors are combined [@neveuDealing2019]. As a result, it is necessary to rely on autonomous data acquisition systems to read out the sensors and store the data, and investigate alternative sensors that do not require manual intervention. Second, commercial data loggers such as CR1000 (Campbell Scientific, USA), DL2e (Delta T Devices, UK), ZL6 (METER Group, USA) or EM50 (ICT International, Australia) are expensive, making it difficult to employ them for hundreds of sensors in large trials. Third, alternative data loggers such as single-board computers or microcontroller boards are cheap, but have limited analogue readout capabilities due to low accuracy analogue-to-digital converters (ADC) (typically in the 10 to 12 bit range) and do not always feature all the necessary hardware to directly interface with various analogue sensors. BeagleBone Black (BeagleBoard.org Foundation, USA) is a single-board computer with analogue capabilities, and Arduino UNO (Arduino, Italy) is a very popular microcontroller board.

From the foregoing, we summarise that there is a need to monitor the dynamic behaviour of plant traits. To this end, a larger amount of sensors needs to be employed that is read out at a higher frequency. This has to occur cost-effectively for various purposes, including basic plant science, breeding, agronomy, and environmental monitoring. Moreover, autonomous data logging systems that do not sacrifice accuracy for cost or ease of use are key to tackle this need. We present an implementation of a data logging system that is designed to address four key needs: sensor scalability, accuracy, cost, and versatility with regards to experimental size and sensor interfaces. First, the system should easily scale to hundreds of sensors without needing a large number of hardware boards. Second, the system should provide accurate sensor readings, with limited influence of noise sources. Third, the system should be low-cost. Fourth, the system should be sufficiently versatile; it should work well in trials that monitor a handful of plants to hundreds of plants. Additionally, most common sensor interfaces should be available to connect various sensors such as light, relative humidity, temperature, and soil moisture sensors. Finally, the system should also have open-source hard- and software, enabling others to build upon this work and tailor it to a specific application.

System Architecture

To meet the key requirements introduced earlier, we selected a distributed sensing architecture. Sensors are connected to nodes. Nearby sensors are connected to the same node, while others are connected to other nodes. Nodes collect the sensing data over a certain area and are interconnected on a linear bus to form a linear sensor array. One of these nodes needs to be connected to a storage device such as a computer or Raspberry Pi. This node aggregates the sensor data and transmits it for storage. The overall architecture is depicted in Figure 1. Different sensor nodes are connected over a Controller Area Network bus (CAN bus) network that operates at a maximum rate of 1 Mbps. The CAN bus is very robust. It is commonly used in automobiles to allow microcontrollers to communicate without the need for a master controller, thanks to a fixed arbitration scheme [@robertboschgmbhCAN1991]. Nodes can easily be added; the only requirement is that the first and last nodes use a 120 Ω resistor to terminate the bus. The bus speed can also be lowered to accommodate longer bus lengths up to 5000 m at 10 kbps [@corriganController2008]. In each node, sensing and control have been split.

To meet the key requirements introduced earlier, we selected a distributed sensing architecture. In this type of architecture, sensing is scattered across different devices, removing the need for a single measurement device that reads all deployed sensors. Nearby analogue or digital sensors are connected to the same measurement device, called node, while others are connected to other nodes. As a result, nodes only have to read out sensors in close proximity, alleviating the need for expensive low-noise cables. To facilitate the readout of data, nodes are interconnected on a linear bus. This bus is used to send measurement and configuration data. Consequently, only a single node on the bus needs to be connected to a storage device such as a computer or Raspberry Pi. The overall architecture is depicted in Figure 1.

archtictural overview archtictural overview legend
System architecture of a typical experiment where multiple plants are monitored using the same set of sensors. Analogue sensors are connected to Planalta and Sylvatica, depending on the type of sensor. For instance, a Linear Variable Displacement Transducer needs a variable input signal, so this type of sensor has to be connected to Planalta because this board can generate analogue signals. Sensors that only require an analogue readout are connected to Sylvatica. Environmental conditions are also measured with digital sensors at node 1 and node N. The data is transmitted and stored on a single board computer (Raspberry Pi) for further analysis. This figure is best viewed in colour. Raspberry Pi illustration by, Lucas Bosch, based on work by Jonathan Rutheiser CC BY-SA.

In challenging conditions such as those present in the field, a robust and high-speed bus has to be selected for the communication between nodes. The Controller Area Network bus (CAN bus) meets these requirements. It operates at a maximum rate of 1 Mbps, while being very robust. It is commonly used in vehicles for communication between microcontrollers without the need for a master controller, thanks to a fixed arbitration scheme [@robertboschgmbhCAN1991]. Nodes can easily be added; the only requirement is that the first and last nodes use a 120 Ω resistor to terminate the bus. The bus speed can also be lowered to accommodate longer bus lengths up to 5000 m at 10 kbps [@corriganController2008].

A Universal Synchronous Receiver-Transmitter connection (UART-connection) is used to interface between the data storage device and the node connected hereto. This computer also provides an interface to configure the nodes.

Each sensor node is comprised of a control board (named Dicio), and one or more sensor boards. These sensor boards vary depending upon the application. Communication between the control and sensor boards is done using the Inter-Integrated Circuit (I2C) protocol. The control node controls this bus and reads each of the sensor boards at predefined time steps [@nxpsemiconductorsUM102042014]. Optionally, a phase-locked signal can be used to achieve cycle synchronisation between sensor boards on a single control node. By default, the I2C protocol operates at 400 kHz, but can also be lowered to 100 kHz for sensors that do not support fast I2C operation.

Two I2C buses are available, one for the readout of the custom sensor boards through a dedicated interface connector, and one for other sensors. Sensors such as digital relative humidity and temperature sensors are connected using screw terminals.

The dual bus system accommodates varying needs. First, the CAN bus interconnects different sensor nodes, thus simplifying the readout system and data storage interconnect. Only a single node needs to be connected to the computer. Furthermore, the CAN bus is robust in noisy environments and well-suited for communication over longer distances (modifiable in the software). The maximum distance is linked to the signalling rate and can be increased at the expense of lower throughput. The system uses a single cable to deliver both signal, power, and synchronisation signals. A synchronisation signal can be used to synchronise sampling between different nodes. Second, the I2C protocol was selected since it is widely used for digital sensors and allows the sensors to be read directly from an Arduino or Raspberry Pi, thus facilitating faster prototyping and stand-alone usage of the sensor boards. Furthermore, most microcontrollers provide a hardware I2C-interface, which reduces cost and complexity. Since the interconnect distance between the sensor boards and the control board is small, there will be less noise and interference. Consequently, the requirements of this bus are a lot more tolerant.

Two boards that cover a wide range of analogue sensors have been designed: Planalta and Sylvatica; Measurement System Design provides more details on the available functionality. These sensor boards cover the most common analogue interfaces found in plant monitoring. Digital sensors should interface directly with the Dicio board. Dicio supports common digital interfaces, including , Serial Peripheral Interface (SPI) and RS-232.

Several Planalta or Sylvatica boards can be connected to the same Dicio control node to increase the number of sensors that can be read while keeping the amount of redundant hardware to a minimum. The two sensor boards, Planalta and Sylvatica, each serve a different purpose. The board named Planalta is designed for sensors that require a variable input voltage such as soil capacitance, LVDT, and impedance sensors. The measurement principle of this board relies on a digital lock-in amplifier (LIA) (see also Digital Signal Processing. The other sensor, named Sylvatica, is designed for sensors that do not require an input signal, like most analogue temperature and relative humidity sensors, for example.

Both sensor boards have generic analogue interfaces and can easily be used with a wide range of sensors that require low input voltages and currents. Sensors should have an operating voltage between 0 and 12 V in the case of Sylvatica and 0 to 3.3 V for Planalta. High-power sensors such as some types of sap-flow sensors require an external power supply. An attenuator is required for both boards if the readout value can be higher than 3.3 V. To keep the hardware low-cost, only single-ended signals are supported. The Planalta board features an option to use a mid-rail referenced signal.

Measurement System Design

Some of the identified design criteria are conflicting, such as low unit cost, and accuracy and versatility. Consequently, a trade-off was made. The main cost in this system is comprised of the necessary components. As a result, the measurement system has been designed to minimise the number of analogue components, thus lowering the overall cost at the expense of more digital processing. Since the software development cost is a one-time investment, most of the signal processing and filtering are done digitally to reduces the hardware cost per unit. Furthermore, the software is easily modifiable due to a simple and generic design. This is ideal for application-specific optimisations.

Integrating all sensing functionality on a single board was, however, deemed too complex and would increase the cost per board while some of the hardware would remain unused. Hence, two sensor boards are designed, Planalta and Sylvatica.

The Planalta board is optimised for sensors that measure modulated signals. A lock-in amplifier is well-suited for this purpose. A digital lock-in amplifier modulates a voltage (or current) based on a reference clock at a specific frequency. This signal is then deformed by the sensor-response, resulting in a new signal that has a different phase and amplitude than the original signal. An ADC whose sample points are synchronised to the same reference clock then digitises the analogue signal. The amplitude and phase can then be determined with high precision, based on the reference signal. Digital LIAs have significant advantages over their analogue counterparts, including better noise performance, phase stability and orthogonality due to the lack of temperature and frequency-dependent drift [@dixonBroadband1989]. Sensors that can take advantage of this measurement principle include LVDT, impedance, and laser-based sensors.

Sylvatica is a board designed for single end-measurements. Sylvatica does not feature an analogue sine wave generator, nor a connection from the input voltage of the sensor to the ADC. Instead, it supports double the number of sensors compared to Planalta (eight sensors vs. four sensors) because more ADC channels are available for sensor outputs.

In what follows, the analogue front-end is discussed first, followed by the digital signal processing of the sampled signals.

Analogue Front-end Design

The front-end design focuses on simplicity and flexibility by using a generic design that is optimised digitally for the application’s needs. The driver uses well-chosen components to limit both non-linearity and noise of the circuit before digitisation.

A general overview of the analogue front-end is depicted in Figure 2. Dashed components are specific to the Planalta board. The left part of the Planalta-only circuit generates a sine wave from a pulse-width modulation (PWM) signal. By varying the duty cycle at a high frequency, a sine wave is formed after low-pass filtering. This signal is then buffered before it passes through the sensor. Buffering ensures that there is no voltage drop due to variations of the load impedance. The buffer can provide up to ±5 mA of current to the sensor. For high-power sensors, an external amplifier has to be used. On the Planalta board, the output voltage signal is also provided as an input to an ADC channel after low-pass filtering.

analogue front end analogue front end legend
Block diagram of the analogue front-end for Sylvatica and Planalta boards. The Planalta-only driver generates a voltage wave that is generated through a pulse-width modulation (PWM) signal, which is low-pass filtered and then buffered before being fed into the sensor. The sensor output on both the Planalta and Sylvatica boards is (optionally) amplified using a Programmable Gain Amplifier (PGA) and low-pass filtered before being sampled.

The circuit that connects the sensor output to the ADC has the same topology on both Planalta and Sylvatica. The output of the sensor is amplified using the PGA113 (Texas Instruments, USA) programmable gain amplifier (PGA), providing a gain between 1 and 200. This PGA can amplify with respect to an offset, on Sylvatica this offset is connected to ground potential, while on Planalta this offset can be soldered to either ground potential or 1.65 V. The signal is then low-pass filtered using a resistor-capacitor filter before digitisation. Digitising the signal as early as possible limits the amount of noise that can enter the system and increases flexibility, since the digital filtering is easy to modify.

The ADC used in this design is ADS8332 (Texas Instruments, USA), and provides a good trade-off between speed, accuracy, and cost. This ADC has a successive approximation register (SAR) architecture, produces 16-bit output data, and has an integrated multiplexer that can rapidly switch between channels, enabling up to 8 analogue signals to be sampled between 0 and 3.3 V. The 16-bit words are well-suited for further processing by the 16-bit microcontroller. The effective ADC resolution at frequencies below 1 kHz is 14.9-bit. Additional details on the analogue front-end are provided in Analogue front-end specifications.

The PGA meets the noise specification for the first design criteria of the selected ADC for the first three gain settings (1, 2 and 5). However, the PGA does just not meet the distortion requirement for a unity gain configuration. The PGA has a Total Harmonic Distortion (THD) of -104 dB, and the ADC -100 dB at 1 kHz. Though, this should not be a limiting factor, since the signal frequencies of interest are typically much lower for Sylvatica. The digital filter on Planalta has a very narrow bandwidth, so this should also be sufficient. A simulation was not performed since there is no simulation model available for this PGA.

For sensors that use a constant supply voltage, there are three possible supply voltages: 3.3 V, 5 V and 12 V. All these supplies are located on Dicio, so all boards that are connected to the same Dicio control board have the same supply voltages available. Each sensor board features one additional power supply of 4 V to power the buffer, PGA and ADC. This is the only supply that is not available to the user from one of the connection ports. Furthermore, this supply voltage generated independently on each sensor board. The 3.3 V supply can deliver up to 1 A.

The 12 V switching power supply enables Sylvatica to interface with common environmental sensors such as HMP50-L and CS616 from Campbell Scientific, USA. The on-board supply can deliver up to 35 mA.

A stable 5 V reference is not always available, for example when the Dicio main power supply is derived from Universal Serial Bus (USB) or when the supply is located at a large distance from the node. For these situations, a power supply (RP604Z501B, Ricoh, Japan) can be soldered on Dicio to stabilise the 5 V supply rail and deliver up to 300 mA.

Digital Signal Processing

After digitisation, the microcontroller further processes the samples to remove unwanted noise and interfering signals such as those originating from the 50 or 60 Hz power grid. The data rate coming from the ADC depends on the board and software configuration. The data rates that are possible in the current design are listed in Table 1.

The maximum sampling frequency of the ADC is 500 ksps. Though, in practice, the upper sampling speed cannot exceed 250 ksps due to limitations of the SPI module of the microcontroller. On the Planalta board, the number of active channels determines the maximum signal frequency for the LIA. The total sample rate (sum of sample rates of every channel) should never exceed 250 ksps.

The sampling frequency is always four times the signal frequency. This simplifies the mixing operation significantly, since the in-phase (I) and quadrature (Q) components of the sine and cosine are simply 0, 1, 0 and -1, and 1, 0, -1 and 0, respectively for 0, \pi/2, \pi, 3\pi/2. These operations involve single-cycle copy and invert instructions. An overview of the whole digital processing cascade is depicted in Figure 3. The whole filter structure is replicated eight times for the four sensors and four reference signals. The decimation factor and filter coefficients of the last stage depends upon the configuration. To obtain an output frequency of 1 Hz, the decimation factor of the last stage has to vary between 2 and 20.

For the Sylvatica board, the signal frequencies of interest are typically much smaller. Biological and environmental sensor responses vary in the range of a few Hertz and below [@inoueNormalized2008; @wallachDevelopment2010; @bakerChlorophyll2008; @maesEstimating2012; @jonesUse1999; @costaThermography2013; @xuEnvirotyping2016]. Consequently, there is no need to maximise the sampling frequency other than to limit the amount of noise. Some plant processes such as the absorption of photons by chlorophyll molecules and chlorophyll fluorescence after photon incidence occur much faster, within 1e-15 s and 1e-9 s respectively [@taizPlant2010]. However, these events cannot be observed by simple sensors in a greenhouse or on the field, so they are not considered here.

To avoid active filters, the sampling frequency is increased to 10 kHz per channel. However, this shifts complexity from the analogue to the digital domain. This is desirable since it is easier to modify and optimise the software for a particular application. An overview of the digital processing is shown in Figure 4.

Overview of the different sampling and signal frequencies employed in the sensor boards Planalta and Sylvatica. N_c is the maximum number of active channels, where R denotes that the output of the signal driver is also sampled (top ADC block in Figure 2). f_c is the analogue cut-off frequency of the anti-aliasing low-pass filter in front of the ADC. f_s is the sampling frequency of a specific ADC channel. The frequency depends upon the number of active channels. f_\text{signal} is the (modulated) signal frequency. The filter cascade reduces the input signal to a 1 Hz signal, as indicated by f_\text{out}.
board N_c [-] f_c [kHz] f_s [kHz] f_\text{signal} [Hz] f_\text{out} [Hz]
Dicio - - - - -
Planalta 1 97.0 200 50\times10^3 1
1+R 39.3 80 20\times10^3 1
2+2R 19.5 40 10\times10^3 1
4+4R 10.3 20 5\times10^3 1
Sylvatica 8 4.8 10 \le0.2 1
schematic representation of the data processing Planalta
Schematic representation of the data processing and flow on the microcontroller of the Planalta board. The data coming from the analogue-to-digital converter ADC contains data from all channels. These samples are copied to new vectors by a demultiplexer (demux) to make the filtering faster. The incoming data is then band-pass filtered before being mixed to remove unwanted signals that the mixer can map to the same frequencies. Afterwards, the data passes through four filter and decimation stages before its final filtering and decimation step. The decimation factor N depends upon the incoming data rate to achieve an output frequency of 1 Hz.
schematic representation of the data processing Sylvatica
Schematic representation of the data processing and flow on the microcontroller of the Sylvatica board. The data coming from the analogue-to-digital converter ADC contains data from all channels. These samples are copied to new vectors by a demultiplexer (demux) to make the filtering faster. The filtering consists of a cascade of low-pass filters, followed by a decimation step of factor 10 to reduce the sample rate.

Digital Interfacing with Sensor Boards

As mentioned, I2C is used to read the sensor data. On this bus, one device has control and can read and write to other devices. In this setup, Dicio is the control device that reads the sensor data from the sensor boards Planalta and Sylvatica. Each device, except for the controller, has a 7-bit address that must be unique on the bus.

To ensure that only valid data are read, a ping-pong buffer system is used at the output. One buffer, buffer A, is written by the software with new data, while the other, buffer B, can be read by the user. When buffer A is full, the roles are revised. Buffer B is written, and buffer A is read. Each buffer stores one sample (this can be increased in the software) and thus has to be read every second. Note that for Sylvatica a sample consists of a single 16-bit value, but for Planalta a single sample consists of up to two or four 16-bit values: two values representing the I- and Q-components of the sensor signal and optionally the driver I- and Q-components.

Results

To validate the system, a small prototype was constructed consisting of one Dicio board, two Sylvatica boards, and one Planalta board. An experiment was conducted during ten days in a growth chamber of 1.45 \times 0.77 \times 1.45 m (height \times depth \times width) at Flanders Research Institute for Agriculture, Fisheries and Food (ILVO), Melle, Belgium with a custom-built frame of 1.00 \times 0.70 \times 1.10 m (height \times depth \times width). All lighting was mounted on this frame, including 32 LED lamps (MAS LED spot VLE D 4.9-50W GU10 927 60D, Koninklijke Philips N.V., The Netherlands) and eleven halogen lights (DECOSTAR 51 PRO 50 W 12 V GU5.3, OSRAM GmbH, Germany). The experiment applied a simulated day-night cycle on two strawberry plants (Fragaria \times ananassa, labelled S1 and S2) and one plum tomato plant (Solanum lycopersicum L., labelled T). The strawberry plants were mature plants, grown in a greenhouse at ILVO during the previous year. The main difference between them was their size, S1 had significantly more leaves and S2. The leaves of S2 were also less green than those of S1. Their pot sizes were the same. The tomato plant was a five-week-old seedling.

The two sensor boards, Sylvatica and Planalta were used to perform the readout of several contact sensors that were connected to the plants at a rate of 1 Hz. We employed four leaf thickness sensors, two leaf length sensors, one soil moisture sensor, one relative humidity and temperature sensor, and one light intensity sensor. The environmental sensors were mounted on a separate board, which was glued to a 3D-printed radiation shield. An overview of the different connections to the sensor boards and the monitored plants is depicted in Table 2. More detailed information about the different sensors is provided in Table 3.

Overview of the sensors connected to each plant. A leaf thickness sensor was connected at the start of the experiment to the tomato plant, but there was an issue with the connection. As a result, this data was not valid and not included in the analysis.
plant sensor type
leaf thickness leaf length soil moisture
strawberry 1 (S1) 2 1 1
strawberry 2 (S2) 2 0 0
tomato (T) 0 1 0
Overview of the different sensors used in the experimental validation of the system.
sensor decription part number and manufacturer interface
leaf thickness AH-303 (AgriHouse, USA) analogue readout
leaf length (LVDT) E100 (Chauvin Arnoux, France) LIA readout
relative humidity and temperature SHT35 (Sensirion AG, Switzerland) digital I2C sensor
light intensity APDS-9301-020 (Broadcom, USA) digital I2C sensor

At the start of the experiment, the plants were watered and left to stabilise for one day in the growth chamber before monitoring started. During the monitoring experiment, the plants were watered twice after visual observation of wilting of S1 just before noon on 24 December and 30 December. These time points are indicated by a dashed green line on Figure 5. A picture of the experimental setting just after the second watering time point (30/12/2019 at 10:45) is depicted in Figure 6. S1 wilted, while the other plants did not show any visual sign of wilting.

Both the leaf thickness and leaf length sensors have a strong temperature dependence, which is eliminated using a simple calibration procedure. A simple temperature calibration was performed. The temperature was gradually increased from 10 °C to 32 °C. For the leaf length sensors, the AgriHouse Calibration Card (AH-300C) was used to calibrate. The leaf length sensors were calibrated using four reference distances.

To calculate the soil water content from the capacitance reading, the following calibration procedure was followed: first, the soil was saturated with water for five days; second, it was left to drip to remove excess moisture for 1 hour; third, the soil was left to dry at ambient temperature conditions for 14 days during which the weight of the pot and sensor readouts were recorded; finally, these were combined with the dry weight and volume of the port to calculate the amount of water in the soil per volumetric unit.

The LVDT sensors, used for leaf elongation measurement, do not provide absolute values of leaf length. Therefore, the first measurement was taken to be the reference distance and was set to zero for both sensors.

A more detailed zoom of a 12 h period is of the grey area is depicted in Figure 7. The captured data contains less noise than expected from Figure 5. Moreover, the oscillatory behaviour of the sensors appeared due to the functioning of the growth chamber. It causes the environmental conditions to oscillate around a predefined setpoint, which in turn are main drivers for the plant response.

experiment overview
Visualisation of the captured data in a growth chamber experiment with a strawberry plant. The dashed green line indicates the watering time point. A detailed figure of the grey shaded line (16:30 on 26 December to 3:30 on 27 December) is depicted in Figure 7.
set-up
(a) Experimental set-up. The environmental sensors are not depicted in this figure, but their cable is. (b) Close-up of the radiation shield that houses the environmental sensors.
zoom of experiment data
This figure is a detailed zoom of the grey shaded in Figure 5 to visualise the accuracy of the system. A further zoom of the grey shaded area in this figure (18:30 to 19:00 on 26 December) is depicted in Figure 8.
double zoomed experimental data
Detail figure of the grey shaded in Figure 7 to visualise of the accuracy of the system. Only the leaf thickness and leaf length of S1 are plotted since there is a significant offset between the leaf thickness of S1 and S2 and leaf length of S1 and T in Figure 7.

Discussion

Evaluation of the Experiment and Future Improvements

Figure 5 depicts the most interesting sensor data throughout the entire experiment. The visual observations of wilting are supported by the readings from the leaf thickness sensors. S1 was wilting (see Figure 6). The leaf thickness, presented in Figure 5(c) (blue), gradually decreased from the previous watering event on 24 December towards this time point of visual observation on 30 December at 10:45. Leaf thickness decreased from around 115 μm to 75 μm in the first drying period with decreasing soil water content down to 200 g/l. In the second drying period, soil water content decreased down to 125 g/l, resulting in a minimum leaf thickness of 53 μm. After re-watering, the leaf thickness quickly recovered. S2 was not wilting, which is supported by a more constant pattern of the leaf thickness towards this time point of observation (Figure 5(c), in green). Once wilting started, there was also a clear decrease in leaf thickness for S1 to approximately 75 μm the first time and 53 μm the second time. The leaf thickness quickly recovered after watering. Furthermore, there was no clear difference before and after the watering time point for S2, indicating that this plant was not perceiving drought stress in this period. S2 was probably not wilting due to its lower leaf area compared to S1, while the pot sizes were the same. As a result, we presume that the water content in the pot of S2 was still sufficient since the water content in both pots was the same when the experiment started.

For S1, the leaf elongation did not show drastic variations in response to the drying conditions, other than a small gradual growth during the experiment. The leaf elongation of the tomato plant T did not show a decrease in response to limited water availability. However, a marked increase in leaf elongation coincided with the two re-watering events, indicating that leaf elongation was slowed down before re-watering.

Besides the effects of drying and re-watering, the leaf thickness and leaf elongation sensors also demonstrated a pattern of shrinkage during the day and increased during the night. Indeed, leaf length and thickness decrease when water loss due to transpiration is not fully compensated by the water uptake, and increase when transpiration decreases and tissues are replenished with water. The pattern of swelling and shrinking was most explicit in the leaf elongation of the tomato leaf T. Additionally, the leaf thickness in plants S1 and S2 did not increase throughout the experiment, as these leaves had already reached their final leaf thickness. The elongation sensors on plants S1 and T, showed a gradual increase in leaf length, indicating that these leaves were still expanding. The elongation of the young tomato leaf was much faster than that of the leaf of the strawberry plant S1 [@deswaefPlant2015; @fabbriAnatomical1986].

From Figure 8, we can conclude that the leaf thickness is strongly influenced by the relative humidity, where the drop in relative humidity corresponds to a similar decrease of the leaf thickness around 18:35. The time offset is probably due to the heterogeneously of the air in the growth chamber. A similar effect is observed between 18:55 and 18:58. When the relative humidity decreases, water will evaporate more quickly, resulting in reduced leaf thickness. Since there is a variation  °C, there is a limited effect of the temperature. The leaf length variation remains limited. Compared to leaf thickness, leaf length features less variation on a short timescale for strawberry plants in this experiment.

As expected, soil water content consistently decreased after watering. However, slight increases at the start of the day are due to the temperature dependence of the sensor. A possible explanation is that the soil temperature was not measured, only air temperature. As a result, there is some over-compensation when the temperature changes drastically due to the slower temperature increase of the soil.

Design Validation and Comparison to Existing Platforms

In the introduction, four key design criteria were identified. A comparison between our design, a commercial data logger, a single board computer with analogue capabilities (BeagleBone Black), and a microcontroller platform (Arduino) is depicted in on a 1 to 5 scale. The higher the scale, the better the performance for this criterion.

platform comparison
Spider chart comparing different data logging approaches. The higher the number, the better the performance for this specific metric. For instance, a score of 5 on the cost scale means this has the best cost performance, i.e. lowest cost.

We evaluate the requirements for two experiments: a smaller experiment where one plant is monitored closely with 15 analogue sensors, and a larger experiment where ten plants are monitored in a greenhouse with 15 analogue sensors connected to each plant. Additionally, the environment will be characterised at each plant in both experiments for temperature, relative humidity, and light intensity. In addition to analogue sensors, a digital sensor is employed that uses I2C to measure the temperature, relative humidity, and light intensity. To simplify the comparison, it is restricted to sensors that do not require an input waveform. All sensors are sampled every 10 seconds.

To evaluate the platforms for each of the criteria, we define the four measures, one for each of the design criteria. First, to assess the sensor scalability, we compare the average number of boards per sensor for each of the trials. Second, to assess the accuracy, we compare the number of bits of the output sample. While this is not the actual accuracy, this estimate provides a first indication thereof. Third, the cost is defined as the average cost per sensor of the entire readout system per platform. Fourth, versatility is assessed by qualitative comparison of the difference between the first three design criteria for the two experimental set-ups aforementioned.

The official Arduino Uno board can interface with up to six sensors on a single board, with a resolution of 10-bit. A single board costs 18.81 € (Mouser, USA, 30 January 2020). Different boards are connected using I2C since a robust protocol is not supported without additional hardware. As a result, three boards are needed for the first experiment and 25 for the second [@Arduino].

The BeagleBone Black Rev. C board can also interface with up to six analogue sensors per board. The resolution is 12-bit per sensors and costs 62.75 € per board (Mouser, USA, 30 January 2020). Different nodes are connected using the CAN bus. As a result, three boards are needed for the first experiment and 25 for the second [@kridnerBeagleBone].

Commercial data loggers are very popular and widely used by researchers in plant and environmental monitoring trials. The cost of such commercial systems is approximately 1000 to 1500 €. A single data logger can typically measure 16 sensors sequentially. And have a resolution of 12-bit to 13-bit. Thus, one data logger is needed in the first experiment and ten in the second. However, in a real set-up researches will typically use a multiplexer to readout all sensors and keep the overall cost of the setup more manageable.

The Gloxinia system discussed here requires the Dicio control board and the Sylvatica sensor board to interface with these sensors, one Sylvatica board can measure up to eight analogue sensors and has a resolution of 16-bit. Different Dicio boards are connected using the CAN bus. A cost overview is depicted in Table 4, where a categorical separation is made. For the first experiment, one Dicio board is needed, while for the second ten boards are required. Per Dicio board, two Sylvatica sensor boards are necessary in both cases.

Cost calculation assuming ten boards are produced of a particular type. The cost calculation is based on the Mouser inventory and component prices of 17 December 2019. Enclosure is optionally available from Hammond Manufacturing for 20.21 €. All prices are excluding VAT.
board cost category (€) total cost (€)
capacitors
resistors
connectors
microcontroller
user interaction
amplifiers
voltage reference
other
ADC
Dicio 3.16 0.74 10.43 6.68 5.09 4.64 1.9 32.64
Sylvatica 3.41 1.25 3.52 6.68 3.26 15.84 4.81 1.35 9.44 49.56
Planalta 6.47 2.03 2.96 6.68 3.26 16.44 4.81 5.25 9.44 57.34
Detailed comparison and summary of the evaluation of the different sensor platforms. The cost values are based on the cost per sensor for the small first experiment. The following scores are given from lowest to highest for a particular metric: −, +/−, +, ++, and +++.
criterion platform
Arduino BeagleBone Black commercial Gloxinia
scalability - + + +++
accuracy [bit] 10 12 12/13 16
cost [per sensor, €] +++ + - ++
versatility +/- ++ +++ ++

Future Improvements and Possibilities

While not tested in the experimental setup, the system should easily scale to large trials that need to monitor sensors over large distances thanks to the CAN bus. Theoretically, there is an upper limit of approximately 400 Dicio boards that are connected to a single CAN bus, based on the differential input resistance and drive capability of the MCP2542FD CAN transmitter. However, we advise that no more than 100 Dicio nodes are connected to the same bus for error-free operation. Each Dicio node supports up to four Planalta and five Sylvatica sensor nodes. A second CAN bus has to be used in case more sensors have to be measured.

To stimulate usage of this data logging tool, both hardware and software have been open-sourced in a GitHub repository. All relevant files can be downloaded from GitHub.

Conclusions

The Gloxinia sensor platform aims to advance monitoring in fundamental and applied plant research from modelling to irrigation and crop management. Four key needs were identified: sensor scalability, accuracy, cost and versatility. The whole platform has been designed to address these needs with an open-source design. The platform comprises of individual sensor nodes that communicate with each other. Each node has a control board Dicio to which sensor nodes are connected. Sylvatica and Planalta are two sensor boards that provide an interface that matches most analogue sensors used in plant research. Digital sensors can also be connected to the control boards. Most of the application-specific optimisations are done in software, making it easier for the user to optimise for a specific application. To validate the accuracy of the system, an experimental trial was set up in a growth chamber. Environmental conditions, leaf length, and leaf elongation were successfully measured at high resolution on one tomato and two strawberry plants to validate the functionality of the system. The overall system scales well due to the multiplexed sampling of up to eight channels on Sylvatica and four on Planalta, accurate 16-bit data acquisition, lost unit cost, and distributed architecture. Consequently, the system strikes a good trade-off between these various requirements, making it well-suited for research, breeding, and precision crop phenotyping.

References