the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The QuantiSlakeTest, dynamic weighting of soil under water to measure soil structural stability
Frédéric M. Vanwindekens
Brieuc F. Hardy
Abstract. We evaluated the performance of a new, simple test to evaluate soil structural stability. The QuantiSlakeTest (QST) consists in a quantitative approach of the slake test, a dynamic weighting of a dried structured soil sample once immersed in water. The objective of this work was threefold: we aimed to (i) derive indicators from QST curves to evaluate soil structural stability regarding the underlying mechanisms of soil disaggregation; (ii) establish the relationship between soil properties and QST indicators; and (iii) assess how QST indicators respond to contrasting soil management practices. To meet these goals, we sampled the soil of 35 plots from three long-term field trials in the silt loam region of Belgium dealing with contrasting organic matter inputs, tillage treatments and P-K fertilisation, respectively. For each plot, QST curves were compared to the three tests of Le Bissonnais, targeting specific mechanisms of soil disaggregation.
Shortly after immersion in water, soil mass increases due to the rapid replacement of air by water in soil porosity. Then soil mass reaches a maximum before decreasing, once mass loss by disaggregation exceeds mass gain by air loss. Our results confirmed that the early mass loss under water is mainly related to slaking, whereas after a longer time period, clay dispersion becomes the dominant process of soil disaggregation. The overall soil structural stability was positively correlated to the soil organic carbon (SOC) content and negatively correlated to the clay content of soil. Accordingly, the SOC:clay ratio was closely related to QST indicators. Nevertheless, for a similar carbon (C) input, green manure and crop residues were more efficient in decreasing clay dispersivity whereas farmyard manure promoted SOC storage and was more efficient against slaking. QST curves had a strong discriminating power between reduced tillage and ploughing regardless of the indicator, as reduced tillage increases both total SOC content and root biomass in the topsoil.
The QST has several advantages. It is rapid to run, doesn't require expensive equipment or consumables and provides a high density of information on both specific mechanisms of soil disaggregation and the overall soil structural stability. As an open access program for QST data management is currently under development, the test has a strong potential for adoption by a widespread community of end-users.
Frédéric M. Vanwindekens and Brieuc F. Hardy
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2022-1092', Nicolas P.A. Saby, 24 Nov 2022
The manuscript presents a very interesting study on a new laboratory method to measure soil structural stability. The subject is well introduced, the statistical analysis is thoroughly executed using a nice and robust dataset.
The manuscript is very well written, well organised, interesting and very clear. There is also a very interesting discussion on the mechanisms controlling soil disaggregation.
the paper addresses a relevant scientific question within the scope of SOIL by providing new soil health indicators determined by a new and simple method. The results are compared with results obtained using a reference method (MWD).
Moreover, the authors paid particular attention to the publication of FAIR data and codes.
I have very little comments
In the summary, l8, it is not clear that you use some indicators based on the curves to the derive information on soil structural stability. Maybe you could add one sentence to specify this point.
L89, ‘in field conditions’, it is not clear to me. SLAKES mobile application should be use in a lab or a room.
L100 there is a repetition of the definition of the QST acronym?
Is Table 1 published on a data repository? If yes, add the DOI in the caption and the text.
L230-240, the choice of the different indicators could appear arbitrary and in particular for the time (ii) or the threshold value (iii). I don’t think it is discussed. More generally, the use of the QST curves could be discussed in terms of statistics. You may use functional data analysis or curve modelling instead. Fajardo et al. use a Gompertz model and interpret the parameters of the curves. Maybe this point could be mentioned for future research. I think there is only one shourt sentence in the conclusion (l495).
Figure 4 and 5 I think we are missing some results on the MWD here to be able to compare.
L286 and Fig 3, yes there is a clear relationship but also there is a residual variability. I think this should be acknowledged in the text.
L344 some QST indicators. It is possible to be more specific?
L386. Add reference to figure 3 for clarity?
L410: something wrong with the reference.
L450 the use of a large soil volume facilitates the soil sample preparation compare to MWD.
L500 this sentence is right only for soils of central Belgium for the moment.
Citation: https://doi.org/10.5194/egusphere-2022-1092-RC1 -
AC1: 'Reply on RC1', Frédéric Vanwindekens, 19 Dec 2022
[N.B. We copied/pasted the comments and questions and added our answers after each points - main questions of referee are in bold in the text below]
The manuscript presents a very interesting study on a new laboratory method to measure soil structural stability. The subject is well introduced, the statistical analysis is thoroughly executed using a nice and robust dataset.
The manuscript is very well written, well organised, interesting and very clear. There is also a very interesting discussion on the mechanisms controlling soil disaggregation.
the paper addresses a relevant scientific question within the scope of SOIL by providing new soil health indicators determined by a new and simple method. The results are compared with results obtained using a reference method (MWD).
Moreover, the authors paid particular attention to the publication of FAIR data and codes.
• We appreciate the kind feedback of Mr Saby and we are very grateful for his time and helpful comments.
I have very little comments
• In the summary, l8, it is not clear that you use some indicators based on the curves to the derive information on soil structural stability. Maybe you could add one sentence to specify this point.
◦ The term "indicators" was mentioned earlier in the paragraph, in the objectives of the paper, but we added details also in l8, for more clearly point that these are calculated from the curve at some typical dynamic step:
◦ "For each plot, QST indicators derived from curves (e.g. relative increase or decrease in soil mass, disaggregation speed, time to reach a threshold value of mass loss …) were compared to the results of the three tests of Le Bissonnais, targeting specific mechanisms of soil disaggregation."
• [X] L89, ‘in field conditions’, it is not clear to me. SLAKES mobile application should be use in a lab or a room.
◦ The strength of SLAKES application relies on the rapid acquisition of data on soil structure with little equipment needed (the measurement can be performed inside but also outside in absence of rain or wind) à we propose to remove “in field conditions”:
◦ “Recently, the SLAKES mobile application provided encouraging results as a tool for rapid data acquisition on soil structure. The test relies on image recognition to measure…”
• [X] L100 there is a repetition of the definition of the QST acronym?
◦ Yes, a mistake, we adapted:
◦ "In this work, we evaluated the performance of a new, simple test to evaluate soil structural stability, named QuantiSlakeTest (QST), a quantitative approach of the slake test."
• [X] Is Table 1 published on a data repository? If yes, add the DOI in the caption and the text.
◦ The data were available via the "gitlab" repository (https://frdvnw.gitlab.io/qst-openscience/). But not in a real "data repository". We added this table in the data repository linked to the "SlakingLab - users, data and codes for running QuantiSlakeTests" community on Zenodo. Reference and doi : Vanwindekens, Frédéric M., & Hardy, Brieuc F. (2022). QST - open data of the article Vanwindekens & Hardy (2022) - table 1 soil properties (1.0) [Data set]. Zenodo. https://doi.org/10.5281/zenodo.7405113
◦ As suggested, we will add the doi link in the text, just after the ref to the table and in the caption of the table 1.
• [X] L230-240, the choice of the different indicators could appear arbitrary and in particular for the time (ii) or the threshold value (iii). I don’t think it is discussed. More generally, the use of the QST curves could be discussed in terms of statistics. You may use functional data analysis or curve modelling instead. Fajardo et al. use a Gompertz model and interpret the parameters of the curves. Maybe this point could be mentioned for future research. I think there is only one short sentence in the conclusion (l495).
◦ This is a very good point, Thank you for this comment. Actually, curve modelling has been tested (we used double exponential models) but the increase in soil mass occurring at the start of the test makes modelling complicated and model parameters correlated poorly with MWDs and soil properties. A perspective that we see is to measure the mass of soil leaving the basket additionally to the mass of soil in the basket. This second measurement would have the advantage to get rid of the interference of air bubbles leaving the basket and therefore limiting the number of factors controlling soil mass evolution under water. At the moment our device doesn’t allow such a measurement. We believe that such a development has the potential to improve greatly curve interpretation regarding the underlying mechanisms of soil disaggregation
• [X] Figure 4 and 5 I think we are missing some results on the MWD here to be able to compare.
◦ We think reviewer point figures: 4 and 6 (same kind of figures, with boxplot). We structured the results of our paper in three main parts: i) comparison QST <-> Le Bissonnais ii) indicators against soil properties (using QST & LeBissonnais) iii) then comparing the treatments of the LTE against soil structural stability, but with a focus on the QST indicators. The two first parts were there to evaluate the QST curves in comparison with recognized ones (Le Bissonnais), while for the comparison of contrasting farming practices, we used only the QST as the central approach of the paper. That's why it did not come into our mind to add boxplot with MWD indicators.
◦ We think also that the 2 graphs are already complete. We already made a "sub-selection" of most relevant QST indicators.
◦ Nevertheless, we found the proposition of the referee and we propose to add 3 graphs in supplementary materials. We will add a reference to these graphs in the text of this part.
• [X] L286 and Fig 3, yes there is a clear relationship but also there is a residual variability. I think this should be acknowledged in the text.
◦ In this part of the results, we focus the presentation of our data on correlation coefficients between variables from soil properties and from indicators. Our aim is not to produce a predictive approach with a model. We observe a standard deviation in the repeated test of a same plot, that could be explained by local soil conditions of the sampling sites (eg. slopes, presence of roots, or earthworms' galleries). We will acknowledge that in the new version of the text.
• [X] L344 some QST indicators. It is possible to be more specific?
◦ We will specify in the text:
◦ "Therefore, we think that that QST indicators Slope~max-60~ might be more specific to slaking than the fast-wetting test of Le Bissonnais."
• [X] L386 Add reference to figure 3 for clarity?
◦ Added:
◦ "It is important to underline that the close linear relationship found between W~max~-W~t0~ and the SOC:clay ratio (see Fig. 3)) has probably no general character and was obtained here because cropland soils of the current study were sampled under standardised conditions of seeding and cover (winter wheat)."
• [X] L410: something wrong with the reference.
◦ Yes, a mistake in the manuscript (omission of a LaTeX command). We adapted:
◦ "… allowed to maintain SOC content to the initial level (about 10 g kg−1; Buysse et al., 2013b). A smaller SOC storage for a similar C input means a higher rate of mineralisation. The formation of water-stable aggregates under the effect of microbial decomposition of root biomass is a known process (Dufey et al., 1986). In the present study, the more microbially active, labile"
• [X] L450 the use of a large soil volume facilitates the soil sample preparation compare to MWD.
◦ We agree that the soil sample preparation of the QST is easier compare to MWD, but we don't think that the sample size is the cause. We will add in the paper, in the beginning of this paragraph:
◦ "The soil sample preparation is far easier than approaches based on MWD as, after removed from the steel cylinder, the sample is dried “as is” and used “as is” in the QST."
◦ If soil sample preparation is shorter, we don’t see it as critical because clods crumbling by hand and soil sieving to 3-5 mm is neither complicated nor time consuming. To our point, the main advantage is the absence of bias compared to sieving that selects aggregates of a certain size, neglecting particles < 3mm that can represent a large fraction of a soil that is frequently tilled
• [X] L500 this sentence is right only for soils of central Belgium for the moment.
◦ We agree (and we are very interested to test the approach in other pedological context). We will adapt this paragraph of the conclusion as following:
◦ "For cropland soils of the silt loam region of central Belgium, we show that the early mass loss under water is mainly related to slaking, whereas after soil saturation with water, clay dispersion becomes the dominant process of soil disaggregation. "Citation: https://doi.org/10.5194/egusphere-2022-1092-AC1
-
AC1: 'Reply on RC1', Frédéric Vanwindekens, 19 Dec 2022
-
RC2: 'Comment on egusphere-2022-1092', Anonymous Referee #2, 18 Feb 2023
First of all, apologies for the delayed response.
This article presents a relevant, novel and worthwhile adaptation of an existing soil structural stability test, and applies it to different case studies, demonstrating its potential for routine evaluations within the context of soil (structural) quality assessments. It certainly deserves to be published. Nevertheless, in my opinion, the introduction needs to be revised, a number of issues have to be clarified (see details below) and additional indicators derived from the QST curves could be tested (without requiring additional measurements).
Introduction
The introduction is quite lengthy. There is a lot of interesting information, but it does not fully serve the purpose of the paper. Part of the difficulty stems from the fact that the authors were not able to decide whether the main focus of the paper is about a new tool or about understanding the effects of various practices on their novel indicators. To me, the paper is first and foremost about a new technique to quantify soil disintegration during immersion, which is then applied to a variety of soil samples from different cropping systems from loess soils in central Belgium. The focus of the introduction should therefore be on existing techniques and why another method like the QST is relevant (L78 and subsequent). L23-49 do not really seem relevant for introducing the QST.
It is also not obvious to me how relevant it is to discuss the hierarchical levels of soil structure, except that the samples should be sufficiently large to encompass all levels so as to be representative of the soil. The fact that the authors do not revert back to these concepts in their discussion clearly indicates that this information is of limited relevance for this paper.
In English, the usual term is ‘aggregate stability’, which in the case of Le Bissonnais is indeed what is being measured (3-5 mm aggregates). It may be worthwhile discussing the distinction the authors make between aggregate stability and structural stability (if any). I would argue that QuantiSlake = structural stability, whereas Le Bissonnais = aggregate stability …
Specific comments
L9 : this is not fully correct. see comments further down on the Le Bissonnais test.
L32 I suppose liming is also an important practice in these systems, that has contributed to changes in pH and base saturation. Liming does not fall under 'organic and mineral fertilizer'. It is a soil amendment.
L35-37 please add reference
L42 add reference for the ‘critical value of SOC content’
L47 ‘…. According to farmers ….’ and scientists !
Starting L46 (till L65), the text becomes somewhat confusion. The paragraph starts about conservation tillage, then switches to erosion, then switches to aggregate stability and the hierarchical structure of aggregates : these are a lot of different ideas for a single paragraph. In general, 1 paragraph = 1 idea.
L57 Ca++ … and Mg++ (as a rule, divalent cations are much more effective than monovalent cations)
L68 The mechanisms listed here are relevant for soil sealing / soil erosion studies. I would argue that in a broader context, mechanical breakdown also occurs as a result of mechanical stresses exerted during tillage operations or traffic.
L86-87 : the method of Koestel is not really about measuring aggregate stability; Likewise, the use of VIS-NIR does not allow to measure stability : it is based on correlations, and therefore will always require reference methods.
L90-93 : there is some confusion here in the text between ‘aggregate size distribution’ (or level of aggregation) and ‘aggregate stability’. Furthermore, the “profil cultural” (and other methods listed here) do not seek to measure aggregate stability.
L95 : not sure I understand why topography and climate would influence the choice of method; please clarify
L108-116 : can be deleted, as it will be repeated in materials and methods / discussion sections
L121 explain abbreviation : CRA-W
L123 : round off to the nearest mm; the decimal really doesn't provide useful information
L148 latin names are provided here but were not given in the previous section when the crops were first named; please correct
L153 “Complete random block with split plot design” : doesn’t sound quite right. ‘split-plot’ is not compatible with ‘completely randomized’ because by definition, at least one treatment (the ‘split’) is not fully randomized. Please clarify.
L165 sampling occurred in April 2019 but no fertilizer was added since 2016 ? Please check
L171 I’m not aware that chlorides can affect soil structure in soils dominated by permanent charge, but potassium (monovalent cation) definitely does!
L198 : did you replicate measurements for each plot ? I believe Le Bissonnais recommends 3-5 replicates per plot (and per test).
L200 : the mechanisms involved in the test are not presented correctly; see comment L211 (below)
L203 : the whole point of rewetting the soil using alcohol prior to shaking in water is also to minimize slaking (and swelling and dispersion)
L203 : No sieving at 50 µm during the immersion phase in ethanol ?
L205 : Start a new paragraph at ‘Two main indicators …’
L209 : how strongly correlated are the MWD and MA indicators ? Is it worth considering both ?
L211-212 : this way of presenting the tests is not strictly correct. rapid wetting involves slaking, dispersion, but also differential swelling! Slow wetting involves differential swelling and dispersion. The mechanical breakdown test seeks to minimize slaking, swelling and dispersion, but it is not obvious that this really mimics drop impact. From the discussion, it appears that the authors are aware of all this, so why present the three test in such a caricatural way ?
Table 1 : please specify the upper and lower limits for the different particle size fractions, as this can be different from one country to another
L225 not very clear. According to L220, the first measurement should be within less than 1 second after plunging the sample into the water. What does ‘right after buoyancy’ refer to ? How do you determine the time to buoyancy? Maybe this could be illustrated graphically.
L229: start new paragraph at ‘Several ….’
L229 If I understand well, this maximum mass is going to depend on the relative rate of disintegration vs. relative rate of wetting, as well as mass of the sample (corrected for buoyancy based on the volume of the solid phase). Doesn’t this also contain information?
L234-236 this has to be better explained. The sentence seems to indicate that it is the slope of the curve at t= 30 sec (with t0 taken as tmax), but this is not what Fig 2 shows. What Fig 2 shows is not the slope of the curve, but the mass loss over a certain time interval (so the slope between t0 and t30, and not the slope at t30).
L240 what about testing (Wmax-Wend) ?
L240 indicators are presented as ‘mass at the end’ and ‘area under the curve’. So units should be ‘g’ and ‘g.sec’. But in later graphs, Wend is presented as unitless and AUC is in 1/sec. So it seems that both indicators have been normalized. This must be explained more clearly.
L241 on line 219, it says the experiment is run for approximately 1000 sec. The AUC is evidently going the depend on the length of the experiment. So is a fixed duration used for all samples ?
The curves are normalized by Wmax, but not by Wend (use Wmax as upper limit and Wenfd as lower limit for normalization): isn’t that introducing some sort of bias ? Also, in the AUC, a large part of the value may come from Wend * duration (the bottom, rectangular area in Fig 2). So a large chunk of the AUC contains the same info as Wend. Normalizing the curves by Wmax and Wend would allow to have an AUC that is independent of Wend.
L243 check this sentence
Caption of Figure 1 : What does 'managing QST laboratory' mean ? Or do you mean 'for managing and graphically displaying QST laboratory data '?
L245 why only the roots remaining in the cage ? Isn’t this a biased estimate, as smaller root fragments may have fallen through the mess during sample breakdown ? Wouldn’t it be more appropriate to pass the entire soil again through a 2-mm mesh sieve to recover all roots?
L256 Although the choice of threshold is to some extent arbitrary, using a 10% probability threshold is somewhat unusual in statistics. Using a higher threshold can sometimes be justified based on the dataset characteristics, but what is the reason for this choice here ? Please justify.
L268 there is no such Fig. 2. Fig. 2 corresponds to the QST and Fig. 3 to another correlation. Perhaps you mean Table 2?
- 3.1 it is also worthwhile noting that correlations are usually stronger with MWD1 than MWD2 (which seems to make sense)
Table 2 : It seems that slope values are expressed in negative values, which was not immediately obvious to me after reading the materials and method section. (I was expecting a higher slope to be negatively correlated with SOC content; the positive correlation stems from the fact that negative values are used). This should be clarified.
L287 : the study of Johannes et al., 2017, did not investigate soil structural stability but soil structural quality (which are two different concepts)! Visual assessment methods don’t allow to assess stability.
In section 3.2, I’m missing information regarding the correlation among indicators from the same test. How are slopes, txx, AUC, … correlated in the QST test ? Is it worth considering all these indicators? Same for MWDs and MAs in Le Bissonnais tests : are they correlated ? This comment also relates to L301-304 : what was the basis for this selection of criteria? Correlation analysis would help justify the selection.
Fig. 4 Wmax – Wt0 are given as unitless, even though in materials and methods the authors talk about mass loss. This is confusing. The same problem arise for slopes, Wend, AUC. The confusion seems to stem from the fact that the indicators are not well explained in the M&M section (one speaks of ‘mass loss’ whereas it should be ‘relative mass loss’).
L336-337 : as mentioned above, it is not correct to attribute single mechanisms to each of the three tests of Le Bissonnais. In the FW-fast wetting test (MWD1), all mechanisms are involved (slaking, dispersion, diff. swelling) except mechanical breakdown. In the SW-slow wetting test (MWD2), dispersion and diff. swelling can be present. So the fact that the QST indicators of L335 are better correlated with FW than SW does not necessarily indicate that slaking is the main mechanisms, but merely that slaking significantly contributes to breakdown in the QST, making the correlation stronger with MWD1. This is a subtle but important difference.
L337-340 : the ‘problem’ with t50, t75, etc… is that these indicators inherently also include what happens during the early times. Wouldn’t it be much more relevant to consider the time it takes for the sample to loose 25% of its mass, then the time from 25 to 50%, then the time from 50 to 75, etc … ? In this way, the various indicators can be expected to be much more independent, and interpretation of processes would also be facilitated. If indeed dispersion and diff. swelling are the main mechanisms at later stages, then one would expect even stronger correlations between MWD2 and, say, the time needed to lose 75-90% of the mass than between MWD2 and t90.
L343 : indeed ! … and from differential swelling. If the authors are aware of this, then why present the Le Bissonnais tests in the wrong way in materials and methods ?
L345 I think the authors misunderstand the FW test of Le Bissonnais. Slaking plays an important role in this test, but it is not the only mechanism at play. It is the difference between the SW and FW tests that predominantly reflects the effect of slaking (though slaking may actually also facilitate dispersion, so the effects are not simply additive). One cannot say that the FW test alone reflects the sole effect of slaking.
L348 indeed, the continuous measurements of QST would allow to better discriminate between mechanism on a single sample, something that cannot be done with the le Bissonnais test. Nevertheless, as suggested above, using time intervals (t0-25, t25-50,t50-75, …) should allow to better discriminate between processes than t0-25, t0-50, t0-75, etc. Furthermore, the difference between MWD1 and MWD2 of Le Bissonnais should better reflect the sole effect of slaking than MWD1 alone, so testing correlations between QST results (early times) and (MWD2-MWD1) might be of interest.
The same reasoning would be applied to the slopes : slope max-60, 60-300, 300-600.
L385 : again, it is not correct to state that VESS assesses structural stability (it is a measures of quality)
L421 : this is rather speculative. It is not certain at all that this would be the case. FYM is very different from green manures or crop residues (much lower C/N ratio right from the start) and therefore will not interact with microbial life in the same way.
L424-425 : I think I sort of understand what the authors mean, but the sentence is awkward and must be clarified
L438-440: Again, I’m not aware that chloride can have such effects on soils with permanent charge. There is no mention of the effect of chloride on stability in the paper of Paradelo. It is related to the monovalent cation (K). This discussion must be revised.
L450 aggregate stability tests such as Le Bissonnais or drop impact tests (Imeson and Vis) where developed in view of relating the results to soil erodibility or even as a way of parameterizing soil detachment during rainfall in erosion models. In that respect, working with aggregates seems fairly relevant (particularly in case of cultivated tilled soils). As a means to assess soil quality, the larger sample size used in the QST seems relevant. Whether QST results could be used to assess soil erodibility remains to be proven, even though the good correlations with Le Bissonnais results are encouraging.
L460 : differential swelling ?
L461 see suggestions above on how to possibly improve the indicators
L463-465 : again, there is a bit of a mixup here between ‘stability’ and ‘quality’
L468 : indeed, some overlap is expected with the current indicators, but this could possibly be improved (see suggestions above)
L484 I may be wrong (those documents are not easily accessible), but I believe that the equation developed in Laon does not relate to detachment by drop impact per se (but rather to the aggregate stability tests of Henin and Monnier as a measure of the soil sensitivity to crusting). To be checked.
L494-497 : future developments are not to be mentioned in the conclusion
L503 : I do not believe this was really demonstrated in this paper (as a matter of fact, linear correlations – as demonstrated in this paper – are not compatible with thresholds) , and it should therefore not appear here.
See additional annotations in the attached pdf file
Frédéric M. Vanwindekens and Brieuc F. Hardy
Data sets
Full git repository - QuantiSlakeTest - qst-openscience Vanwindekens F. M. and Hardy B. F. https://gitlab.com/FrdVnW/qst-openscience
Model code and software
R-package - slaker - Analysing the data of QuantiSlakeTest approach. R-package and Web Application Frédéric M. Vanwindekens and Roisin C. https://gitlab.com/FrdVnW/slaker
Notebook with codes, figures and tables - qst-openscience Vanwindekens F. M. and Hardy B. F. https://frdvnw.gitlab.io/qst-openscience/
Video supplement
A visualisation of the QuantiSlakeTest, comparing two contrasted samples Frédéric M. Vanwindekens https://youtu.be/G9UweThvHYI
Frédéric M. Vanwindekens and Brieuc F. Hardy
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
395 | 145 | 12 | 552 | 4 | 5 |
- HTML: 395
- PDF: 145
- XML: 12
- Total: 552
- BibTeX: 4
- EndNote: 5
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1