Electronic Journal of Polish Agricultural Universities (EJPAU) founded by all Polish Agriculture Universities presents original papers and review articles relevant to all aspects of agricultural sciences. It is target for persons working both in science and industry,regulatory agencies or teaching in agricultural sector. Covered by IFIS Publishing (Food Science and Technology Abstracts), ELSEVIER Science - Food Science and Technology Program, CAS USA (Chemical Abstracts), CABI Publishing UK and ALPSP (Association of Learned and Professional Society Publisher - full membership). Presented in the Master List of Thomson ISI.
2005
Volume 8
Issue 4
Topic:
Forestry
ELECTRONIC
JOURNAL OF
POLISH
AGRICULTURAL
UNIVERSITIES
Kormanek M. , Walczykova M. 2005. CHARACTERISTICS OF HORIZONTAL DEFORMATION OF SOME FOREST SOILS, EJPAU 8(4), #39.
Available Online: http://www.ejpau.media.pl/volume8/issue4/art-39.html

CHARACTERISTICS OF HORIZONTAL DEFORMATION OF SOME FOREST SOILS

Mariusz Kormanek1, Maria Walczykova2
1 Department of Forest Work Mechanization, Faculty of Forestry, University of Agriculture in Krakow, Poland
2 Department of Machinery Management, Ergonomics and Production Processes, Faculty of Production and Power Engineering, University of Agriculture in Krakow, Poland

 

ABSTRACT

The paper presents results of measurements of horizontal shearing of soil, made by the help of wheel bevameter on four types of forest sites, i.e. fresh mixed coniferous forest (fMCF), boggy mixed pine forest (bMPF), moist mixed broadleaved forest (mMBF) and moist broadleaved forest (mBF), with consideration of different moisture variants. On the basis of obtained data two relationships were determined: between soil horizontal deformation and the circumferential (driving) force as well as between slip and the coefficient of adhesion. Grouping of the investigated soil according to their shearing stresses was made by Ward’s method of agglomerative cluster analysis. For each of the obtained agglomeration new data of slip vs. coefficient of adhesion were calculated.

Key words: forest soils, coefficient of adhesion, slip, shearing stresses, cluster analysis.

INTRODUCTION

Environmentally friendly forest development is considered to be inevitable with regard to restraining soil erosion processes, preservation of water resources as well as protection of the landscape and natural environment [12]. Timber harvesting that enables using wood resources has at the same time negative influence on the environment. As it requires great physical effort of men it must be mechanized [14, 16, 19] and mechanized felling, bucking and skidding directly affect the forest stand and soil. Relocation of machines under forest conditions takes place both on skidding roads and forest sites [6, 9, 18], and it can be considered in two aspects. The first, it is effect of the vehicle wheels that are of prime importance as far as soil damage is concerned [1, 8, 11, 21]. The second is related to the machines traction properties. The wheel ability to transmit driving force on a given soil is decisive with regard to the vehicles traction properties such as pulling force and slip properties, and in fact determines the vehicle mobility under given soil and terrain conditions [1, 7, 17]. Knowledge of traction parameters in relation to the forest soils makes possible to determine the effects that the vehicle used under given conditions can have on the performed job as well as to foresee them [1].

On the basis of available literature review dealing with forest machinery management, it can be stated that only very limited number of research attempts were made towards improving knowledge of forest soil traction problems.

In the presented paper problems of traction properties of four forest sites were undertaken and its scope included:

MATERIALS AND METHODS

Variables that determine mutual influence of pneumatic wheel and the ground are basically traction forces making the vehicle movement possible, as well as slips generated by these forces [1]. In classical theory of mobility the coefficient of adhesion is applied in assessment of mechanical vehicles traction properties. This approach is commonly used in approximate calculations of agricultural and forest tractors operating parameters [7, 10].

In method based on coefficient of adhesion µ, the calculated force of adhesion Fµ is understood as a maximum force that wheel can transmit on the ground. The force is proportional to the wheel loading Qk and its value depends on coefficient of adhesion µ characteristic of the given ground (substrate) (eq. 1) [1, 2, 5, 15, 20].

Fµ = Qk µ     (1)

The driving wheel will move, when its circumferential (driving) force will be less or equal to the adhesion force (eq. 2).

Fk ≤ Fμ     (2)

Otherwise slipping of the wheel occurs. Values of coefficient of adhesion of a pneumatic tyre depend on type of ground and its mechanical and physical properties as well tyre properties. The above mentioned brings about considerable difficulties in obtaining accurate values of coefficient of adhesion. In calculation only approximate values obtained in measurements are used [5]. Those available in literature are basically related to agricultural soils. Appropriate data for forest soils are in fact non existing.

Layout of measurements

Terrain measurements of horizontal shearing of forest soils were made by the help of wheel bevameter [18] on four types of forest sites, i.e. fresh mixed coniferous forest (fMCF), boggy mixed pine forest (bMPF), moist mixed broadleaved forest (mmBF) and moist broadleaved forest (mBF). With consideration of different soil moisture contents altogether 8 variants were investigated.

Occurrence of features characteristic of the given forest site was used as a criterion of its choice. Open canopy of the forest stands that would not restrain relocation of the measuring devices, and as far as possible flat surface, were taken into consideration as additional factors during laying out the experimental site.

On each site four blocks of ca 30 m long were laid out. Each block consisted of 18 repetitions, which gave 72 replicates for one variant of measurement. Each wheel rut, formed as a result of the soil shearing process, was subject to measurement of its dimensions.

Physical properties of soils were characterized by dry volumetric bulk density, moisture content by weight and penetration resistance. Measurements of the latter were made by cone penetrometer with an electronic data logger, in 120 replicates for each experimental site. Volumetric cylinders were used for density and moisture determination, collecting 5 samples from 3 depth, up to 0.3 m.

The bevameter soil shearing wheel was equipped with tyre of 165/80-R13 dimensions, inflated to pi = 240 kPa and loaded by Qk = 4777 N, according to the catalogue data.

Processing of measurement data and calculation of the soil horizontal deformation

Quantity obtained in terrain measurements was torque Mk, measured on the bevameter soil shearing wheel and recorded at every 1.5° of the wheel turning angle.

In order to get the relationship between horizontal displacement of soil j and the circumferential shearing (driving) force Fk, the following calculations were made. Values of the turning torque Mk were divided by the wheel dynamic radius Rd, determined on the site during measurements (eq. 3).

     (3)

Conversion of the wheel turning angle into soil horizontal displacement was done by determination of distance (arch) O, along which wheel circumferential points were relocated during the wheel turn by angleb (eq. 4).

     (4)

In order to get coefficients of adhesion µ, measured circumferential shearing forces Fk were divided by the wheel vertical loading (eq. 5).

     (5)

The next step consisted of determination of slip value s for each value of coefficient of adhesion µ. According to the literature [3] it was assumed that the slip reaches its maximum (i.e 100%), when the wheel turns by an angle corresponding to relocation of its circumferential points along the full length of the soil-wheel contact area. Thus, division of the 100% value by the wheel turning angle embracing the length of the contact area resulted in slip value per 1° of the wheel turn, for which corresponding coefficient of adhesion was determined (eq. 5).

It is known that on deformable surfaces the ability of driving wheels to transmit circumferential forces depends on relation between shearing stresses τ and the soil horizontal displacement j. The above mentioned relations Fk = f(j) were than converted into relations τ = f(j) through dividing the shearing (driving) force by the wheel-soil contact area (eq. 6).

     (6)

Dimensions for calculation of the contact area (eq. 7) were determined under terrain measurements.

     (7)

Width of the wheel rut b corresponded in fact with the width of tyre, the length was calculated from eq. 8 [13]. Depth of the rut z and tyre deformation D were measured for the given terrain conditions.

L= [[D (z + Δ) – (z + Δ)2 ] + (D Δ – Δ2 )]2     (8)

RESULTS

Characteristics of the forest site properties

On both sites of fresh mixed coniferous forest (fMCF) (TABLE 1) podsolic soils occured from sandy soils textural group [4]. Upper organic layer present on both compartments was made up of fresh organic matter. On compartment No 188 it was raw humus richly covered by reedgrass, while on compartment No 244 the organic matter was covered by bilberry.

Soil on boggy mixed pine forest (bMPF) was classified as boggy rot on light loamy sand [4]. This place differed from the rest of investigated sites by its low bulk density and high moisture content (TABLE 1). Average penetration resistance there was lowest from among all investigated sites (FIGURE 1).

Table 1. Moisture and density of soils on the investigated forest sites (w1, w2, w3 – denotations of moisture contents; higher number corresponds with higher moisture)

Sites

Depth
(m·10 –2)

Moisture by weight

Dry bulk density

Average
(%)

Coeff. of var.
(%)

Average
(Mg/m3)

Coeff. of var.
(%)

fMCF _188_w1

0-10

17.2

54.5

1.05

24.9

11-20

7.6

81.2

1.40

13.6

21-30

6.2

21.0

1.40

4.5

fMCF _244_w2

0-10

21.9

51.8

0.61

37.2

11-20

5.9

50.0

1.17

12.7

21-30

3.3

44.1

1.28

9.8

bMPF

0-10

160.3

21.4

0.61

35.4

10-20

202.4

21.4

0.65

22.1

20-30

177.3

40.9

0.79

28.8

mMBF _244_w1

0-10

93.0

68.1

0.55

41.8

10-20

44.3

43.9

1.00

26.0

20-30

39.8

82.8

1.05

31.7

mMBF_291_w2

0-10

158.9

52.5

0.22

28.7

10-20

47.6

66.7

1.01

20.7

20-30

28.8

38.4

1.19

28.1

mMBF_291_w3

0-10

268.4

9.24

0.22

28.7

10-20

45.2

35.8

1.01

20.7

20-30

27.3

45.2

1.19

28.1

mBF_w1

0-10

34.8

11.4

1.14

6.8

10-20

29.9

10.7

1.36

4.3

20-30

30.3

6.7

1.37

5.9

mBF_w2

0-10

52.0

12.5

1.14

6.8

10-20

40.8

8.0

1.36

4.3

20-30

41.4

6.0

1.37

5.9

Gley soils characterized the moist mixed broadleaved forest (mMBF) site. Texturally they were numbered among strong loamy sands (compartment No 244), and sands layered by non skeleton light loamy sand (compartment No 291) [4]. On both compartments the layer up to 0.1 m was made of raw organic matter.

Soil physical properties measured on the moist broadleaved forest (mBF) were characterized by lowest coefficient of variability (TABLE 1). This was mainly caused by lack of soil cover. Up to 0.35 m the soil on this site was classified as heavy silty loam. Under dry conditions it formed dense structure, while at higher moisture it became plastic.

Average penetration resistance concerning all sites and their moisture variants were presented in FIGURE 1.

Figure 1. Average penetration resistance on the investigated forest sites

Shearing forces and slips

Horizontal shearing of soils was carried out on the above described sites and for each research variant results were drew up in form of charts. For each considered situation two charts were presented – relation between shearing strength (wheel driving force) vs. soil horizontal displacement and slip vs. coefficient of adhesion.

On the soil of fresh mixed coniferous forest site (fMCF) maximum value of shearing strength occurred at 0.11m and 0.125 m soil displacement, which corresponds with lower and higher moisture respectively (FIGURE 2). Charts in FIGURE 3 and data included in TABLE 2 show that higher coefficient of adhesion amounting to 0.73 was obtained on weaker soil of fMCF_244_w2 (TABLE 1). This phenomenon was caused by presence of the above-ground part of bilberry and its root system on that site. Maximum driving force that can be transmitted on the considered soils of fresh mixed coniferous forest site (fMCF) would cause 40.2% and 55.3% of slippage on denser (188_w1) and weaker (244_w2) variant respectively (TABLE 2).

Figure 2. Shearing strength in relation to soil horizontal displacement – fresh mixed coniferous forest (fMCF)

Figure 3. Relationship between slip and coefficient of adhesion – fresh mixed coniferous forest (fMCF)

In the case of boggy mixed pine forest (fMCF) (FIGURE 4) the soil horizontal displacement at maximum shearing strength, corresponding with maximum driving force, came to 0.081m and beyond that point its value was slowly decreasing along the wheel rut. Boggy rot soil of low bulk density and high moisture content (TABLE 1) as well as lack of living plants cover contributed to forming the longest wheel rut (0.29 m) from among all investigated sites and their moisture variants. Curve of slip dependence on coefficient of adhesion (FIGURE 5) showed characteristic bending at the moment when maximum value of coefficient of adhesion was reached. It happened at 30% slip (TABLE 2), after which values of the driving force on that soil were dropping.

Figure 4. Shearing strength in relation to soil horizontal displacement – boggy mixed pine forest (fMCF)

Figure 5. Relationship between slip and coefficient of adhesion – boggy mixed pine forest (fMCF)

Results of measurements on moist mixed broadleaved forest (mMBF) in three moisture variants depict charts in FIGURE 6. Shape of the curves is characteristic for soils revealing both cohesion and friction. Values of coefficient of adhesion were decreasing with increasing moisture content (FIGURE 7) and were located between 0.73-0.77 (TABLE 2).

Figure 6. Shearing strength in relation to soil horizontal displacement – moist mixed broadleaved forest (mMBF)

Figure 7. Relationship between slip and coefficient of adhesion – moist mixed broadleaved forest (fMCF)

Curves of shearing stresses (wheel driving forces) on soil of moist broadleaved forest (mBF) showed asymptotic character without revealing maximum value. It happened at both moisture variants. At lower moisture content lugs hardly sunk into the hard soil and the shearing tyre gave worse grip. At higher moisture soil became plastic and yielded without showing clear shearing maximum. Lack of plant cover considerably contributed to the above described phenomenon. Obtained coefficient of adhesion amounted to 0.52 and 0.60 for dry and moist soil of moist broadleaved forest (mBF) respectively (TABLE 2).

Figure 8. Shearing strength in relation to soil horizontal displacement – moist broadleaved forest (mBF)

Figure 9. Relationship between slip and coefficient of adhesion – moist broadleaved forest (mBF)

Table 2. Values of slip obtained for different coefficient of adhesion for all investigated forest soils

Coefficient of adhesion
µ

Forest sites

fMCF

bMPF

mMBF

mBF

w1

w2

w1

w2

w3

w1

w2

0.05

2.2

3.3

2.8

2.3

2.4

2.2

4.5

4.01

0.10

3.9

5.1

4.1

4.1

4.1

3.9

7.5

6.32

0.15

5.3

6.8

5.4

5.4

5.4

5.0

9.5

7.60

0.20

6.6

7.9

6.4

6.7

6.7

6.2

11.3

8.88

0.25

7.7

9.0

7.5

7.9

7.8

7.3

12.9

10.14

0.30

8.8

10.2

8.6

9.1

8.9

8.3

14.6

11.41

0.35

10.0

11.3

9.8

10.4

9.9

9.4

16.4

12.87

0.40

11.3

12.4

11.0

11.6

11.2

10.7

18.3

14.69

0.45

12.7

13.8

12.4

12.9

12.5

12.0

21.0

17.53

0.50

14.6

15.6

13.4

14.5

14.0

13.5

24.3

32.51

0.52

15.0

16.3

14.6

15.2

14.8

14.4

25.7

41.5

0.54

16.3

17.4

15.8

15.9

15.5

15.2

28.2

 

0.55

17.0

18.0

16.3

16.3

15.9

15.6

29.4

 

0.60

20.7

21.2

19.6

18.7

18.6

18.7

40.3

 

0.65

26.7

26.2

29.6

22.0

21.8

24.3

   

0.70

40.2

37.1

 

27.4

31.3

33.6

   

0.73

 

55.3

 

32.6

39.4

57.6

   

0.74

     

34.8

45.9

     

0.75

     

37.4

       

0.77

     

60.4

       

Grouping of the investigated substrates

Shearing strengths (FIGURES 2, 4, 6, 8) were converted into shearing stresses (FIGURE 10) according to equation 6 and agglomerative cluster analysis by Ward’s method was carried out. Results of the procedure were shown in FIGURE 11.

From among 8 research variants differing from each other in type of soil, moisture and plant cover, four groups (clusters) emerged.

Cluster 1: fMCF_w2; mMBF_w1
Cluster 2: fMCF_w1; mMBF_w2; mMBF_w3
Cluster 3: mBF_w1
Cluster 4: bMPF; mBF_w

Order in which groups of substrates were listed corresponds with decrease of their resistance to shearing.

Figure 10. Shearing stresses vs. soil displacement for the investigated forest soils

Figure 11. Chart of cluster analysis: fresh mixed coniferous forest (fMCF), boggy mixed pine forest (bMPF), moist mixed broadleaved forest (mMBF); moist broadleaved forest (mBF); w1; w2 ; w3 – moisture contents

Cluster 1 includes forest soils characterized by ability to transmit high shearing stresses. Here can be found sites with sandy soils at low and medium moisture content. Nevertheless, their resistance to shearing results mainly from the plant covers. That’s why substrate of fresh mixed coniferous forest (fMCF_w1) at lower moisture (w2) turned up here, in spite of less favorable soil physical properties than of its equivalent fMCF_w2.

Cluster 2 includes more resistant sandy soil substrates. Moist mixed broadleaved forest (mMBF) at the considered moisture contents showed best traction ability.

In the third group (cluster 3) substrate of moist broadleaved forest (mBF_w1) was found. As it was already mentioned, hard surface, bad grip of tyre and lack of plant cover can be named as reasons for reaching rather small shearing stresses.

The last group (cluster 4) is made up of two weakest substrates – moist broadleaved forest (mBF) at plastic conditions and boggy mixed pine forest (bMPF).

For each group of substrates agglomerative curves of shearing stresses vs. soil horizontal displacement (FIGURE 12) and slip vs. coefficient of adhesion (FIGURE 13) were drawn up. In TABLE 3 values of slip for determined clusters were included.

Figure 12. Shearing stresses vs. soil displacement for obtained clusters

Figure 13. Relationship between slip and coefficient of adhesion for obtained clusters

Table 3. Values of slip depending on coefficient of adhesion for determined groups of forest soils

Coefficient of adhesion
µ

Slip s (%)

Cluster 1

Cluster 2

Cluster 3

Cluster 4

0.000

0.0

0.0

0.0

0.0

0.050

2.7

2.3

4.5

3.2

0.100

4.6

4.0

7.5

4.8

0.150

6.1

5.2

9.5

6.2

0.200

7.5

6.5

11.3

7.3

0.250

8.7

7.6

12.9

8.4

0.300

9.9

8.7

14.6

9.6

0.350

11.1

9.7

16.4

10.6

0.400

12.3

11.1

18.5

12.2

0.450

13.6

12.4

21.0

13.7

0.500

15.4

14.1

24.3

16.3

0.550

17.3

16.0

29.4

22.

0.580

18.9

17.9

35.3

35.7

0.600

20.1

19.2

40.3

 

0.649

24.3

24.3

   

0.650

24.4

24.5

   

0.700

31.5

34.3

   

0.723

37.3

48.1

   

0.750

50.7

     

0.752

54.0

     

CONCLUSIONS

  1. Relationships between slips and coefficients of adhesion determined for some soils make possible rough assessment of the vehicle traction properties working under comparable conditions.

  2. Agglomeration of the investigated soils with consideration of their moisture variants showed that type of the forest site does not condition to the full extent the soil resistance to horizontal shear. Heavy soil of moist broadleaved forest (mBF) at plastic state behaved with regards to traction similarly as weak soil of boggy mixed pine forest (bMPF).

  3. Plant cover of the forest site soils had considerable and positive influence on their traction properties.


REFERENCES

  1. Ala-Ilomäki J. 1993. A terramechanical approach for evaluating mobility and ground disturbance during skidding and forwarding: preliminary trials. Special Report SR-86, 1-14.

  2. Arczyński S. 1984. Teoria ruchu samochodu. [Theory of car mobility]. Wydawnictwo Politechniki Warszawskiej, s. 253 [in Polish].

  3. Bekker M. G. 1969. Introduction to terrain – vehicle systems. Ann Arbour. University of Michigan Press, s. 820.

  4. Brożek S., Zwydak M. 2003. Atlas gleb lesnych Polski. [Atlas of forest soils in Poland]. Centrum Informacyjne Lasów Państwowych, Warszawa [in Polish].

  5. Dajniak H. 1985. Ciągniki. Teoria ruchu i konstruowanie. [Tractors. Theory of mobility and construction]. Wydawnictwo Komunikacji i Łącznosci, Warszawa, s. 741 [in Polish].

  6. Froehlich H.A. 1986. Causes and effects of soil degradation due to timber harvesting. Canadian Forestry Service Victoria B.C. Timber harvesting session, 3-12.

  7. Grecenko A. 1994. Vlastnosti terénních vozidel. [Properties of off-road vehicles]. Vysoká škola zemedelská v Praze, Technická fakulta, s. 118 [in Czech].

  8. Krag R. K., Holmsen S. D., Wong T. B. 1993. Effects of log-skidder travel on soil bulk density. Forest Engineering Research Institute of Canada. Special Report, SR–87, Timber Harvesting Session, 1-18.

  9. Laurow Z. 1996. Szlaki technologiczne w procesie pozyskania drewna. Czesc II. Szlak a srodowisko. [Technological roads in wood harvesting. Part II. Skidding road and the environment]. Przegląd Techniki Rolniczej i Leśnej, 6, 23-25 [in Polish].

  10. Lukáč T. 1997. Lesné dopravníctvo. [Forest transport]. Technická univerzita vo Zvolene, s. 227 [in Slovak].

  11. Marshall V.G. 2000. Impacts of forest harvesting on biological processes in northern forest soils. Forest Ecology and Management, 133, 43-60.

  12. Raport o stanie lasów. [Report on forests condition]. 2002. http://www.lasypanstwowe.gov.pl [in Polish].

  13. Schwanghart H. 1991. Measurement of contact area, contact pressure and compaction under tires in soft soil. Journal of Terramechanics, Vol. 28, 4, 309-318.

  14. Sosnowski J. 2002. Szlaki zrywkowe w proekologicznym gospodarstwie lesnym. Czesc I. Szlaki a przygotowanie stanowisk roboczych. [Skidding roads in environmentally friendly forestry enterprise. Part I. Strip roads and preparation of working stands]. Sylwan, 8, 73-80 [in Polish].

  15. Studziński K. 1973. Teoria, konstrukcja i obliczenia samochodu. [Theory, construction and calculation of a car]. Wydawnictwo naukowo-techniczne, Warszawa, s. 834 [in Polish].

  16. Walczyk J., Walczykova M. 2002. Trends in development of multifunctional skidding machines. Mobilné energetické prostriedky – Hydraulika – Ekologické oleje – Zivotné prostredie. Fakulta enviromentálnej a výrobnej techniky. Technical University in Zvolen. ISBN 80-228-1171-8, 126-132.

  17. Walczyk J., Walczykova M. 2002. Traction properties of some skidding tractors and the possibilities of their use in the Carpathian Forests. Logistics of Wood Technical Production in the Carpathian Mountains. Faculty of Forest of Technical University in Zvolen. ISBN 80-228-1170-X, 291- 299.

  18. Walczykova M., Walczyk J., Kormanek M. 2002. DETERMINATION OF FOREST SOIL TRACTION PARAMETERS. EJPAU, 5(2), #6

  19. Więsik J. 1998. Czynniki decydujące o wyborze maszyn do pozyskania drewna w Polsce. [Factors decisive in selection of machines for timber harvesting]. Przeglad Techniki Rolniczej i Lesnej, 6, 6-9 [in Polish].

  20. Wong J. Y. 1993. Theory of ground vehicle. A Wiley – Interscience Publication, s. 435.

  21. Wronski E. B., Murphy G. 1994. Responses of forest crops to soil compaction. In: Soil compaction in crop production. Elsevier Science B. V., 317-342.

List of symbols

A – soil-wheel contact area (m2)
b – width of the tyre (m)
D – wheel diameter (m)
Fk – driving force (kN)
Fµ – adhesion force (kN)
j – soil horizontal displacement (m)
L – length of soil-wheel contact area (m)
Mk – wheel driving torque (N·m, kN·m)
O – length of arc (m)
pi – inflation pressure (kPa)
Qk – wheel loading (kN)
Rd – dynamic radius (m)
s – slip (%)
z – depth of wheel rut (m)
β – wheel turning angle (°)
Δ – tyre deformation (m)
µ –coefficient of adhesion (–)
µmax –maximum coefficient of adhesion (–)
τ – shearing stresses (kPa)


Mariusz Kormanek
Department of Forest Work Mechanization, Faculty of Forestry, University of Agriculture in Krakow, Poland
al. 29 Listopada 46
31-423 Kraków
Poland
Phone: +48 12 662 5024
email: rlkorma@cyf-kr.edu.pl

Maria Walczykova
Department of Machinery Management, Ergonomics and Production Processes, Faculty of Production and Power Engineering, University of Agriculture in Krakow, Poland
Balicka 104, 30-149 Cracow, Poland
Phone: (012) 662 4634
email: rtwalczy@cyf-kr.edu.pl

Responses to this article, comments are invited and should be submitted within three months of the publication of the article. If accepted for publication, they will be published in the chapter headed 'Discussions' and hyperlinked to the article.