2019, Vol. 56
表1(Table 1)
2. 邢台学院资源与环境学院，河北邢台 054001
水耕人为土（水稻土）是我国最重要的土壤资源之一。作为重要的粮食生产基地[1]与有机碳（SOC）库[2-3]及温室气体——CH4的重要排放源[4-5]，其发生和演变过程对土壤管理有重要的意义，因而受到了广泛的关注[6-10]。成土母质作为土壤形成的物质基础对土壤的理化性质具有重要的影响[11]。虽然不同母质起源的土壤经过一定时间的水耕熟化后均可以发育为水耕人为土，但母质不同很可能会导致其发育速率和演变特征存在一定的差异。因此母质如何影响水耕人为土的发育速率和演变特征是人为土研究中一个值得探讨的问题。中国土壤系统分类作为定量化的分类系统，采用可测定的指标定义了水耕人为土并与其起源土壤进行区分。这从根本上解决了水耕人为土与起源土壤之间的界线问题，对研究水耕人为土发生有利，但如何度量水耕人为土的发育程度和发育速率仍是土壤发生学的一个难题。在最初的研究中，主要利用水耕熟化程度来反映水耕人为土的发育程度[12]。实际上水耕人为土发生包括了两个相互矛盾的过程，即水耕熟化和水耕淋溶过程。水耕熟化使土壤肥力正向发展，而强烈的水耕淋溶使土壤肥力衰退。水耕人为土的发育是这两个过程共同作用的结果（特定条件下还有一些附加过程，如潜育化过程等）。因此，从土壤发生学角度而言，水耕熟化过程并不完全等同于水耕人为土的发育过程。当前亟需引入一些定量指标来衡量水耕人为土的发育速率和发育程度。相比较而言，基于土壤理化性质的属性距离[13]和基于形态特征的土壤发育指数[14-17]等指标在表征土壤发育程度上有独特的优势。这些指标能够综合反映土壤理化性质和土壤形态特征的变化，并排除了单个土壤属性用于指示土壤发育程度时的不确定性。为此，本研究选择南方丘陵区3种常见母质发育的水耕人为土时间序列为研究对象，尝试利用属性距离和土壤发育指数来表征水耕人为土的发育程度，揭示母质对水耕人为土发育速率的影响，以期为今后的水耕人为土发生学研究提供参考。
1 材料与方法 1.1 研究区概况本研究选取我国亚热带地区3种母质发育的水耕人为土从十年到百年尺度的时间序列（表 1），分别为发育于广西龙胜各族自治县（龙脊梯田）紫色砂页岩坡积物的水耕人为土时间序列（PS10，0 a；PS11，30 a；PS12，100~300 a；PS13，300 a）、发育于江西省进贤县第四纪红黏土的水耕人为土时间序列（RC10，0 a；RC11，100 a；RC12，100~300 a；RC13，300 a）和发育于江西省余江县红砂岩坡积物的水耕人为土时间序列（RS10，0 a；RS11，30 a；RS12，60~200 a；RS13，200 a）。研究区属亚热带季风气候，年均气温分别为18.2、17.3和17.6℃，年均降水量分别为1 557、1 549和1 789 mm。3个研究区的年均气温和降水量相差不大，本研究认为其气候条件基本一致，不再单独讨论。3种母质发育的水耕人为土时间序列建立和判定的依据分别见文献[18]，采样点信息见表 1，剖面照片见图 1。土壤类型的判别主要基于文献[19]。
|
表 1 剖面基本信息 Table 1 Descriptions of the studied soil profiles |
|
图 1 不同母质水耕人为土时间序列的土壤剖面照片 Fig. 1 Soil profiles of Hydragric Anthrosols different in parent material and chronosequence |
结合史料分析与实地考察，选取不同母质发育的水耕人为土时间序列典型土壤剖面，按照土壤发生层采样法采集土样并描述土壤形态（包括土壤颜色、土壤结构和结持性等）[20]，土壤颜色根据《中国标准土壤色卡》[21]确定。所有土样均在水稻收割排水后采集，土壤剖面选在稻田的中心部位，坡底剖面避开容易积水的洼地。土壤样品采集后在室内自然风干，挑出枯枝落叶、根系和大2 mm的非土壤物质，四等分法取土，先后过10、60和100目的尼龙筛，装好备用。
1.3 测定项目与分析方法SOC采用重铬酸钾-外加热法测定；pH（水土比2.5:1）采用电位法测定；土壤颗粒组成采用激光粒度仪测定；游离铁（Fed）采用连二亚硫酸钠-柠檬酸钠-重碳酸钠浸提（DCB），原子吸收光谱仪测定[22]。磁化率（MS）的测定：称取适量（以装满为准）过10目筛的风干土样（精确度0.01 g），装入10 cm3无磁性塑料盒中，用MS-2B型磁化率仪（英国Bartington公司）分别在低频（0.47 kHz，MSlf）和高频（4.7 kHz，MShf）磁场中进行测定，每个样品测定3次，最后取其平均值。
1.4 属性距离的计算基于Carré和Jacobson[13]提出的理论计算了水耕人为土剖面（剖面PS11和RS1为雏形土，下同）与起源土壤剖面间的发生层距离（Horizon distance，HD）。在计算HD时采取自上而下层次对应的方法[13]。起源土壤与水耕人为土对应层次间的HD计算如式（1）和式（2）。
| $ HD = {D_h}\left( {{h_a}, {h_b}} \right)/N $ | (1) |
| $ {D_h}\left( {{h_a}, {h_b}} \right) = \sum\limits_{f = 1}^N {\left| {{v_{a, f}} - {v_{b, f}}} \right|} $ | (2) |
式中，HD为发生层距离，Dh（ha，hb）为发生层ha与hb的Manhattan距离，N为参与计算的土壤属性的个数，va, f与vb, f分别为发生层ha与发生层hb的某种土壤属性数值。本研究用于HD计算的土壤属性主要为在水耕人为土发育过程中有规律变化的黏粒、SOC、pH、Fed和MS共5个土壤属性，并在计算HD之前对土壤属性数据进行标准化，使所有的数据介于0~1之间。
考虑到土壤厚度对土壤肥力供应有一定的影响。因此本研究在计算剖面距离（Pedological distance，PDD）时，将土壤厚度考虑进去（剖面厚度取120 cm，不足120 cm的取实际值）。各土壤剖面PDD由该剖面各个发生层的HD值乘以该发生层的厚度后求和而得。
1.5 剖面发育指数的计算基于Harden[14]、Harden和Taylor[15]提出的理论，利用剖面描述特征计算土壤发育指数。首先对土壤各个发生层的描述数据进行赋值，接着计算水耕人为土发生层各种赋值数据与起源土壤的差值（自上而下层次对应），然后标准化使所有数据介于0~1之间。每个发生层所有标准化后的形态属性取平均值得到该发生层的发生层指数（Horizion index，HI）。剖面每个发生层的HI值乘以该发生层的厚度后求和得到该剖面的发育指数（Profile development index，PDI）。标准化时各属性的最大值采用了Harden[14]提出的最大理论数值，这样可以统一基准值，便于不同研究间的相互比较。选取红度指数、黑化指数、总质地、结构和湿结持性[14]来计算土壤发育指数。由于采样时土壤均处于湿润状态，土壤的干结持性没有描述。已往的研究也表明由于形态属性内在的相关性，个别属性的缺失对结果的影响并不大[14]。pH在计算属性距离时已经用到，为避免重复就未用于剖面发育指数的计算。
由于水耕人为土的发育过程和方向同自然土壤不同，各项描述数据在数值化时有各自的独特性，如：（1）种稻后土壤颜色的色调变黄，明度降低，土壤红度的赋值标准和自然土壤相反；（2）水耕人为土耕作层经常受到翻耕扰动及水耕熟化的影响，土壤结构和土壤湿结持性的发育方向与其他层次不同，因此采用了单独的标准进行赋值（*标记）。详细数值化过程如表 2。
|
|
表 2 剖面形态描述数据的数值化与形态属性的计算 Table 2 Quantification of profile descriptions and calculation of morphological properties |
水耕人为土剖面较自然土壤剖面内部开始有不同程度的分异。主要表现为水耕表层（包括耕作层和犁底层）和氧化还原层的出现，坡底种稻年限长的水耕人为土剖面内部分异更明显。随着种稻年限的增加，土壤颜色的色调变黄，彩度降低。耕作层由于有机质（SOM）含量的提高土壤颜色的明度与彩度值降低，根系较多；犁底层变得愈加紧实（容重介于1.21~1.83 g·cm-3）；氧化还原层有较多的新生铁、锰斑纹或结核形成。这表明淹水种稻后土壤形态学特征发生一系列有规律的变化。PS母质起源土壤的色调为7.5 YR（干态，下同），种稻后逐渐向2.5 Y演化。该母质的土壤有轻度的可塑性和黏结性。RC母质起源土壤的色调为10 R，种稻后逐渐向5 YR演化。该母质土壤的可塑性和黏结性中等。RS母质起源土壤的色调为10 YR，种稻后逐渐向2.5 YR演化。该母质土壤的可塑性和黏结性极差。
同样，土壤黏粒、pH、SOC、Fed和MS等随着种稻年限的增加也呈现了明显的变化。更详细的数据分别参见文献[18]和[24]。RC和RS母质土壤种稻后黏粒含量逐步降低，基本呈现出随着种稻时间的增加而逐渐降低的趋势（图 2）。PS母质土壤种稻后黏粒未有明显降低，剖面下部甚至有较明显的增加。PS、RC和RS母质水耕人为土的pH分别介于4.46~5.41、5.01~6.25和5.2~6.19，略高于起源土壤（分别介于4.65~4.91、4.55~5.35和4.56~5.54）。3种母质土壤种稻后pH均趋于中性，这和已有的研究结果一致[1]。PS和RC母质土壤种稻后SOC的剖面分布规律相似：SOC含量均随土壤深度的增加而下降，水耕人为土耕作层的SOC含量均明显高于起源土壤；种稻年限较短的土壤与起源土壤相比，耕作层SOC含量增加明显，但耕作层以下增加不明显；种稻时间较长的土壤与起源土壤相比，SOC在下层也有较明显的增长。说明这两种母质土壤在种稻初期SOC的增加主要集中在耕作层，到一定阶段后，下层也表现出明显增加。相比之下，RS母质的起源土壤种稻后SOC没有明显增加，种稻年限较短的土壤的SOC含量甚至低于起源土壤，而种稻约200 a的土壤也未表现出较起源土壤有大幅度的增加。PS和RC母质起源土壤Fed的剖面分布相对均一，而RS母质起源土壤的剖面内部变异较大。种稻后Fed含量均逐渐降低，也基本表现出随着种稻时间的增加而逐渐降低的趋势。PS和RS母质起源土壤的MS值较低，种稻后演变规律类似，土壤上、下层呈明显分异（图 2）。种稻后上层土壤的MS值明显降低，在较短的时间内（30 a内）就会降低到一个极低的数值，之后没有明显变化。种稻后下层土壤的MS值降低没有上层明显，种稻年限较短的土壤（PS11和RS11）有明显的增加（RS11的增加尤为明显）。RC母质起源土壤的MS值较高，种稻之后土壤上、下层也呈明显分异。上层土壤MS的演变规律与PS和RS母质上层土壤的演变规律类似。种稻后下层土壤的MS值虽然也下降，但下降幅度没有上层那么大。总体而言，随着种稻年限的加长，3种母质的水耕人为土最终发育成结构相似的MS剖面构型（MS值低且剖面均一）。
|
图 2 不同母质水耕人为土时间序列基本理化性质的剖面分布 Fig. 2 Vertical distribution of basic soil physico-chemical properties in soil profiles of Hydragric Anthrosols different in chronosequence and parent material |
发生层距离（HD）反映的是土壤各发生层理化性质的综合变化程度，而发生层指数（HI）反映的是土壤各发生层形态特征的综合变化程度，值越大代表变化程度越大。两者从不同方面反映土壤发生层的发育程度。3种母质水耕人为土的HD与HI值大致呈现了随着种稻年限的增加而增加的趋势，这和水耕人为土的发育时间相符（图 3）。说明水耕人为土发生层的发育程度可以通过HD与HI两个定量指标表现出来。上层土壤的HD与HI值通常较下部土层更高。说明种稻后上层土壤受人为水耕的影响较大，剖面内部逐渐分异。将HD或HI值和种稻年限的比值近似看作发生层的平均发育速率。通过不同母质种稻年限相近的水耕人为土之间的相互比较可知，RC母质水耕人为土的HD与HI值和种稻年限的比值通常较大，这说明RC母质水耕人为土发生层的发育速率要较PS母质与RS母质水耕人为土快。PS母质水耕人为土发生层的发育速率次之，RS母质最慢。RS母质种稻年限较长剖面（RS12和RS13）耕作层的HI值仍很低，说明RS母质不利于水耕人为土耕作层的发育。
|
图 3 不同母质水耕人为土时间序列属性距离和剖面发育指数的剖面分布 Fig. 3 Distribution of soil attribute distance and soil developmental indices in profiles of Hydragric Anthrosols different in chronosequence and parent material |
PDD或PDI值是剖面理化性质或形态特征变化程度综合的度量，值越大代表该剖面理化性质或形态特征总体的变化程度越大。本研究中PDD和PDI值与已判定的土壤相对种稻年限十分吻合，可以用PDD和PDI值来表示水耕人为土整体的发育程度。将PDD或PDI值和种稻年限的比值近似看作整个剖面的平均发育速率。通过不同母质种稻年限相近的水耕人为土之间的相互比较可知，RC母质土壤的PDD和PDI值与种稻年限的比值最高，PS母质土壤次之，RS母质土壤最低。这表明RC母质水耕人为土的平均发育速率最高，PS母质水耕人为土的平均发育速率次之，RS母质水耕人为土的平均发育速率最慢。种稻年限短的土壤的HD和HI值与种稻年限的比值明显高于种稻年限比较长的土壤。这说明在种稻初期，土壤的发育速率较快（土壤的理化性质与形态特征的变化快），之后随着发育程度的加深，土壤发育速率变慢。
2.3 水耕人为土发育与土壤剖面形态特征和基本理化性质的相关关系在过去的几十年里，大量土壤时间函数被提出，时间跨度几百年至几百万年不等，涉及了不同母质、不同气候下土壤的发育过程[25-27]。其中一些研究认为土壤性质发育与时间成对数或指数关系[25-26]。本研究表明土壤在种稻初期发育速率较快，以后随着发育程度的加深，发育速率变慢，水耕人为土发育速率和成土年龄的关系大致符合对数关系。
水耕人为土理化性质和形态特征的演变并不是完全孤立的两个过程，而是存在着千丝万缕的联系。水耕人为淹水状况对土壤铁锰氧化物的迁移转化非常重要。土壤淹水后，在SOC等因素的作用下，氧化铁发生溶解、螯合溶解与还原溶解等，以Fe2+进入土壤溶液。水溶性Fe2+扩散或下渗至犁底层下面的氧化区，分别以斑状、管状和结核状沉淀下来。铁的这种淋溶和淀积被认为是水耕人为土形成的重要特征之一[1]。这一过程同时伴随了土壤理化性质的规律性变化，如Fed在剖面内部的迁移和再分配、SOC的积累、土壤酸碱度趋于中性和耕作层磁性矿物的衰竭等。将水耕人为土各发生层的HD与HI值做相关分析，发现二者呈极显著线性正相关（r = 0.759，n = 49，P < 0.01），这进一步证实了水耕人为土理化性质与形态特征的演变之间存在密切联系。这也表明水耕人为土剖面形态特征经过定量后和理化性质（属性距离）一样能够表征水耕人为土发育的时间趋势，在指示水耕人为土的发育程度和判定土壤相对种稻年限上有重要意义。
从表 3可见Fed、SOC、pH、MS均与HD值显著相关，表明它们在计算HD值中起着重要作用。红度指数与HI值的相关性最高，可以解释其变异的83.1%。其次是黑化指数和湿结持性。这表明土壤颜色和湿结持性在计算HI值中起着重要作用。种稻后上述理化性质和形态学特征相对容易变化的水耕人为土，其发育速率通常就快，反之，就慢。RC母质的土壤，土层较厚，颜色较红，黏粒含量较高，保肥和保水状况好，土粒的黏结性和可塑性较好，SOC容易累积，土壤结构体和水耕表层（耕作层和犁底层）容易发育。PS母质的土壤，土层较厚，黏粒含量虽然不高，但细颗粒（黏粒+细粉粒，< 0.01mm）含量高（PS母质种稻的细颗粒平均含量高达56.7%），保肥和保水状况好，土粒有一定的黏结性和可塑性，SOC容易累积，土壤结构体和水耕表层也容易发育。土壤发育速率低的RS母质的土壤，黏粒和细颗粒含量均很低（RS母质种稻的细颗粒平均含量为24.6%），土层薄，保肥和保水状况以及黏结性和可塑性差，SOC很难累积，土壤结构体和水耕表层难发育。母质对水耕人为土发育过程的影响可以通过影响其发育速率表现出来。与已有的研究结果[14-17, 25-27]相比，水耕人为土的平均发育速率远较自然土壤快，并在人为作用下快速定向发育。这种新的定量方法可以提高对水耕人为土成土过程的理解并能实现不同地区水耕人为土发育速率的定量比较，有较好的使用前景。但这些方法仍处于试验探讨阶段，有待更多的研究进行验证和进一步改进。
|
|
表 3 土壤的发育程度与土壤剖面形态学特征和基本理化性质的相关关系 Table 3 Relationships of soil development degree with soil profile morphology and basic soil properties |
水耕人为土剖面形态特征经定量后和属性距离一样能够表征水耕人为土发育的时间趋势，在指示水耕人为土的发育程度和估算发育速率上有重要意义。SOC、pH、黏粒含量和MS在计算HD中贡献很大，而土壤颜色和湿结持性在计算HI中起着重要作用。RC母质的土壤，土层较厚，颜色红，黏粒含量较高，保肥和保水状况好，土粒的黏结性和可塑性较好，SOC容易累积，发育速率最快。PS母质的土壤，土层较厚，细颗粒（黏粒+细粉粒，< 0.01mm）含量高，保肥和保水状况好，土粒有一定的黏结性和可塑性，SOC容易累积，发育速率次之。RS母质的土壤，黏粒和细颗粒含量均很低，土层薄，保肥和保水状况以及黏结性和可塑性差，SOC很难累积，发育速率最慢。母质对水耕人为土发育过程的影响可以通过影响其发育速率表现出来。研究同时显示，水耕人为土平均发育速率远较自然土壤快，并在人为培育下快速定向发育。但随着种稻年限的增加，水耕人为土的发育速率普遍下降，水耕人为土发育速率和成土年龄的关系大致符合对数关系。
| [1] |
李庆逵. 中国水稻土. 北京: 科学出版社, 1992. Li Q K. Paddy soils of China (In Chinese). Beijing: Science Press, 1992. (  3) |
| [2] |
Pan G X, Li L Q, Wu L S, et al. Storage and sequestration potential of topsoil organic carbon in China's paddy soils . Global Change Biology, 2003, 10(1): 79-92.
( 1) |
| [3] |
Pan G X, Li L Q, Zhang Q, et al. Organic carbon stock in topsoil of Jiangsu Province, China, and the recent trend of carbon sequestration . Journal of Environmental Sciences(China), 2005, 17(1): 1-7.
( 1) |
| [4] |
Li C S. Modeling trace gas emissions from agricultural ecosystems . Nutrient Cycling in Agroecosystems, 2000, 58(1/3): 259-276. DOI:10.1023/A:1009859006242
( 1) |
| [5] |
Pathak H, Li C S, Wassmann R. Greenhouse gas emissions from India rice fields:Calibration and upscaling using the DNDC model . Biogeosciences, 2005, 2(2): 113-123.
( 1) |
| [6] |
Zhang G L, Gong Z T. Pedogenic evolution of paddy soils in different soil landscapes . Geoderma, 2003, 115(1/2): 15-29.
( 1) |
| [7] |
Cheng Y Q, Yang L Z, Cao Z H, et al. Chronosequential changes of selected pedogenic properties in paddy soils as compared with non-paddy soils . Geoderma, 2009, 151(1/2): 31-41.
( 0) |
| [8] |
Chen L M, Zhang G L, Effland W R. Soil characteristic response times and pedogenic thresholds during the 1000-year evolution of a paddy soil chronosequence . Soil Science Society of America Journal, 2011, 75(5): 1807-1820. DOI:10.2136/sssaj2011.0006
( 0) |
| [9] |
韩光中, 黄来明, 李山泉, 等. 水耕人为土磁性矿物的生成转化机制研究回顾与展望. 土壤学报, 2017, 54(2): 309-318. Han G Z, Huang L M, Li S Q, et al. Review and prospect of researches on production and transformation of magnetic minerals in paddy soils during pedogenesis (In Chinese). Acta Pedologica Sinica, 2017, 54(2): 309-318. ( 0) |
| [10] |
黄来明, 邵明安, 陈留美, 等. 水耕人为土时间序列铁氧化物与磁化率演变特征. 土壤学报, 2017, 54(1): 1-13. Huang L M, Shao M A, Chen L M, et al. Evolution characteristics of iron oxides and magnetic susceptibility in Stagnic Anthrosols along Chronosequences (In Chinese). Acta Pedologica Sinica, 2017, 54(1): 1-13. ( 1) |
| [11] |
Schaetzl R J, Anderson S. Soils:Genesis and geomorphology . New York: Cambridge University Press, 2005.
( 1) |
| [12] |
龚子同. 水稻土的耕性及其在土壤分类上的意义. 土壤学报, 1979, 16(2): 85-93. Gong Z T. The tillage of paddy soil and its significance on soil classification (In Chinese). Acta Pedologica Sinica, 1979, 16(2): 85-93. ( 1) |
| [13] |
Carré F, Jacobson M. Numerical classification of soil profile data using distance metrics . Geoderma, 2009, 148(3/4): 336-345.
( 3) |
| [14] |
Harden J W. A quantitative index of soil development from field descriptions:Examples from a soil chronosequence in central California . Geoderma, 1982, 28(1): 1-28. DOI:10.1016/0016-7061(82)90037-4
( 6) |
| [15] |
Harden J W, Taylor E M. A quantitative comparison of soil development in four climatic regimes . Quaternary Research, 1983, 20(3): 342-359. DOI:10.1016/0033-5894(83)90017-0
( 1) |
| [16] |
Vidic N J, Lobnik F. Rates of soils development of the chronosequence in the Ljubljana Basin, Slovenia . Geoderma, 1997, 76(1/2): 35-64.
( 0) |
| [17] |
Tsai C C, Tsai H, Hseu Z Y, et al. Soil genesis along a chronosequence on marine terraces in eastern Taiwan . Catena, 2007, 71(3): 394-405. DOI:10.1016/j.catena.2007.03.011
( 2) |
| [18] |
Han G Z, Zhang G L. Changes in magnetic properties and their pedogenetic implications for paddy soil chronosequences from different parent materials in South China . European Journal of Soil Science, 2013, 64(4): 435-444.
( 2) |
| [19] |
Cooperative Research Group on Chinese Soil Taxonomy. Chinese Soil Taxonomy . Beijing: Science Press, 2003.
( 1) |
| [20] |
Schoeneberger P J, Wysocki D A, Benham E C, et al. Field book for describing and sampling soils . 2nd ed. Lincoln: National Soil Survey Center, 2002.
( 1) |
| [21] |
中国科学院南京土壤研究所, 中国科学院西安光学精密机械研究所.中国标准土壤色卡.南京: 南京出版社, 1989 Institute of Soil Science, Chinese Academy of Sciences, Xi'an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences. Chinese standard soil colour charts(In Chinese). Nanjing: Nanjing Press, 1989 ( 1) |
| [22] |
Mehra O P, Jackson M L. Iron oxide removal from soils and clays by a dithionite-system buffered with sodium bicarbonate . Clays and Clay Minerals, 1960, 7(5): 317-327.
( 1) |
| [23] |
南京农业大学. 土壤农化分析. 北京: 中国农业出版社, 1986: 153-154. Nanjing Agricultural University. Soil and agricultural chemistry analysis (In Chinese). Beijing: China Agriculture Press, 1986: 153-154. ( 0) |
| [24] |
韩光中, 张甘霖. 母质对南方丘陵区水耕人为土理化性质演变的影响. 土壤学报, 2014, 51(4): 98-106. Han G Z, Zhang G L. Influence of parent material on evolution of physico-chemical properties of Hydragric Anthrosols in hilly regions of South China (In Chinese). Acta Pedologica Sinica, 2014, 51(4): 98-106. ( 1) |
| [25] |
Colman S M. Rock weathering rates functions of time . Quaternary Research, 1981, 15(3): 250-264. DOI:10.1016/0033-5894(81)90029-6
( 3) |
| [26] |
Righi D., Huber K., Keller C. Clay formation and podzol development from postglacier moraines in Switzerland . Clay Minerals, 1999, 34(2): 319-332. DOI:10.1180/000985599546253
( 1) |
| [27] |
Huang W S, Tsai H, Tsai C C, et al. Subtropical soil chronosequence on Holocenemarine terraces in eastern Taiwan . Soil Science Society of America Journal, 2010, 74(4): 1271-1283. DOI:10.2136/sssaj2009.0276
( 2) |
2. College of Resources and Environment, Xingtai University, Xingtai, Hebei 054001, China