Each reference sample was composed of a mixture of 180 g of glass beads (0.250 mm <x<0.500 mm), 20 g of pure air dry Portland cement and 35 g of tap water. Once homogenized, the mixture was filled into a 100 cm3 (5 cm height/5 cm diameter) stainless steel ring by frequently tapping it on a table to ensure that it was uniformly packed. The ring was closed at the bottom by a lid. Any excess material on the top was removed with a spatula. Each sample was allowed to cure for 72 h at room temperature. The bottom lid was subsequently replaced with a Eijkelkamp nylon cloth supported by a rubber band. The sample, with the ring, the cloth and the rubber, was weighted. The empty ring, the cloth and the rubber were previously weighted separately.

The ring test was organized into 14 soil physics laboratories. An example of a reference sample was sent to each laboratory alongside the material needed to construct five other samples. A total of 84 reference samples were constructed. The ring test consisted of three successive rounds of SWRC measurements. At the beginning of each round, each sample was initially saturated for 48 h. The mass of each sample was then measured after equilibration at different matric potential (or suction) values: ψ=10 hPa (log⁡10(ψ)=1.00), ψ=50 hPa (log⁡10(ψ)=1.70), ψ=100 hPa (log⁡10(ψ)=2.00) and ψ=300 hPa (log⁡10(ψ)=2.48). Equilibration times were 5 d at 10 hPa, 7 d at 50 hPa, 10 d at 100 hPa and 15 d at 300 hPa. Finally, samples were weighed after drying for 72 h in an oven drying at 60 ∘C. Gravimetric water content (wc in g g-1) was calculated by the ratio of water masses over dry masses (water content=fresh mass-dry massdry mass). The six samples from each laboratory were divided into three exchange modalities (Fig. ); two samples were kept by the same laboratory throughout the three rounds of measurements (“stay”), two samples were sent to different laboratories between rounds (“move”) and the last two samples were sent to a different laboratory for the second round but were sent back to the original laboratory for the third round (“back”). This scheme was designed to estimate intralaboratory and interlaboratory variability as well as the effect of sample transfer between laboratories.

The final dataset consisted of 250 SWRCs. Two curves were missing. Since each SWRC was composed of four successive measurement points whose relative value may depend on the previous point, the data were not independent. Statistical analyses were then based on parameter values of fitted functions. To model our dataset, a linear function with log⁡10(ψ) as the independent variable was adjusted to the measured wet part of SWRCs (Eq. ). 1wci=β0+(β1⋅(log⁡10(ψ)i-1))+ϵi Water content (wci) was linearly expressed as a function of log⁡10(ψ) values (log⁡10(ψ)i=1.00, 1.70, 2.00, 2.48). The log⁡10(ψ)i-1 was intended to set the first retention point at the intercept. The index i represented the ith data point. β0 and β1 represented the mean intercept and the slope over all data. The term ϵi represented the residuals. The next step was to adjust a single linear model to each SWRC (Eq. ). 2wcin=β0+z0n︸varying intercept+(β1+z1n︸varying slope)⋅(log⁡10(ψ)i-1)+ϵin The n index represented the nth SWRC. Depending on the modelled SWRC, intercepts (z0n) and slopes (z1n) were allowed to vary around a general intercept (β0) and slope (β1). This type of model refers to a linear mixed (effect) model.

The purpose of this study was also to investigate the interlaboratory variability as well as the differences between samples. Another linear mixed (effects) model was used to consider the by-sample and by-laboratory variability using adjustment terms called “random effects”. The first random effect, u0j and u1j, respectively adjusted β0 and β1 depending on the analysing laboratory j (j ∈ [1, …, 14]). The other random effect, v0k and v1k, respectively adjusted β0 and β1 depending on the sample k (k ∈ [1, …, 84]). This mixed-effect model was described by Eq. (). 3wcijk=β0+u0j+v0k︸varying intercept+(β1+u1j+v1k︸varying slope)⋅(log⁡10(ψ)i-1)+ϵijk

Finally, the effect of sample changes between round 1 and round 3 on the intercept (w0) and the slope (w1) was modelled through a “fixed-effect” covariate. The covariate depended on a dummy variable associated with the round number; for the first round, roundi=-0.5 and for the third round, roundi=0.5. This later model (Eq. ) was applied only to data associated with the back samples and stay samples to avoid laboratory effects. The results were compared to determine whether the differences in measurements between rounds 1 and 3 were due to transport or differences caused by wear of the samples not related to transport. 4wcik=β0+v0k+w0⋅roundi︸varying intercept+(β1+v1k+w1⋅roundi︸varying slope)⋅(log⁡10(ψ)i-1)+ϵik All parameters from each models were estimated using Bayesian statistics. Posterior distributions were sampled with a Markov chain–Monte Carlo (MCMC) algorithm implemented in C++ through an R package called “RStan” . Priors are noninformative, i.e. centred normal distributions for random effects parameters and uniform distributions for general intercept, general slope, fixed effect covariates and variance parameters. Sensitivity analyses of priors and validations of models were also carried out. Inference was based on Bayes factors and Bayesian credible intervals of posterior distributions. More details are available in the Supplement.

Each laboratory was given the same procedure to measure the SWRC. However, it allowed some freedom, and some laboratories did not perfectly implement it. For instance, laboratory 8 dried the samples at 100 ∘C instead of 60 ∘C. Hence, laboratories used slightly different procedures to those shown in Table . Laboratories mostly used the sand box (SB) at 10, 50 and 100 hPa. At 300 hPa, the suction plate (SP) and the pressure plate (PP) were dominating. Lab 7 was the only one to use the sand/kaolinite box (SKB).

The simple linear regression (Eq. ) with the log⁡10(ψ) as predictor was used to model the dataset. The posterior probability distribution of the general intercept (β0), the slope (β1), and the standard deviation of the residuals (σϵ) are shown in Fig. . The mean value of σϵ was fairly high (0.0126 g g-1). Indeed, as shown in Fig. , the variability of measured SWRCs was large (spreading of the curves). The following steps were devoted to explaining the origin of this variability.

The next step was to model a single linear regression for each SWRC (Eq. ). The posterior probability distribution of the general intercept (β0), the slope (β1), the individual intercept (z0n) and slope (z1n), and the standard deviation of the residuals (σϵ) are shown in Fig. . The intercept (z0n) and slope (z1n) parameters are different for each individual SWRC. These individual parameters explain the variability that exists between all SWRC. Hence, the mean value of σϵ presented in Fig.  decreased by approximately 60 % compared to the previous model (Fig. ) and now only represents a fitting error introduced by the choice of modelling SWRCs by linear regressions.

From these results, one can determine the standard deviation of z0n and z1n for each of the stay samples (between the three rounds). As each laboratory measured two stay samples, an estimate of the intralaboratory variability of each laboratory can be made by pooling the density estimates of the standard deviation of the two samples (Fig. ). As the intralaboratory variability associated with each of the two stay sample is determined separately and then merged together, this allows us to minimize the effect of possible variations between the two stay samples. Intralaboratory variability is therefore defined as the variability between retention curves, modelled by linear regressions, measured on a similar (stay) sample within a same laboratory that uses a given measurement procedure.

The estimate of the mean intralaboratory standard deviation of all laboratories pooled together (Fig.  bottom row) was 0.00533 g g-1 (95 % credible interval (CrI) 0.00018–0.01138 g g-1) for the intercept and 0.00519 g g-1log⁡10 (hPa)-1 (95 % CrI 0.00038–0.01068 g g-1log⁡10 (hPa)-1) for the slope (Table  in the Appendix). Figure  (top row) also shows that the intralaboratory variability was quite different depending on the laboratory. Some laboratories succeeded in repeating similar SWRCs results on a same sample, while others failed.

Although all laboratories were given the same procedure to build the reference samples, the conditions under which they were constructed differed between laboratories. Hence, the bulk density of samples at the beginning of the experiment was variable depending on the laboratory that constructed the sample (Table ). Indeed, the difference between the mean bulk density of samples constructed by the lab 1 (highest bulk density) and lab 14 (lowest bulk density) was 0.1573 g cm-3. Hence, the later linear mixed (effect) model was used to estimate the interlaboratory variability on the SWRC considering the differences between samples (Eq. ). Densities of the posterior probability distribution of the general intercept (β0) and slope (β1), the random effect of laboratory on the intercept (u0j) and slope (u1j), the random effect of sample on the intercept (v0k) and slope (v1k) and the standard deviation of the residuals (σϵ) are shown in Fig. . The mean value of σϵ presented in Fig.  decreased by approximately 40 % compared to the simple linear model (Fig. 1). Indeed, a part of the variability between SWRCs has been explained by sample and laboratory random effects. Parameter values of the laboratory random effects (u0j and u1j) show how SWRCs systematically deviate depending on the analysing laboratory. Differences between samples were also estimated with parameter values of the samples random effects (v0j and v1j). The wider dispersion of the laboratory random effect parameters indicates that the analysing laboratory explained a larger proportion of the overall variance than the analysed sample. Indeed, on the intercept, the mean laboratory random effect standard deviation (σu0) was 0.00872 g g-1, while it was 0.00350 g g-1 for the sample random effect (σv0). The same observation applies to the slope with a mean standard deviation of 0.00602 g g-1log⁡10 (hPa)-1 for the laboratory random effect (σu1) and 0.00451 g g-1log⁡10 (hPa)-1 for the sample random effect (σv1). The mean laboratory random effect standard deviations on the intercept and slope values (σu0 and σu1) represent an estimation of the interlaboratory variability.

It should be noted that results presented here are only representative of the reference samples that were measured in this particular case. Thus, estimates of interlaboratory (and intralaboratory) variability values are not directly transferable to other samples of a different nature with different retention characteristics.

In order to assess the effect of possible sample changes on the SWRC measurements, a last model was separately fitted to the data from back and stay samples (Eq. ). The Bayes factor indicated that the predicted data of back samples were 46.60 times more probable under the model that takes the round effect into account than the model without the round effect. Moreover, the 95 % credible interval of the posterior probability distribution of the “round” effect (Fig. ) laid outside 0 for the intercept (w0b) and the slope (w1b). Therefore, back samples changed between round 1 and round 3.

On the other hand, the Bayes factor indicated that the predicted data of stay samples were 4348 times less probable under the model that takes the round effect into account than the model without the round effect. This suggests that the changes in the retention properties of the stay samples between the first and third rounds were negligible. In contrast, the 95 % credible interval of the posterior probability distribution of the round effect (Fig. ) suggests that the stay samples changed between round 1 and round 3 as it laid outside 0 for the intercept (w0s) and slope (w1s). Thus, for the stay samples, the Bayes factor and 95 % credible interval yielded two opposite conclusions. Also, the dry mass of stay samples increased between round 2 and 3 (p=0.016*), which was not the case with back samples (p=0.199 ns). The round effect on stay samples is discussed in Sect. .

In total, 57 of the 250 measured SWRCs showed an increase in water content between at least two increasing suction steps (Fig. ). Whatever the origin, this increase in water content is physically impossible. It appeared that the occurrence of these anomalies depended on the analysing laboratory, with some having no anomalies and others having a large number of occurrences. Indeed, laboratories 3, 11 and 14 together accounted for 35 of the 57 anomalies recorded (Table ). Moreover, this anomaly has happened more than once for some samples such as for the samples 1, 2, 11, 15, 56, 61, 62, 63, 65 and 66. The samples 15, 63, 65 and 66 showed this anomaly in two different laboratories. Also, it occurred more often the drier the sample was. There were 55 occurrences between 100 and 300 hPa, while there were only 9 between 50 and 100 hPa and 3 between 10 and 50 hPa. For some SWRCs, there was more than one occurrence.

This study confirms that there are systematic differences in the measurement of SWRCs depending on the laboratory (Fig. ). This is true even for laboratories using similar devices (e.g. lab 6 vs 9). These systematic differences in the measurement of SWRCs attributed to laboratories resulted in a large interlaboratory variability. The portion of variability attributed to differences between laboratories was larger than the portion of variability attributed to intrinsic differences between samples. This is concerning since it was shown, through the comparison of the bulk densities (Table ), that the samples were different even at the very beginning of the experiment. From saturation to drying, all laboratories used slightly different procedures (Table ). Similarly to the argument of on the measurement of macropore volume, total pore volume and saturated hydraulic conductivity, we believe that procedural differences between laboratories can be at the origin of this interlaboratory variability. The identification of the aspects of the procedures that influence SWRC measurements is challenging since these were not studied in isolation. This is a multidimensional problem that remains beyond the scope of this article. Nonetheless, an attempt is made to hypothesize the effect of some of these procedural aspects. It should also be mentioned that the true value of the SWRC was unknown. Laboratories were compared according to their relative position with respect to the others and not against a fixed target value.

Differences between laboratories are unlikely to be associated with differences in the analysed samples. The intrinsic differences between samples were considered by the model (Eq. ) using the sample random effect. Indeed, there was no correlation between the intercept parameter of each laboratory and the average bulk density of the samples analysed by each laboratory (r=-0.0078).

A first possible source of this interlaboratory variability could be attributed to the different devices used. To our knowledge, no study has yet attempted to compare SWRCs obtained with SB, SKB, SP or PP. Nonetheless, found that SWRCs measured with SP were less reproducible (wider spread) than those measured with the evaporation method in the ψ range of 0–300 hPa. This wider spread could be associated with contact issues with the plate and/or with the smaller sample size used with the SP method which could be smaller than the representative elementary volume of the soil. They also found that, for sandy soils, water contents are systematically smaller for SWRCs obtained with SP than with the evaporation method. These differences could be due to smaller initial saturation degrees of the sandy samples measured with the SP compared to those measured with the evaporation method. Temperature effects and dynamic non-equilibrium effects as found by could also have played a role in these differences. also compared macropore volume (as equal to the difference between the total pore volume or water content at saturation and water content at a matric potential of 60 hPa) obtained with SB and with pressure cells on real soil samples. It can be observed that the volume of macropores determined by laboratories using the SB is generally greater with a greater variability than when determined by laboratories using the pressure cell. However, it is not clear from their results whether this difference is due to the method of obtaining the water content at 60 hPa (SB or pressure cell) or to the method of estimating the total pore volume (from bulk density or weighting at saturation, etc.) which was also different between laboratories.

Moreover, all laboratories used two different devices between 10 and 300 hPa, except labs 5, 10 and 12 that kept the same device for each pressure step. Switching from a suction to pressure system may affect the measurement of the SWRC. In a suction device, the suction is applied by a hanging water column via a continuum of water. In a pressure device, the pressure is applied via the air to the soil and plate water, while the plate bottom is at atmospheric air pressure. So, the propagation of the applied tension/pressure might be different.

The procedures for the dry mass measurement of the samples may also have played a role in the observed differences. Indeed, the estimation of the intercept parameter of laboratory 8, which dried the samples at 100 ∘C, was higher than the ones from the other laboratories, which dried the samples at 60 ∘C (Fig. ). This suggests that the dry masses measured by laboratory 8 were lower than those measured by the other laboratories. This may be explained by the fact that the assumed “zero” water content corresponds to two different water content values associated with two different water potentials. Indeed, the equilibrium potential of the water contained in a sample after drying depends on the drying temperature and the relative humidity of the air inside the oven . At liquid–vapour equilibrium and at constant relative humidity, the same sample dried at 60 ∘C will have a higher water content than if it was dried at 100 ∘C. This difference in water content will also depend directly on the shape of the dry part of the retention curve. Furthermore, for the same sample at the same initial condition, the drying time required to remove the same amount of water increases as the drying temperature decreases. If liquid–vapour equilibrium has not been reached after the prescribed time of drying, it is possible that the amount of water released is lower when drying at 60 ∘C than at 100 ∘C.

Another possibility to explain differences is the way laboratories maintained hydraulic contact between the draining porous media and the sample, enabling water to be released from the sample until hydrostatic equilibrium is reached. When the draining porous media were rigid (e.g. ceramic), some laboratories used a “contact material” to improve the hydraulic contact (Table ). From this study, it cannot be concluded that the use of contact materials did or did not improve the hydraulic contact between the samples and the porous plates. The results of laboratory 11 suggest that contact issues may have occurred even if filter paper was used as the contact material.

Nevertheless, the use of contact materials may sometimes be useful when considering laboratories using the same devices. Hence, it appears that the use of filter paper by laboratory 5 resulted in more water being released (more negative slope) than laboratory 10, which did not use any contact material but used the same devices (Fig. ). also found that filter paper allowed more water to be released than polyester fabric and synthetic knitwear at 5000 and 15 000 hPa with the PP. The use of kaolinite by lab 6 and loamy soil by lab 9 as contact material seems to cause more water to be released between 100 and 300 hPa than laboratories 3 and 14 that did not used any contact material but used the same devices (results not shown). However, when looking at the entire domain of suction (from 10 to 300 hPa), the effect of kaolinite or loamy soil was negligible (Fig. ). also found kaolinite ineffective in speeding equilibrium (or increasing hydraulic conductance), with inconsistent effects, at 15 000 hPa. Further work should be done to determine which contact materials are useful depending on the specific situation.

An option to check if the hydrostatic equilibrium is achieved is to connect the porous drainage medium to a graduated cylinder and monitor the water the amount of water drained from the sample. Once no more water flowing out of the sample is observed, hydrostatic equilibrium is considered to be achieved. This setup has to be sealed in order to ensure that there is no evaporation. The advantage of such a system is that one does not need to assume the equilibration time a priori. To our knowledge, this setup was used by laboratory number 8. However, it is still possible, with this setup, that the hydraulic contact is broken and the flow of water is stopped before hydrostatic equilibrium is reached, which is referred to as an apparent hydrostatic equilibrium.

It should also be mentioned that with devices using a hanging water column as suction regulation system, the applied suction is usually expressed in centimetres of water column. Units of hectopascals (hPa) and centimetres are commonly considered equivalent, but in fact 1 cm of vertical water column corresponds to 0.98 hPa. This is usually overlooked when units are transformed (see Table ). This bias may constitute a small part of the variability between laboratories and calls for harmonization of units.

In addition, the reference level compared to the sample at which the suction is applied varies between laboratories and devices used. Some laboratories applied the prescribed suction to the bottom of the samples, while others applied it to the middle (see Table ). Laboratories that applied suction to the bottom of the samples systematically applied 2.5 cm more suction than those that applied it to the middle. This difference could have easily been corrected for after the measurement. However, in practice, this information is never transmitted from the laboratory to the end user. Hence, when pooling data from different laboratories using different reference levels to apply suction, this unknown difference introduces variability.

There might be other procedural aspects that could be responsible for these differences between laboratories (saturation procedure, porous plate maintenance during the experiment, means of preventing air leakages and evaporation, maintenance of the ceramics and the sandboxes, weighting procedure, or maintenance of the scales, etc.) (Table ). Big errors can be avoided by a quality check of the results. To our knowledge, only one laboratory used a reference sample to control the quality of their SWRC measurements as a standard operating procedure.

Some laboratories successfully reproduced SWRCs of the same (stay) sample, while others failed (Fig. ). Nevertheless, the estimate of the mean intralaboratory variability over all laboratories was smaller than the mean interlaboratory variability but was more uncertain since it was drawn from fewer samples and since the intralaboratory variability was quite different between laboratories. Obviously, this variability can partly be attributed to the different methods and procedures that existed between laboratories that were discussed above. Some procedures ensured fairly good repeatability of results, while others did not.

The two laboratories with the greatest intralaboratory variability on the slope (see Figs.  and : light pink and dark purple curves) were also among those with the most anomalies (see Table : lab 11 and 14). Concerning laboratory 8 (see Fig. : cyan curve), the bimodal shape of the intralaboratory variability on the intercept shows that for one sample the variability was high, while it was low for the second. This bimodality clearly indicates that the estimation of individual intralaboratory variabilities is rather uncertain as they are only based on two samples measured with only three repetitions. This calls for further trials on reference samples to obtain a more reliable estimate of intralaboratory variability. Nevertheless, this provides an insight into the way forward to improve data quality management in soil physics laboratories.

It appears that there was a slight effect of the transfer between laboratories on the back samples. The values of w0b and w1b indicate that SWRCs of back samples globally have a smaller intercept with a flatter slope after being transported (Fig. ). This pattern might indicate a shift to more small pores. A possible explanation for these changes in porosity is the calcium carbonation of the cement. This reaction (Ca(OH)2+CO2atm⇌CaCO3+H2O) forms CaCO3 precipitates inside the pore network, inducing a shift of the pore size distribution towards smaller pores, a decrease of the total porosity, pore clogging, and a loss of pore connectivity in cement-based materials . This hypothesis is also motivated by the fact that the dry masses of the samples increased significantly between rounds 2 and 3 for the samples stay and the dry masses did not decrease significantly for the samples back even if losses of materials were reported by the laboratories. Indeed, estimated that the carbonation-induced increase in bulk density (due to CO2 fixation) from a non-carbonated to a fully carbonated cement paste was 1.60 to 2.03 g cm3. However, the actual contribution of this phenomenon to changes in the retention properties of each reference sample is difficult to estimate, as the degree of carbonation was influenced by environmental factors (CO2 concentration, air humidity, or water content of the cement, etc.) which have not been controlled.

Nevertheless, this significant transport effect could have led to an overestimation of the interlaboratory variability, as a part of the variability of the SWRC measurements can be attributed to sample changes. The use of cement to construct such reference samples should certainly be avoided in the future.

Although it was not significant in general for stay samples, some laboratories still seem to report sample changes between rounds. The changes followed the same patterns as for back samples, which were significant it that case (Fig. ). Indeed, for stay samples of laboratories 2, 3, 4, 6, 9, 10, 11, 12 and 13, it seems that the water content at 10 hPa (log⁡10(ψ) = 1) systematically decreased with time. This was particularly visible for lab 9 (Fig. ). The origin of the changes of the stay samples can partly be attributed to the same origin as for the back samples. The degree of changes may therefore be influenced by the way the samples were handled and stored, resulting in less wear for non-transported samples. Nevertheless, the wear of the stay samples implies that the estimation of the intralaboratory variability was certainly inflated as it included the variability attributed to sample changes between rounds. It should be mentioned that for some laboratories with the highest intralaboratory variability (laboratories 8 and 14), this trend was not visible, indicating that for these laboratories the variability attributed to procedures was probably more important than the variability attributed by sample changes.

Many reasons might be elicited to explain the fact that some SWRCs showed an increase in water content between at least two increasing suction steps (Fig. ). Obviously, this happened depending on a combination of reasons related to the laboratory but also to the sample being analysed (Table ).

A possible reason can be the lack of hydraulic contact between the draining porous media and the sample, preventing water from being released in time from the sample. This is supported by the higher frequency of outliers when the sample was drier, as hydraulic conductivity decreases as the sample dries . Indeed, there is a possible scenario in which samples may absorb water but may not be able to release it according to the driving (higher) pressure. Measurements in a pressure chamber typically involve placing samples on pre-wetted ceramic plates. However, especially when a wet contact material is used, an unsaturated sample may start absorbing water (from the plate and the contact material) and re-saturate before the chamber is pressurized. Once the chamber is pressurized the excess of water may not be drained if the hydraulic contact is not well established. Hydraulic contact could have been hampered by the rigid nature and non-flat bottom topography of reference samples which did not fit the porous plate or by the use of shrinkable contact materials. When using the pressure plate, it is also possible that a “backflow” of water from the ceramic to the sample may occur between the release of the pressurized air and the disconnection of the sample from the plate . Nevertheless, increasing SWRC also occurred with sandboxes and sand/kaolinite boxes, where the applied suction was not released when the sample was disconnected.

The results presented in Sect. 3.3 and 3.4 show that interlaboratory and intralaboratory variability exists in the measurement of the SWRC. As discussed in Sect. 4.1, we suspect that some part of this variability can be attributed to the different methods and procedures used by each laboratory. Therefore, this variability could potentially be reduced by improving procedures and methods. Ideally, these should be adapted in such a manner that they allow the closest possible estimation of the “actual” SWRC of any soil sample. The prerequisite for this to be feasible is to have a fixed point of comparison, i.e. a reference sample whose retention properties are well known and remain relatively stable over time. However, these two requirements were not fulfilled in this study. Therefore, the way in which procedures and methods should be adapted remains an open question. Nevertheless, based on the results of this study, we can suggest some general future directions to further improve the quality of the laboratory SWRC measurements.

We believe that the reproducibility of SWRC measurements within a same laboratory would be improved if a reference sample was used by each laboratory as an internal quality control. This implies that the SWRC of the reference sample must remain relatively stable over time, but its true value should not necessarily be known a priori. The reference samples used in this study are already used by one laboratory as an internal quality control. According to their experience, the samples remain fairly stable for about 10 measurements as long as they are not oven-dried.

Also, an important element that would allow us to move forward would be to have a point of comparison. Thus, the emphasis should be placed on the development of reference samples for interlaboratory comparisons. Internal trials are underway to build reference samples with clays or sintered glass beads. Also, reference samples based on parallel capillary bundles with adjustable diameters could allow a “theoretical” reference retention curve to be calculated from Jurin's law. The feasibility of developing such a reference sample is also being explored.

The interlaboratory and intralaboratory variability could be reduced by improving and standardizing procedures and harmonizing methods and data. Ideally, all laboratories should endorse a unique standard operational procedure for the same method, and methods should be harmonized between each other. Improving and standardizing procedures requires a full assessment of the effect of each step of the procedures on the final SWRC measurement on a reference sample that fulfils the above-mentioned conditions (a priori known and relatively stable). Harmonization of methods can be achieved with interlaboratory comparisons on the same types of reference samples.

Finally, since procedures and methods could have an impact on the final measurement of the SWRC, the transparency of the procedures and methods used in SWRC datasets should be ensured.