Introduction
작물의 생산성은 기상 환경, 품종 등 여러 요인에 의해 결정되며, 특히 토양의 특성과 양분 순환에 밀접한 관련이 있다 (Yamagishi et al., 2003). 작물의 생산성 증대를 위해 다양한 무기질 비료가 개발되었으나 지속적으로 과도한 무기질 비료를 시비할 경우 토양과 지하수가 오염되고, 이로 인한 작물의 생리적 장애 등 여러 가지 피해가 발생한다고 보고되었다 (Chang and Dregne, 1955; Bernstein, 1975; Gimeno-García et al., 1996; Vitousek et al., 2009; Sun et al., 2012). Ca, Mg, Na, SO42- 등의 양분 집적은 삼투압 구배에 영향을 끼쳐 작물의 수분 흡수를 어렵게 한다. 또한 타 양분과의 길항작용으로 이온 흡수를 저해해 생리장해를 유발하므로, 적절한 토양관리를 통한 토양 환경의 개선이 중요하다 (Chang and Dregne, 1955). 토양 환경 개선은 단순히 토양 내 양분을 많이 보유하는 것이 아니라, 작물의 건강한 생육과 수량 증대를 위해 토양 이화학적 특성과 미생물 활동 등의 생물학적 성질이 조화롭게 어우러진 상태를 의미한다 (Weaver et al., 1994; Yang et al., 2001; Kim et al, 2003; Hong et al., 2009; Kim et al., 2023).
토양 이화학적 특성 및 미생물 활성 등은 재배환경 및 토양 관리에 따라 영향을 받으며, 특히 유기농업은 농약과 무기비료의 사용을 금지하기 때문에 토양 내 양분의 형태가 관행 농가와는 다른 특성을 나타낸다 (Barak et al., 1997; Maeder et al., 2002; Pimentel et al., 2005; Liang et al., 2015). 고추 재배지의 유기농업 토양은 pH, 유기물 함량, 탄질율 (C/N ratio), 유효인산이 관행 대비 높으며, 전기전도도 (EC)는 낮은 수치를 보이는 등 농업 형태에 따라 토양의 화학적 특성의 차이가 보고된 바 있다 (Kim et al., 2020). 또한 국내 유기 농경지 토양 조사 결과, 관행 농경지보다 대부분 용적밀도가 낮고 pH, 전기전도도 및 urease와 같은 효소 활성이 비교적 높은 것으로 나타났다 (Kim et al., 2023).
한편, 유기농업은 유기농 육성 정책 등을 통해 급속도로 성장하고 있으며 (Gomiero et al., 2008), 특히 충청북도의 유기농 재배 면적은 1,634 ha로 2015년 이후로 유기재배 인증 면적이 지속적으로 증가하고 있다 (https://www.enviagro.go.kr). 하지만 관행농가의 경우 1999년부터 주기적으로 각 지역의 농업기술원과 공동으로 토양 실태 조사가 실시되고 있는 것과 달리 유기 농경지의 지역별 토양 조사는 미비한 실정이다 (Kim et al., 2019). 따라서 본 연구는 충북 지역의 유기 농경지에서 토양의 이화학적 및 생물학적 특성을 분석하여, 유기 농경지의 특성을 파악하고 토양 건강성 평가를 위한 기초 자료로 활용하고자 한다.
Materials and Methods
조사지역 선정 및 시료 채취
본 연구는 2023년 기준 국립농산물품질관리원에서 경지유형이 ‘밭’이면서 5년 이상 유기 인증을 받은 충북지역의 토양 시료를 116점 채취하여 토양의 화학적 및 생물학적 특성을 분석하였으며 (Table 1), 채취된 필지에 대한 정보는 Table 2에 제시하였다. 토양 시료는 퇴비가 포함되지 않도록 유기물층을 제거한 후 표토층에서 채취하였으며, 시료의 일관성을 확보하기 위해 필지 내 서로 다른 5곳의 시료를 혼합하여 실험에 사용하였다. 토양 화학적 특성 분석에 사용된 토양은 풍건한 후 2 mm 표준체를 통과한 시료를 이용하였으며, 생물학적 특성에 사용된 시료는 4°C 냉장 보관하여 습토 상태로 사용하였다.
토양 화학적 특성 분석
토양 화학적 특성은 NAAS (2010)의 종합검정실 분석 매뉴얼에 따라 수행되었다. 토양의 pH 및 전기전도도는 토양과 증류수의 혼합물을 1:5 비율로 진탕한 후, pH 미터 (Star A211, Orion, USA)와 전도도계 (VersaStar, Orion, USA)를 사용하여 측정하였으며, 전기 전도도는 측정값에 희석배수를 곱하여 값을 산출하였다. 유기물 (OM) 함량과 총 질소 함량은 원소분석기 (Vario max Cube, Elemental, Germany)를 이용하여 측정하였으며, 유기물 함량은 분석한 총 탄소 (total carbon, T-C) 함량에 변환 계수인 1.724를 곱하여 계산하였다. 유효인산 농도는 Lancaster 방법을 이용하여 비색 분석 (Cary3500, Agilent, USA)하였다. 교환성 양이온은 1 M 암모늄 아세테이트 (pH 7.0)로 추출한 후, 유도 결합 플라즈마 분광광도계 (5100 ICP-OES, Agilent, USA)로 정량하였다.
Table 1.
Description of sampling site.
Region
|
Section
No.
|
Year of Initial
Organic Certification
|
Location
|
Soil Texture
|
Sampling date
|
Latitude
|
Longtitude
|
Northern
|
Chungju
|
|
Cu1
|
2016
|
37°06
|
127°53
|
sandy loam
|
2022-10-25
|
|
Cu2
|
2016
|
37°06
|
127°53
|
sandy loam
|
2022-10-25
|
|
Cu3
|
2016
|
37°06
|
127°53
|
loamy sand
|
2022-11-09
|
|
Cu4
|
2009
|
37°08
|
127°49
|
sandy loam
|
2022-11-09
|
|
Cu5
|
2009
|
37°08
|
127°49
|
sandy loam
|
2022-11-09
|
|
Cu6
|
2009
|
37°09
|
127°49
|
sandy loam
|
2022-11-09
|
|
Cu7
|
2009
|
37°09
|
127°49
|
sandy loam
|
2022-11-09
|
|
Cu8
|
2009
|
37°08
|
127°49
|
sandy loam
|
2022-11-09
|
|
Cu9
|
2009
|
37°09
|
127°49
|
sandy loam
|
2022-11-09
|
Jecheon
|
|
Je1
|
2010
|
37°06
|
128°24
|
loamy sand
|
2022-10-25
|
|
Je2
|
2010
|
37°06
|
128°24
|
sandy loam
|
2022-10-25
|
|
Je3
|
2015
|
37°06
|
128°24
|
loam
|
2022-10-25
|
|
Je4
|
2015
|
37°06
|
128°24
|
loam
|
2022-10-25
|
|
Je5
|
2017
|
37°12
|
128°06
|
sandy loam
|
2023-09-25
|
|
Je6
|
2010
|
37°03
|
128°07
|
sandy loam
|
2023-09-25
|
Danyang
|
|
D1
|
2008
|
37°08
|
128°19
|
silt loam
|
2022-10-18
|
|
D2
|
2008
|
37°08
|
128°19
|
loam
|
2022-10-18
|
|
D3
|
2008
|
37°07
|
128°19
|
loam
|
2022-10-18
|
|
D4
|
2016
|
37°06
|
128°24
|
silt loam
|
2022-10-18
|
|
D5
|
2016
|
37°06
|
128°24
|
loam
|
2022-10-18
|
|
D6
|
2016
|
37°06
|
128°24
|
loam
|
2022-10-18
|
|
D7
|
2016
|
37°06
|
128°24
|
loam
|
2022-10-18
|
|
D8
|
2016
|
37°06
|
128°24
|
loam
|
2022-10-18
|
|
D9
|
2016
|
37°06
|
128°24
|
loam
|
2022-10-18
|
|
D10
|
2016
|
37°05
|
128°24
|
loam
|
2022-10-19
|
|
D11
|
2010
|
37°06
|
128°24
|
loam
|
2022-10-19
|
|
D12
|
2016
|
37°05
|
128°23
|
silt loam
|
2022-10-19
|
|
D13
|
2016
|
37°05
|
128°24
|
loam
|
2022-10-19
|
|
D14
|
2016
|
37°05
|
128°24
|
loam
|
2022-10-19
|
|
D15
|
2016
|
37°08
|
128°23
|
loam
|
2023-09-19
|
|
D16
|
2016
|
37°08
|
128°23
|
loam
|
2023-09-19
|
|
D17
|
2016
|
37°08
|
128°23
|
loam
|
2023-09-19
|
|
D18
|
2016
|
37°06
|
128°24
|
loam
|
2023-09-19
|
|
D19
|
2016
|
37°05
|
128°25
|
loam
|
2023-09-19
|
|
D20
|
2016
|
37°06
|
128°24
|
loam
|
2023-09-19
|
|
D21
|
2016
|
37°06
|
128°24
|
loam
|
2023-09-19
|
|
D22
|
2016
|
37°06
|
128°24
|
loam
|
2023-09-19
|
|
D23
|
2016
|
37°05
|
128°24
|
silt loam
|
2023-09-19
|
|
D24
|
2016
|
37°06
|
128°24
|
silt loam
|
2023-09-19
|
|
D25
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-19
|
|
D26
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-19
|
|
D27
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-19
|
|
D28
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-19
|
|
D29
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-19
|
|
D30
|
2015
|
37°08
|
128°20
|
silty clay
|
2023-09-20
|
|
D31
|
2016
|
37°05
|
128°22
|
silt loam
|
2023-09-20
|
|
D32
|
2015
|
37°07
|
128°20
|
silt loam
|
2023-09-20
|
|
D33
|
2015
|
37°08
|
128°19
|
loam
|
2023-09-20
|
|
D34
|
2015
|
37°08
|
128°18
|
loam
|
2023-09-20
|
|
D35
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-20
|
|
D36
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-20
|
|
D37
|
2016
|
37°05
|
128°24
|
loam
|
2023-09-20
|
|
D38
|
2016
|
37°06
|
128°23
|
loam
|
2023-09-20
|
|
D39
|
2016
|
37°06
|
128°23
|
loam
|
2023-09-20
|
|
D40
|
2016
|
37°06
|
128°23
|
loam
|
2023-09-20
|
|
D41
|
2016
|
37°06
|
128°23
|
silty clay loam
|
2023-09-20
|
|
D42
|
2016
|
37°06
|
128°24
|
loam
|
2023-09-20
|
|
D43
|
2016
|
37°05
|
128°23
|
loam
|
2023-09-20
|
|
D44
|
2016
|
37°06
|
128°24
|
silty clay loam
|
2023-09-20
|
|
D45
|
2016
|
37°05
|
128°23
|
loam
|
2023-09-20
|
|
D46
|
2016
|
37°05
|
128°23
|
loam
|
2023-09-20
|
|
D47
|
2016
|
37°05
|
128°23
|
loam
|
2023-09-20
|
|
D48
|
2016
|
37°06
|
128°23
|
loam
|
2023-09-20
|
|
D49
|
2016
|
37°05
|
128°23
|
loam
|
2023-09-20
|
|
D50
|
2009
|
37°05
|
128°22
|
loam
|
2023-09-20
|
Central
|
Cheongju
|
|
CJ1
|
2015
|
36°42
|
127°37
|
sandy loam
|
2022-11-01
|
|
CJ2
|
2015
|
36°41
|
127°38
|
loam
|
2022-11-01
|
|
CJ3
|
2016
|
36°42
|
127°37
|
loam
|
2022-11-01
|
|
CJ4
|
2016
|
36°41
|
127°37
|
sandy loam
|
2023-09-18
|
|
CJ5
|
2016
|
36°41
|
127°36
|
sandy loam
|
2023-09-18
|
|
CJ6
|
2016
|
36°41
|
127°35
|
sandy loam
|
2023-09-18
|
|
CJ7
|
2015
|
36°41
|
127°37
|
loam
|
2023-09-18
|
Jeungpyeong
|
|
JP1
|
2011
|
36°48
|
127°37
|
sandy loam
|
2023-10-05
|
|
JP2
|
2011
|
36°49
|
127°37
|
sandy loam
|
2023-10-05
|
Jincheon
|
|
JC1
|
2017
|
36°53
|
127°18
|
sandy clay loam
|
2023-10-05
|
|
JC2
|
2011
|
36°56
|
127°26
|
sandy clay loam
|
2023-10-05
|
|
JC3
|
2011
|
36°56
|
127°26
|
sandy clay loam
|
2023-10-05
|
|
JC4
|
2011
|
36°56
|
127°26
|
sandy clay loam
|
2023-10-05
|
Geosan
|
|
GS1
|
2016
|
36°47
|
127°53
|
silt loam
|
2022-11-01
|
|
GS2
|
2016
|
36°47
|
127°53
|
loam
|
2022-11-01
|
|
GS3
|
2015
|
36°39
|
127°54
|
sandy loam
|
2022-11-01
|
|
GS4
|
2011
|
36°52
|
127°56
|
silt loam
|
2022-11-09
|
|
GS5
|
2014
|
36°39
|
127°51
|
sandy loam
|
2023-09-18
|
|
GS6
|
2014
|
36°38
|
127°51
|
sandy loam
|
2023-09-18
|
|
GS7
|
2007
|
36°38
|
127°46
|
silt loam
|
2023-09-18
|
|
GS8
|
2007
|
36°38
|
127°46
|
loam
|
2023-09-18
|
|
GS9
|
2012
|
36°46
|
127°42
|
loam
|
2023-10-05
|
|
GS10
|
2012
|
36°46
|
127°42
|
loam
|
2023-10-05
|
Eumseong
|
|
ES1
|
2008
|
36°53
|
127°38
|
loam
|
2023-10-05
|
|
ES2
|
2008
|
36°54
|
127°37
|
sandy clay loam
|
2023-10-05
|
|
ES3
|
2008
|
36°55
|
127°37
|
sandy clay loam
|
2023-10-05
|
Southern
|
Boeun
|
|
BE1
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE2
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE3
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE4
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE5
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE6
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE7
|
2017
|
36°23
|
127°52
|
loam
|
2022-11-01
|
|
BE8
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE9
|
2017
|
36°23
|
127°52
|
loam
|
2022-11-01
|
|
BE10
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE11
|
2017
|
36°23
|
127°52
|
sandy loam
|
2022-11-01
|
|
BE12
|
2017
|
36°23
|
127°51
|
sandy loam
|
2022-11-01
|
|
BE13
|
2017
|
36°23
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE14
|
2017
|
36°23
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE15
|
2017
|
36°23
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE16
|
2017
|
36°23
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE17
|
2017
|
36°22
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE18
|
2017
|
36°22
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE19
|
2017
|
36°23
|
127°51
|
sandy loam
|
2023-09-25
|
|
BE20
|
2017
|
36°22
|
127°52
|
sandy loam
|
2023-09-25
|
|
BE21
|
2017
|
36°23
|
127°51
|
sandy loam
|
2023-09-25
|
|
BE22
|
2017
|
36°23
|
127°49
|
sandy loam
|
2023-09-25
|
|
BE23
|
2017
|
36°23
|
127°48
|
sandy loam
|
2023-09-25
|
Okcheon
|
|
OC1
|
2017
|
36°21
|
127°50
|
sandy loam
|
2023-09-25
|
|
OC2
|
2017
|
36°21
|
127°50
|
sandy loam
|
2023-09-25
|
Table 2.
List of total number of certified organic soils for more than 5 years and sampling sites in this study.
Region
|
Sampling number of soils
|
Total number of certified
organic soils1 |
Ratio (%)
|
Northern
|
Chungju
|
9
|
37
|
24.3
|
Jecheon
|
6
|
13
|
46.2
|
Danyang
|
50
|
160
|
31.3
|
Central
|
Cheongju
|
7
|
49
|
14.3
|
Jeungpyeong
|
2
|
5
|
40.0
|
Jincheon
|
4
|
11
|
36.4
|
Geosan
|
10
|
29
|
34.5
|
Eumseong
|
3
|
16
|
18.8
|
Southern
|
Boeun
|
23
|
161
|
14.3
|
Yeongdong
|
-
|
5
|
0
|
Okcheon
|
2
|
17
|
11.8
|
Total
|
116
|
503
| |
토양의 생물학적 특성 분석
토양의 생물학적 특성은 토양 호흡량과 효소 활성도를 측정하였다. 토양 호흡량 분석은 300 mL 용기에 습토 100 g을 넣고, 0.5 M NaOH를 담은 vial을 용기 중앙에 고정한 후 밀폐하여 25°C에서 72시간 동안 배양하였다 (Bottomley et al., 2020). 배양이 끝난 후 0.5 M NaOH 용액에 페놀프탈레인 지시약을 추가한 뒤, 1 M HCl을 사용하여 용액이 무색으로 변할 때까지 적정하였다. 이후, 무색으로 변한 용액에 브로로크레졸 그린 지시약을 이용해 용액이 푸른색에서 노란색으로 변화할 때까지 0.2 M HCl로 적정하였다. 이 과정에서 소모된 0.2 M HCl의 양을 다음 식을 사용하여 CO2 발생량으로 변환하여 토양 호흡량을 산출하였다.
M : Hydrochloric acid concentration (0.2 mol L-1)
T : 0.2 M HCl Titrate volume (mL)
44 : CO2 molecular weigth (g mol-1)
S : Soil weight (kg)
day : Days of incubation
토양 효소 활성도는 총 4가지 항목 (dehydrogenase, Urease, β-Glucosidase, Phosphatase)에 대해 측정하였다. Dehydrogenase는 2,3,5-Triphenyl tetrazolium chloride (TTC)가 Triphenyl formazan으로 환원된 양을 정량하였다 (Casida et al., 1964). 습토에 CaCO3를 혼합한 후, TTC와 증류수를 섞어 입구를 밀봉하고 37°C에서 24시간 배양하였다. 배양한 용액은 탈지면으로 막은 깔때기를 이용해 여과하였으며, 이때 탈지면의 붉은색이 없어질 때까지 지속적으로 메탄올을 추가하였다. 여과액은 최종부피가 100 mL가 되도록 메탄올을 mass up한 후, UV/Vis Spectrophotometer를 이용해 485 nm 파장에서 비색 정량하였다. Urease는 100 mL 삼각플라스크에 5 g의 습토와 2.5 mL의 79.9 mM urea 기질을 섞어 37°C에서 2시간 배양하였다. 배양 후 2.5 mL의 증류수와 30 mL의 KCl을 첨가하여 30분간 교반한 후 여과하고, 여과액 1 mL와 증류수 4 mL, 0.3 M NaOH, sodium salicylate solution, 39.1 mM dichlorisocyanurate solution 2 mL를 30분간 혼합교반하여 690 nm 파장에서 비색 정량하였다 (Kandeler and Gerber, 1988). β-Glucosidase와 Phosphatase는 각각 ρ-nitrophenyl-β-D-glucopyranoside (PNG)와 sodium-p-nitrophenyl phosphate (pNPP) 기질을 첨가하여 배양한 후, 전환된 생성물인 p-nitrophenol (PNP)을 400 nm 파장에서 비색 정량하였다 (Tabatabai and Bremner, 1969; Eivazi and Tabatabai, 1988). 즉, β-Glucosidase는 습토에 톨루엔, 구연산-인산완충액, 증류수 및 50 mM PNG 기질을 첨가하고, Phosphatase는 톨루엔, 0.5 M Tris (pH 6.5), 증류수 및 50 mM pNPP 기질을 첨가하여 30°C에서 1시간 배양하였다. 배양 후 Ethanol을 넣고 15초간 교반 및 여과하여, 여과액에 2 M Tris를 넣어 흡광도를 측정하였다.
통계 분석
토양 화학적, 생물학적 특성 등 지역 간 평균 비교를 위해 ANOVA 검정을 하였고, 유의성이 인정되는 처리구는 95% 수준에서 Duncan’s 사후 검정으로 유의성 정도를 분석하였다. 또한, 항목 간 연관성 분석을 위해 Pearson’s correlation 분석으로 상관관계를 평가하였으며, 모든 통계는 IBM SPSS statistics 21 (statistical package for social science, IBM) 통계프로그램을 활용하여 분석하였다.
Results and Discussion
토양의 화학적 특성 결과
토양의 화학적 특성을 분석한 결과, 전기전도도를 제외한 모든 항목에서 지역 간의 유의한 차이가 나타났다 (Table 3). 북부지역은 pH 7.4 (4.6 - 8.4), 토양 유기물 함량 45.1 g kg-1 (7.0 - 141.9 g kg-1), CEC 13.9 cmolc kg-1 (4.2 - 26.4 cmolc kg-1)로 타지역에 비해 통계적으로 유의하게 높은 수치로 조사되었고, 중부지역은 유효인산 함량에서 849 mg kg-1 (49 - 1887 mg kg-1)로 타지역보다 유의하게 높았다.
농촌진흥청의 밭토양 화학성 적정 기준 (NAS, 2017)에 따라 충북 권역별 미달/적정/초과 농가 비율을 조사하였을 때, 북부지역은 pH의 적정 범위를 벗어난 토양 비율이 74%로 가장 높게 나타났으며, 칼슘 이온 (Ca2+) 함량에서도 적정 범위보다 높은 평균 10.2 cmolc kg-1 (3.1 - 20.4 cmolc kg-1)로 다른 지역보다 높은 수치를 나타냈다 (Table 3, Fig. 1). 이는 주로 강원도 남부의 영월지역과 충북 북부의 제천, 단양 지역에 분포하는 석회암 지대의 특성에서 기인된 것으로, 본 연구에서 채취한 북부지역의 시료 중 81.5%가 석회암 지대에 포함된다 (data not shown). 일반적으로 석회암에서 유래된 토양은 pH가 높고 점토와 칼슘 함량이 풍부한 적색토 특성을 나타낸다 (ATI, 1992; Hur et al., 2009).
Table 3.
Chemical properties of upland soils in Chungbuk province.
Region
|
pH1
(1:5 H2O)
|
EC
(dS m-1)
|
OM
(g kg-1)
|
Av. P2O5
(mg kg-1)
|
Exch. cations (cmolc kg-1)
|
Ca
|
K
|
Mg
|
Northern
|
7.4
a2 |
1.10
a
|
45.1
a
|
635
b
|
10.17
a
|
0.82
a
|
2.87
a
|
Central
|
6.4
b
|
1.18
a
|
27.3
b
|
849
a
|
7.57
b
|
1.03
a
|
2.21
ab
|
Southern
|
6.3
b
|
0.77
a
|
28.2
b
|
485
b
|
7.28
b
|
0.43
b
|
2.07
b
|
Optimum range
|
6.0 - 7.0
|
≤2.00
|
20 - 30
|
300 - 550
|
5.0 - 6.0
|
0.5 - 0.8
|
1.5 - 2.0
|
한편, 유기물 함량의 적정 범위를 초과한 토양 비율은 북부지역에서 66%로 가장 높았고, 남부지역 40%, 북부지역 31% 순이었다 (Fig. 1). 또한 적정 범위를 나타낸 토양 비율은 북부지역, 중부지역, 남부지역이 각각 26%, 38%, 44%였고, 적정 범위에 미달하는 비율은 8%, 31%, 16%로 나타나, 충북 권역의 유기 농경지에서는 유기물 함량의 과다한 비율이 가장 높았다. 유효인산 함량의 경우 남부지역을 제외한 모든 지역에서 조사 필지의 60% 이상이 적정 범위를 초과하였으며, 교환성 양이온 중 칼슘 이온 (Ca2+)은 모든 지역에서 60% 이상, 마그네슘 이온 (Mg2+)은 40% 이상에서 영양 과다 현상이 나타났다. 이러한 영양불균형 현상은 지역과 관계없이 다수의 필지에서 나타났으며 이는 유기 농경지에서 주로 가축분 퇴비를 통해 토양 양분 관리가 이루어지기 때문으로 판단된다.

Fig. 1.
Distribution status of soil chemical properties by region of soil managed as organic farming for more than 5 years in Chungbuk province.
선행 연구에 따르면, 유기농가는 절반 이상 (57%)이 가축분이 포함된 부숙 퇴비를 사용하고 있으며, 시비하는 가축분의 종류에 따라 토양 내 유효인산 및 칼슘 이온 (Ca2+) 등의 금속 양이온 농도가 다르게 증가하고, 이들이 지속적으로 토양에 축적된다고 보고된 바 있다 (Lee et al., 2004; Yao et al., 2007; Ge et al., 2011; Lee et al., 2017; Lee et al., 2018). 또한 유기재배 기간이 길어질수록 토양의 전반적인 양분 함량이 증가하는 경향을 보이며 (Lee et al., 2004), 특히 관행 재배지의 밭토양에서 상대적으로 경미한 문제로 여겨지는 염류집적 현상이 유기농업에서는 빈번히 나타나고 있다 (Lee et al., 2017). 토양 내 유기물 함량은 작물 생산의 중요한 요소로 유기물 함량과 작물 생산량의 양의 상관관계가 보고된 바 있어 적절한 유기물 시비관리는 필수적이다 (Loveland and Webb, 2003; Schjonning et al., 2018; Hong et al., 2024). 그러나 본 연구에서도 충북의 다수 지역에서 영양불균형 현상이 나타났으며, 현재 충북지역의 유기재배 농가에서는 토양 내 양분을 고려하지 않고 유기물 투입을 위해 퇴비를 과다하게 사용하고 있는 것으로 판단된다. 따라서 토양 검정에 의한 시비 처방으로 적정량의 퇴비를 살포하고 물리성 개선과 과잉 집적된 염류 제거에 효과적인 녹비작물, 두과작물 및 심근성작물을 윤작함으로써 토양의 근본적인 건전성 개선이 필요할 것으로 보인다 (Willer, 1998).
토양의 생물학적 특성 결과
토양 호흡량은 미생물이 유기물 분해 시 에너지 획득 과정에서 방출하는 이산화탄소량을 산출한 것으로, 북부지역은 평균 55.4 CO2 mg kg-1 day-1 (3.6 - 204.4 CO2 mg kg-1 day-1), 중부지역은 46.9 CO2 mg kg-1 day-1 (0 - 202.7 CO2 mg kg-1 day-1)으로 나타났다 (Fig. 2). 이는 이전의 Kim et al. (2023)이 보고한 유기 농경지 토양 호흡량 평균 44.0 CO2 mg kg-1 day-1보다 더 높은 경향을 보여주었으며 남부지역의 농경지는 평균 26.6 CO2 mg kg-1 day-1 (0 - 99.4 CO2 mg kg-1 day-1)로 다소 낮은 활성을 나타냈다.
효소 활성은 토양 내 유기물을 무기태로 전환하여 식물의 양분 가용화에 직접적인 영향을 미치는 중요한 요소로, 토양의 변화에 신속하게 반응하여 토양의 질을 평가하는데 주로 사용된다 (Benitez et al., 2006). 이때, 효소는 기질 특이적이기 때문에 단일 효소 반응으로 토양의 질을 평가하는 데 한계가 있으므로 (Nannipieri et al., 2012), C, N, P 순환에서 중요한 역할을 하는 4가지 효소 (Dehydrogenase, β-Glucosidase, Urease, Phosphatase)를 측정하였다. 그 결과, 토양 내 효소 활성의 모든 항목에서 북부지역과 다른 지역 간의 통계적으로 유의한 차이를 보였다 (Fig. 2). 그 중, C 순환에 관여하는 Dehydrogenase와 β-Glucosidase, N 순환에 관여하는 Urease는 북부지역이 다른 지역보다 통계적으로 유의하게 높았다. 일반적으로 유기물이 풍부할 경우, 토양 내 미생물에게 충분한 기질을 제공하여 바이오매스 함량이 높게 유지되고 높은 효소 활성이 나타나 토양의 효소 활성과 유기물 함량 간의 상관관계가 있다고 보고된 바 있다 (Lagomarsino et al., 2009; Moeskops et al., 2010; Zhao et al., 2010; Mohammadi, 2011). 본 결과에서도 유기물 함량이 다른 지역에 비해 약 1.6배 높은 북부지역에서 Dehydrogenase, β-Glucosidas, Urease 활성이 각각 평균 149.3 TPF µg g-1 24 h-1(9.4 - 458.7 TPF µg g-1 24 h-1), 111.1 PNP µg g-1h-1 (24.2 - 261.2 PNP µg g-1 h-1), 210.0 NH4-N µg g-1 24 h-1 (0 - 430.9 NH4-N µg g-1 24 h-1)로 다른 지역에 비해 유의하게 높았으며 이는 유기물 함량에 따라 효소 활성이 증가한다는 결과와 일치한다.

Fig. 2.
Box-plots of soil biological properties of upland soils in Chungbuk province. zThe same letter in each column indicated no significant difference in soil parameters from the three different regions determined by Duncan’s multi-range test at p ≤ 0.05.
P 순환에 중요한 역할하는 Phosphatase는 토양 내 유기태로 존재하는 인 (P)이나 유기성 잔해물에서 유래한 인 (P)을 가수분해하는 효소로, 유기물 함량이 증가하거나 pH가 감소할 때 Phosphatase 활성이 증가한다는 여러 연구 결과가 있다 (Joner and Jakobsen, 1995; Sarapatka et al., 2004; Van Aarle and Passard, 2010). 본 결과에서는 북부지역이 평균 453.5 PNP µg g-1 h-1 (77.0 - 1596.5 PNP µg g-1 h-1)로, 중부지역 880.7 PNP µg g-1 h-1, 남부지역 718.4 PNP µg g-1 h-1에 비해 가장 낮은 값을 나타났으며 이는 pH의 영향으로 판단된다 (p < 0.05). 토양의 pH는 북부 (7.4) > 중부 (6.4) > 남부 (6.3)순으로, Phosphatase 활성과 역순의 결과가 나타났으며 pH가 낮아짐에 따라 Phosphatase 활성이 증가한다는 결과와 일치한다.
토양의 화학적 특성과 생물학적 특성 간 상관관계
Pearson 상관계수 분석 결과 (Table 4), 토양 화학적 특성과 생물학적 특성의 다양한 항목에서 유의한 상관관계가 나타났다. pH는 토양 호흡량 (r = 0.37), Dehydrogenase (r = 0.56), β-Glucosidase (r = 0.34), Urease (r = 0.51) 등과 양의 상관관계를, Phosphatase (r = -0.64)과는 음의 상관관계를 나타냈다. pH는 토양 화학적 특성 중 중요한 지표로, 양분의 가용성과 미생물 군집, 미생물 활성 등 다양한 생태적 요소에 영향을 미친다 (Frankenberger and Johanson, 1982; Dick and Wang, 2000). 대부분의 생물학적 특성은 pH와의 상관관계가 기존의 선행연구 결과와 유사한 경향을 보였다 (Eivazi and Tabatabai, 1990; Puissant et al., 2019).
Table 4.
Correlation coefficient between chemical and biologocal properties of soil examined in this study.
|
pH
|
EC1 |
OM
|
Av. P2O5 |
T-N
|
Ex. Ca
|
Ex. K
|
Ex. Mg
|
SR
|
DHA
|
GLU
|
URE
|
AP
|
pH
|
1
| | | | | | | | | | | | |
EC
|
-0.21
*
|
1
| | | | | | | | | | | |
OM
|
0.20
*
|
0.03
|
1
| | | | | | | | | | |
Av. P2O5 |
-0.01
|
0.14
|
0.09
|
1
| | | | | | | | | |
T-N
|
-0.01
|
0.13
|
0.46
**
|
0.28
**
|
1
| | | | | | | | |
Ex. Ca
|
0.47
**
|
-0.10
|
0.41
**
|
0.05
|
0.10
|
1
| | | | | | | |
Ex. K
|
0.04
|
0.09
|
0.12
|
0.45
**
|
0.43
**
|
0.06
|
1
| | | | | | |
Ex. Mg
|
0.33
**
|
-0.05
|
0.41
**
|
0.28
**
|
0.18
|
0.60
**
|
0.27
**
|
1
| | | | | |
SR
|
0.37
**
|
-0.09
|
0.44
**
|
0.03
|
0.29
**
|
0.34
**
|
0.16
|
0.27
**
|
1
| | | | |
DHA
|
0.56
**
|
-0.16
|
0.32
**
|
0.06
|
0.24
**
|
0.23
*
|
-0.01
|
0.33
**
|
0.45
**
|
1
| | | |
GLU
|
0.34
**
|
0.05
|
0.21
*
|
0.08
|
0.35
**
|
0.27
**
|
0.12
|
0.24
**
|
0.29
**
|
0.61
**
|
1
| | |
URE
|
0.51
**
|
-0.12
|
0.35
**
|
0.24
*
|
0.35
**
|
0.30
**
|
0.30
**
|
0.37
**
|
0.50
**
|
0.67
**
|
0.55
**
|
1
| |
AP
|
-0.64
**
|
0.30
**
|
-0.10
|
0.10
|
0.17
|
-0.30
**
|
0.06
|
-0.26
**
|
-0.23
*
|
-0.33
**
|
0.01
|
-0.39
**
|
1
|
유기물 함량은 토양 호흡량 (r = 0.44), Dehydrogenase (r = 0.32), β-Glucosidase (r = 0.20), Urease (r = 0.35)와 양의 상관관계를 보였다. 또한 전질소 함량, 칼슘 이온 (Ca2+) 및 마그네슘 이온 (Mg2+)도 이들 생물학적 특성과 강한 상관관계를 나타내며, 이는 유기물 함량의 영향에서 파생된 것으로 사료된다. 이러한 결과는 유기물 함량이 증가할수록 생물학적 특성이 더욱 활성화됨을 시사하며, 가축분 퇴비를 장기간 사용하면 토양 내 탄소 (C), 질소 (N) 순환에 관련된 효소의 활성이 증가한다는 Li et al. (2015)의 연구 결과와 일치한다. 그러나 유기물 함량에 따라 칼슘 이온 (Ca2+) 및 마그네슘 이온 (Mg2+) 또한 높은 양의 상관관계 (r = 0.41)를 보임에 따라, 현재 충북지역 유기 밭토양 농경지에서 문제가 되는 양분 불균형을 초래할 수 있으므로, 토양 내 유기물 투입에 주의가 필요할 것으로 보인다. 또한 토양 화학적 및 생물학적 특성과 농업생산성과 환경 건강에 미치는 영향을 규명하기 위해 추후 토양 관리법 및 수확량 검정 등 환경 요인에 대한 조사가 필요할 것으로 판단된다.
Conclusions
본 연구에서는 충북지역의 5년 이상 유기 인증이 유지된 농가를 중심으로 토양의 화학적 및 생물학적 특성을 분석하여 지역별 비교 평가하였다. 분석 결과, 전기전도도를 제외한 모든 항목에서 지역 간 유의한 차이가 나타났으며, 북부지역의 pH, 유기물 함량 등이 가장 높게 나타났다. 특히, 북부지역에서는 채취한 필지의 74%가 알칼리성 토양으로 나타났으며, 이러한 차이는 주로 석회암 지대의 영향을 받는 것으로 판단된다.
유기물 함량이 30 g kg-1 이상인 필지 또한 북부지역에서 가장 많이 발견되었으며, 북부지역의 높은 유기물 함량은 토양 호흡량과 Phosphatase를 제외한 효소 활성 증가에 기여하며, 높은 상관관계를 보여주었다. 그러나 북부지역뿐만 아니라 충북의 다수 유기 농경지 필지에서 유기물 함량 및 유효인산 등의 과다로 영양 불균형 문제가 발생하고 있다. 이러한 결과는 유기물의 증가가 생물학적 특성을 활성화시키는 장점이 있지만, 동시에 양분 불균형 문제를 초래할 수 있어 적정한 유기물 투입과 영양 균형을 고려한 관리가 필요하다는 점을 시사한다. 본 연구를 통해 충청북도 내 5년 이상 유기인증을 받은 밭토양의 화학적 및 생물학적 특성을 분석하였으며, 이를 바탕으로 유기 농경지의 토양 건강성를 위한 기초 자료로 활용할 수 있을 것으로 판단된다.
Funding
This work was carried out with the support of Project No.RS-2022-RD010404 from the Rural Development Administration of the Republic of Korea.
Conflict of Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Author Contribution
Do JW: Conceptualization, Data curation, Writing-original draft, Visualization, Kim SK: Formal analysis, Data curation, Lee SH: Formal analysis, Data curation, Lee CR: Funding acquisition, Writing-review & editing, Park JH: Supervision, Conceptualization, Writing-review & editing, Kim JH: Supervision,Conceptualization, Writing-review & editing.
Data Availability
Data will be provided on reasonable request.
Acknowledgements
The authors thanks Chungcheongbuk-do Agricultural Research and Extension Services for providing the research sites and experimental materials.
References
ATI (Agricultural Technique Institute). 1992. Korean Soil Introduction. Agricultural Technique Institute, Suwon, Korea.
Barak P, Jobe BO, Krueger AR, Peterson LA, Laird DA. 1997. Effects of long-term soil acidification due to nitrogen fertilizer inputs in Wisconsin. Plant Soil 197:61-69.
10.1023/A:1004297607070Benitez E, Nogales R, Campos M, Ruano F. 2006. Biochemical variability of olive-orchard soils under different management systems. Appl. Soil Ecol. 32:221-231. https://doi.org/10.1016/j.apsoil.2005.06.002
10.1016/j.apsoil.2005.06.002Bottomley PJ, Angle JS, Weaver R. 2020. Methods of soil analysis, Part 2: Microbiological and biochemical properties. John Wiley & Sons, Inc., Hoboken, NJ, USA.
Casida LE, Klein DA, Santoro T. 1964. Soil dehydrogenase activity. Soil Sci. Soc. Am. J. 47:599-603.
10.1097/00010694-196412000-00004Eivazi F, Tabatabai M. 1990. Factors affecting glucosidase and galactosidase activities in soils. Soil Biol. Biochem. 22:891-897.
10.1016/0038-0717(90)90126-KEivazi F, Tabatabai MA. 1988. Glucosidases and galactosidases in soils. Soil Biol. Biochem. 20:601-606.
10.1016/0038-0717(88)90141-1Ge T, Nie S, Wu J, Shen J, Xiao H, Tong C, Huang D, Hong Y, Iwasaki K. 2011. Chemical properties, microbial biomass, and activity differ between soils of organic and conventional horticultural systems under greenhouse and open field management: A case study. J. Soils Sediments 11:25-36. https://doi.org/10.1007/s11368-010-0293-4
10.1007/s11368-010-0293-4Gimeno-García E, Andreu V, Boluda R. 1996. Heavy metals incidence in the application of inorganic fertilizers and pesticides to rice farming soils. Environmental Pollution 92:19-25. https://doi.org/10.1016/0269-7491(95)00090-9
10.1016/0269-7491(95)00090-915091407Gomiero T, Paoletti MG, Pimentel D. 2008. Energy and environmental issues in organic and conventional agriculture. Crit. Rev. Plant Sci. 27:239-254.
10.1080/07352680802225456Hong SD, Lee CS, Hwang SW. 2009. Soil Fertility. Korean J. Soil Sci. Fert. 42:178-194.
Hong YK, Kim JW, Lee KM, Lee DJ, Lee CR, Kim SC. 2024. Selecting indicator of soil properties to correlate with crop productivity in organic farming system. Korean J. Soil Sci. Fert. 57:35-47. https://doi.org/10.7745/KJSSF.2024.57.1.035
10.7745/KJSSF.2024.57.1.035Hur SO, Jung KH, Sonn YK, Ha SK, Kim JG. 2009. Determination of pedo-transfer function using the relation between soil particle distribution, organic matter and water movement in soil originated from limestone. Korean J. Soil Sci. Fert. 42:132-138.
Joner EJ, Jakobsen I. 1995. Growth and extracellular phosphatase activity of arbuscular mycorrhizal hyphae as influenced by soil organic matter. Soil Biol. Biochem. 27:1153-1159. https://doi.org/10.1016/0038-0717(95)00047-I
10.1016/0038-0717(95)00047-IKandeler E, Gerber H. 1988. Short-term assay of soil urease activity using colorimetric determination of ammonium. Biol. Fertil. Soils 6:68-72.
10.1007/BF00257924Kim JW, Hong YK, Lee CR, Kim SC. 2023. Comparison of physicochemical and biological soil properties in organic and conventional upland fields. Korean J. Soil Sci. Fert. 56:77-89. https://doi.org/10.7745/KJSSF.2023.56.1.077
10.7745/KJSSF.2023.56.1.077Kim LY, Cho HJ, Han KH. 2003. Effects of tile drain on physicochemical properties and crop productivity of soils under newly constructed plastic film house. Korean J. Soil Sci. Fert. 36:154-162.
Kim YH., Kong MS, Lee EJ, Lee TG, Jung GB. 2019. Status and changes in chemical properties of upland soil from 2001 to 2017 in Korea. Korean J. Environ. Agric. 38:213-218. https://doi.org/10.5338/KJEA.2019.38.3.28
10.5338/KJEA.2019.38.3.28Kim YS, Lee YM, Weon HY, Sang MK, Song JK. 2020. Comparative analysis of soil microbial communities between conventional and organic farming systems in pepper cultivation. Korean J. Organic Agric. 28:235-250. http://doi.org/10.11625/KJOA.2020.28.2.235
Lagomarsino A, Moscatelli MC, Di Tizio A, Mancinelli R, Grego S, Marinari S. 2009. Soil biochemical indicators as a tool to assess the short-term impact of agricultural management on changes in organic C in a Mediterranean environment. Ecol. Indic. 9:518-527. https://doi.org/10.1016/j.ecolind.2008.07.003
10.1016/j.ecolind.2008.07.003Lee CR, Kim PJ, Oh Y, Park CB, Park KL, Nam HS, Park GC. 2018. Biomass, nitrogen, and phosphorus productivities of green manure by barley and hairy vetch mixtures. Korean J. Organic Agric. 26:719-729. https://doi.org/10.11625/KJOA.2018.26.4.719
10.11625/KJOA.2018.26.4.719Lee CR, Ok JH, An MS, Lee SB, Park KL, Hong SG, Kim MG, Park CB. 2017. Soil chemical properties of long-term organic cultivation upland. Korean J. Organic Agric. 25:161-170. https://doi.org/10.11625/KJOA.2017.25.1.161
10.11625/KJOA.2017.25.1.161Lee YJ, Choe DH, Kim SH, Lee SM, Lee YH, Lee BM, Kim TW. 2004. Long-term changes in soil chemical properties in organic arable farming systems in Korea. Korean J. Soil Sci. Fert. 37:228-234.
Li J, Cooper JM, Li Y, Yang X, Zhao B. 2015. Soil microbial community structure and function are significantly affected by long-term organic and mineral fertilization regimes in the North China Plain. Appl. Soil Ecol. 96: 75-87.
10.1016/j.apsoil.2015.07.001Liang Y, Wu L, Clark IM, Xue K, Yang Y, Van Nostrand JD. 2015. Over 150 years of long-term fertilization alters spatial scaling of microbial biodiversity. MBio 6:e240-e215. https://doi.org/10.1128/mBio.00240-15
10.1128/mBio.00240-1525852161PMC4453543Loveland P, Webb J. 2003. Is there a critical level of organic matter in the agricultural soils of temperate regions: A review. Soil Tillage Res. 70:1-18. https://doi.org/10.1016/S0167-1987(02)00139-3
10.1016/S0167-1987(02)00139-3Maeder P, Fliessbach A, Dubois D, Gunst L, Fried P, Niggli U. 2002. Soil fertility and biodiversity in organic farming. Science 296:1694-1697. https://doi.org/10.1126/science.1071148
10.1126/science.107114812040197Moeskops B, Buchan D, Sleutel S, Herawaty L, Husen E, Saraswati R, Setyorini D, DeNeve S. 2010. Soil microbial communities and activities under intensive organic and conventional vegetable farming in West Java, Indonesia. Appl. Soil Ecol. 45:112-120.
10.1016/j.apsoil.2010.03.005Mohammadi K. 2011. Soil microbial activity and biomass as influenced by tillage and fertilization in wheat production. American-Eurasian J. Agric. Environ. Sci. 10:330-337.
Nannipieri P, Giagnoni L, Renella G, Puglisi E, Ceccanti B, Masciandaro G, Fornasier F, Moscatelli MC, Marinari S. 2012. Soil enzymology: Classical and molecular approaches. Biol. Fertil. Soils 48:743-762. https://doi.org/10.1007/s00374-012-0723-0
10.1007/s00374-012-0723-0NAAS (National Academy of Agricultural Science). 2010. Methods of soil analysis. RDA, Suwon, Korea.
NAS (National Institute of Agricultural Science). 2017. Fertilizer recommendation for crops. Third edition. Rural Development Administration, Wanju, Korea.
Pimentel D, Hepperly P, Hanson J, Seidel R, Douds D. 2005. Organic and conventional farming systems: Environmental and economic issues. Report 05-1. The Rodale Institute, Kutztown, PA, USA.
Puissant J, Jones B, Goodall T, Mang D, Blaud A, Gweon HS, Malik A, Jones DL, Clark IM, Hirsch PR, et al. 2019. The pH optimum of soil exoenzymes adapts to long-term changes in soil pH. Soil Biol. Biochem. 138:107601. https://doi.org/10.1016/j.soilbio.2019.107601
10.1016/j.soilbio.2019.107601Sarapatka B, Dudová L, Krsková M. 2004. Effect of pH and phosphate supply on acid phosphatase activity in cereal roots. Biol. Bratislava 59:127-131.
Schjonning P, Jensen JL, Bruun S, Jensen LS, Christensen BT, Munkholm LJ, Oelofse M, Baby S, Knudsen L. 2018. The role of soil organic matter for maintaining crop yields: Evidence for a renewed conceptual basis. Adv. Agron. 150:35-79. https://doi.org/10.1016/bs.agron.2018.03.001
10.1016/bs.agron.2018.03.001Sun B, Zhang L, Yang L, Zhang F, Norse D, Zhu Z. 2012. Agricultural non-point source pollution in China: Causes and mitigation measures. Ambio 41:370-379. https://doi.org/10.1016/j.envpol.2023.121953
10.1016/j.envpol.2023.12195337307861Tabatabai MA, Bremner JM. 1969. Use of p-nitrophenyl phosphate for assay of soil phosphatase activity. Soil Biol. Biochem. 1:301-307.
10.1016/0038-0717(69)90012-1Van Aarle IM, Plassard C. 2010. Spatial distribution of phosphatase activity associated with ectomycorrhizal plants is related to soil type. Soil Biol. Biochem. 42:324-330. https://doi.org/10.1016/j.soilbio.2009.11.011
10.1016/j.soilbio.2009.11.011Vitousek PM, Naylor R, Crews T, David MB, Drinkwater L, Holland E, Johnes P, Katzenberger J, Martinelli L, Matson P. 2009. Nutrient imbalances in agricultural development. Science 324:1519-1520. https://doi.org/10.1126/science.1170261
10.1126/science.117026119541981Weaver RW, Angle S, Bottomley P, Bezdicek D, Smith S, Tabatabai A, Wollum A. 1994. Methods of Soil Analysis, Part 2: Microbiological and Biochemical Properties. SSSA, Madison, WI, USA.
10.2136/sssabookser5.2Willer H. 1998. Okologischer landbau in Europa. SEOL, Germany. pp. 400.
Yamagishi J, Nakamoto T, Richner W. 2003. Stability of spatial variability of wheat and maize biomass in a small field managed under two contrasting tillage. Field Crops Res. 81:95-108. https://doi.org/10.1016/S0378-4290(02)00213-7
10.1016/S0378-4290(02)00213-7Yang JE, Cho BO, Shin YO, Kim JJ. 2001. Fertility status in northeastern alpine soils of South Korea with cultivation of vegetable crops. Korean J. Soil Sci. Fert. 34:1-7.
Yao LX, Li GL, Tu SH, Sulewski G, He ZH. 2007. Salinity of animal manure and potential risk of secondary soil salinization through successive manure application. Sci. Total Environ. 383:106-114. https://doi.org/10.1016/j.scitotenv.2007.05.027
10.1016/j.scitotenv.2007.05.02717572477Zhao B, Chen J, Zhang J, Qin S. 2010. Soil microbial biomass and activity response to repeated drying-rewetting cycles along a soil fertility gradient modified by long-term fertilization management practices. Geoderma 160: 218-224.
10.1016/j.geoderma.2010.09.024