Veteen dispergoituvan saven aggregaattijakauma, stabiilisuus ja vapautuminen kahdella subtrooppisella Oxisolilla

SOILS AND PLANT NUTRITION

Aggregaattijakauma, stability and release of water dispersible clay for two subtropical Oxisols

Distribuição de agregados, estabilidade e liberação de argila dispersa em água para dois Latossolos subtropicais

Antonio Carlos de AzevedoI, *; Darrel Gene SchulzeII

I UFSM – Depto. de Solos, Prédio 42 – 97105-900 – Santa Maria, RS – Brasilia
IIPURDUE UNIVERSITY – Agronomy Dept, Lilly Hall of Life Sciences, 915 W. State Street, West Lafayette, IN 47907-2054, USA.

ABSTRACT

Maaperän käyttäytyminen sen hiukkasten hajaantumisen ja kasautumisen suhteen on erittäin tärkeää ympäristön ja maatalouden maaperän toimintojen kehittämisen kannalta. Tässä tutkimuksessa selvitettiin, miten maankäyttö vaikuttaa aggregaattien jakautumiseen ja stabiilisuuteen ja miten vesihajoavan saven (Water Dispersible Clay, WDC) vapautuminen liittyy disaggregoitumiseen Oxisoleissa subtrooppisesta Brasiliasta. Näytteitä kahdesta Oxisolista, jotka oli kerätty kolmelta syvyydeltä alueilta, joilla ei ollut maanmuokkausta (NT), tavanomaista maanmuokkausta (CT) ja kotoperäistä kasvillisuutta (NV), ravisteltiin 250 ml:n muovipulloissa enintään 27 tunnin ajan. Kiviaineksen massa mitattiin viidessä kokoluokassa, jotka vaihtelivat 53-2000 µm:n välillä. Suurin osa yli 500 mm:n kokoisista aggregaateista hävisi ensimmäisten 7,5 tunnin ravistelun aikana samanaikaisesti WDC:n vapautumisen lisääntymisen kanssa ilman, että maaperän suspensioiden pH-arvo ja sähkönjohtavuus muuttuivat, eikä pienempien aggregaattien määrä lisääntynyt. Näissä maaperissä ei siis ole aggregaattihierarkiaa, ja WDC:n vapautuminen johtui 500-2000 mm:n aggregaattien rikkoutumisesta. Maankäyttö vaikuttaa kunkin kokoluokan aggregaattien massaan, mutta aggregaattien vakaus riippuu niiden koosta, ei maankäytöstä.

Key words: land use, disaggregation, soil tillage, soil management

RESUMO

ResumO

Soil behavior in relation to its aggregate and dispersion state is of high agricultural and environmental relevance. Tässä tutkimuksessa selvitettiin, miten eri käyttötarkoitukset vaikuttavat aggregaatin jakautumiseen ja vakauteen ja miten veteen dispergoituneen saven (WDC) vapautuminen käyttäytyy suhteessa hajoamiseen brasilialaisissa subtrooppisissa latosoleissa. Näytteitä kahdesta latosolista, jotka kerättiin kolmelta syvyydeltä paikoilta, joilla ei ollut kylvöä (NT), tavanomaista (CT) ja alkuperäistä kasvillisuutta (NV), sekoitettiin 250 ml:n muovipulloissa enintään 27 tunnin ajan. Kiviaineksen massa mitattiin viidessä kokoluokassa 53-2000 µm. Suurin osa yli 500 mm:n kokoisista aggregaateista hävisi ensimmäisten 7,5 tunnin sekoittamisen aikana, mikä johti WDC:n vapautumisen lisääntymiseen ilman, että suspension pH- ja sähkönjohtavuusarvot olisivat muuttuneet, eikä pienempien aggregaattien massa olisi kasvanut. Näin ollen näissä maaperissä ei ole aggregaattihierarkiaa, ja WDC-päästöt johtuivat 500-2000 mm:n kokoisten aggregaattien hajoamisesta. Käyttötyyppi vaikutti kunkin kokoluokan kiviaineksen massaan, mutta kiviaineksen stabiilius on riippuvainen kiviaineksen koosta, ei käyttötyypistä.

Avainsanat: maankäyttö, disaggregaatio, maanmuokkaus, maaperän hoito

TULOKSET

Koska maaperän aggregaatit ovat dynaamisia ja reagoivat nopeasti ympäristömuutoksiin, kiinnostus niitä kohtaan maaperän laatuindikaattoreina on kasvussa (Caravaca et al., 2004; Boix-Fayos et al., 2001). Lisäksi maaperän hiukkasten kyky suojata orgaanista ainesta hapettumiselta antaa mahdollisuuden lähestyä maaperää hiilinieluna ilmakehään (Torn et al., 1997). Maankäytön muutoksilla, jotka johtuvat puunkorjuusta (Yanai et al., 2003) ja/tai metsästä maatalouteen siirtymisestä, on laaja vaikutus aggregaatioon (Carpenedo & Mielniczuk, 1990; Perin et al., 2003) ja hiilidynamiikkaan trooppisissa, happamissa maaperissä (Leite et al., 2004; Zinn et al., 2005). Aggregaattihierarkian puuttuminen (Oades & Waters, 1991) ja rauta- ja alumiinioksidien suuri määrä aiheuttavat sen, että säätyneistä Oxisoleista peräisin olevien aggregaattien stabiilisuus ja koko vaativat erityistutkimuksia. Brasiliassa tällaisilla mailla on yleensä hyvin stabiileja aggregaatteja, joiden halkaisija on alle 2 mm. Tämän seurauksena saviset Oxisolit käyttäytyvät kuin keskirakeiset maat, mikä mahdollistaa maataloustoimien, kuten maanmuokkauksen tai sadonkorjuun, suorittamisen pian sateen jälkeen (Buol & Eswaran, 2000). Näiden maaperien osalta ei tunneta hyvin mahdollisuutta vapauttaa veteen dispergoituvaa savea (Water Dispersible Clay, WDC) disaggregoitumisen seurauksena (Azevedo & Bonumá, 2004), ja tällainen tieto voi auttaa hallitsemaan näitä maaperiä paremmin maankäytön muutosten aikana, ja se on tärkeää arvioitaessa maataloudessa käytettävien rikkakasvien torjunta-aineiden, tuholaismyrkkyjen ja muiden ksenobioottisten yhdisteiden liikkuvuutta ympäristössä (Seta & Karathanasis, 1996; Bertsch & Seaman, 1999).

Tämän tutkimuksen tavoitteena oli arvioida hoitokäytäntöjen pitkän aikavälin vaikutusta halkaisijaltaan alle 2 mm:n kokoisten aggregaattien jakautumiseen ja pysyvyyteen sekä WDC:n vapautumiseen hajoamisen aikana kahden Etelä-Brasiliasta peräisin olevan Oxisolin osalta.

MATERIAALI JA MENETELMÄT

Kaksi valittua kohdetta sijaitsivat Sul Riograndensen ylätasangolla, joka on Serra Geral -muodostuman laavavirtojen päälle kehittynyt geomorfinen alue (kuva 1). Santo Ângelon alueen maaperä oli basalttikivistä kehittynyttä Typic Haplorthoxia, joka kattaa 7,26 prosenttia Rio Grande do Sulin osavaltion pinta-alasta ja esiintyy 200-400 metrin korkeudella merenpinnasta. Ilmasto on Köppenin järjestelmän mukaan Cfa-ilmastoa, jonka vuotuinen keskilämpötila on 19,5 ºC ja keskimääräinen sademäärä 1 850 mm vuodessa (BRASIL, 1973). Näytteet tavanomaisesta (CT) ja muokkaamattomasta (NT) järjestelmästä kerättiin Santo Ângelossa, Rio Grande do Sulin osavaltiossa (28º16′ eteläistä leveyttä, 54º13′ läntistä pituutta, noin 280 m korkeudella merenpinnasta). Maanmuokkauskoe perustettiin vuonna 1979 pellolle, joka oli ollut vehnä-soija-vuoroviljelyssä vuodesta 1964 (Dalla Rosa, 1981). Näytteet häiriintymättömästä, metsäisestä maaperästä kerättiin lähimmältä alkuperäisen metsäkasvillisuuden suojellulta alueelta (28º12′ eteläistä leveyttä, 54º13′ läntistä pituutta) noin 15 kilometrin päässä koealueelta. Passo Fundon alueen maaperä oli tyypillistä Haplohumox-maata, joka on kehittynyt basaltin ja hiekkakiven sekoituksesta ja jota esiintyy 460-700 metrin korkeudella merenpinnasta. Ilmasto on Cfa1-ilmastoa, vuotuinen keskilämpötila 18ºC ja keskimääräinen sademäärä 1 750 mm vuodessa. Näytteet kerättiin Passo Fundossa, Rio Grande do Sulin osavaltiossa (28º14′ eteläistä leveyttä ja 52º24′ läntistä pituutta) vuonna 1983 perustetussa kokeessa (Kochham & Denardim, 1997). Alkuperäisen kasvillisuuden alaiset maaperänäytteet otettiin koeaseman metsäiseltä suojelualueelta 0,5 kilometrin etäisyydeltä koealueesta.

Kummastakin koealasta kerättiin satunnaisesti kolme näytettä käsittelyistä NT ja CT koekentältä, ja kolme näytettä alkuperäisestä kasvillisuudesta (NV) kerättiin myös satunnaisesti metsistä. Koekentällä ja metsissä olevien keräyspisteiden väliset etäisyydet pyrittiin pitämään samankaltaisina (5-10 m). Tässä tutkimuksessa NT:stä, CT:stä ja NV:stä käytettiin nimitystä ”maankäyttö”. Näytteet 0-5 ja 10-15 cm:n syvyydestä kerättiin pienistä, noin 0,3 × 0,3 × 0,20 m:n syvyisistä kuopista, ja näytteet 40-60 cm:n syvyydestä kerättiin kauhakuopalla. Ilmakuivatut näytteet murskattiin kevyesti ja läpäistiin 2 mm:n seulan läpi (hieno maa-aines). Kunkin maankäyttömuodon kolme kenttäkertaa kustakin paikasta yhdistettiin, sekoitettiin perusteellisesti ja varastoitiin sitten suljettuihin muovipusseihin. Näin ollen laboratorioanalyysit tehtiin yhdestä yhdistelmänäytteestä kustakin kolmesta syvyydestä, kolmesta maankäyttömuodosta ja kahdesta maa-aineksesta, eli yhteensä kahdeksastatoista näytteestä.

Orgaaninen hiili määritettiin mädättämällä K2Cr2O7:ssa ja titraamalla Fe(NH4)2(SO4) 2.6H2O:lla, ja hiukkaskokojakauma määritettiin pipettimenetelmällä sen jälkeen, kun se oli dispergoitu 6 %:n NaOH:lla (EMBRAPA, 1997). Molemmat analyysit tehtiin kahdesti kullekin yhdistelmänäytteelle.

Sekä WDC:n että aggregaatin stabiilisuuden määritys komposiittinäytteistä tehtiin kolme kertaa (kolme ajoa) WDC:n määrityksen standardimenetelmällä (USDA, 1996). Lyhyesti viisi 10 g:n alikvoottia hienoa maa-ainesta punnittiin ja sijoitettiin 250 ml:n muovipulloihin. Pullot täytettiin 175 ml:lla deionisoitua (DI) vettä ja ravisteltiin (120 liikettä minuutissa, 4 cm:n vaakasuora siirtymä) 0, 3,75, 7,5, 15 ja 27 tunnin ajan (yksi aliquot jokaisella ajanjaksolla). Tämän jälkeen suspensiot kaadettiin viidestä seulasta koostuvan pesän läpi, jotka olivat 1000 µm, 500 µm, 250 µm, 106 µm ja 53 µm. Hajonnut savi ja siltti pestiin varovasti seuloissa olevasta maaperästä DI-vedellä ja kerättiin 1 litran sylintereihin WDC:n määrittämistä varten. Seuloihin jäänyt maa-aines kuivattiin 110º:ssa 24 tuntia ja punnittiin. Seulojen pesän läpi mennyt maa-aines kerättiin yhden litran sylintereihin ja suspendoitiin uudelleen WDC-mittausta varten pipettimenetelmällä (USDA, 1996). Liukenemisen oletettiin olevan vähäistä, koska näytteet, jotka olivat altistuneet vain nopealle kostutukselle, olivat lähes täysin aggregoituneita (0 tunnin ravisteluaika kuvissa 2 ja 3, aggregaattien summa).

Kunkin kokoluokan aggregaattien massa, Ai, laskettiin seuraavasti:

jossa Bi on maa-aineksen kuiva massa kokoluokassa i, Ci on hiekan kuiva massa kokoluokassa i ja D on maa-aineksen alkuperäinen uunikuiva massa. Ci oli kunkin yhdistelmänäytteen hiekkapitoisuuden kolmen edellisen määrityksen keskiarvo. Hiekka vähennettiin yhtälön osoittajasta, jotta yksittäisiä hiekanjyviä ei laskettaisi aggregaateiksi. Seuraavassa keskustelussa käytetään lyhyyden vuoksi termiä ”aggregaatit” viittaamaan Ai, , ja ”näyte” viittaamaan yhdistettyyn näytteeseen.

Kiviainesjakauma analysoitiin jaettuna koealasuunnitelmana, jossa oli kolme käsittelyä (NT, CT, NV), joissa oli kolme lohkoa (kukin kolmesta ajosta) toistettuna ajan kuluessa. Kunkin maalajin, syvyyden ja kokoluokan tiedot analysoitiin erikseen. Ravisteluaikaa pidettiin jaetun lohkokokeen kokonaisyksikkönä ja maankäyttöä alayksikkönä. Maamassan prosenttiosuutta koskeville tiedoille tehtiin neliöjuurimuunnos ennen varianssianalyysia, jotta saavutettiin virheen varianssin homogeenisuus. Virhe (a) yhdistettiin virheeseen (b), koska se ei ollut useimmissa tapauksissa merkitsevä (P = 0,25). Ajan ja maankäytön ja ajan väliset vaikutukset jaettiin ortogonaalisiksi polynomikontrasteiksi. Riippuvien muuttujien keskiarvojen regressioita maankäytön ja ajan funktiona seurasi varianssianalyysi (ANOVA), jossa regressiomalli määritettiin merkittävien käsittelyvaikutusten ja kontrastien perusteella. Taulukossa 2 esitetään ANOVA:n avulla tehtyjen regressiomallien vertailu kunkin maankäytön (ravisteluajankohtien) erittelyn välillä.

TULOKSET JA JOHTOPÄÄTÖKSET

Koska hiekkapitoisuus vähennettiin kussakin kokoluokassa pidättyneen maa-aineksen massasta (yhtälö ), hiekkaisemmassa tyypillisessä Haaplohumoksessa (taulukko 1) oli vähemmän maa-aineksen aggregaatteja kussakin kokoluokassa kuin tyypillisessä Haaplohumoksessa.

Ennen ravistelua (ravisteluaika 0 h; kuvat 2 ja 3) Typic Haplorthoxista 73-91 % ja Typic Haplohumoxista 54-69 % esiintyi aggregaatteina. Vain kahta poikkeusta lukuun ottamatta regressoidut erittelymallit erosivat (P = 0,01) kolmen maankäyttömuodon välillä 0-5 ja 5-10 cm:n näytteenottosyvyyksillä, mutta maankäytöllä oli vähemmän vaikutusta erittelymalleihin 40-60 cm:n syvyydellä (taulukko 2).

Joitakin yhteisiä käyttäytymismalleja voitiin havaita molemmissa maaperissä (kuvat 2 ja 3): 106-53 µm:n ja 250-105 µm:n kokoisia aggregaatteja oli hyvin vähän maaperässä alkuperäisen kasvillisuuden alla 0-5 cm:n ja 10-15 cm:n syvyyksillä; halkaisijaltaan yli 500 mm:n kokoiset aggregaatit hallitsivat sekä 0-5 cm:n että 10-15 cm:n syvyyksillä, mutta 40-60 cm:n syvyydellä halkaisijaltaan 2000-1000 µm:n kokoisia aggregaatteja esiintyi vähiten; ja vaikka käppyröiden keskinäinen pinoamisjärjestys yksittäisessä kuvaajassa muuttuikin kokoluokasta ja syvyydestä riippuen, havaittiin, että käppyränmuodot olivat yleisesti ottaen samankaltaisia kunkin kokoluokan osalta.

WDC:tä vapautui vähän kaikissa näytteissä alkukostutuksen yhteydessä (0 tunnin ravistelu), mikä tukee olettamusta, että lahoaminen ei ollut merkittävä hajoamisprosessi koeolosuhteissa.

Suurimmat aggregaatit (2000-1000 µm ja 1000-500 µm) hajosivat nopeasti ensimmäisten 7,5 tunnin ravistelun aikana, mutta pienempien aggregaattien massa ei kasvanut samanaikaisesti, mikä osoitti, että tällaiset aggregaatit hajosivat pääasiassa primäärisiksi maahiukkasiksi eikä pieniksi aggregaateiksi (kuvat 2 ja 3). Näin ollen näillä mailla ei ollut aggregaattihierarkiaa (Oades & Waters, 1991). Koska WDC lisääntyi, kun yli 500 µm:n kokoiset aggregaatit hajosivat ilman pH:n muutosta suspensiossa, hajoaminen näyttää olevan tärkein WDC:n tuotantomekanismi.

Aggregaattihierarkian puuttuminen mahdollistaa disaggregoitumisen kuvaamisen ensimmäisen kertaluvun prosessimallilla (Beare & Bruce, 1993; Parkin & Robinson, 1992; Olson, 1963):

jossa Ai,t on kokoluokkaan i kuuluvien aggregaattien massa ajankohtana t, Ai,0 on aggregaattien massa ajankohtana 0 ja k on kaarevuusparametri. Tässä mallissa ei oteta huomioon seulaan i lisättyjä kiviaineksia, jotka syntyvät sen yläpuolella olevassa seulassa i+1 olevien kiviainesten hajoamisesta, joten sitä voidaan käyttää vain maaperässä, jossa ei ole kiviaineshierarkiaa. Kaarevuusparametrin (k) oletettiin olevan hajoamisnopeusindeksi, ja mitä suurempi sen absoluuttinen arvo on, sitä pienempi on aggregaatin vakaus. Mitään tunnistettavaa mallia ei havaittu, kun k-arvot järjestettiin maankäytön mukaan. Kun k:n vaihteluväli piirrettiin kuitenkin aggregaatin kokoluokan mukaan, ilmeni kohtuullisen selvä suuntaus (kuvat 4 ja 5), joka osoittaa, että k-arvot ryhmittyivät peräkkäin suurempiin arvoihin kokoluokan kasvaessa. Tämä viittaa siihen, että maankäytöllä oli suurempi vaikutus kunkin kokoluokan maa-aineksen määrään (kiviainesjakauma, kuvat 2 ja 3) kuin kiviaineksen vakauteen.

Maankäytön vaikutus k-arvoihin voidaan päätellä arvojen hajonnasta kullakin syvyydellä ja kokoluokalla (kunkin viivan varrella olevat pisteet kuvissa 4 ja 5). Vaikka maankäytön vaikutus ei ole hallitseva, se on suurin (suurempi vaihteluväli) yli 500 mm:n kokoluokissa ja pintahorisontissa (kuvat 4 ja 5). Tämä on sopusoinnussa Oades & Watersin (1991) ehdottaman aggregaatiomallin kanssa, jossa suuret aggregaatit ovat enemmän riippuvaisia sienihyfoista ja hienoista juurista ja siten maankäytöstä, kun taas pienet aggregaatit riippuvat enemmän maaperän kolloidisista ominaisuuksista ja kemiasta.

Kun otetaan huomioon, että WDC:n vapautuminen liittyi läheisesti yli 500 mm:n kokoisten aggregaattien hajoamiseen, että tavanomainen ja muokkaamaton maanmuokkausjärjestelmä edistivät suurten aggregaattien määrän vähenemistä ja että aggregaattihierarkiaa ei ole olemassa, huomattavan suuri määrä savea voi mahdollisesti hävitä siirryttäessä metsänkäytöstä maatalouskäyttöön. Näiden tulosten mukaan suurten kiviainesten säilyttäminen on kaikissa tapauksissa olennaista WDC:n kasvun välttämiseksi. Tämän lisäksi savimineraalien ja orgaanisten kolloidien muodostamien kompleksien, jotka lisäävät hajoamispotentiaalia, pitäisi olla suurempia pintahorisontissa (Tombácz et al., 2004).

LÄHTEET

CAPES-Brasilialle (BEX 1316-96/6), Purdue Agricultural Research Programs -lehden artikkelinumero 2006-17943, ja A. Dallarosalle ja J. Beckerille COTRISAsta; J. E. Denardinille ja R. A. Kochhamille EMBRAPAsta.

AZEVEDO, A.C.; BONUMÁ, A.S. Partículas coloidais, dispersão e agregação em Latossolos. Ciência Rural, v.34, s.609-617, 2004.

BEARE, M.H.; BRUCE, R.R. Vesistabiloitujen aggregaattien mittausmenetelmien vertailu: vaikutukset ympäristövaikutusten määrittämiseen maaperän rakenteeseen. Geoderma, v.56, s.87-104, 1993.

BERTSCH, P.M.; SEAMAN, J.C. Characterization of complex mineral assemblages: implications for contaminant transport and environmental remediation. Proceedings of National Academy of Science USA, v.96, s.3350-3357, 1999.

BOIX-FAYOS, C.; CALVO-CASES, A.; IMESON, A.C.; SORIANO-SOTO, M. D. Maaperän ominaisuuksien vaikutus joidenkin Välimeren maiden aggregaatioon ja aggregaattien koon ja stabiilisuuden käyttö maan huonontumisen indikaattoreina. Catena, v.4, s.47-67, 2001.

BRASIL. Levantamento de reconhecimento dos solos do Estado do Rio Grande do Sul. Brasília: Ministério da Agricultura, SNLCS, 1973. 431p.

BUOL, S.W.; ESWARAN, H. Oxisolit. Advances in Agronomy, v.68, s.151-195, 2000.

CARAVACA, F.; LAX, A.; ALBALADEJO, J. Aggregaatin stabiilisuus ja hiiliominaisuudet hiukkaskokojakaumissa Espanjan puolikuivien alueiden viljely- ja metsämaissa. Soil and Tillage Research. v.78, s.83-90, 2004.

CARPENEDO, V.; MIELNICZUK, J. Aggregaattitilanne ja aggregaatin laatu purppuranvärisillä latosoleilla, joihin on sovellettu eri hoitojärjestelmiä. Brazilian Journal of Soil Science , v.14, s.99-105, 1990.

DALLA ROSA, A. Mekaaniset ja kulttuuriset käytännöt viljelyn vaurioittaman maaperän fysikaalisten ominaisuuksien palauttamisessa – Santo Ângelon maaperä (dystrofinen Latossolo Roxo). Porto Alegre: UFRGS/Faculdade de Agronomia, 1981. 136p. (Valtiotieteellinen väitöskirja).

EMBRAPA. Manual de métodos de análise de solos. 2. Rio de Janeiro: EMBRAPA, 1997. 212pp.

KOCHHAM, R.A.; DENARDIM, J.E. Vehnä-, soija- ja maissiviljelmien käyttäytyminen fosfaattilannoituksen suhteen no-till- ja tavanomaisissa maanmuokkausjärjestelmissä. In: CONGRESSO BRASILIERO DE CIÊNCIA DO SOLO, 26., Rio de Janeiro, 1997. SBCS, 1997.

LEITE, L.F.C.; MENDONÇA, E.S.; MACHADO, P.L.O.A.; FERNANDES FILHO, I.E.;NEVES, J.C.L.. Maaperän orgaanisen hiilen kehityssuuntausten simulointi Acrisolin maaperän orgaanisen hiilen osalta ilman maanmuokkausta ja kiekkokyntöä Century-mallin avulla. Geoderma, v.12, s.283-295, 2004.

OADES, J..M.; WATERS, A.G. Aggregaattihierarkia maaperässä. Australian Journal of Soil Research, v.29, s.815-828, 1991.

OLSON, J.S. Energian varastointi sekä tuottajien ja hajottajien tasapaino ekologisissa järjestelmissä. Ecology, v.44, s.322-331, 1963.

PARKIN, T.B.; ROBINSON, J.A. Log-normaalin datan analysointi. Advances in Soil Science, v.20, s.193-235, 1992.

PERIN, E.; CERETTA, C.A.; KLAMT, E. Tempo de uso agrícola e propriedades químicas de dois Latossolos do Planalto Médio do Rio Grande do Sul. Revista Brasileira de Ciencia do Solo, v.27, s.665-674, 2003.

SETA, A.K.; KARATHANASIS, A.D. Water dispersible colloids and factors influencing their dispersibility from soil aggregates. Geoderma, v. 74, s.255-266, 1996.

TOMBACZ, E.; LIBOR, Z.; ILLES, E.; MAJZIK, A.; KLUMPP, E.. Reaktiivisten pintakohtien ja humushappojen kompleksoitumisen merkitys savimineraalin ja rautaoksidipartikkelien vuorovaikutuksessa. Organic Geochemistry. v.35, s.257-267, 2004.

TORN, M.S.; TRUMBORE, S.E.; CHADWICK, O.A.; VITOUSEK, P.M.; HENDRICKS, D..M. Mineral control of soil organic carbon storage and turnover. Nature, v.389, s.170-173, 1997.

USDA. Soil survey laboratory methods manual. 3.0. Washington: NRCS, NSSC. 1996. 693p.

YANAI, R.D.; CURRIE, W.S.; GOODALE, C.L. Soil Carbon Dynamics after Forest Harvest: An Ecosystem Paradigm Reconsidered. Ecosystems, v.6, s.197-212, 2003.

ZINN, Y.; LAL, R.; RESCK, D.V.S. Texture and organic carbon relations described by a profile pedotrasnfer function for Brazilian Cerrado soils. Geoderma, v.127, s.168-173, 2005.