Geografisk Tidsskrift, Bind 74 (1975)AggregatstabilitetsanalyseHenrik Madsen Side 70
Madsen, Henrik,
1975: Aggregatstabilitetsanalyse. Geografisk
A method to describe the stability and size distribution of soil aggregates by dry sifting. The profit and disadvantages of the method are discussed. A methematical adaption is given and further the influence of humus on the soil aggregates is discussed, as well as the possibility of the method for giving significant values in relation to soil erosion. Henrik Madsen, M.
Se. Geographical Institute, University of
IndledningDenne artikel er et forsøg på at finde en metode til beskrivelse af en jords aggregatstabilitet og -fordeling ved tørsigtning, samt at påpege fordele og ulemper ved den valgte procedure. Under den matematiske bearbejdelse af materialet vil tesen: »mængden (gram) af aggregater, der destrueres pr. tidsenhed er proportional med den tilstedeværende intakte mængde af aggregater« blive udsat for en nøjere vurdering. Det vil endvidere blive forsøgt at belyse humusindholdets betydning for aggregatstørrelsesfordelingen og -stabiliteten. Til sidst vil analysens muligheder inden for muldflugtstudierne blive berørt. StrukturstabilitetsbegrebetVed en jords aggregatstabilitet forstås aggregaternes modstand mod mekanisk nedbrydning eller destruktion ved vands indvirken (sprængning ved fugtning eller regndråbers destruktive influens). Aggregaternes nedbrydningshastighed er en funktion af kationbelægningen på kolloiderne, jordvandets saltkoncentration, sesquioxidindholdet, rodudviklingen i jorden og ler- og humusindholdet, der ved komplexforbindelser danner broer mellem de større primærpartikler. Endvidere spiller aggregaternes fugtighed en uhyre stor rolle, og aggregater, der ved et bestemt vandindhold er stabile, vil under andre fugtighedsforhold være ustabile. Til bestemmelse
af aggregatstabiliteten kan benyttes a) tør og våd
sigtning, hvor især den mekaniske påvirknings
b) dråbemetoden:
ved denne metode udsættes aggregaterne c) udvaskning med
NaCl: ved denne metode bliver
d) behandling af
aggregaterne med alkohol og organiske Aggregatstabilitetens betydningStrukturen er af stor betydning for jordens vand- og varmebalance, udluftning samt for planternes rodudvikling 0.1. En ødelæggelse af strukturen vil ændre hele det fysiske miljø i jorden og nedsætte den biologiske aktivitet. Det er derfor vigtigt, at have en god og stabil struktur, der kan modstå nedbørens destruktive virkning (især på marker uden vegetation), vindens erosion og mekanisk ødelæggelse ved landbrugsmaskiner, der kan medføre pløjesåler. Det er især gennem forbedring af infiltrationskapaciteten og den dermed forbundne mindre run-off, at strukturen har betydning for den vandinfluerede soilerosion, hvorimod det er aggregaternes størrelse og deres stabilitet, der har betydning ved vinderosion. I regioner, hvor vinderosion kan finde sted, vil en god strukturstabilitet hindre, at aggregaterne bliver slået i stykker af eventuelt saltationsmateriale og derved selv blive omdannet til dette. Hvis en sådan situation opstår, vil resultatet blive en afblæsning af en del af overjorden ved suspension, saltation og surface creep. Dette medfører en forringelse af markens dyrkningsværdi som vil tage år at genoprette, om nogensinde. Den her anvendte
tørsigtningsmetode vil især kunne Prøveudtagning og foreløbig analyseDe til undersøgelsen anvendte aggregater blev udtaget fra en brunjord (cambisol) ved Karlstrup Kalkbrud nær Køge. Vegetationen på prøveudtagningsstedet var græs. Jordbundsprofilet, tilhørende texturanalyse (Andreasens pipettemetode for fraktionen mindre end 63^1) og pHcaCia ses i fig. i. Man bør bemærke
det store lerindhold i jorden med et Side 71
Aslyngske system må jorden betegnes som en lerjord. pH er konstant gennem profilet til dybden 135 cm, hvorefter der er en jævn stigning ned til kalken. Skønt det rene danienkalk fandtes i 190 cm, sporedes frit calciumcarbonatførst i dybden 135 cm, og da kun i yderst ringe mængde. Selv i dybden 175 cm var calciumcarbonatindholdetså ringe, at Scheiblers apparat ikke kunne anvendes.Humusprocenten blev bestemt ved Ter Meulens metodepå prøverne i A-horisonten. Resultatet blev: i o—20 cm 4,2 %
humus Hvad strukturen angik, viste A-horisonten en svag brokket struktur, hvorimod polyeledstruktur dominerede i Bhorisonten. Der var rodudvikling i hele profilet med hovedparten af rødderne koncentreret i de øverste 10 cm. Arbejdshypoteser og analyseprogramUndersøgelsen kan
deles i to dele: i. en
undersøgelse af aggregaterne separat i
størrelsesfraktioner 2. en
aggregatanalyse, hvor alle fraktioner deltager i
analysen Undersøgelserne
under del i koncentreres om følgende a. en bestemmelse
af halveringstiden for de forskellige b. en undersøgelse af tesen: mængden (vægt) af aggregater, der destrueres pr. tidsenhed, er proportional med den tilstedeværende intakte mængde af aggregater dvs.: dx/dt = H-kx=}lnx = -4-kt + c=) log x^x,, = -r- kt/2,303 i) hvor x er mængden af intakte aggregater. x0 betegner intakte aggregater ved begyndelsen (tiden o), og xt er intakte aggregater til tiden t. k betegnes aggregatkonstanten og varierer mellem o og uendeligt. Aflæses log Xt/x0
som funktion af t, fås, hvis teorien Undersøgelserne
under del 2 koncentreres om følgende c. de fine
fraktioners modstand mod vinderosion. Da de to analysemetoder er forskellige vil det være rimeligt at behandle dem særskilt. Analysemetodik og resultater fra undersøgelsen omhandlende aggregatstabiliteten indenfor de enkelte størrelsesfraktioner vil derfor først blive omtalt. DEL 1Analysens udførelseDa fremgangsmåden
er den samme for horisontal og vertikalsigtning, Jordprøverne tørres et par døgn i fri luft og sigtes derpå (i hele phi-intervaller) i tredive sekunder i hånden, hvorved en grovsortering opnås. Derefter påføres hver sigte en bund og sigtes femten sekunder i sigtemaskinen, hvorved mindre aggregater og støv, der er på vej ned gennem sigterne, fjernes. Der begyndes med den groveste sigte, og bundmaterialet hældes over i den næste sigte (en phi mindre), der så bliver behandlet. Derved sikres, at en aggregatstørrelsesfordeling kan blive udtegnet. (Fejlen ved dobbeltsigtningen af nogle af aggregaterne omtales senere). Ved vertikalsigtningen førtes bundmaterialet ikke over i den næste sigte. Aggregaterne i de forskellige størrelsesfraktioner vejes (x gram) og sigtes derefter i bestemte tidsrum, og det destruerede materiale Dt (materialet i bunden) vejes efter hver sigtning. Aggregaterne sigtes i den sigte, de under forproceduren blev holdt tilbage i. Når den sidste sigtning er tilendebragt, udvaskes aggregaterne, og residualet (primærpartikler større end hullerne i sigten) vejes og subtraheres fra totalmængden, x0 = x -f- p, hvor p er mængden af primærartikler. I de fraktioner, hvor aggregatmængderne er så store, at de vil hæmme sigtningen Side 72
ved ikke at kunne
bevæge sig rimeligt frit, udtages repræsentative
Resultater diskussiona. Som mål for aggregatstabiliteten i de forskellige størrelsesfraktioner, kan man anvende aggregaternes halveringstid, altså det tidspunkt hvor X(/x0 = Vz. Ofte vil man på grund af aggregaternes store stabilitet ikke sigte til t*/2 men beregne den, hvilket vil blive diskuteret senere. Halveringstiden er ikke et entydigt mål for aggregatstabiliteten, da alle aggregater inden for en størrelsesfraktion ikke behøver at have samme stabilitet. Som eksempel på dette kan fraktionerne 16-8 mm, 8-4 mm og 4-2 mm fra dybden 90-110 cm fremdrages. Det ses, at halveringstiderne næsten er identiske, medens andre X^/XQ værdier giver t-værdier, der divergerer en del. Tages Xt/Xo = 0,3, fås en forskel på 6 min, mellem 16-8 mm prøven og 4-2 mm prøven. Side 73
Ved aflæsning af
kurverne fås følgende halveringstider Ud af tabellen
ses, at tJ/2tJ/2 generelt vokser med faldende
For visse fraktioner må evt. begge værdier anvendes, hvis det det ene øjeblik er flyvende materiale (vertikalstabiliteten) og i det næste fastliggende (horisontalstabiliteten). Det bør bemærkes, at de i artiklen behandlede fraktioner til halveringstidsbestemmelse alle er for store til at være saltationsmateriale af betydning. Man bør derfor
medtage endnu mindre størrelsesfraktioner, Side 74
b. De opnåede kurver, se fig. 3abc, er ikke rette linier, som man kunne forvente udfra de teoretiske beregninger side 71 (b). Det bemærkes derimod, at næsten alle kurverne er krumme og på en sådan måde, at hældningen på tangenten er faldende gennem kurven. Dette kan skyldes
flere ting, hvoraf skal nævnes: ikke samme
stabilitet, det vil sige, at de ustabile let bukker
2. aggregatstørrelsesfordelingen inden for samme fraktion begunstiger den aftagende hældning på kurverne ved at have mange små og få store aggregater. Det omvendte tilfælde kan også tænkes, og da vil den modarbejde krumningen af kurven. 3. aggregatmængden og dennes ændring i sigten. En for stor aggregatmængde i sigten vil hæmme sigtepåvirkningen i begyndelsen, da aggregaterne vil være låst fast. Denne effekt vil nærmest modarbejde formen på kurverne, da den formindsker destruktionen i begyndelsen. Ved færre aggregater i sigten kan der opstå en gnidningseffekt, hvor aggregaterne foruden afslibning mod sigten også sliber 4. store
aggregater vil ved sammenstød med mindre 5. aggregaterne
bliver under sigtningen afrundede. tagende hældning gennem kurverne, og samtidigt analysens største svaghed. Ved sigtningens begyndelse er aggregaterne ofte kantede, og disse spidser og hjørner vil først blive angrebet, da kraften her angriber på et meget lille areal. Resultatet er, at aggregaterne hurtigt bliver gformede dermed mere resistente. Den aftagende hældning gennem kurven bliver derfor mere et udtryk for en formændring, end for at aggregaterne, der er tilbage, var mere stabile fra begyndelsen. Denne afrundingseffekt er størst ved den horisontale sigtning, da aggregaterne der slides ved rulning og glidning, hvorimod vertikalsigtningen foruden horisontalsigtningseffekterne også deler og flækker aggregaterne, men dog i så ringe mængde, at afrundingen heller ikke der undgås. Da kurverne systematisk buer, vil det være umuligt at beskrive aggregatsøndefdelingen ved hjælp af en ret linie. Man må derfor forkaste tesen: »Mængden af aggregater der destrueres pr. tidsenhed er proportional med den tilstedeværende intakte aggregatmængde« og gå over til en matematisk analyse af problemet. Matematisk analyseTil beskrivelse af de fremkomne kurver, fig. 2 abc, er valgt polynomiet y= axb. Konstanterne aogber blevet bestemt ved mindste kvadraters metode, y er i dette tilfælde mængden af destrueret materiale målt i procent af totalaggregatmængden, og ikke som i fig. 2 den intakte aggregatmængde. Resultatet af analysen ses i nedenstående Af fit kvotienternen ses, at tilpasningen med polynomiet er brugbar, dog er der en tendens til ensidig afvigelse, og formlen er derfor ikke velegnet til extrapolation. Udtrykket y=a+ bcx, hvor (o, i) ikke medtages i beregningen, har vist sig velegnet til beskrivelse af en del af kurverne uden ensidig afvigelse, og kan da bruges til at extrapolere tl/2. AggregatstørrelsesfordelingenFor at klarlægge de to prøvers aggregatstørrelsesfordeling er denne udtegnet som histogrammer og kumulative kurver,fig. 5 (a, b, c). Som ordinatakse i histogrammerne haves sigtefraktionens procentvise vægtandel af totalaggregatmængdenstørre end 125u og som abscisseakse sigtefraktionerneafsat i phi. Det samme er tilfældet for de kumulative kurver, hvor ordinataksen er den summerede vægt. Af tabel i ses, at textursammensætningen af jorden i de to prøveudtagningshøjder minder en del om hinanden,og Side 75
den,ogisær bør man bemærke det høje lerindhold i beggejorde. Derfor vil det være rimeligt at tolke forskellene i størrelsesfordelingen som et resultat af de to prøvers forskellige humusindhold. Som mål for
aggregatstørrelsen kan man anvende dso dgo ( o- 15) =
-1,2 phi Det ses, at den
humusholdige prøves aggregater gennemsnitligt Som et mål for
jordens aggregering kan man anvende I dette tilfælde
fås: o- 15 = 6,1 phi
Det ses, at underjorden er mere aggregeret end overjorden. Man må gøre sig klart, at metoden er behæftet med den fejl, at der er en nedre begrænsning på aggregaterne (i dette tilfælde 125^), medens texturanalysen medtager alle partikler uanset størrelse. Man må derfor regne med andre aggregeringsværdier, hvis aggregater ned til 37// medtages. Det er i det hele taget et problem, hvor små aggregater man skal medtage i sine sigtninger. De fleste laboratorier har normalt kun sigter ned til størrelsen 63^ eller 37^. Ved aggregatanalyser, hvor hovedparten af elementerne er mindre end 3 y//, må specialapparatur udvikles. I det aktuelle tilfælde vil en 37^ sigte være tilstrækkelig, da kun få procent af prøverne er mindre, og endda kun en del af det er aggregater. Da de gunstigste dyrkningsvilkår opnås med et såbed, hvor der er mange aggregater mellem 3 mm og 0,5 mm, ses af de kumulative kurver, at humusindholdet har en gunstig indvirkning på størrelsesfordelingen. Dette emne vil blive videreført under del 2, hvor endnu to humusholdige prøver vil blive inddraget i analysen. DEL 2I dette afsnit
vil en samlet aggregatanalyse blive beskrevet, a. den samlede
aggregatstabilitet målt ud fra ændringen b.
humusindholdets indflydelse på aggregatstørrelsesforde-
lingen. c. de fine
fraktioners modstand mod vinderosion. Analysens udførelseDer bliver ved denne analyse arbejdet med størrelsesklasserne i phi og med mindste sigte 4,75 phi (37^1). Fremgangsmåden ved forbehandlingen af prøven er den samme som under del i. Efter at aggregaterne er udskilt i de forskellige størrelsesklasser, halveres aggregatmængden i hver af fraktionerne. Den ene halvdel bruges som nulprøve og skal vise aggregatstørrelsesfordelingen ved tiden o. Den anden del behandles ved horisontalsigtning i femten minutter. Sigterne er stablet med faldende sigtediameter nedefter, hvorved opnås, at de destruerede aggregater ikke bliver til bundmateriale som under del i, men til aggregater i mindre fraktioner. Vægten af materialet i de forskellige fraktioner bestemmes for både nul prøven og 15 min. prøven, og efter udvaskning af aggregaterne ligeledes vægten af enkeltkornene. Derved sikres, at aggregatmængden i prøverne kan bestemmes. Til analysen
anvendtes to prøver fra Karlstrup-profilet Side 76
På baggrund af sigtningerne er kurverne i fig. 6 udtegnet. Udregningen af de anvendte procenter til kurverne (minus kurve 15 min. b fig. 6c) er sket på baggrund af den aktuelle aggregatmængde, som er større end 37/4, dvs. at de kumulative kurver vil vise 100% ved 37|x uanset aggregatmængdens størrelse. For kurverne 15 min. b i fig. 6c er procenterne regnet ud på baggrund af nulprøvens aggregatmængde større end 37^. Slutprocenten på den kumulative kurve bliver da den procentdel af aggregatmængden gatmængdenstørre end 37^, der stadig er større end den værdi efter 15 minutters sigtning. Den NØ-SV gående skravering viser det procentvise nettotab ved sigtningen og den NV-SØ gående skravering den procentvise nettogevinst. Det ses, at for begge prøver har der været tab i de store fraktioner, hvilket er naturligt, da de i fraktionen destruerede aggregater ikke kan erstattes »ovenfra«. Det bemærkelsesværdige er, at prøven 40-55 cm har et nettotab helt ned i 4-2 mm klassen, dvs. tab i fire fraktioner, hvorimod 25-35 cm prøven kun har tab i to fraktioner. Ud fra histogrammerne kan man endvidere se, at fordelingskurvens form i 25-35 cm Prøy en er uforandret efter sigtningen, blot er maximumværdien (prøvens modalklasse) forskudt en phi til højre. Derimod forandrer formen sig for 40-55 prøvens vedkommende og bliver nærmest udjævnet. Dette kan skyldes det lavere ler- og humusindhold i den sidstnævnte prøve. a. Som mål for den samlede aggregatstabilitet kan man anvende den procentdel af aggregaterne, der hverken er nettotab eller -gevinst. Rent visuelt er den samlede aggregatstabilitet det uskraverede område i histogrammerne. I det aktuelle tilfælde bliver resultatet: 25-35 : 66,9
4°-55 :
5J>B5J>8 Et andet mål for den samlede aggregatstabilitet er den relative forskydning i dso eller forskydningen målt i phi. Denne forskydning måles mellem kurverne o min. og 15 min a i fig. 6c. Opgives denne værdi sammen med slutværdien på kurve 15 min. b fig. 6c fås et mål for ndringen aggregatstørrelsesfordelingen større end 37^ samt hvor stor en del af aggregaterne, der »overlevede behandlingen«. I det aktuelle tilfælde måles forskydningen i phi og resultatet er: 25-35 cm: 1,0 phi
80,5% Resultater diskussionJo større phi-forskydningsværdier, jo dårligere er den samlede aggregatstabilitet. Den sidstnævnte metode har den fordel, at værdien giver udtryk for, hvor mange gange aggregaterne gennemsnitligt er halveret, da en phi svarer til en halvering. En svaghed ved metoden er, at den ikke tager højde for lokale ændringer i yderfraktionerne. Af det tænkte eksempel fig. 7 ses, at en prøve kan have samme dsophi-værdi før og efter en behandling, men aggregatstabiliteten har ikke været 100 %, dvs. at i dette tilfælde vil dsophi-metoden overvurdere den samlede aggregatstabilitet. b. Det vil i dette afsnit blive forsøgt at finde en sammenhæng mellem prøvens humusindhold, aggregaternes modalklasse og mængden af aggregater mellem 3 mm og 0,5 mm. For at kunne sammenligne de fire prøvers aggregatstørrelsesfordeling må disse være defineret inden for det samme størrelsesområde, dvs. at 125^ i det aktuelle tilfælde må blive den nedre grænse, da det var den mindste sigte under del i. Der må altså udregnes nye procenter for 25-35 °S 4°~55 prøvernes størrelsesfordelinger. Ved at gøre størrelsesfordelingerne sammenlignelige bliver også aggregatmængden mellem 3 mm og 0,5 mm sammenlignelig. Det vil derfor være muligt ud fra de forbehandlede prøvers aggregatstørrelsesfordelinger at finde sammenhængen mellem prøvens humusindhold, prøvens modalklasse og mængden af aggregater mellem 3 mm og 0,5 mm udregnet i procent af totalaggregatmængden (Tabels). Det kan altså konstateres, at selv hvor der er texturforskelle (90-110 2i % ler, 40-55 12 % ler) har humusindholdet en kraftig indvirkning på aggregatstørrelsesfordelingen., således at modalklassens størrelse bliver omvendt proportional med humusindholdet. Det ses endvidere, at med øget humusindhold stiger antallet af aggregater i de fraktioner, der befordrer en god plantekultur. De ovennævnte betragtninger gælder kun inden for de i tabellen opgivne humusprocenter, dvs. at man ikke kan slutte ud fra det ovennævnte, at aggregatstørrelsesfordelingen vil blive endnu mere gunstig ved at øge humusprocenten til f.eks. 8. c. Ud fra resultaterne i tabel 5 har man en mulighed for at vurdere de fine fraktioners vinderodibilitet, hvis prøverne som overfladelag var exponerede for vindens virke. Fra studier om muldflugt på danske jorde i begyndelsen af 6oerne ved man, at vinderosion ikke vil forekomme i nævneværdigt omfang på marker, hvis primærpartikler er bundet i aggregater større end 0,5 mm. Det vil derfor være gunstigt, om de vindexponerede og derved truede jorde har deres primærpartikler mindre end 500^ bundet i aggregater større end 500^. Det er vigtigt at beskytte jordens fine fraktioner mod afblæsning, da det er på kolloidernes adsorptionskomplexer, næringsionerne er beliggende, samt at de finere fraktioners elementer ofte danner rammerne om de porer, der tilbageholder det tilgængelige vand. For en jord med en god struktur mod vinderosion vil brøken (vægt% textur { 0,5 mm — vægt% aggregater (0,5 mm) / vægt % textur (0,5 mm være nær i. Størrelsen betegnes VE (0,5. Man bør bemærke, at værdien ikke beskriver den absolutte mængde, der er truet, men er et relativt mål for hvor stor en del af fraktionerne (textur) mindre end 0,5 mm, der er i fare for afblæsning. Man kan teoretisk tænke sig, at der vil være en større afblæsning mængdemæssigt fra en mark, der har en højere VE (0,5 end en anden, hvis den sidstnævnte f.eks. er en meget mager jord. Hvis man er interesseret i at bestemme, hvor stor en del af lerfraktionen, der er bundet i aggregater større end 0,5 mm, må man lave en texturanalyse separat på aggregaterne større og mindre end denne værdi. Man kan også ud fra aggregatstørrelsesfordelingskurverne danne sig et indtryk af, hvor stor en del af aggregatmængden, der vil være truet under en storm. Det er af uvurderlig betydning, at aggregaternes stabilitet er stor, således at aggregaterne større end 0,5 mm ikke under bombardement af saltationsmaterialet bliver gjort »egnede« til afblæsning. Denne stabilitet kan bestemmes som ændringen iVE {0,5 værdien ved mekanisk påvirkning. I tabel 6 er VE
(0,5 udregnet for 0-15, 25-35 °g 4°~55 Side 78
Det ses, at 25—35 prøven er lidt mere resistent mod vinderosion end 40-55 prøven i begyndelsen, men senere under saltationsmaterialets bombardement vil forskellen blive forstørret kraftigt. Man må altså i det aktuelle tilfælde forvente en større procentvis afblæsning af finfraktionen fra 40-55 prøven end fra 25-35 prøven, hvis de andre faktorer, der betinger soil erosion, er ens for de to jorde. SUMMARYA method is discussed to describe the stability of soil aggregates by dry sifting. The analysis is divided into two parts: i) a total measurement including all size fractions for stability and 2) a separate investigation of the stability of the individual size fractions. The soil analysed
is a cambisol taken at Karlstrup Kalkbrud At the analysis of the separate size fractions the following hypothesis was put forward: the amount (weight) of aggregates destroyed per unit of time is proportional to the present, intact amount of aggregates investigated. This was not true, as the aggregates by sifting passed through a change of shape and became more resistant. The size distribution of aggregates, the different stability within the size fractions, and the change of quantity of material in the sieve through time can all have an influence on the progress of the destruction of aggregates. The mathematical analysis indicated that the polynomium Y= axb described the produced diagrams to satisfaction. However, it is not possible to use the formula for extrapolation because of a certain trend to systematic deviation. The bisecting time for the different size fractions can be used as a unit of measurement for the stability of the aggregates. However this value is not synonymous, i.e 2 samples may have the same bisecting time (tVz value), but not the same stability. The influence of the humus content on the distribution of size of aggregates is discussed by making a comparison between the unit content of the four samples and their size distribution of aggregates. At increasing contents of humus (up to 5 %} the size - distributions modal-class 32-16 mm was falling to 4-2 mm, at the same time would the amount of aggregates between 3 and ilt mm rise considerably. The last mentioned development is favourable for a seed bed. At the total analysis the stability of the aggregates was measured, as the amount of aggregates from which there were neither netloss nor netgain for the size fractions at the sifting. Another method is to determine the displacement of dso (measured in phi) at the sifting, and at the same time indicate, which percentage of the aggregates is still above the minimum The fine fraction's resistence to wind erosion can be evaluated on the basis of size distribution and stability of the aggregates after a dry sifting. As particles above 0.5 mm are only slightly movable by wind it is possible to estimate to which extent the fine fractions are exposed to erosion. The formular: shows the share
of particles less than 0.5 mm that is aggregated LITTERATURAslyng, H. C.
(1968): Klima, jord og vandbalance i jordbruget. Aslyng, H.
C.: Øvelsesvejledning i kulturteknik. KVL København.
Baver,
Gardner & Gardner (1972): Soil Physics. S. 130-229,
Koenigs, F.
F. R. & Schuffeien, A. C.: Mathematical treatment
of
aggregate-stability determination. 7. int. cong. soil
sei. vol. i. marker.
Geografisk tidsskrift bd. 59.
|