Connecting the dots between computational tools to analyse soil-root water relations

In the recent years, many computational tools, such as image analysis, data management, process-based simulation and upscaling tools, were developed to help quantify and understand water flow in the soil-root system, at multiple scales (tissue, organ, plant and population). Several of these tools work together or, at least, are compatible. However, for the un-informed researcher, they might seem disconnected, forming a unclear and disorganised succession of tools. In this article, we present how different pieces of work can be further developed by connecting them to analyse soil-root-water relations in a comprehensive and structured network. This “explicit network of soil-root computational tools” informs the reader about existing tools and help them understand how their data (past and future) might fit within the network. We also demonstrate the novel possibilities of scale-consistent parameterizations made possible by the network with a set of case studies from the literature. Finally, we discuss existing gaps in the network and how we can move forward to fill them. Highlights Many computational tools exist to quantify water flow in the soil-root system. These tools can be arranged in a comprehensive network that can be leveraged to better interpret experimental data.

Relative distribution of root water uptake between root segments when water is equally available in space (units: %) (Couvreur et al., 2012) High pressure flow meter Device designed to measure the root system conductance by perfusing pressurized water into a root system opposite from the natural direction of the transpiration stream (Tsuda and Tyree, 2000b) Root pressure probe Device designed to measure the hydraulic conductance of a single root through variations of water pressure and flow at the cut end of a root (Steudle and Jeschke, 1983) Cell pressure probe Device designed to measure the hydraulic conductivity of the membranes of a single plant cell by observing the relaxation time of water pressure pulses applied to the cell (Steudle, 1980) RSML Root System Markup Language: File format for the storage of root system data (Lobet et al., 2015) MTG Multi-Tree Graph: File format for the storage of multidimensional tree information, designed for plant models.
( Godin et al., 1999) FSPM Functional-Structural Plant Model: computer model of plant that combine a detailed representation of the plant 3D architecture with functional properties (radial hydraulic conductivity, solute permeability, etc.) (Godin and Sinoquet, 2005) Parameter Fixed input of the models that characterizes a specific property of the system, within the scope of study. For instance, branching rates are parameters in root architectural models State variable Variable that characterizes the state of the system at any moment of the simulation. For instance, the water potential within the plant is a state variable in water FSPM Boundary conditions Variables constraining the model at its external boundaries for the entire duration of the simulation. For instance, in model of water flow within the root system, evaporative demand or soil water potentials are the plant boundary conditions.

Introduction
Water deficit is one of the most dramatic abiotic stresses in agriculture (Cattivelli et al., 2008). It occurs when leaf water supply is limited by either the low potential of soil water, and/or by the high hydraulic resistance of the soil-plant system (Sperry et al., 2002). At this point, the atmospheric demand for water is hardly met, leading to a decrease in water potential within the leaf tissues. As a result, stomata will close, reducing the plant transpiration and photosynthesis.
To investigate when such limitation occurs, the complex plant-soil-atmosphere system is often conceptualised as a multidimensional hydraulic network, in which both soil and root hydraulic properties may substantially control shoot water supply (Draye et al., 2010;Schoppach et al., 2013;Lobet et al., 2014a).
The structural properties of the roots compose the first dimension of the soil-root hydraulic network. Structural properties refer to the physical position and arrangement of the objects of interest. They can be conceptualised at the tissue/organ (transverse anatomy, fig. 1A), plant (root architecture, fig. 1B), or population scale (rooting density profile, fig. 1C).
A second dimension, overlaying structural properties, encompasses the system functional properties. When studying water movement, functional properties often refer to hydraulic conductivities or reflection coefficients. Like in the structural layer, these properties can be defined at different scales. Local radial and axial hydraulic conductivities can be defined at the organ scale ( fig. 1D) while the entire root system of a single plant can be characterized by its conductance ( fig. 1E) and would relate to plant water status (Alsina et al., 2011). An extension of this property to the population scale is the plant population hydraulic conductance per unit horizontal area ( fig. 1F), common in canopy models (Cox et al., 1998), and recently integrated in root models (Cai et al., 2017).
Finally, a third dimension describes the plant environment. In this contribution, we focus on the soil compartment, which includes the rhizosphere and the bulk soil ( fig. 1G). Their respective spatial domains are concentric around individual roots, and their properties differ substantially, so that the rhizosphere is often considered to critically affect plant water availability under water deficit (Carminati et al., 2011;Ahmed et al., 2014). The bulk soil and rhizosphere hydraulic A c c e p t e d M a n u s c r i p t 6 properties may be described by their water retention and hydraulic conductivity curves. The former defines the pressure needed to extract water from the porous media, and the latter the relation between water flux and water potential gradient in space (van Genuchten, 1980). The water potential that defines the energy level of water is a critical environmental variable, driving the flow of water in the soil-plant system. Similarly to the plant, the soil could be divided into functional and structural components and described according to the studied scale. However, we did not explicit this separation in the following as we rather focus on the plant property description in this study.
Each element of the network is dynamic and heterogeneous. Root systems grow, develop and take up water, while soil water content continuously changes in response to root water uptake and climatic conditions, potentially resulting in complex system behaviour. In addition, some key variables and parameters are hard (if not impossible) to quantify experimentally. As a result, the whole system is difficult to apprehend, and novel approaches might prove useful to study it.
In the recent years, many computational tools (image analysis, data analysis, process based modelling and upscaling tools) were developed to help quantify and understand water dynamics in the soil-plant system. Some of these tools were developed to work together, or at least be compatible. However, for the uninformed researcher, they might seem disconnected, forming a collection of tools with, at best, a common target (plant-water relation exploration) but unrelated to each other.
The overall objective of this paper is to draw and discuss the role of a functional landscape of interconnected experiments and models for the study of soil-plant water relations. It is articulated as 3 sub-objectives: (i) to inform readers about existing procedures and tools used for the quantification of water flow at the organ, root system and plant population scales, as well as their interconnections forming a comprehensive, though non-exhaustive, network, (ii) to provide examples of studies combining experiments, analytical and modelling tools in this network, motivating the use of such approaches to enhance interpretations of available and future data, and (iii) to identify gaps in the network and argue for a better integration of future tools in this workflow with appropriate experiment and model design. A web interface was developed to help researchers use the network: it is available at https://plantmodelling.shinyapps.io/water_network/.
Through four examples, we illustrate how the dots, consisting in apparently scattered data and tools, can be connected together in a comprehensive network. These examples span over the different scales mentioned above (organ, root system and population) and for all of them, we  et al., 2015) are fully automated tools that can quantify anatomical properties (such as the number of cells or the mean size of each cell type) but do not provide a digitized cell network.
As a part of the functional layer, cell hydraulic properties are hard to estimate as water fluxes are difficult to measure at this scale. The cell pressure probe enables this estimation from measurements of water pressure relaxation times of individual cells, at a high time cost (Steudle, 1980;Hose et al., 2000). Osmotic pressure can be measured using nanoliter osmometer (Malone et al., 1989) or scanning electron microscopy (Enns, 1998). However, the latter is expensive and generally not part of the standard equipment of a plant physiology laboratory. At the organ scale, the root pressure probe enables the measurement of axial and radial conductivities of root segments (Steudle and Jeschke, 1983) and junctions to the stem (Meunier et al., 2018b). Some properties of the system can hardly be determined experimentally such as the partitioning of water pathways across cell layers (apoplastic or cell-to-cell) (Barzana et al, 2012).
Detailed root cross-section anatomical descriptions and a minimal set of empirical cell hydraulic properties (e.g. permeability of cell walls and membranes) enables the simulation of water flow across root cross sections. Like at other scales, water flow in the system is solved using transfer equation with boundary water pressures and conductances as input parameters. Such a model can estimate the equivalent hydraulic conductivity of the root cylinder as well as the partitioning of water flow between apoplastic and symplastic compartments of the system. For instance, by combining measurements of cell and root permeability with a hydraulic model, (Bramley et al., 2009)demonstrated that water flow is primarily apoplastic in lupin roots. A c c e p t e d M a n u s c r i p t 8 A recent study took advantage of these computational tools to estimate the contribution of pearl millet root types to water uptake. Five types were identified based on cross-section anatomy: primary roots, crown roots and 3 types of lateral roots (Passot et al., 2016). A cross-section was thoroughly digitized for each type of root using CellSet ( 2C). The model was used to estimate the radial conductivity of a typical segment of each root type. In this example, different tools (image analysis and modelling) were combined to estimate radial and axial conductivities, based on easy-to-acquire experimental data (cross section images). While complementary measurements of root hydraulic properties will always remain an asset (e.g. in order to cross-validate the estimated properties), this method opens the way to high-throughput estimations of root hydraulic properties.

Root system architecture
Unlike plant cells, the root system of annual crops has a convenient macroscopic scale and all elements (roots) are visible to the human eye. However, the main difficulties faced when retrieving the root system architecture are the hidden nature of this part of the plant, the large number of elements that can possibly overlay and the fragility of the smallest roots, making the full excavation of a complete intact root system particularly difficult. Direct manual methods exist to measure single root architectural traits such as the angle of crown roots with a protractor (Trachsel et al., 2010) or with the basket method (Oyanagi et al., 1993;Uga et al., 2011), the length of individual roots with a ruler (Pritchard et al., 1990;Trachsel et al., 2010) or a combination of several root architectural traits (Trachsel et al., 2010). However, these manual methods do not give access to the full root architecture.
Several digital tools have been developed and are now widely used to access root architectural traits, mostly from images of root system grown in specific experimental setups (see (Paez-Garcia et al., 2015) for a review of existing root phenotyping strategies). These image analysis tools are listed in www.plant-image-analysis.org and will not be detailed here (Lobet et al. 2013, Lobet 2017. The only point to underline is that each tool generally corresponds to a specific growth medium and image capture technique (eg: RooTrak applies to root system growing in 3D and imaged with X-ray Micro Computed Tomography (Mairhofer et al., 2012). While some of these tools have been designed to retrace a full root system architecture (often with an important manual input), many of them only extract some root architectural traits (eg. mean lateral root length, number of seminal roots, crown root emergence angle…). Furthermore, even with the use of specifically designed image analysis tools, whole root system digitization   et al., 2018), are designed to simulate root systems from a limited number of traits, given as input parameters. The major interest of root architectural models is to generate a large number of contrasted root system architectures. Root system modeling enables the exploration of several variants for the same mean traits and the simulation of contrasted architectures, even from synthetic datasets. These contrasting architectures can then be tested in different scenarios, to identify traits that would be beneficial in challenging environments.
An example that illustrates how root architecture models can be applied to interpret experimental data of other root zone processes is given by Schnepf et al. (2016) Water flow in the root system At the root system scale, understanding which root traits positively influence plant water uptake dynamics for a given pedo-climatic situation remains an important research question. Ideotypes have been proposed, but are always tied to a specific environment (Dencic, 1994;Guswa, 2010;Pagès, 2011;Wasson et al., 2012;Lynch, 2013). Different traits, either functional or structural, have been found to maximize the final crop yield depending on the environment (Tardieu, 2012). Ultimately, we need more than single traits or final yield to have a better understanding of plantenvironment interactions. We need to understand how water flow within the plant is A c c e p t e d M a n u s c r i p t 10 dynamically regulated, both spatially and temporally. Unfortunately, accurately measuring water flow is often the limiting step of the experimental pipeline. Several techniques exist to dynamically measure changes in soil water content, such as X-ray computed tomography (Hainsworth and Aylmore, 1983), Electrical Resistivity Tomography (Garré et al., 2013), neutron tomography (Carminati et al., 2010;Esser et al., 2010;Zarebanadkouki et al., 2014a;Tötzke et al., 2017), light transmission imaging (Garrigues et al., 2006) or Magnetic Resonance Imaging (Jahnke et al., 2009;Pohlmeier, 2010;Rascher et al., 2011). These techniques can be deployed for a relatively high number of plants. However, due to water capillary flow within the soil domain, observed changes in soil water content are rarely (if never) a direct indication of the location of root water uptake. Water uptake rate itself can be estimated using more advanced but timeconsuming lab techniques that use tracers, such as deuterated water that is monitored using neutron radiography (Warren et al., 2013;Zarebanadkouki et al., 2014b).
Functional-structural plant models (FSPM) are often used to decipher plant-environment relationships (Godin and Sinoquet, 2005). FSPMs couple a complete representation of the root system architecture (or whole plant or shoot) with functional properties. Their input parameters are both functional and structural. For FSPMs simulating soil-root interactions, hydraulic parameters can be obtained using a root pressure probe (Steudle et al., 1987) or the outputs of organ scale models but, as stated earlier, are generally difficult to acquire. Thus they are frequently adapted from the available literature. Rhizosphere hydraulic properties can also be coupled to FSPM (Schneider et al., 2009) and would constitute a central component of plant water availability (Carminati et al., 2011). Rhizosphere properties are however difficult to parametrize, and would display complex temporal dynamics (Carminati and Vetterlein, 2013). The FSPM structural input consists of an explicit representation of the root architecture (see Root system architecture section and the related previous case study). Together with the root system geometry, hydraulic properties define the root system hydraulic architecture (Lobet et al., 2014a) and are critical for water stress determination (Leitner et al., 2014b;Vadez, 2014).
Water-related FSPMs provide an exhaustive description of the root water relations (uptake rates, water potentials, etc.) in both space and time. Thus, they constitute an important way to integrate different types of information about properties of the root system and soil state variables in the root zone, which can be obtained experimentally, and to translate this information into a distributed pattern of water flows and local state variables (e.g. water potentials at the soil-root interface) within the root zone. The latter type of information is, as of today, hardly accessible experimentally. An exhaustive review of FSPMs related to water flow can be found in Ndour et al. (Ndour et al., 2017). FSPMs can also be used in so-called inverse modelling studies. In such case, the output of the model is known and the model is used to estimate one of the input parameters. For instance, the most likely distribution of root hydraulic properties (that are usually assumed to be age-and order-dependent) can be estimated using a soil-root water flow model and laboratory measurements (Zarebanadkouki et al., 2016). In this study, measurements of local water fluxes were obtained from neutron radiography at different locations in the root system (Zarebanadkouki et al., 2014a). As the experiment took place in a rhizotron, the root system could be fully digitized using an appropriate image analysis tool, which provided accurate information on the root system topology and positions in space ( fig. 4A and see also Root system architecture section). Water uptake patterns and axial flows within the root system could then be modelled by applying existing water flow equation resolution algorithms (Doussan et al., 1998;Meunier et al., 2017b) to the segmented root system (fig. 4B). The water flow model requires boundary conditions that need to be estimated or measured. In this case, the water potential at the root collar was measured using a pressure probe and root-soil interface water potentials were estimated from water content distribution.
Such a coupling allowed the authors to estimate the parameters of the root hydraulic conductivity function that best fitted the water flow measurements ( fig. 4C and D). These parameters then allowed for novel predictions including water uptake and axial flow distributions everywhere in the root system and not only at observed segments, in homogeneous and heterogeneous soil conditions or under various evaporative demands.

Water flow at the population scale
The population level is a pivotal scale. It interfaces with general circulation models that represent, among others, the circulation of the atmosphere and its interaction with land surface for climate forecasting (Shrestha et al., 2014). It also introduces variables of critical agronomic interest like crop yield per acre (Holzworth et al., 2018). Like at other scales, robust predictions of the system behaviour require the ability to quantify system properties and a proper validation, here involving field scale observations of water fluxes. These fluxes can be estimated with heavy instrumentation and data analytics, for instance using eddy covariance flux towers (McElrone et al., 2013)  A c c e p t e d M a n u s c r i p t 12 system, but rather extract synthetic metrics such as a root length density profile, or root crown data (angles, numbers etc.). Some root architectural traits can be derived from data obtained with these techniques using root architecture models and inverse modeling, as stated above (Garré et al., 2012;Vansteenkiste et al., 2013). Functional plant properties, such as root system and stomatal conductances, which coordinate shoot water supply (Alsina et al., 2011) and underground water uptake distribution (Meunier et al., 2018a), can be characterized on individual plants with low-throughput instruments, such as the high-pressure flowmeter (Tsuda and Tyree, 2000a) or the porometer (Rodrigues et al., 2008), then scaled to the population level using the planting density. The main limitation of plant measurement in the field is often the limited sample size, that might not reflect the general behaviour of the system. The same critique can be made about soil hydraulic properties estimated on small and (un)disturbed samples as they may not be representative of the hydrological behaviour at the population level (Hopmans et al., 2002;Vrugt et al., 2004).
These limitations motivate the use of effective descriptions of population water relations, tailored for this specific scale, such as the transpiration correction for "soil water stress" (Verhoef and Egea, 2014) or one-dimensional soil water and nutrient transfer principles (Baram et al., 2016). Two major methodologies address the parametrization of effective field water relations. First, the artificial neural network approach takes advantage of the availability of large amounts of data to train a model. It was used to predict canopy water fluxes from state variables such as the vapor pressure deficit and soil moisture (Whitley et al., 2012;Couvreur et al., 2016). Second, the inverse modelling approach (as described previously) builds on state-of-the-art models to simulate spatio-temporal series of the system state. The model parameter values producing simulated series that best match field observations are considered optimal and representative of the system behaviour. This approach was used to connect models of soil and plant water flow to observations of soil moisture and transpiration in an almond orchard, in order to estimate soil and plant properties, as well as the hardly measurable leaching of water below the root zone (Vrugt et al., 2001). Numerous variables can be used for inverse modeling, such as soil water content, isotopes distributions or root length density profiles.
Going one step further, simplistic macroscale models can be derived from equations of water flow at a lower scale, offering an interesting trade-off between functional simplicity and realism. This type of model involves scale-consistent properties and processes. A cross-validation is thus possible between parameter values estimated directly at the macroscale of interest (e.g. plant population hydraulic conductance per surface area) and derived from the lower scale (e.g. upscaled values derived from root architectural and hydraulic properties). In order to parametrize such a macroscale model of water dynamics in the soil-wheat system, Cai et al. Regarding soil properties, soil water retention curves were fitted on simultaneous soil water content and pressure head measurements (Cai et al., 2016) with the software RETC (van Genuchten et al., 1991) (fig. 5C). However, the parameters of the soil hydraulic conductivity curve were not experimentally determined ( fig. 5D). The vertical distribution of roots in the field (root length density profiles over time) was extracted from in situ rhizotube pictures, with the software Rootfly (Zeng et al., 2008) (fig. 5E). Its relative distribution is typically used as proxy for the water uptake distribution in uniformly wet conditions (Feddes et al., 1978;Šimůnek and Hopmans, 2009) (fig. 5F). Plant hydraulic properties could not be observed in situ for the wheat population ( fig. 5F-G).
An inverse modelling strategy was therefore used to find the "optimal" soil and plant hydraulic properties ( fig. 5H) that best fitted the observed soil water status dynamics. The optimized plant hydraulic parameters were cross-validated with properties at the individual plant scale. Conductance parameters obtained for winter wheat at the same stage of maturity using the hydraulic architecture approach turned out to be consistent with the inversely modeled properties at the population scale (Cai et al., 2017).
In order to limit the number of parameters, this approach requires the assumption that system properties are time invariant (e.g. soil hydraulic conductivity curve). Because root system conductance tends to scale with root length, the root conductance per unit root length was assumed invariant in order to accommodate for root growth. Such a constraint also matters when accounting for the spatial heterogeneity of root development under different soil/microclimate environments in macroscale simulations (Cai et al., 2018).

Discussions and perspectives
Many computational tools exist to better understand water dynamics in soil-plant systems. These tools span different scales (organ, plant and population), types (image analysis, data storage, simulation models) and computational languages (Python, Fortran, C++, C#, Java, ...). For the average user, this multitude might seem overly complex and hard to understand. Yet, most of the tools could work together and form a continuous network. Using this network, experimental data can be transferred from scale to scale and generate new insights ( fig. 6, and case studies developed above). Modelling tools currently present in the network are listed in table 1. This list is non-exhaustive as the objective of this paper is less to review existing tools than to encourage their integration in order to enhance our understanding soil-plant water relations. For image analysis tools, we refer the reader to the www.plant-image-analysis.org database (Lobet et al., 2013;Lobet, 2017). We created an interactive online visualization of the network, which contains links to the different tools and makes the understanding of Figure 6 easier. We also added a submission form, such as the community could update the network with new (or missing) tools. The online visualisation tool (open-source) is available at https://plantmodelling.shinyapps.io/water_network/.
Identifying gaps in the tool network Analogies between tool connectivity patterns at each scale in figure 6 reveal the existence of network gaps (represented by the dashed lines in fig.6). Yet, these gaps should not be filled just for the sake of symmetry. Here, we analyze what the function of filling these gaps could be.
The plant "structure" row in Fig. 6 has the most striking pattern, with imaging techniques systematically feeding image-analysis tools. These tools extract two types of data: (i) explicit spatial structures (e.g. RSML, MTG), and (ii) structural pattern properties (e.g. growth rates, branching rates). At the population and plant scales, root development models (Dupuy et al., 2010;Pagès et al., 2013) offer the possibility to convert root pattern properties into predicted root structures. While root anatomical patterns can be automatically characterized by image analysis tools such as PHIV-RootCell (Lartaud et al., 2014), no root anatomy development model exists at that scale. In the perspective of generating a mechanistic model of a whole plant from the cell scale (Band et al., 2014), a root development model would become essential. It would fulfill two main functions: (i) conducting predictions and test hypotheses related to root anatomical development, and (ii) allowing the spatial and temporal interpolation of root anatomies between experimental observations. Models using explicit root anatomical structures to test hypotheses about hormone signalling (Stoma et al., 2008), tropisms (Robbins and Dinneny, 2018) or radial water flow (Couvreur et al., 2017) have emerged lately. However, models of axial water flow remain largely underexplored.
In the broadly used Poiseuille-Hagen model, only the quantity and diameter of xylem vessels are accounted for. Yet, multiple studies have demonstrated that xylem porous plates, pit membranes and persistent primary cross-walls limit root axial conductivity (Sanderson et al., 1988;Shane et al., 2000;Choat et al., 2008Choat et al., , 2012Brodersen et al., 2013;Tixier et al., 2014) and affect the profile of water uptake along a single root (Meunier et al., 2017a) or its partitioning among root types (Meunier et al., 2018b;Ahmed et al., 2018). Furthermore, cavitation may substantially alter xylem conductivity, particularly in leaves and stems (Rodriguez-Dominguez et al., 2018). X-ray micro-computed tomography now enables the study of this phenomenon in vivo in all parts of the plant (root, stem and leaf) (Skelton et al., 2017) and also the in vivo embolism removal by root pressure (Brodersen et al., 2018). This process is already simulated by models working at the tissue scale, with experimental functions of percentage loss of conductivity (Mirfenderesgi et al., 2016;Couvreur et al., 2018). However, modelling tools are still missing at the vessel scale and A c c e p t e d M a n u s c r i p t 15 their development may reveal a complexity that is neglected at higher scales. Such numerical representations of xylem vessel structure and hydraulic properties in the axial direction would allow the use of alternatives to Hagen-Poiseuille law (Giordano et al., 1978;Lewis and Boose, 1995). Indeed, Hagen-Poiseuille law describes water flow at equilibrium state but cavitation is a process happening fundamentally out of equilibrium. Lewis and Boose point out that "Ideally, the exact solutions should be used to calculate [volume flow rate] in xylem conduits, but the equations are difficult to solve without the aid of computer" (Lewis and Boose, 1995). Computer availability is no longer an issue and explicit models of xylem flow are now emerging (Mrad et al., 2018). We expect that filling this gap will shed light on the role of cross-walls in the generation of root hydraulic types, and in root-leaf preferential connectivity (David et al., 2012).
Similarly, at the soil-root interface, imaging tools are now available to precisely observe processes at the scale of the soil particle and root hair (Keyes et al., 2013). Connecting such soilroot interface geometrical descriptions to root hydraulic anatomical models would open new avenues to understand how root hairs enhance plant water availability in dry soils (Carminati et al., 2017).
Soil water fluxes were only explicitly considered in the last case study (population scale). In other case studies, soil was either neglected (organ and plant scales) or included as static boundary conditions. However, in all cases, a model of water flow in the soil domain can be coupled to the plant water flow. Such analyses were for instance carried out at the plant level (Javaux et al., 2008;Huber et al., 2014Huber et al., , 2015Postma et al., 2017b) or the population level (Gijsman et al., 2002;Hack-ten Broeke et al., 2016). Such models may incorporate multiple soil characteristics such as macropores (Landl et al., 2016), solute convection-dispersion (Schroeder et al., 2012) or specific rhizosphere properties (Schwartz et al., 2016;Roose et al., 2016). For an extensive review of existing soil models, we refer the reader to Vereecken et al. (Vereecken et al., 2016).

Limitations and future developments
Simulating water fluxes in roots with this collection of tools can either help understanding plant water relations as a main goal or be a tool for further application. These tools could also be used as a side usage of a dataset obtained for other purposes. The advent of imaging in plant sciences A c c e p t e d M a n u s c r i p t 16 and the huge progress made in image analysis allow generating high quality quantitative data of plant structure, suitable for model parameterization. Many models exist at different scales and we highlighted the many possibilities to combine these tools. This set of tools greatly increases the potential of interpretation of experimental data. Yet many authors still publish rich datasets without using modelling tools to interpret them. Using models is not trivial for a large part of the plant science community. Coupling several tools and models together is still rarely applied. Several requirements seem essential to facilitate the use of this pipeline whenever it may add value to the data.
A lot of image analysis tools and models flourish in the water transport domain. Thanks to the wide breadth of scientific literature available (scientific papers, reviews, websites…), developers are usually aware of already existing tools and keep in mind to justify the interest of their tool in this landscape. However, we suggest that further efforts should be made to render the tools compatible with existing ones. In this context, the existence of several tools at the same place of the network (eg. RootTyp and CRootBox for root architecture simulation) is not conflicting. Each one can best suit one scale or specific situation. In our opinion, special attention should be paid to the data format. Indeed, the output data format of upstream tools must be compatible with the input format required by downstream ones. If this is not the case, easy-to-handle tools must exist to convert these data. The multiplication of formats and the need to convert data from one type to another may discourage the use of some of the models. In this respect, the existence of standard formats, such as the Root System Markup Language (Lobet et al., 2015) for root architecture, smooths the interconnection between tools.
When a new tool is created, the documentation of its potential connections with existing ones (e.g. in the user guide) would benefit the whole network. It is indeed expected that the knowledge and the use of all modelling tools will increase. It also underlines the need to keep the models and their documentation updated. Pioneering tools sometimes get outdated by new ones that do similar tasks but that are more user-friendly, faster, use the latest formalisms, or are better connected with newly existing tools. Therefore, either the interconnection between tools needs to be part of a huge maintenance effort for already existing tools, or the acceptance that pioneering tools are doomed to sink into oblivion.
Making the different tools freely available to the community is also a key aspect in their long term maintenance (Lobet, 2017). Many different repositories and licences exist so that everyone should be able to find a combination that suits their (and their institution's) needs. Free access to the tools' source codes would indeed greatly facilitate their evolution, reproducibility of the in silico experiments and allow future developers to interconnect them more easily.     Details of the connected dots for estimate soil and plant-scale conductivities through inverse modelling. Colored parts are the tools, models and inputs used in the approach. Specific tools names were added where relevant. The red arrow indicate the optimisation step used in the inverse modelling. The red boxes highlight the variables evaluated using the inverse modelling. See text for references to letters A-H.

Figure 6
Full network of tools and data used to quantify water flow in the soil plant system. The network connects experimental procedures, computational tools and data related to water flow in the soil-plant system. It is organised by scales (organ, plant and population) and by the types of information (structural or functional, see text for details).A. Tools to quantify the water flow at the organ scale (root section). B. Tools to quantify the structure of root organ. C. Tools to quantify the water flow at the plant scale (root system). D. Tools to quantify the structure of root system. E. Tools to quantify the water flow at the population scale. F. Tools to quantify the structure of root profiles. G. Tools to quantify water flow in the soil.