矿物岩石地球化学通报  2016, Vol. 35 Issue (2): 353-359   PDF    
山东平阴黄土剖面常量元素地球化学特征
徐树建1, 倪志超1,2, 丁新潮1,2     
1. 山东省水土保持与环境保育重点实验室, 临沂大学水土保持与环境保育研究所, 资源环境学院, 山东 临沂 276000;
2. 山东师范大学 地理与环境学院, 济南 250014
摘要: 为探讨末次冰期以来平阴地区的化学风化过程及其古环境,对山东地区平阴黄土剖面进行了常量元素含量、粒度及光释光年代测定。结果显示,平阴黄土常量元素含量呈现出SiO2 >Al2O3 >CaO >Fe2O3 >Na2O >K2O >MgO的变化特征;CIA指数、Na/K值及A-CN-K图解等表明平阴黄土风化程度较弱,处于未受化学风化-初等化学风化阶段;CIA指数、残积系数和风化淋溶系数指示其化学风化强度从剖面下部往上逐步变弱,黄土序列的常量元素含量及化学风化强度与粒度有很强的相关性。
关键词: 平阴黄土     常量元素     化学风化     地球化学     CIA    
Geochemical Characteristics of Major Elements of the Pingyin Loess in Shandong Province
XU Shu-jian1, NI Zhi-chao1,2, DING Xin-chao1,2     
1. Shandong Provincial Key Laboratory of Water and Soil Conservation and Environmental Protection, Institute of Water and Soil Conservation and Environmental Protection, College of Resources and Environment, Linyi University, Linyi, Shandong 276000, China;
2. School of Geography and Environment, Shandong Normal University, Jinan 250014, China
Abstract: In order to discuss the chemical weathering process and paleo-environment of the Pingyin area since the last glacial period, detailed analysis of major elements and grain size and optically stimulated luminescence(OSL) dating of the loess samples from the Longqiao profile in the Pingyin area, Shandong Province have been carried out in this study. Results show that contents of major elements of the Pingyin loess samples are varied in an order of SiO2 >Al2O3 >CaO >Fe2O3 >Na2O >K2O >MgO. CIA values, Na/K ratios, and features in the A-CN-K diagram suggest that the Pingyin loess was relative weakly weathered in a primary chemical weathering stage. CIA values, Eluvial Coefficients and Coefficients of Weathering and Eluviation suggest that the chemical weathering intensity is gradually decreased from the bottom to the top in the Pingyin loess profile. There are strong correlations among contents of major elements, intensities of chemical weathering, and grain sizes of the Pingyin loess.
Key words: Pingyin loess     major elements     chemical weathering     geochemistry     CIA(Chemical Index of Alteration)    

黄土是第四纪环境变化研究的重要载体,元素地球化学是其中的主要替代性指标。风尘堆积物中元素的地球化学行为及其组合特征较完整地记录了沉积物形成演变过程中的气候环境变化(刘东生,1985文启忠,1989陈骏和李高军,2011张玉柱等,2012郭媛媛等,2013Song et al.,2014Tang et al.,2015),通过研究黄土序列中元素迁移、富集程度与过程,可深入理解黄土沉积的古气候意义,有效揭示黄土-古土壤化学风化强度及古气候环境演变过程(程燕等,2003宋友桂和史正涛,2010汪海斌等,2011)。

地球化学指标在洛川黄土-古土壤(陈骏等,19972001)、西峰红黏土(陈旸等,2001)、西宁黄土、西峰黄土、吉县黄土(Jahn et al.,2001)、灵台黄土(Yang et al.,2006)、镇江下蜀土(李徐生等,1999Chen et al.,2008)、巫山黄土(张玉芬等,2014)、辽南周家沟黄土(张威等,2013)的研究中取得了丰富的结果。众多学者从粒度、磁化率、古地磁、黏土矿物、碳酸钙等物理指标(曹家欣等,1987张祖陆等,2004Xu and Wang,2011Peng et al.,2011Zhang et al.,2012)方面对山东地区的黄土进行了研究,而元素地球化学方面的研究虽有一些初步工作,但有待进一步深入研究。

不同气候条件与地质营力下存在着与之相应的沉积物,这些沉积物与它所处环境之间的关系是通过沉积矿物及元素组成的迁移或富集及元素重新组合产生新矿物来实现的,因此它们不但具有原岩的组成特征,而且记录着它们形成时的气候环境,这是利用元素地球化学的研究来探讨气候变化历史的理论基础(李徐生等,1999Chen et al.,2008)。沉积物中主要的造岩元素Si、Al、Fe、Ca等主要以氧化物形式存在。根据这些氧化物的表生地球化学性质,研究其百分含量在剖面中的变化规律,可以对古环境进行指示(李长安等,2013)。常量元素分析的基础是沉积物对母岩的主元素组合和重矿物组合的继承性。一般在短距离搬运和化学风化很弱的条件下,二者具有较好的可比性。在进行元素组合分析时,要考虑到搬运过程中的稀释作用,即应注意相对含量而非绝对含量。基于地球化学组成在不同粒级的差异,广泛深入研究不同地区的风尘物质对于揭示黄土物源以及环境变迁具有重要意义(Xiong et al.,2010Feng et al.,2011谢远云等,2012Liang et al.,2013)。山东地区的风尘堆积受海洋和大陆的双重控制,既有典型黄土的一般特征,又有独特的区域特点,其物质来源广泛,岩性特征与其他地区黄土有明显区别。本文重点对山东地区平阴黄土的粒度和常量元素地球化学特征进行研究,并与中国典型风成堆积物样品进行对比,以此探讨山东平阴黄土的常量元素地球化学特征及其与典型风成堆积物样品的异同以及这些特征对物质来源及其气候环境演化的指示意义。

1 样品与方法 1.1 研究区概况

研究对象是位于山东省济南市平阴县龙桥村的龙桥剖面(36°17′33.6″N,116°21′28.8″E,海拔约76.6 m)(图 1)。该剖面位于鲁中南山地丘陵的西北麓,属于温带大陆性季风气候,年平均气温为13.6℃,年平均降水量在500~900 mm。平阴黄土直接覆盖在黄河济南平阴段的二级阶地上,出露良好,自山麓向上逐渐变薄,沉积连续稳定,层序清楚。剖面上部厚700 cm,黄色粉砂,岩性均一,无层理和流水作用痕迹,垂直节理发育;下部为发育很好的钙结核层,结核直径为2~50 mm,质地坚硬,厚度150 cm左右,未见底。

图 1 山东平阴黄土剖面的地理位置 Fig. 1 Location of the Pingyin loess profile, Shangdong Province
1.2 研究方法与年代测定

野外样品采集从剖面下往上,以5 cm间隔连续取样,共采集样品141个。样品粒度在 Mastersizer 2000激光粒度仪上测量,仪器测量范围0.02~2000 μm,重复检测相对误差小于1%。元素含量由Panalytical Magix PW2403型X荧光光谱仪测得,以10 cm为间隔共分析70个样品。首先将自然风干样品研磨至200目(75 μm)以下,称取4.0 g样品加压制饼,然后上机测量,共获得28种元素和7种氧化物,重复测量误差在2%左右。粒度测试在山东省水土保持与环境保育重点实验室完成,元素测试在兰州大学西部环境教育部重点实验室进行。

平阴龙桥剖面共采集7个光释光年代样品,在青岛海洋地质研究所完成。

2 实验结果与分析

测量结果(图 2)表明,平阴黄土主要为末次冰期时期沉积的黄土,与黄土高原马兰期黄土相对应。

图 2 平阴黄土剖面光释光年代 Fig. 2 OSL dating results of the Pingyin loess profile

对于常量元素,XRF仪以氧化物形式给出,未测定元素以CO3给出。龙桥剖面CO3含量为1.17%~27.02%,平均为7.97%。含量校正公式为:

E校正值=E测量值/(1-CO3/100)

校正后,样品中含量在1%以上的氧化物有SiO2、Al2O3、CaO、Fe2O3、MgO、K2O、Na2O。其中SiO2含量最高,为42.92%~78.15%,平均69.28%;Al2O3含量为9.59%~13.57%,平均11.94%;CaO为1.38%~37.02%,平均7.76%;Fe2O3为3.34%~4.83%,平均4.00%;K2O为1.93%~2.84%,平均2.68%;Na2O为1.46%~3.43%,平均2.27%;MgO为0.95%~1.87%,平均1.44%。氧化物含量平均值的大小顺序为:SiO2>Al2O3>CaO>Fe2O3>K2O>Na2O>MgO,其中SiO2、Al2O3、CaO为最主要化学成分,元素含量总和平均为88.98%。样品中含量在0.1%~1%之间的氧化物有TiO2、P2O5和MnO。

表 1图 3显示,平阴龙桥黄土剖面中SiO2含量与UCC较接近,CaO的含量相对UCC显著富集,Al2O3、Fe2O3、MgO、K2O相对UCC亏损;Na2O则相对UCC(Taylor and McLennan,1985)显著亏损。化学风化过程中碱金属元素易淋溶迁移,稳定元素相对富集(乔彦松等,2010)。平阴黄土的常量元素特征与中国典型风成沉积物相似,但仍然存在一定的差异。平阴黄土Fe2O3含量与辽南黄土的马兰黄土、洛川黄土及伊犁黄土相似,相对UCC亏损,远小于宣城风成红土。CaO含量与伊犁黄土相似,比上陆壳及其他典型黄土沉积物明显偏高,这是因为测量时未进行醋酸前处理,样品中含有较多碳酸盐矿物。碳酸盐矿物易发生风化淋失再淀积(曹家欣等,1987),平阴龙桥剖面下部CaO含量较低,表明黄土堆积期,降雨量丰富,碳酸盐矿物淋溶迁移,在剖面底部形成坚硬钙结核层。Na2O含量与其他风成沉积物相似,明显低于UCC,但其含量明显高于洛川黄土(Xu and Wang,2011)及其他典型风成沉积物,表明其形成于较干冷的气候条件,成土作用弱,Na淋失较弱。以硅酸盐形式存在的CaO与Na2O化学性质相似,所以硅酸盐中的CaO也存在明显的亏损(曹家欣等,1987)。

表 1 平阴黄土剖面常量元素含量及其与其他典型沉积物的对比 Table 1 Major element compositions of the Pingyin loess and other deposits

图 3 平阴黄土剖面与中国各地风尘堆积常量元素及UCC比较 Fig. 3 Comparison among major element contents of the Pingyin loess, other aeolian deposits, and the UCC

结合镇江下蜀土(Chen et al.,2008)、宣城风成红土(李徐生等,1999)、辽南周家沟黄土等中国东部黄土分布区黄土发现SiO2含量明显比洛川黄土、UCC含量高,这与中国东部黄土分布区黄土粒度较粗,形成过程中近物质来源的加入有关;与SiO2相反,MgO、K2O、Al2O3含量比洛川黄土低,这与SiO2、Na2O在粗粒级中富集,MgO、K2O、Al2O3在细粒级中富集的粒度效应(陈骏等,1997陈旸等,2001Jahn et al.,2001Yang et al.,2006Chen et al.,2008)相吻合。因此,Na2O含量较高原因除气候干冷外,可能还与平阴龙桥剖面粒度较粗有关。常量元素含量的粒度效应对于中国东部黄土分布区黄土近物质来源有一定的指示意义。

常量元素氧化物含量在剖面中的变化记录了成壤过程中元素的迁移、富集情况(乔彦松等,2010)。在整个沉积序列中,CaO元素的变化较大;其他元素的变异系数较小。这是由于Ca离子半径较小,容易随土壤溶液或地表水迁移而发生淋溶与淀积。同时CaO较大的含量变化,对其他元素起到明显的稀释或富集作用(曹家欣等,1987),这种影响比元素本身的变化还有显著,因此单个元素绝对含量变化不能真实反映地层的化学风化程度。

3 化学风化程度判别指标

元素含量变化是探讨古气候变化的理想指标(Buggle et al.,2011)。由于元素淋滤迁移再淀积等因素的影响,单个元素绝对含量变化不能真实反映化学风化强度的演化特征,常采用元素含量或元素比值来衡量沉积物形成时的化学风化程度。

3.1 CIA指数和Na/K值

化学蚀变指数(CIA)可以指示长石风化成黏土矿物的程度,是衡量沉积物来源物质的风化作用和沉积后风化作用的代用指标(Nesbitt and Young,1982张玉柱等,2012郭媛媛等,2013Song et al.,2014),计算公式为:

CIA=[Al2O3/(Al2O3+CaO*+Na2O+K2O)]×100

式中,氧化物均为分子摩尔数,CaO*为硅酸盐相,CaO*的计算采用McLennan(1993)的方法。因而,CIA主要反映的是硅酸盐(主要是长石矿物)的风化,不存在元素迁移后再淀积的情况,能很好地反映沉积物形成时的化学风化情况。

研究表明(图 4),典型风成沉积物及UCC的CIA值为:上陆壳(Taylor and McLennan,1985)UCC平均值47.92,处于基本未受化学风化的阶段;伊犁盆地黄土(Song et al.,2014)平均值53.84(46~61),西宁(Jahn et al.,2001)黄土平均值59.39(58~61),西峰黄土(Jahn et al.,2001)平均值61.45(60~62),吉县(Jahn et al.,2001)黄土平均值62.11(60~65),洛川黄土(陈骏等,2001)平均值63.73(62~66),处于寒冷干燥气候条件下低等化学风化阶段;辽南周家沟黄土(张威等,2013)平均值66.09(60~80),灵台黄土(Yang et al.,2006)平均值66.56(60~70),洛川古土壤(陈骏等,2001)平均值67.36(65~69),西峰红黏土(陈旸等,2001)平均值69.11(67~72),陆源页岩CAS(Taylor and McLennan,1985)平均值70.36,镇江下蜀土(Chen et al.,2008)70.45(66~76),九江风成红土(Hong et al.,2013)84.85(75~89),处于温暖、湿润条件下中等的化学风化阶段;宣城风成红土(Hong et al.,2013)88.78(85~90),处于炎热、潮湿气候条件下强烈的高等化学风化阶段。各地区风成沉积物的化学风化程度变化幅度不大,基本随着温度升高降水增多呈增大趋势。

图 4 平阴黄土剖面与其他黄土剖面化学风化参数CIA与Na/K关系散点图 Fig. 4 The diagram of CIA values vs. Na/K molar ratios of the Pingyin loess and other aeolian deposits

平阴黄土CIA值为42.94~63.26,处于未受化学风化-初等化学风化阶段,变幅较大,这与伊犁盆地黄土46~61、辽南周家沟黄土60~80、九江风成红土75~89相似。平阴黄土CIA平均值(53.67)比UCC略高,接近伊犁黄土,低于洛川黄土、辽南周家沟黄土、宣城风成红土。UCC、平阴龙桥、伊犁盆地黄土、洛川黄土、辽南周家沟黄土、洛川古土壤、陆源页岩、镇江下蜀黄土、九江风成红土和宣城风成红土,CIA值依次升高,化学风化强度逐渐加强。

Na/K值可以衡量风化剖面中斜长石的分解程度(张玉柱等,2012郭媛媛等,2013)。平阴黄土样品Na/K值与CIA指数呈明显负相关(相关系数R=-0.95,n=70),UCC、平阴龙桥、伊犁盆地黄土、洛川黄土、周家沟黄土、洛川古土壤、陆源页岩、镇江下蜀黄土、九江风成红土、宣城风成红土,Na/K值依次逐渐减小,化学风化作用强度逐渐加强。Na/K比所反映的化学风化程度与CIA结果相吻合。从 图 4可以看出,不同风化程度的Na/K值趋势不同。同时随着CIA数值的增大,Na/K值减小趋势逐渐变缓,一方面说明斜长石数量减少,一方面说明钾长石开始风化。

3.2 A-CN-K三角模型

Nesbitt和Young(1982)提出了大陆化学风化趋势预测的A-CN-K(即Al2O3-CaO*+Na2O-K2O)三角模型图,用来反映沉积物化学风化趋势及主矿物成分变化(McLennan,1993Buggle et al.,2011)。

平阴黄土与洛川黄土、镇江下蜀土、辽南周家沟黄土等典型风尘堆积物的数据点分布于去Na、Ca的同一风化路径,与典型大陆风化趋势一致,说明平阴黄土风尘物质同样起源于广泛的上部陆壳,在风力搬运沉积过程中经过高度混合。平阴黄土数据点大致与CN-A连线平行,Na、Ca含量减少,K及Al含量基本不变,说明斜长石开始风化分解,Ca、Na淋溶迁移,钾长石含量相对稳定。与其他典型风尘堆积物相比,平阴黄土明显远离CAS,更靠近UCC,说明其化学风化强度远比其他典型沉积物及CAS要弱,处于脱Ca、Na的低等化学风化阶段。值得注意的是,平阴黄土化学风化强度最大的数据点与洛川黄土、辽南周家沟马兰期黄土极为相似,其化学风化程度相似,表明平阴黄土化学风化强度最强时与洛川黄土、辽南周家沟马兰期黄土想当。

图 4图 5可见,有3个数据点的CIA值明显比上陆壳UCC低,其CaO含量在25%以上,表明碳酸盐矿物含量较高可能会影响CIA、Na/K和A-CN-K三角模型图数据点的分布,偏离样品实际情况。

图 5 平阴黄土的A-CN-K三角图解 Fig. 5 Triangle diagram of A-CN-K of the Pingyin loess
3.3 化学风化程度与环境变化

残积系数反映沉积物成壤强度的大小变化,风化淋溶系数主要反映盐基的淋溶状况(张玉芬等,2014)。随着气候温湿程度加剧、风化程度加深,化学风化作用加强,CaO、Na2O、K2O淋溶迁移,含量减少,而Al2O3、Fe2O3含量相对稳定,故残积系数增大,化学风化淋溶系数减小。从 图 6可以看出,淋溶系数与Na/K值随着深度同步变化,呈现良好的相关性,与化学蚀变指数呈明显的反相关;残积系数与化学蚀变指数CIA呈现出高度的正相关。

图 6 平阴龙桥黄土的地球化学特征 Fig. 6 Geochemical characteristics of the Pingyin loess

剖面底部CIA值较高,残积系数较大,化学风化淋溶系数较小,说明风尘沉积物堆积初期,气候较湿润,化学风化强度较强,碳酸盐矿物被大量淋溶。与剖面CaO含量较低相吻合。剖面自下至上CIA值、残积系数呈现一种减弱的趋势,Na/K、退碱系数、淋溶系数呈现增加的趋势,表明气候逐渐变冷、源区化学风化强度逐渐减弱,这与末次冰期气候总体上变冷的大方向相吻合。

4 讨论

图 4图 5可见,CIA值越大,化学风化越强。从整体来看,中国黄土的化学风化强度基本上自西向东逐渐增加,同时由纬度引起的温度降雨量的增加,使得化学风化程度由北向南逐渐增强。

山东地区黄土位于中国黄土分布区的东部边缘,北连辽南黄土、河北与北京地区黄土,西邻内陆黄土高原,南有长江下游地区的下蜀黄土与之相对应(张祖陆等,2004)。然而山东平阴黄土的化学蚀变指数CIA平均值(53.67)远比辽南周家沟黄土(66.19)、洛川黄土(63.73)、镇江下蜀土(70.56)要低,同时平阴马兰黄土比辽南周家沟黄土、洛川黄土、镇江下蜀土同期马兰黄土变化范围要大。气候(特别是降雨量)是影响化学风化程度的主导因子(乔彦松等,2010),平阴黄土属于温带大陆性季风气候,与属于北亚热带季风气候的镇江下蜀土相比温度较低,降雨量较少,化学风化程度要低。然而与属于暖温带大陆性季风气候的洛川相比,平阴黄土降水量多,化学风化强度却比洛川黄土弱,这主要是因为近源物质颗粒较粗的粒度效应的结果。辽南地区黄土与平阴黄土同属于温带大陆性季风气候,在末次冰期时期气候条件相似,而两者化学风化程度相差较大,可能是因为末次冰期时期莱州湾裸露地层和黄泛平原提供了更多的近源物质的缘故。这与常量元素的粒度效应结论相吻合。

平阴黄土序列早期化学风化强度较高,与洛川黄土、辽南周家沟黄土重叠,说明末次冰期以前平阴黄土与洛川黄土、辽南周家沟黄土化学风化强度近似,受到近源物质的影响较小,其物质来源主要西北气流携带来的内陆粉尘物质,化学风化强度较高;随着末次冰期气候变冷,降雨量减少,风力加强,莱州湾裸露地层及黄泛平原松散沉积物为平阴地区提供了丰富的近物质来源,黄土砂粒含量明显增加,粒度变粗,化学风化强度减弱。7.7万年以来CIA指数逐渐下降,造成这种现象的原因除了末次冰期干冷的气候原因外,更主要的可能与大量近物质来源的加入有关。平阴黄土序列经历了较弱的化学风化,成土变化的影响远小于沉积分选和源区变化,CIA指数明显受到粒度变化的影响,不能有效指示剖面元素迁移富集规律及化学风化强度。

5 结论

(1)平阴黄土以SiO2、Al2O3、CaO为主,常量元素含量呈现出SiO2>Al2O3>CaO>Fe2O3>Na2O>K2O>MgO的变化特征。

(2)山东平阴黄土风化程度较弱,处于未受化学风化-初等化学风化阶段。CIA指数、残积系数和风化淋溶系数表明从剖面底部往上随着时间的推移化学风化强度逐步变弱。

(3)随着末次冰期气候变冷,风力加强,降雨量减少,近源松散沉积物为平阴地区提供了丰富的近物质来源,黄土砂粒含量明显增加,粒度变粗,化学风化强度减弱。

(4)CIA黄土序列的常量元素含量及化学风化强度明显受到粒度变化的控制。对于中山东地区的黄土,常量元素含量及CIA不能有效指示剖面元素迁移富集规律及化学风化强度。

参考文献
[1] Buggle B, Glaser B, Hambach U, Gerasimenko N, Markovi Ac'1 S. 2011. An evaluation of geochemical weathering indices in loess-paleosol studies. Quaternary International, 240(1-2):12-21
[2] Chen Y Y, Li X S, Han Z Y, Yang S Y, Wang Y B, Yang D Y. 2008. Chemical weathering intensity and element migration features of the Xiashu loess profile in Zhenjiang, Jiangsu Province. Journal of Geographical Sciences, 18(3):341-352
[3] Feng J L, Hu Z G, Ju J T, Zhu L P. 2011. Variations in trace element(including rare earth element) concentrations with grain sizes in loess and their implications for tracing the provenance of eolian deposits. Quaternary International, 236(1-2):116-126
[4] Hong H L, Wang C W, Zeng K F, Gu Y S, Wu Y B, Yin K, Li Z H. 2013. Geochemical constraints on provenance of the mid-Pleistocene red earth sediments in subtropical China. Sedimentary Geology, 290:97-108
[5] Jahn B M, Gallet S, Han J M. 2001. Geochemistry of the Xining, Xifeng and Jixian sections, Loess Plateau of China:Eolian dust provenance and paleosol evolution during the last 140 ka. Chemical Geology, 178(1-4):71-94
[6] Liang L J, Sun Y B, Beets C J, Prins A M, Wu F, Vandenberghe J. 2013. Impacts of grain size sorting and chemical weathering on the geochemistry of Jingyuan loess in the northwestern Chinese Loess Plateau. Journal of Asian Earth Sciences, 69:177-184.
[7] McLennan S M. 1993. Weathering and global denudation. The Journal of Geology, 101(2):295-303
[8] Nesbitt H W, Young G M. 1982. Early Proterozoic climates and plate motions inferred from major element chemistry of lutites. Nature, 299(5885):715-717
[9] Peng S Z, Zhu L J, Xiao G Q, Qiao Y S, Gao Z D, Chen D D. 2011. Magnetostratigraphy and provenance of the Qingzhou loess in Shandong province. Journal of Arid Land, 3(3):184-190
[10] Song Y G, Chen X L, Qian L B, Li C X, Li Y, Li X X, Chang H, An Z S. 2014. Distribution and composition of loess sediments in the Ili Basin, Central Asia. Quaternary International, 334-335:61-73
[11] Tang L, Wang X S, Zhang S Q, Chu G Q, Chen Y, Pei J L, Sheng M, Yang Z Y. 2015. High-resolution magnetic and palynological records of the last deglaciation and Holocene from Lake Xiarinur in the Hunshandake Sandy Land, Inner Mongolia. The Holocene, DOI:10.1177/0959683615571426
[12] Taylor S R, McLennan S M. 1985. The continental crust:Its composition and evolution:An examination of geochemical record preserved in sedimentary rocks. London:Blackwell Scientific Publications
[13] Xiong S F, Ding Z L, Zhu Y J, Zhou R, Lu H J. 2010. A-6 Ma chemical weathering history, the grain size dependence of chemical weathering intensity, and its implications for provenance change of the Chinese loess-red clay deposit. Quaternary Science Reviews, 29(15-16):1911-1922
[14] Xu S J, Wang T. 2011. Comparative study on the grain size characteristics of loess deposit both on Miaodao Islands and on the Laizhou Bay Plain and its implications for provenance. Procedia Environmental Sciences, 10:1869-1875
[15] Yang S L, Ding F, Ding Z L. 2006. Pleistocene chemical weathering history of Asian arid and semi-arid regions recorded in loess deposits of China and Tajikistan. Geochimica et Cosmochimica Acta, 70(7):1695-1709
[16] Zhang W G, Dong C Y, Ye L P, Ma H L, Yu L Z. 2012. Magnetic properties of coastal loess on the Midao islands, northern China:Implications for provenance and weathering intensity. Palaeogeography, Palaeoclimatology, Palaeoecology, 333-334:160-167
[17] 曹家欣, 李培英, 石宁. 1987. 山东庙岛群岛的黄土. 中国科学(B辑),(10):1116-1122
[18] 陈骏, 季峻峰, 仇纲, 鹿化煜, 朱洪兵. 1997. 陕西洛川黄土化学风化程度的地球化学研究. 中国科学(D辑), 27(6):531-536
[19] 陈骏, 安芷生, 刘连文, 季峻峰, 杨杰东, 陈旸. 2001. 最近2.5 Ma以来黄土高原风尘化学组成的变化与亚洲内陆的化学风化. 中国科学(D辑), 31(2):136-145
[20] 陈骏, 李高军. 2011. 亚洲风尘系统地球化学示踪研究. 中国科学(D辑), 41(9):1211-1232
[21] 陈旸, 陈骏, 刘连文. 2001. 甘肃西峰晚第三纪红黏土的化学组成及化学风化特征. 地质力学学报, 7(2):167-175
[22] 程燕, 张小曳, 鹿化煜, 王丹, 王旭龙. 2003. 最近140 ka以来黄土元素地球化学演化及其古气候意义. 海洋地质与第四纪地质, 23(3):103-108
[23] 郭媛媛, 莫多闻, 毛龙江, 郭伟民, 顾海滨. 2013. 澧阳平原岩板垱剖面地球化学特征与风化强度研究. 地理科学, 33(3):335-341
[24] 李徐生, 杨达源, 鹿化煜. 1999. 皖南风尘堆积序列氧化物地球化学特征与古气候记录. 海洋地质与第四纪地质, 19(4):75-82
[25] 李长安, 张玉芬, 熊德强, 周耀, 郭洁. 2013. "巫山黄土"常量元素地球化学特征. 地球科学, 38(5):916-922
[26] 刘东生. 1985. 黄土与环境. 北京:科学出版社
[27] 乔彦松, 赵志中, 王燕, 傅建利, 王书兵, 蒋复初. 2010. 川西甘孜黄土-古土壤序列的地球化学演化特征及其古气候意义. 科学通报, 55(3):255-260
[28] 宋友桂, 史正涛. 2010. 伊犁盆地黄土分布与组成特征. 地理科学, 30(2):267-272
[29] 汪海斌, 于英鹏, 刘现彬, 冯兆东. 2011. 黄土高原S1古土壤的地球化学特征及其对物源的指示. 第四纪研究, 31(2):338-346
[30] 文启忠. 1989. 中国黄土地球化学. 北京:科学出版社
[31] 谢远云, 孟杰, 郭令芬, 何葵, 康春国. 2012. 地球化学组成在不同粒级中的分布特征:以哈尔滨城市道路表土为例. 地理科学, 32(11):1397-1403
[32] 张祖陆, 辛良杰, 聂晓红. 2004. 山东地区黄土研究综述. 地理科学, 24(6):746-752
[33] 张玉柱, 黄春长, 庞奖励, 查小春, 周亚利, 周自翔, 王娟, 顾洪亮. 2012. 泾河下游古洪水滞流沉积物地球化学特征研究. 沉积学报, 30(5):900-908
[34] 张威, 董应巍, 于洋, 刘蓓蓓, 李永化, 李媛媛, 王美霞. 2013. 辽南黄土化学风化特点及其环境意义. 海洋地质与第四纪地质, 33(5):163-171
[35] 张玉芬, 邵磊, 熊德强. 2014. "巫山黄土"元素地球化学特征及成因和物源意义. 沉积学报, 32(1):78-84