Bayesian Simultaneous Estimation of Unsaturated Flow and Solute Transport Parameters from a Laboratory Infiltration Experiment

Numerical modeling has become an irreplaceable tool for the investigation of water flow and solute transport in the unsaturated zone. The use of this tool for real situations is often faced with lack of knowledge of hydraulic and soil transport parameters. In this study, advanced experimental and numerical techniques are developed for an accurate estimation of the soil parameters. A laboratory unsaturated flow and solute transport experiment is conducted on a large undisturbed soil column of around 40 cm length. Bromide, used as a nonreactive contaminant, is injected at the surface of the undisturbed soil, followed by a leaching phase. The pressure measurements at different locations along the soil column as well as the outflow bromide concentration are collected during the experiment and used for the statistical calibration of flow and solute transport. The Richards equation, combined with constitutive relations for water content and permeability, is used to describe unsaturated flow. Both linear and non-equilibrium mobile–immobile transport models are investigated for the solute transport. All hydraulic and mass transport parameters are inferred using a one-step Bayesian estimation with the Markov chain Monte Carlo sampler. The results prove that the pressure and concentration measurements are able to identify almost all hydraulic and mass transport parameters. The mobile–immobile transport model better reproduces the infiltration experiment. It produces narrower uncertainty intervals for soil parameters and predictive output concentrations.


Introduction
Numerical modeling of water flow and solute transport through the vadoze zone is essential to estimate the amount of recharge and/or to prevent soil and subsurface water contamination. It is becoming increasingly important for several applications as water quality management, water resources planning in urban zones, and for mitigating groundwater pollution. The processes of flow and mass transport in the unsaturated zone depend on the driving forces (matric and gravitational potential) and on the soil properties. Thus, to perform real field simulations, the numerical models require an accurate prior knowledge of the soil parameters. These parameters are usually determined using laboratory unsaturated flow and solute transport experiments performed using soil samples collected from the field. Inversion methods can then be applied to fit the model outputs to the laboratory measurements and to assess the soil parameters.
Different laboratory experiments have been conducted for the estimation of the soil parameters e.g., [1][2][3][4][5][6][7][8][9][10][11][12][13]. Small scale laboratory column experiments were often used for the characterization 2 of 16 of hydraulic and mass transport properties [1,2,11,14,15]. Mishra and Parker [16] showed that the simultaneous estimation of hydraulic and mass transport parameters yields more accurate results than the two-steps estimation in which hydraulic and transport parameters are determined sequentially. Inoue et al. [17] used electrical conductivity and matric pressure head measurements at different depths for the flow-transport inversion with the local search Levenberg-Marquardt algorithm. Recently, several Bayesian approaches, where measurements are combined with prior parameter information to provide posterior parameter distributions, have been investigated for the estimation of the unsaturated hydraulic soil parameters, among others [10,12,[18][19][20][21]. The term Bayesian is used to describe statistical inversion by considering [22]: (i) That model variables are random, (ii) that randomness describes the degree of information for their realization, and (iii) that the solution of the estimation problem is the posterior probability distribution from which several statistics can be obtained. Bayesian hydraulic soil parameter estimation was investigated for a laboratory drainage experiment in [10]. In [12], it is shown that the infiltration laboratory experiment yields more accurate hydraulic parameters than the percolation-drainage experiment. Bayesian inference of both hydraulic and transport unsaturated soil parameters has been investigated for hypothetical experiments in the case of linear transport [21] and in the case of pesticide transport [23] in unsaturated soils. Moreira et al. [11] considered the estimation of hydraulic and solute transport parameters based on a Bayesian inversion of a laboratory infiltration experiment performed on a small column of 7.3 cm length. One tensiometric sensor was used to monitor the pressure head in the middle of the column. The pressure head and outflow concentrations measurements were employed for parameters estimation. Because of the insufficient collected data, some hydraulic soil parameters were not included in the inversion procedure to prevent identification difficulties. These parameters were considered as known (measured) and they were assigned informative Gaussian priors. Furthermore, Moreira et al. [11] used a sequential procedure in which the hydraulic parameters were firstly estimated using the pressure head measurements. Then, in a second step, these estimates were employed as prior information with Gaussian distributions to determine the mass transport parameters from the concentration measurements.
In the present study, a laboratory infiltration experiment is conducted in order to estimate soil hydraulic and mass transport parameters. Bromide is used as a nonreactive contaminant and injected at the surface of the soil column. The current research differs from the literature by (i) considering a large column of around 40 cm length of an undisturbed sandy soil, (ii) performing a simultaneous estimation of all hydraulic and mass transport soil parameters using the outflow breakthrough curve of bromide and the pressure head measurements at three different locations, (iii) assuming no prior information for all hydraulic and solute transport parameters for which we attribute uniform prior distributions with large intervals to reflect our poor prior knowledge of their values, and (iv) estimating all parameters using a one-step Bayesian procedure via the Markov chain Monte Carlo (MCMC) method [24].
The direct problem is based on the Richards equation coupled with the Darcy's law and the mass conservation equation [25]. The Mualem and van Genuchten models are used as constitutive relations [26,27]. The linear advection-dispersion and the non-equilibrium mobile-immobile models are investigated for bromide transport. Bayesian inversion results are discussed in terms of fitting between simulations and measurements and in terms of uncertainty ranges associated with the estimated parameter values and to the model outputs.

Laboratory Experiment
A cylindrical PVC tube of 60 cm length, beveled at its end to form a cutting edge, was sunk into the soil in the region of Bekalta in the east center of Tunisia. A homogeneous soil column of 38.4 cm length and 15 cm diameter has been carefully extracted. The soil, essentially formed by sand, was placed over a porous plate of 6 mm high with a saturated water content of 0.42 cm 3 cm −3 . The soil column (sand + porous plate) is set on a container filled with gravel ( Figure 1). An orifice, located face to the soil bottom, allows effluent recovery and fixing the pressure head h [L] at the lower boundary to zero. The vertical axis Z is downward oriented with the origin located at the soil surface. Three tensiometric sensors, connected to a data logger and a computer, allow monitoring the pressure head near the top (Z 1 = 2.5 cm), middle (Z 2 = 17.5 cm), and bottom of the column (Z 3 = 32.5 cm). The outflow bromide concentration is measured using an ion selective electrode (HI4102 from HANNA instruments). Initially, hydrostatic equilibrium is obtained by verifying that each tensiometer indicates a suction equal to its distance to the bottom of the column. A constant water flow rate of 14 cm 3 /min was applied at the top surface of the soil using a high-accuracy Masterflex ® L/S peristaltic pump. Cotton fiber wicks were deposited at the surface of the soil to obtain a uniform distribution of the flow rate over the soil area. The tensiometric data and the cumulative outflow were monitored until t = 660 min. Infiltration occurs during a period of 562 min, it starts at t = 24 min and was stopped at t = 586 min. During a first period (24 min < t < 220 min), the pump injects a potassium bromide (KBr) solution from a reservoir with a concentration of 0.08 g/L. The period of injection ensures around one pore volume of injected KBr solution. Then, clean water was injected at the surface of the soil with the same flow rate for the period 220 min < t < 586 min. This ensures around two pore volumes of injected clean water. To limit evaporation, the upper end of the column was sealed with the exception of a small orifice, which ensures equilibrium with atmospheric pressure. face to the soil bottom, allows effluent recovery and fixing the pressure head h [L] at the lower boundary to zero. The vertical axis Z is downward oriented with the origin located at the soil surface. Three tensiometric sensors, connected to a data logger and a computer, allow monitoring the pressure head near the top (Z1 = 2.5 cm), middle (Z2 = 17.5 cm), and bottom of the column (Z3 = 32.5 cm). The outflow bromide concentration is measured using an ion selective electrode (HI4102 from HANNA instruments). Initially, hydrostatic equilibrium is obtained by verifying that each tensiometer indicates a suction equal to its distance to the bottom of the column. A constant water flow rate of 14 cm 3 /min was applied at the top surface of the soil using a high-accuracy Masterflex ® L/S peristaltic pump. Cotton fiber wicks were deposited at the surface of the soil to obtain a uniform distribution of the flow rate over the soil area. The tensiometric data and the cumulative outflow were monitored until t = 660 min. Infiltration occurs during a period of 562 min, it starts at t = 24 min and was stopped at t = 586 min. During a first period (24 min < t < 220 min), the pump injects a potassium bromide (KBr) solution from a reservoir with a concentration of 0.08 g/L. The period of injection ensures around one pore volume of injected KBr solution. Then, clean water was injected at the surface of the soil with the same flow rate for the period 220 min < t < 586 min. This ensures around two pore volumes of injected clean water. To limit evaporation, the upper end of the column was sealed with the exception of a small orifice, which ensures equilibrium with atmospheric pressure.

Numerical Model
The 1D unsaturated flow through the column can be described using Richards equation (RE):

Numerical Model
The 1D unsaturated flow through the column can be described using Richards equation (RE): where θ is the water content, t is time, C s (h) represents the specific moisture capacity, S s is the specific storage, θ s is the saturated water content, q d is the water velocity, and K(h) is the hydraulic conductivity. Equation (1) is closed using the constitutive relations of Mualem [26] and van Genuchten [27] expressing the hydraulic conductivity in terms of saturation and the saturation in terms of pressure head, respectively.
where θ r is the residual water content, S e is the effective saturation, K s is the hydraulic conductivity at saturation, m = 1 − 1/n, and α and n are the shape parameters specific to the soil.
The flow system of Equations (1) and (2) is subject to the following boundary conditions. At the upper boundary, we impose a constant flux (q inj = 0.08 cm/min) during the injection period (24 min < t < 586 min); otherwise, the upper flux is zero. A zero-pressure head (Dirichlet boundary condition) is imposed at the column lower boundary. Initial conditions correspond to a hydrostatic pressure distribution (h = z − L).
Two transport models are used to investigate bromide transport through the studied soil column. The first model is the well-known linear transport model, based on the advection-dispersion equation: where C is the bromide concentration, D = a l q d + d m is the dispersion coefficient, d m is the molecular diffusion coefficient fixed to a small value of 3.10 −4 cm 2 /min as common in literature, and a l is the soil dispersion coefficient. The transport and flow equations are coupled together through the water content θ and the Darcy velocity q d . The second transport model is based on the non-equilibrium approach. It is usually called mobile-immobile model. This model proceeds by considering that the liquid phase in the porous medium can exist as mobile and immobile phases. No advective solute transport is assumed in the immobile zone and a first-order kinetic diffusion transport process occurs between the mobile and immobile zones. The mobile-immobile model in variably saturated porous media writes [28,29]: where subscripts (m) and (im) represent the mobile and immobile phases, respectively, ω is the first-order exchange coefficient between phases, and θ is the total water content. We assume that the mobile and immobile water contents can vary in space and time but the ratio of mobile water content to the total water content ( f = θ m /θ) is assumed to be constant. As a consequence, the fraction of immobile water θ im is also constant during the transient flow (see Equation (6)). This assumption is less restrictive and more representative of reality than assuming a constant immobile water content θ im during the unsaturated flow [11]. For both transport models, initial conditions correspond to a zero-solute concentration inside the whole column. Boundary upper condition (z = 0) corresponds to a Dirichlet condition with the concentration fixed to that of the bromide solution during pollutant injection and to zero during clean water injection. The lower boundary condition (z = L) corresponds to a zero concentration gradient.
Numerical solution of flow and transport equations in unsaturated porous media is a challenging problem [30][31][32]. The systems of equations in both linear transport and mobile-immobile transport cases are solved with the standard finite volume method. A uniform mesh of 192 cells of 0.2 cm is employed for the discretization of the soil column and porous plate. This spatial discretization yields mesh independent solutions. Temporal discretization is performed with a variable order ODE (ordinary differential equation) solver via the method of lines (MOL), as in Fahs et al. [33,34]. MOL has been successfully used for solving flow and/or solute transport in unsaturated porous media in many studies e.g., [10,12,21,35,36]. The ODE solver used in this work (DLSODIS: Double-precision Livermore Solver for ordinary differential equations-implicit form and sparse matrix) allows variable time-step and variable order integration with time error control [37,38].

Bayesian Parameter Inference
The vector of observations y mes is formed by the measurements of the pressure head at three different locations inside the column and the breakthrough curve of bromide at the outflow. The specific storage is fixed to a very small value (S s = 10 −8 cm −1 ). All hydraulic and solute transport soil parameters are assumed unknowns and sought by inverse modeling. In the case of linear transport model, the estimated parameters are: K S , θ S , θ r , a l , α, and n for the soil and the saturated conductivity for the porous plate (K p S ). Hence, the vector of unknown parameters is ξ = K S , K p S , θ S , θ r , a l , α, n . In the case of mobile-immobile model, the vector of unknown parameters contains, in addition, the fraction of mobile water content (f ) and the exchange coefficient (ω). Hence, in this case, it writes The soil parameters for both linear and mobile-immobile situations are estimated using a Bayesian approach where the measurements are combined with prior parameter information to determine the posterior probability distribution function (pdf). In the following, the pdfs are inferred using DREAM (ZS) software [39] based on the Markov chain Monte Carlo (MCMC) method. The MCMC method has been used by several authors in hydrogeology, e.g., [11,[40][41][42][43][44]. With MCMC, random sequences of parameter sets are generated and converge asymptotically toward the target distribution [45]. The latter can then be used to evaluate mean estimated parameter values and confidence intervals to characterize parameter uncertainty. Due to the absence of prior information on the parameter values, independent uniform prior distributions are assumed for all soil parameters. With these non-informative distributions, all possible parameter values are equally likely. Furthermore, the parameter intervals are chosen sufficiently large (Table 1) to reflect our usual poor prior knowledge of the model parameter values.
Error measurements of the tensiometric pressure sensors and of the ion selective electrode are assumed to be normally and independently distributed with unknown standard deviations. The latter are considered as hyperparameters and are simultaneously estimated with the physical parameters. We denote by σ = (σ h , σ c ) the vector of hyperparameters.
Setting the calibration problem in a Bayesian framework yields the following posterior pdf, where SS h (ξ) and SS C (ξ) are defined as the sums of the square of differences between the observed and predicted data of the pressure head and concentration, respectively.
mod ) pressure heads at time t k for the number of pressure head observations (N h ).
With the MCMC sampler, a new candidate  The last parameter sets that adequately fit the model onto observations are used to appraise the quality of parameter estimate (such as the mean and the standard deviation of the estimated parameters). The MCMC sampler is used hereafter with three parallel chains. Results are considered stationary if the criterion, defined by Gelman and Rubin [47], is verified (R stat ≤ 1.2) and if the chains are not autocorrelated. Table 1. Ranges of the soil parameters.

Parameter
Lower Bound Upper Bound

Results and Discussion
The surface injection induces an increase of the pressure head at the tensiometer locations. When the surface injection is stopped, the measured pressure heads decrease as the water front moves below the tensiometers (Figure 2). For each tensiometer, the pressure head is almost constant during the injection period. Note that since the tensiometers are separated by the same distance (15 cm between each other), the jump in the pressure levels during the injection period between the first and second tensiometers is similar to that between the second and third ones (Figure 2). The results of Figure 2 show that only the tensiometer located at Z 1 = 2.5 cm from the soil surface continuously measures negative pressure heads. Hence, it is the only tensiometer that remains under unsaturated conditions after the arrival of the water front. The other tensiometers (located at Z 2 = 17.5 cm and Z 3 = 32.5 cm from the soil surface) measure positive pressure heads indicating that the latter are under the water table. Note that after the end of the infiltration period (i.e., for t > 586 min), desaturation occurs and the pressure head at each tensiometer decreases until reaching its initial hydrostatic negative value corresponding to −35.9 cm, −20.9 cm, and −5.8 cm at, respectively, Z 1 = 2.5 cm, Z 2 = 17.5 cm, and Z 3 = 32.5 cm from the soil surface. The pressure measurements during this desaturation period ( Figure 2) are prone to contain useful information for the characterization of the unsaturated soil parameters. The measured concentration at the column outflow is given in Figure 3. This figure depicts the logarithm of the bromide concentration in the effluent since the potential measured by the ion-selective electrode is directly proportional to the logarithm of the ion concentration.
Parameter identifiability is assessed for the linear transport model using MCMC with 30,000 model runs. Figure 4 shows that the chains have converged after 14,000 model runs. The values of the mean estimated parameters and the corresponding confidence intervals, for the linear transport model, are depicted in Table 2. The results of this table show that, except the residual water content, all the parameters are well estimated since the posterior intervals of the estimated parameters ( Table 2) are strongly narrowed, as compared to their prior intervals (given in Table 1). The significant difference between the prior and posterior distributions of a parameter testifies the high model sensitivity of this parameter. Furthermore, the marginal distributions of all soil parameters (except θ r ) exhibit almost bell-shaped posterior distributions ( Figure 5).
Water 2019, 11, x FOR PEER REVIEW 7 of 16 the mean estimated parameters and the corresponding confidence intervals, for the linear transport model, are depicted in Table 2. The results of this table show that, except the residual water content, all the parameters are well estimated since the posterior intervals of the estimated parameters ( Table 2) are strongly narrowed, as compared to their prior intervals (given in Table 1). The significant difference between the prior and posterior distributions of a parameter testifies the high model sensitivity of this parameter.        The saturated water content of the soil is also well estimated with an associated uncertainty less than 5%. The residual water content has no influence on the pressure head nor on the output concentration. Indeed, the posterior interval of r θ is similar to its prior interval, which indicate that this parameter is not influent and, hence, cannot be well estimated. This non identifiability is due to the fact that dry conditions were not attained during the experiment. The results of the Figure 5 show that the parameters α and n are quite well estimated, although they are strongly correlated (

r =
). Finally, the dispersivity transport parameter l a is quite well estimated with a high mean estimated value of 2.33 cm and a large confidence interval ( Figure 5). Figures 2 and 3 depict the pressure head and the outflow concentration evaluated at the measurement points using results of the MCMC inversion. A good agreement is observed between the predicted and simulated pressure heads using the mean estimated parameter values. However, a less satisfactory matching is observed between the predicted and measured concentrations. The disagreement between concentration profiles is especially notable at the leaching phase. Indeed, the measured concentrations show an asymmetric and long-tailed breakthrough curve, which cannot be reproduced by the symmetrical simulated concentration profile. The 95% uncertainty bands corresponding to the 97.5 and 2.5 percentiles of the distributions are shown in Figures 2 and 3. In these figures, the green regions represent the total predictive uncertainty, which takes into account both parametric uncertainty and measurement errors. All pressure and concentration observations fit within the 95% total uncertainty range. Small uncertainties were obtained for pressure head The saturated conductivities of the soil K S and that of the porous plate K p S are accurately estimated and show very narrow uncertainty intervals. The sand soil is around 100 times more permeable than the porous plate. Note that K S and K p S are moderately correlated (r = −0.85). The saturated water content of the soil is also well estimated with an associated uncertainty less than 5%. The residual water content has no influence on the pressure head nor on the output concentration. Indeed, the posterior interval of θ r is similar to its prior interval, which indicate that this parameter is not influent and, hence, cannot be well estimated. This non identifiability is due to the fact that dry conditions were not attained during the experiment. The results of the Figure 5 show that the parameters α and n are quite well estimated, although they are strongly correlated (r = 0.99). Finally, the dispersivity transport parameter a l is quite well estimated with a high mean estimated value of 2.33 cm and a large confidence interval ( Figure 5). Figures 2 and 3 depict the pressure head and the outflow concentration evaluated at the measurement points using results of the MCMC inversion. A good agreement is observed between the predicted and simulated pressure heads using the mean estimated parameter values. However, a less satisfactory matching is observed between the predicted and measured concentrations. The disagreement between concentration profiles is especially notable at the leaching phase. Indeed, the measured concentrations show an asymmetric and long-tailed breakthrough curve, which cannot be reproduced by the symmetrical simulated concentration profile. The 95% uncertainty bands corresponding to the 97.5 and 2.5 percentiles of the distributions are shown in Figures 2 and 3. In these figures, the green regions represent the total predictive uncertainty, which takes into account both parametric uncertainty and measurement errors. All pressure and concentration observations fit within the 95% total uncertainty range. Small uncertainties were obtained for pressure head responses, while very large uncertainties were assigned to the concentration outputs. This phenomenon is related to the bad matching between measured and simulated concentrations. This bad matching is responsible of the large dispersivity coefficient found by the MCMC sampler in order to reproduce the long-tailed breakthrough curve of the observed concentrations.
The Bayesian inversion of the mobile-immobile model, which includes nine unknown parameters, was performed with the MCMC sampler using 40,000 model runs. The convergence of the MCMC chains requires more simulations to reach the convergence than with the linear transport model. As reported in Figure 6, the overall chains have converged after 24,000 model runs. The mean values of the estimated parameters and the associated confidence intervals are summarized in Table 3. As previously, the posterior confidence intervals of all parameters, except θ r , are much narrower than their prior intervals. The histogram of marginal distribution of θ r shows almost a uniform distribution similar to its prior one. The histograms of marginal distributions of the rest of parameters show almost bell-shaped distributions (Figure 7) significantly different from their prior uniform distributions which indicate the high sensitivity of the parameters.
The results of Table 3 and Figure 7 show that the new two parameters f and ω are well estimated. Note that the fraction of the mobile water content to the total water content f is very accurately estimated with an associated uncertainty less than 2%. A moderate correlation (r = −0.73) is observed between the parameters f and ω. The comparison between the results of the linear transport model ( Table 2) and the mobile-immobile one (Table 3) show that the two models yield almost similar mean parameter estimated values and similar confidence intervals for the saturated hydraulic conductivities (K S and K p S ) and the van-Genuchten parameters (α and n). As previously, the parameters K S and K p S are moderately correlated (r = −0.87), whereas the parameters α and n are very strongly correlated (r = 0.99). The mean estimated values of the saturated water content of the two models are almost identical. However, the associate uncertainty intervals are significantly different. Indeed, the inversion of the mobile-immobile model yields an uncertainty interval three times narrower than with the linear transport model. Hence, θ S is much more accurately estimated with the mobile-immobile transport model, although both models yield the same mean estimated parameter value. As with the linear transport model, the residual water content is irrelevant and, hence, cannot be identified. Among the parameters estimated by the inversion of the two transport models, the dispersivity coefficient shows the most significant differences both in terms of mean estimated value and in term of associate uncertainty range. Indeed, the inversion of the mobile-immobile transport model yields a mean estimated a l value that is 3.3 times less than the one obtained with the linear transport model. The confidence interval is also more than 3 times narrower with the mobile-immobile transport model.  Because the two transport models yield almost the same mean estimated values and uncertainty ranges for all hydraulic parameters ( S K , p S K , S θ , r θ , α , n ), the matching between measured and calibrated pressure heads by inverting the mobile-immobile transport model (not shown) is similar to the matching depicted in Figure 2 obtained by inverting the linear transport model. Likewise, the total predictive uncertainty associated with the pressure head output with the mobile-immobile model (not shown) is similar to the one observed with the linear transport model (grey bands in Figure 2). Measured and simulated output concentrations using the mean estimated parameter values obtained from the inversion of the mobile-immobile transport model are plotted in Figure 8. Contrary to the results of the linear transport model (Figure 3), a very good matching is observed between measured and simulated concentrations during the whole experiment (including the leaching phase). Furthermore, the total uncertainty (grey bands) characterizing the concentration output is much narrower with the mobile-immobile transport model ( Figure 8) than with the linear one ( Figure 3). Hence, for the investigated undisturbed soil, the mobile-immobile transport model better reproduces bromide infiltration and yields much more accurate output concentrations than the linear one. Because the two transport models yield almost the same mean estimated values and uncertainty ranges for all hydraulic parameters (K S ,K p S ,θ S ,θ r ,α,n), the matching between measured and calibrated pressure heads by inverting the mobile-immobile transport model (not shown) is similar to the matching depicted in Figure 2 obtained by inverting the linear transport model. Likewise, the total predictive uncertainty associated with the pressure head output with the mobile-immobile model (not shown) is similar to the one observed with the linear transport model (grey bands in Figure 2).
Measured and simulated output concentrations using the mean estimated parameter values obtained from the inversion of the mobile-immobile transport model are plotted in Figure 8. Contrary to the results of the linear transport model (Figure 3), a very good matching is observed between measured and simulated concentrations during the whole experiment (including the leaching phase). Furthermore, the total uncertainty (grey bands) characterizing the concentration output is much narrower with the mobile-immobile transport model ( Figure 8) than with the linear one ( Figure 3). Hence, for the investigated undisturbed soil, the mobile-immobile transport model better reproduces bromide infiltration and yields much more accurate output concentrations than the linear one.

Conclusions
In this study, a laboratory infiltration experiment was conducted by injecting bromide solution at the surface of an undisturbed soil column. The local pressure head measurements inside the column and the outflow concentration measurements were used to simultaneously estimate all hydraulic and transport soil parameters. Unsaturated water flow is simulated using the Richards equation coupled with the Mualem/van Genuchten models. The linear advection-dispersion and the non-equilibrium mobile-immobile models were investigated for modeling the bromide propagation. A Bayesian inversion was performed using the pressure head and concentration measurements. In the case of linear transport model, seven soil parameters were estimated: The saturated conductivity S K , the saturated water content S  , the residual water content r  , the dispersivity coefficient l a , and the van Genuchten parameters (  and n ) of the soil and the saturated conductivity p S K of the porous plate. In the case of mobile-immobile model, the vector of unknown parameters contains, in addition, the fraction of mobile water content f and the exchange coefficient  . All hydraulic and transport soil parameters were estimated in a one-step procedure using the MCMC sampler starting from non-informative uniform distributions for all parameters. The obtained results are summarized hereafter: - The flow through the investigated soil can be well reproduced by the Richards equation combined with the Mualem/van-Genuchten models. - The pressure head and outflow concentration measurements are able to identify all the unsaturated hydraulic and solute transport soil parameters, except the residual water content r  which was not identified with either of the models since dry conditions were not attained during the laboratory experiment. -A significant difference was observed between the two models for the estimation of the dispersivity coefficient. Indeed, with the mobile-immobile transport model, the estimated value of l a was three times lower than the one obtained with the linear transport model. The associated uncertainty was alpressure head and outflow concentration measurements are able to identify all the unsaturated hydraulic and solute transport soil parameters, except the residual water content

Conclusions
In this study, a laboratory infiltration experiment was conducted by injecting bromide solution at the surface of an undisturbed soil column. The local pressure head measurements inside the column and the outflow concentration measurements were used to simultaneously estimate all hydraulic and transport soil parameters. Unsaturated water flow is simulated using the Richards equation coupled with the Mualem/van Genuchten models. The linear advection-dispersion and the non-equilibrium mobile-immobile models were investigated for modeling the bromide propagation. A Bayesian inversion was performed using the pressure head and concentration measurements. In the case of linear transport model, seven soil parameters were estimated: The saturated conductivity K S , the saturated water content θ S , the residual water content θ r , the dispersivity coefficient a l , and the van Genuchten parameters (α and n) of the soil and the saturated conductivity K p S of the porous plate. In the case of mobile-immobile model, the vector of unknown parameters contains, in addition, the fraction of mobile water content f and the exchange coefficient ω. All hydraulic and transport soil parameters were estimated in a one-step procedure using the MCMC sampler starting from non-informative uniform distributions for all parameters. The obtained results are summarized hereafter: - The flow through the investigated soil can be well reproduced by the Richards equation combined with the Mualem/van-Genuchten models. - The pressure head and outflow concentration measurements are able to identify all the unsaturated hydraulic and solute transport soil parameters, except the residual water content θ r which was not identified with either of the models since dry conditions were not attained during the laboratory experiment. -significant difference was observed between the two models for the estimation of the dispersivity coefficient. Indeed, with the mobile-immobile transport model, the estimated value of a l was three times lower than the one obtained with the linear transport model. The associated uncertainty was also three times smaller with the mobile-immobile transport model.

-
The linear transport model yielded a good agreement between the measured and simulated pressure heads, but a less satisfactory matching was observed between measured and simulated concentrations. Furthermore, small uncertainties were obtained for pressure head responses, while very large uncertainties were assigned to the output concentrations. - The mobile-immobile transport model better reproduced the infiltration experiment. Indeed, a very good matching was obtained between measured and simulated concentrations. Furthermore, the concentration uncertainty region was much narrower than with the linear transport model.
Finally, the results of this work point out the necessity of testing different models for interpreting real experiments and to employ accurate Bayesian methods for model calibration. These methods provide accurate information concerning the estimated parameters, their confidence intervals, and model output uncertainties.
Author Contributions: J.Z. and S.K. worked on the laboratory experiment; F.L. worked on the numerical simulations; A.Y. and M.F. supervised the research and finalized the manuscript.
Funding: This research was partially supported by the French National Research Agency (ANR Hydrocrizsto) and the Tunisian-French joint international laboratory NAILA (http://www.lmi-naila.com/).