Laboratory Experiments of Drainage, Imbibition and Infiltration under Artificial Rainfall Characterized by Image Analysis Method and Numerical Simulations

: Two laboratory experiments consisting of drainage/imbibition and rainfall were carried out to study flow in variably saturated porous media and to test the ability of a new measurement method. 2D maps of water content are obtained through a non-invasive image analysis method based on photographs. This method requires classical image analysis steps, i.e., normalization, filtering, background subtraction, scaling and calibration. The procedure was applied and validated for a large experimental tank of internal dimensions 180 cm long, 120 cm wide and 4 cm deep that had been homogenously packed with monodisperse quartz sand. The calibration curve relating water content and reflected light intensities was established during the main monitoring phase of each experiment, making this procedure very advantageous. Direct measurements carried out during the water flow experiments correspond to water content, pressure head, temperature, and cumulative outflow. Additionally, a great advantage of the proposed method is that it does not require any tracer or dye to be injected into the flow tank. The accuracy and other benefits of our approach were also assessed using numerical simulations with state-of-the-art computational code that solves Richards’ equation.


Introduction
The vadose zone (VZ), situated between the soil surface and the groundwater table, supports many processes linked to the hydrological cycle. The dynamics of water in the VZ govern interrelationships between precipitation, surface runoff, infiltration, groundwater recharge and evapotranspiration [1][2][3][4]. The disciplines of soil physics, hydrology and agronomy study the VZ mainly through the lens of water flow and distribution. Furthermore, biological activities and transport phenomena are also affected by the specificities of this unsaturated part of the soil [5,6]. In other words, the VZ takes a central role within the critical zone extending from the vegetation canopy to the bottom of the groundwater aquifer, and as such, it is involved in many environmental and socio-economic issues.
Hence, the quantification of water dynamics in porous media is important to assess the effects of water management actions dealing with agricultural irrigation or water supply plans. The water content, and by extension the saturation of the media, is a key variable involved in the characterization and modelling of flow in the VZ. The relationships between pressure head, water content, and hydraulic conductivity have been widely studied with respect to the water retention curve or soil water characteristic curve [7,8]. Water content is also involved in the transport of solutes and pollutants (evaluation of pollution of groundwater reservoir, salinization, etc.), or in chemical reactions and interactions, because not all the void spaces are connected by water, contrary to what happens in groundwater flow where the total porosity is taken into account [9]. Biological activity and development in porous media also depend on water content due to mobility-attachment effects, amount of dissolved oxygen (aerobic vs anaerobic zones), etc. [5,10].
Consequently, the measurement of water content is of great importance in many applications, both for research investigations and practical guidance. In recent decades, various technical solutions (neutron scattering techniques, gamma attenuation techniques, dielectric techniques, ground penetrating radar, optical techniques, etc.) and logistical supports (data-logger, geographic information system, etc.) have been developed, increasing our observation capacity [11][12][13][14][15][16][17][18]. The chosen and implemented devices can be distinguished in terms of accuracy, robustness, coverage, and volume of measurements, also taking into account their life cycle costs. In this context, the definition of scale for both the measuring and modelling issues remains a crucial point. According to Romano (2014) [14], the spatial scale of soil moisture observations/computations should be defined by the triplet of "support", "spacing", and "extent". As described by this author, the "support" is the volume in which the average value of soil moisture is obtained, "spacing" is the distance between the sensors, and "extent" represents the entire domain over which the measurements are carried out. The chosen technical solution determines the spatial resolution of the sensing probe and, consequently, the support of the measuring scale. In modelling issues, the spatial grid where the simulated outputs are computed represents the extent, which is the main information mentioned that deals with scale. Unfortunately, a mismatch in scale between observations and simulations often occurs (see [14]). The study depicted in this article focuses on an experimental system with external dimensions 203 cm long, 148 cm wide and 8 cm deep. The imaging technique applied allows the combination of local measurements with full photographs to provide 2D maps of water content that can be more efficiently compared to the simulated results.
Hence, in the present contribution, the domain of interest for water content measurement is the local scale. Usually, this means that the characteristic length of the soil volume investigated by the sensor probe ranges from about 0.10 m to 0.50 m. Among the various methods available for measuring water content at the local spatial scale (see [14]), a distinction can be made between direct and indirect techniques. The thermo-gravimetric method, based on the difference between weight measurements before and after drying, is a direct, simple and inexpensive approach to access the water content. Nonetheless, it entails the destruction of the soil sample, and hence prevents the repetition of the measurement at the same location. With respect to indirect measurement methods, the most used techniques rely on neutron thermalization (gamma-ray attenuation or neutron scattering effects), dielectric-based sensing techniques, or other geophysical methods. When considering a substance subject to an electromagnetic field, its dielectric properties depend on the polarization of its molecules and are described by the apparent relative dielectric permittivity. Within the frequency range of roughly 50-100 MHz to approximately 2 GHz, the water content of the soil sample can be correlated with the measured apparent relative dielectric permittivity of soil, which is principally affected by the apparent relative dielectric permittivity of water. Various techniques are traditionally used to access dielectric permittivity, including time domain transmissometry (TDT) or time domain reflectometry (TDR). Also, the capacitance and impedance techniques are used to estimate soil moisture [19][20][21]. Electromagnetic-based devices for the measurement of water content require calibration relationships that can be more or less suited depending on soil properties (compaction, organic matter, …) (see [22]). Nowadays, capacitance and impedance sensors are receiving more attention than TDR devices because of their lower price and required resources. Nonetheless, their accuracy and precision are much more prone to be affected by soil temperature and salinity variations. In the present article, impedance theta probes have been used; previous investigations [23] have given satisfactory results. Their size and technical specifications are appropriate to the dimensions of the domain and the local measurements expected from the experimental setup developed.
To get a profile or a map of water content without disturbing the flow domain, non-destructive measurement options have to be investigated. Non-invasive imaging techniques represent the most promising approach; they can carry out global observations leading to quantitative measurements without inserting (m)any probes in the porous media. Oostrom et al. (2007) [24], in their review, summarized different flow experiments and their associated qualitative/quantitative monitoring devices. In the context of hydrology, they focused particularly on (i) photon-attenuation-based methods, such as the gamma radiation method and the X-ray attenuation method, and (ii) photographic methods, such as the light reflection method (LRM), the light transmission method (LTM), and the multispectral image analysis method (MIAM). For our study, LRM was used because it can be conducted with non-transparent porous media and wide flow tanks; LTM is dedicated rather to transparent media of small dimensions, while MIAM is much more expensive to access saturation of different phases in porous media. Image analysis based on reflected light has largely been used for intermediate-scale experiments (e.g., [23,[25][26][27][28][29][30][31][32][33][34][35][36][37][38]). Application of quantitative photometric procedures for unsaturated flow characterization is less widespread, and even less so with the LRM. The first application began with the LTM when Hoa (1981) [39] proposed a linear relationship between water content and the degree of illumination coming from an emergent light passing through a sand sample and collected by a photosensitive surface. Application of the LRM for characterizing flow dynamics in porous media was mainly restricted to qualitative observations, focusing for instance on the location of the water table. Then, Persson (2005) [40] and Yoshimoto et al. (2011) [41] directly measured the degree of saturation of an unsaturated porous medium with the LRM by relating variations in the color of the medium to water content changes. A feasibility study with small samples of different soil materials was proposed in the technical note of Persson (2005) [40]. Otherwise, Yoshimoto et al. (2011) [41] used a bright red food dye and disconnected calibration and experiments to achieve this goal. Belfort et al. (2017) [23] combined the advantages of these two approaches to determine water content, and Cremer et al. (2017) [42] have proposed an automatic image analysis that detects fingers in infiltration patterns and evaluates their width, velocity, and frequency. Hence, it is important to notice that the efficiency of optical measurement techniques to produce an accurate estimation of the variables of interest (i.e., saturation, concentration, etc.) relies on the quality of image processing. Therefore, the developed procedure has to account for data conversion, evaluation of the resolution size, selection of a sensitive channel, correction of brightness values for the LRM, and definition of a calibration curve. Concerning color specification, it could be defined either with traditional grey-scale values [34,43], or with red, green and blue (RGB) mean intensities [41]. Gerhard and Kueper (2003) [44] found that the best correlation with fluid saturation was obtained for the saturation component using the hue-saturation-brightness (HSB) system. Persson (2005) [40] adopted a linear combination of the saturation and value component of the hue-saturation-value (HSV) system; he also recommended the red channel for the RGB system. Finally, according to   [27] and McNeil et al. (2006) [29], the green channel is the most sensitive for determining concentration; Belfort et al. (2017) [23] found the same result for water content. The majority of these studies used tracers or dyes to improve the visualization of processes using imaging techniques. For instance, Schincariol et al. (1990) [33] used Rhodamine WT liquid in their experiment on saltwater plumes, Konz et al. (2008; [26,27] used Cochineal Red A (E124) and Cremer et al. (2017) [42] used Brilliant Blue dye. In the present study, we adopt the procedure proposed by [23] that characterizes the water content directly from the reflected light intensity of the media without any dye or tracer. The methodology includes, among other things, a normalization to avoid the effects of light fluctuations in contrast to Persson (2005) [40], and also a selection of the sensitive channel. The construction of the calibration curve is also directly derived from the experimental setup.
Several experimental studies dealing with variably saturated flow conditions have focused on issues such as water table recharge, the effects of capillary fringes, heterogeneity and anisotropy on the bench-or field scale (e.g., [38,[45][46][47][48][49][50][51][52]). Other issues concern the estimation of unsaturated hydrodynamic parameters from column experiments, fingering and unstable flows observed during infiltration in dry unsaturated media, rainfall-induced slope failure, device calibration, flux quantification in VZ, etc. [42,[53][54][55][56][57]. As previously stated, photography is currently used for qualitative explanations and observations because of evident visualization benefits; nonetheless, an advanced procedure might be required [42]. Our article focuses on the quantification of water content by photometric procedure, precisely where the neutron methods encountered difficulties due to the large size of the tank and where other methods could be much more invasive or disturbing.
Among the various motivations for developing numerical models in the domain of subsurface and groundwater flow, it should be noted that they are widely considered (i) as interpretative tools for investigating/understanding different involved processes and flow patterns, (ii) as predictive tools for testing different factors/conditions and assessing their impacts, (iii) as supporting tools for identifying practical solutions and planning the collection of field data, and (iv) as visualization tools for communicating key results. Checking simplification hypothesis and validating numerical models are of particular importance for the modelling community. This means that modelers have to verify if their numerical results are in agreement with the physics of the problem being examined, and subsequently if the mathematical and numerical approaches chosen are appropriate. To this end, analytical-numerical or experimental-benchmarking problems are performed and inter-comparisons carried out (for instance, [58,59]). Restricted access to the underground, punctual location of measuring devices with their own limited control volume, along with difficulties in porous media characterization, are significant drawbacks that induce the development of non-invasive procedures, among which the photometric procedure. Even if the goal of this paper is not to validate a new model for water flow in the VZ, we are convinced that numerical results have to be included to underline the benefit of the image analysis method.
Based on a previous article explaining how to construct 2D maps of water content from photometric images and the originalities of the method (see [23]), this article aims (i) to apply and test this procedure on a larger experimental tank with different experimental boundary conditions including drainage, infiltration under artificial rainfall, and imbibition, and (ii) to compare observation-deduced maps with numerical results based on Richards' model for VZ flow.

The 2D Flow Tank and Porous Media
The original flow visualization experiments are conducted in a rectangular flow tank of external dimensions 203 cm long, 148 cm wide, and 8 cm deep. The back, lateral and bottom walls consist of 2-cm-thick polycarbonate sheets to give more rigidity and the front side is fabricated in a 2-cm-thick transparent sheet of Plexiglas for facilitating visual observation of moisture changes. The flow tank is divided into three distinct chambers: a central flow chamber containing the porous medium and two lateral chambers (20 cm 2 × 120 cm) containing freshwater and connected to outflow systems for controlling the application of boundary conditions. The internal dimensions of the porous media region are approximatively 180 cm × 120 cm × 4 cm (L × W × D). The partitioning between chambers is realized with a device of polycarbonate blocks which are variably perforated to simulate different levels of water imposed at the boundaries depending on the performed experiments. O-rings embedded in the side rails provide a seal with forces applied by clamping screws distributed along the length and height of the chambers. Two external aluminum stiffening bars are placed vertically in the middle of the flow chamber in order to minimize bending of the Plexiglas wall under the pressures exerted by sand and water (see Figure 1). At the top, two removable shells are added to maintain a constant spacing of the flow chamber without impeding the filling of the tank. The back wall is drilled with many sensor ports that can be used depending on the experimental setup. The tank has also five openings at the bottom, which are regulated by valves. An opening placed under each lateral chamber is connected to an overflow outlet whose height can be adjusted as required. All openings, at the bottom side and on the lateral separation bars, are covered with a nylon filter mesh to avoid sand discharge as water leaves the tank. The tank is homogenously packed, as far as possible, with monodisperse quartz sand (reference K30, from Kaltenhouse, France). The filling of the flow chamber is achieved in layers of a few millimeters at a time by using a system composed of a funnel and PVC tubes. Overflow outlets are progressively moved higher up to allow more water than sand in the tank, in order to avoid the trapping of air. A rod is also used to ensure that the medium is well placed and is tapped gently with a piece of PVC to produce the most uniform compaction.
Hence, the material surrounding the flow tank is comprised of a digital balance with a weighing capacity of 34 kg and a readability of 0.1 g, two light devices placed horizontally at the bottom front of the tank, and an opaque black curtain that is arranged around the tank and the camera. Notice that two white cards (with dimensions of around 7.5 cm × 14.5 cm) have been attached to the front face of the tank for the purposes of normalization, i.e., to correct possible fluctuations in lighting.

Image Acquisition and Instruments Used to Take Direct Measurements
A Nikon digital camera D80, placed on a fixed metal tripod at a distance of approximately 3 meters from the tank, automatically records the images in the raw format (unprocessed files ".nef"). All camera parameters (ISO100, shutter speed 1/6 s, aperture F11.0) are set manually to take large-sized images with dimensions of 3872 pixels wide by 2592 pixels high. The acquisition device is kept as uniform as possible throughout the duration of the experiment because of the requirements of image processing, primarily pixel-by-pixel intensity tracking at different times.
Direct measurements carried out during the water flow experiments consist of water content, pressure head, temperature, and cumulative outflow. Pressure transmitters used (PR-41X/81955, Keller Druckmesstechnick) are floating piezo resistive transducers, and measure pressure head from −0.1 to 0.1 bar with ±0.2% FS (full scale) accuracy. To ensure good transmission of pressure and suction, a small porous ceramic cup (of 6 mm diameter and 28 mm length) is inserted into the porous media and attached to a Plexiglas tube filled with deionized de-aired water which is connected to each pressure device. After having completely removed air bubbles, the device with ceramic cups may remain saturated up to a pressure of less than 0.1 MPa (or 1 atm).
Water content is measured using three theta probes, Type ML2x (Delta-T Devices Ltd.). This kind of sensor is based on medium impedance measurements [60]. The soil probe consists of four inert metal electrodes that are 6 cm in length and are inserted vertically into the surface of the porous media. The soil sampling volume indicated by the manufacturer consists of a cylinder 4 cm in diameter by 6 cm in length (for a volume of approximately 75 cm 3 ) surrounding the central rod. For a range of measurements between 0.05 and 0.6 m 3 m −3 , the accuracy reaches ± 0.01 m 3 .m −3 . The theta probe was calibrated through column experiments that used different sand/water mixtures prepared over the whole range of water content values (from 0 to θs with a step of 0.05).
Details on the positioning of the pressure and water content sensors, with plans and photographs, are provided in Figures S1 and S2 of the supplementary material.

Method Dealing with Photometric Procedure
The photometric procedure for obtaining water content maps is fully described in a previous article [23], and the present section summarizes the main steps automatically performed by the program:  Importing the images, converting to 16-bit .tiff images and extracting their information.
proc I x z t through background subtraction: It through the insertion of a measurement zone.
A rectangular zone is inserted that corresponds to the projected soil sampling volume of the Theta probe sensor on the front pane; its physical dimensions, 4 cm × 6 cm, represent 87 pixels × 129 pixels on the picture. Circles of 11 px diameter (almost 0.5 cm, i.e., the dimension matching the influence diameter of a ceramic porous cup) can also be inserted at the positions where pressure heads are monitored. Average intensities are automatically computed and stored in the results output file.
 Obtaining    ,, img x z t through image calibration: img proc  Converting processed intensities to water contents to construct moisture maps.  Computing the error of the estimated water content: In Equation (4) s b a is the estimated covariance between â and b (see [23]).

Numerical Method for Variably Saturated Flow Modelling
Richards' equation (RE, see Equation (8)) is used to model variably saturated flow in porous media. This means that the flow velocity is estimated using Darcy-Buckingham's law, i.e., assuming that air remains at atmospheric pressure.
K h H f tt (8) In the mixed form of RE given above, H (L) and h (L) are the hydraulic piezometric and pressure heads, respectively, such that  H h z ; z (L) is the elevation, which is taken to be positive upward;   (10) in which Se is the effective saturation (−),  r is the residual water content (L 3 L −3 ),  is a parameter related to the mean pore size, n is a parameter reflecting the uniformity of the pore size distribution, s K is the saturated conductivity and r K is the relative conductivity, and  1 n .
Parameters can be obtained by inverse modelling (see [53]). In addition to the hydraulic parameters characterizing the porous medium  9) and (10)) requires the definition of initial and boundary conditions (at the surface, bottom and lateral parts of the flow chamber) characterizing the problem studied. In terms of initial conditions, the fully saturated condition is required for the calibration procedure. Full details of the numerical methods can be found in previous research articles [62,63]. Concerning the discretization, an unstructured mesh composed of 11,154 triangles is used for both experiments; this means that the domain is represented by triangular cells of about 2 cm 2 .

First Experiment: Simple Drainage-Imbibition Cycle
The first flow experiment consisted of a single drainage-imbibition cycle. It was carried out over 3.75 days. The flow chamber was initially fully saturated, and the water level was adjusted with overflow outlets so that no ponding conditions were encountered at the surface. To avoid modification of the porous media arrangement in the experimental tank during the experiment, which could affect the results mainly due to background subtraction, three complete imbibitiondrainage cycles were carried out over the previous few weeks. Initial conditions were maintained for 3 hours after the beginning of image acquisition and sensor monitoring. The drainage phase consisted of ten successive steps, followed by ten successive steps of imbibition; each step was characterized by a displacement of the right overflow outlet by 10 cm (see Table 1) and a duration of 3 to 6 hours. The evolution of the lateral boundary condition during the experiment is imposed with the right overflow outlet (relative to the clear Plexiglas front pane), whereas the left outlet is kept in place to conserve the reference of the top water level. Two pressure sensors were connected to valves above the lateral chambers in order to record the level of water on each side of the flow chamber. In this manner, the real boundary condition encountered in our laboratory experiment could be set in the numerical model. The evolution of water in the right and left lateral chambers as a consequence of the overflow outlet displacement is depicted in Figure 2. The correction of fluctuations in brightness is of great importance, and Figure 3 illustrates the necessity of this step. In fact, Figure 3a shows the evolution of raw intensities over the course of the experiment. The standard deviations of the two series corresponding to the right and the left white cards are greater than 300. The fluctuations occurring could be linked to the stability of the electrical power supply which is affected differently during nights and days because of activities in the building. In Figure 3b, Equation (1) is used to normalize the raw intensities. The images are on the same intensity level and standard deviations are limited to 4.93.10 −4 , suggesting that the correction method works well.  In basic image processing steps, filtering is generally applied to reduce random superimposed variations and impulse noise, for instance coming from grain size in a photograph. In this view, "Kernel" operations are generally employed to accomplish neighborhood averaging or ranking on pixel values (see [64]). The median filter consists of replacing a pixel value by the median value of a selected neighborhood pattern. The median smoothing technique is used in studies dealing with contaminant hydrology or porous media characterization, mainly with radii from 1 to 3 pixels (i.e., from 3 × 3-to 7 × 7-pixel squares) [29,34,65,66]. Nonetheless, filtering techniques are sometimes counter indicated, because they can have a negative effect on the sharp edges, i.e., they introduce diffusion, which affects the resolution and the identification of edges [67,68]. In the section devoted to the large flow tank, filtering steps have not been included, mainly because of the larger pixel size. The correction of imaging defects can be linked to the issue of measurement area, which has been considered to a moderate extent in the literature. In [26], the concentration at a point is given by the median of 100 pixels, which corresponds to a surface of 31 mm 2 . [41] considered horizontal samples to determine that the luminance value is stable for a surface of 200 × 200 pixels, i.e., a zone measuring about 1 cm 2 . Hence, for our experiment, we determined the measurement area following the same approach as [26]. Squares of 108 pixels in length (around 5 cm × 5cm) at different locations were considered, and the processed intensities of the 11,664 pixels were stored. Non-corrected pixel values and median smoothing values of radius 1 to 19 px were investigated in order to identify the reasonable definition of the maps that will be compared to modelling results.
In the case of vertical unsaturated water flow, it is difficult to isolate the impact of the optical heterogeneity of porous media. This might be done easily at the beginning or at the end of the experiment with a fully saturated condition in the tank, or when the upper part of the central flow chamber is not submitted to transient flow. Tests were carried out with raw images, normalized and processed intensities at different locations of the tank, and at different times to analyze the optical heterogeneity of porous media. Therefore, Figure 4a illustrates the results in the left upper part of the flow chamber, in a square of 25 cm 2 (11,664 pixels), after a long time of drainage (72 h), allowing consideration of a stabilized condition in the tank. Note that imbibition steps begin without affecting the observation point. The effect of median smoothing clearly appears in Figure 4a, which depicts the 11,664 pixel values derived from 1 × 1 pixel and the median over 5 × 5, 11 × 11 and 39 × 39 pixels. The fluctuation of intensities increases with increasing resolution. The comparison of 95% confidence intervals in Figure 4b shows that the dispersion of results decreases with increasing radius of median smoothing. The purpose of this picture is to identify possible dimensions for the size of the measurement area. A length of around 0.7 cm (radius = 7, corresponding to an area of about 48 mm 2 and a 95% confidential interval less than 1.0 × 10-5) should be the critical minimal size to respect. According to the tests performed at different times and locations, the resolution for the water content maps is based on a unit length of 1 cm; this means that each image is decomposed into a grid with a square mesh size of 1 cm × 1 cm and each cell returns the median value of its interior pixels. These results are in good agreement with previous recommendations by [26] and [41]. Consequently, for both kinds of sensor, pressure transmitters, and water probes, the soil sampling volumes are above this critical size. The calibration curve between light intensity and water content was constructed with 21 images referring to the state at the end of each flow step in order to avoid transient flow (cf. Table 1). Since three theta probes were monitored, a total of  63 m N points are depicted in Figure 5. For each point of the calibration curve, the error in the processed intensities was computed by considering the 10 images recorded during the 10 minutes before each flow step. Hence, error bars have been depicted in Figure 5. As prescribed in the final step of the procedure, a linear regression is performed and the fitting coefficients corresponding to Equation (3), computed for the entire range of points is presented in Table 2.   Table 2. The 68% confidence region, computed with Equation (4), is also shown in Figure 5. As can be also noticed with the high value of the r-squared coefficient ( 2 r in Table 2), these elements indicate a good fit between measurements and confirm the linear regression model used.
To conclude this part, we have to indicate that the calibration component,  We look at the error in water content at specific locations. Figure 6 depicts the evolution of water content during the drainage imbibition experiment for the three theta probe sensors. The photometric procedure is applied in the corresponding zones and the evolution of computed water content is also shown. Then, absolute and relative errors over time are calculated as follows: where j refers to desired printing time (between 1 and jmax = 5399; total duration is 90 h), i corresponds to the different dates of the photographs (in fact each minute),  mes is the measured water content, and  cal is the value obtained with the photometric procedure on the 4 × 6 cm 2 zone corresponding to the sensor location. Absolute and relative cumulative errors are shown in Figure 7. The linear trend of the cumulative error curves until 2.3 days indicates a relatively constant error; after 200,000 seconds, i.e., during the imbibition phase, the error of sensor 3 (and to a lesser extent that of sensor 1) increases. An explanation can be found in the water supply, which is achieved with a single overflow outlet giving a true 2D process in the tank. Figure 7 shows that the water content obtained from photographs is overestimated during drainage and underestimated during imbibition, which is in concordance with previous results obtained with the photometric method [23].   (11)) and relative (dashed lines, computed with Equation (12)) errors between measurements and computation.
The measurement of cumulative outflow mass was performed with a digital balance, even if, for this large flow tank experiment, the measurement was subject to more technical difficulties (presence of lateral chambers, length of pipes, etc.). Figure 8 depicts the mass balance measurements during the drainage steps, the cumulative mass obtained from photographs throughout the complete experiment, and the theoretical calculation at the end of each step. The contributions of the lateral chambers and the fixed overflow outlet are taken into account to correct the results of the theoretical and photographic methods. A good agreement was obtained with the knowledge that the error is amplified by the 3D computation of the cumulative mass. As previously stated, the water content coming from the photometric procedure is overestimated during drainage and the error on the cumulative mass does not exceed 3 kg, which is less than 12%. Figure 9 depicts water content profiles at different times during drainage for both the photometric procedure and numerical simulation. Notice that for the first experiment, we considered that the soil material may be described by the following parameter values: Ks = 5.10 −3 cm.s −1 ,  s = 0.375,  r = 0.099,  = 0.04 cm −1 , n = 7.45 and S S = 1.10 −10 cm −1 (values corresponding to Sim1 in [23]). The general trends are similar, except that the photometric procedure is able to capture the sharpness of the infiltration front, while numerical modelling would require better parameter estimation. The evolution of pressure heads simulated and measured at two specific control points located near the surface of the flow tank is represented in Figure 2. Again, it shows that differences between local measurements and numerical results can be observed, but that the trend is globally well reproduced. Hence, comparing Figures 2 and 10 underlines that photographs are of great interest in following the elevation of the water table and especially the dynamic of the experiment in the vadose zone from a spatial and temporal point of view.

Second Experiment: Drainage-Rainfall Infiltration-Drainage-Imbibition Cycle
The second experiment was carried out as follows (cf. Table 1): 1) the flow tank was fully saturated and observed for 90 minutes; 2) then, the drainage phase was realized in two stages: the overflow outlets were moved downward by 50 cm and 48.5 cm. An equilibrium interval of 5 hours was included between these two steps. 3) After an equilibrium time of 17 hours, a constant rainfall was applied to a central upper zone of about 60 cm for 7 hours. 4) A free drainage step and an equilibrium phase were carried out for 41 hours. 5) Finally, the imbibition phase was achieved in 3 stages with successive increases of the overflow outlets level: an increase of 48.5 cm was carried out, followed by a break of 5 hours, a second rise of 49 cm with a break of 1 hour and a final rise of 1 cm to reach the fully saturated state.
Concerning the rain simulation device, it was manually adjusted with an overall hourly injection rate of about 22.74 liters, and, despite the use of a filter at the pump inlet, we did not operate in a closed circuit to avoid clogging the injection nozzle. On the other hand, this setup allowed continuous monitoring of the cumulative mass of water flowing from the overflow outlets with a digital mass balance. For this second experiment, several soil material properties were compared (sim1 to sim4 provided in [23]) and the values corresponding to sim4 were adopted: Ks = 2.74 10 −2 cm.s −1 ,  s = 0.376,  r = 0.088,  = 0.054 cm −1 , n = 7.37 and S S = 1.10 −10 cm −1 . Due to the 2D nature of the flow, which may be explained with the simulation of rain on part of the porous medium at the surface of the tank, the choice of parameters had a greater impact on the numerical results. Figure 10 displays the cumulative mass of water collected at the outlets of the overflow during the parts of the experiment composed of the two-step drainage and the rainfall and stabilization periods. The comparison between measurements and simulated results shows very good correlation and the relative error in the volume of water injected during rainfall does not exceed 3%. Figure 10. Evolution of the cumulative mass of water leaving the tank during the drainage and rainfall part of the 2nd experiment; the red line depicts the results of simulation and the black line represents the measurements.
As shown in Figure 11, our methodology for image analysis works very well for creating water content maps in a more complex flow situation than drainage and imbibition. The simulation results presented at different times of the experiment are in close agreement with the water content maps. Even if the medium is homogeneous, and while the simulated maps are naturally very regular, the photographs make it possible to illustrate certain irregularities in the water flow that could be due to disparities in the compaction of the porous medium during filling. Compared to the results obtained with a small tank (L: 40 cm × W: 14 cm in [23]), Figure 11 illustrates the effectiveness of the method, which does not generate too many discrepancies in working with a larger tank; it could be noticed that the mesh size of the picture (1 cm 2 ) was not modified, even if this would have allowed the water content maps to be smoothed. The presence of the reinforcement bars on the experimental device cannot be completely removed during image processing, but the various steps make it possible to limit the impact of the shadows generated by the two projectors.  Figure 11. Results obtained for the second experiment including rainfall: presentation of the original photographs, maps resulting from image analysis method (right side) and numerical simulation (middle).

Conclusion
The purpose of the present study was to test the methodology proposed by Belfort et al. (2017) [23] in a more complex flow situation and a much larger experimental tank. Indeed, the internal dimensions of the flow chamber are 120 cm × 180 cm × 4 cm (L × W × D), which corresponds to a viewing area 90 times larger than for the validation experiment [23]. As for the experimental conditions, in addition to the drainage and imbibition steps, a rain simulation was carried out on the central upper part of the tank, introducing true 2D flows.
First, we demonstrated the interest of the different steps of our methodology; mainly, normalization and then filtering according to a chosen radius. The calibration line was constructed according to a well-established protocol based on different drainage and imbibition levels. The confidence interval and error bars have been depicted to show the quality of the calibration. For the first experiment, the absolute and relative cumulative errors are very good, with global amplitude that does not exceed 8%. The monitoring of the cumulative mass of water leaving (or entering the interior of) the flow chamber also shows a good match between measurements and results from the photometric method. For the vertical line located in the middle of the tank, the representation of simulated water content profiles and those obtained from the photos are very close together. In the second experiment carried out, which included a step of constant rainfall simulation, the cumulative mass was depicted to validate the numerical simulation compared to the measurement. The constant intensity of rain was clearly observed with a steady state occurring in the tank for both the experiment and the simulation. Then, water content maps obtained from image analysis methods and numerical models were depicted. It appears that the photometric method did not lead to regular maps compared to those resulting from simulations, and that instabilities were present. Nevertheless, the trend of the infiltration front, even during typical 2D flow related to rain infiltration, and the values of the water content in the main zones of the domain show great similitudes. Finally, the photometric method gave satisfactory results even in more complex conditions, and this article confirmed the potential of our previously published methodology. Again, we underline the interest of applying the methodology directly during the experiment because of the specific structure of the porous medium and environmental conditions. Our previous study concluded by expressing our desire to develop a form of benchmark problem for unsaturated flow inspired by Vauclin's experiment [51]. With these two experiments, we can say that this objective has been achieved, since we have a wide range of water content and pressure measurements, as well as complete water content maps. Now, all this data could serve to develop and test parameter estimation methodologies directly in 2D without having to limit themselves to spot measurements. Ongoing work would also consider heterogeneous domain and/or more complex geometries.

Supplementary Materials:
The following are available online at www.mdpi.com/xxx/s1, Figure S1: Plans of the tank, Figure S2: Pictures of the tank. Videos of both experiments are available at https://zenodo.org/record/3462133 Author Contributions: B.B., S.W. and F.L. took part in the design conception and the experimental setup; they achieved measurements and developed the methodology of image analysis. Marwan Fahs was involved in the numerical developments. Benjamin Belfort worked on the writing of this paper with the support of all the co-authors.

Funding:
The study was partially financed by Strasbourg University and ENGEES; we acknowledge their scientific councils for their confidence.