颜色是土壤、沉积物最直观的物理特征之一^{[1, 2, 3]}。传统的Munsell色系是基于肉眼对样品与标准色卡图的对比，带有较大的主观性，且得到的结果是定性的，难以区别细微的色差，不适宜数学分析^{[4, 5, 6]}。1976年CIE提出了一个新的三维色系(CIE/L^*a^* b^*)^[7]，它用便于仪器测量的3个参数L^*a^* b^* (亮度、红度和黄度)定量地表示物体颜色，避免了肉眼判别颜色的主观误差，使得物体颜色的描述更加客观、精准和便捷^[3]。鉴于CIE/L^*a^* b^*对颜色表征的这些优点，它被逐渐引入到深海沉积^{[8, 9, 10, 11]}、黄土—古土壤序列^{[1, 12, 13, 14, 15, 16, 17, 18, 19]}、湖泊沉积^{[20, 21, 22, 23, 24]}、南方红土^{[3, 25]}、白垩纪地层^[26]和现代表土^{[1, 27]}的颜色定量描述，以探讨海相沉积和陆相沉积颜色变化可能蕴含的不同尺度的气候信息。

风成沙—古土壤沉积序列广泛分布于中国沙漠东部^{[28, 29, 30]}，如呼伦贝尔^[31]、科尔沁^[32]、浑善达克^[33]和毛乌素^[34]沙地的固定半固定沙丘上都普遍保留有风成沙—古土壤沉积序列，它们可能记录了千年尺度的气候变化信息。然而，到目前为止，对这类沉积序列在沙丘剖面呈现的颜色变化所蕴含的古气候信息报道还很少。本研究以浑善达克沙地为研究区，选择典型风成沙—古土壤沉积序列，用CIE/L^*a^* b^*色系定量表征颜色，探讨风成沙—古土壤沉积序列的颜色差异，影响颜色变化的主要因素以及颜色与古气候变化的可能关系。

浑善达克沙地(42°30′N~44°30′N，112°30′E~118°E)位于内蒙古高原东南部锡林郭勒大草原上，东西长约360 km，南北宽10~30 km，面积约2.14×10⁴ km² (图 1)。它是一个典型的固定、半固定沙地，沙丘相对高度为15~20 m。年降水量150~ 400 mm，年平均气温2℃~5℃，全年盛行西风、西北风^[35]。研究区全新世沉积以风成沉积为主^{[36, 37, 38]}，广泛分布于沙地内部及其周边地区。沙丘剖面普遍可见4~6层风成沙和3~5层灰—黑色砂质古土壤^{[36, 37, 38]}。

图 1 浑善达克沙地概况与剖面位置 Fig. 1 Schematic map of Otindag sandy land and the location of profile sites

图选项

基于研究区大量的野外考察，选择该沙地东、南、西、北和中部7个天然保留较好的沙丘剖面(图 1)：Profile 1 (42°53′ N，116°21′ E)、Profile 2 (43°15′ N，116°08′ E)、Profile 3 (42°40′ N，115°56′ E)、Profile 4 (42°55′ N，115°25′ E)、Profile 5 (42°52′ N，115°24′ E)、Profile 6 (42°56′ N，114°48′ E)、Profile 7 (42°43′ N，116°10′ E)。其中Profile 1~Profile 3和Profile 6为固定沙丘人工开挖后的自然剖面，其余3个剖面为固定沙丘被风蚀后残存的露头剖面。这7个剖面最显著的特征是风成砂层与灰/黑色砂质土壤相互叠覆交替出现(图 2)，砂质土壤最多可达4层；在灰/黑色砂质土壤中可见虫孔、白色钙质胶膜、土壤团粒状结构和明显的成壤特征，且砂质土壤层的机械组成明显细于砂层的颗粒组成。根据每个沙丘剖面风成砂层与灰/黑色砂质土壤层的颜色变化，从其底部依次向上采样，记录取样深度(图 2)。

图 2 浑善达克沙地沙丘剖面特征与样品采集 Fig. 2 Field paleosol-sand sedimentary profiles and sampling sites across Otindag sandy land

图选项

将沙丘剖面中采集的风成沙、砂质古土壤样品在自然条件下风干，随机分为4等份，分别测定颜色、粒度、有机碳和总氮。样品颜色的测定使用日本Minolta-CM2002光谱光度计测量，将样品粉碎均匀后，平铺在一个无色透明的玻璃板上，仪器可定量输出颜色数值，结果用L^* (亮度)、a^* (红度)、b^* (黄度)表示。粒度分析使用英国Malvern公司的Mastersizer 2000激光粒度仪测定，实验的前处理按照肖举乐等^[39]提出的步骤进行，其测量误差为±2%。样品有机碳的测定使用重铬酸钾—外加热法，按照GB7857—87标准进行测定，分析误差为±0.1%^[40]。总氮的测定使用开氏定氮法，按照GB7173—87标准进行测定，分析误差为±0.1%^[40]。

所有颜色参数的统计计算、相关分析、回归分析和均值的非参数检验(Cruskal-Wallis秩和检验)，都在IBM SPSS Statistics 20.0软件上运行，统计检验显著水平P＜0.05。

浑善达克沙地7个天然风成沙—古土壤沉积序列的颜色测定结果显示：L^*(亮度)值介于32.41~65.89，其中风成沙的L^*值为45.74~65.89，灰色古土壤的L^*值为41.27~49.81，黑色古土壤的L^*值为32.41~41.07；a^*(红度)值介于2.73~7.52，其中风成沙的a^*值为3.41~7.52，灰色古土壤的a^*值为3.31~7.50，黑色古土壤的a^*值为2.73 ~5.25；b^*(黄度)值介于-1.68~17.16，其中风成沙的b^*值为7.46~17.16，灰色古土壤的b^*值为2.75~8.84，黑色古土壤的b^*值为-1.68~1.84。整体上L^*值变化较大，而a^*和b^*值变化较小。

所有风成沙—古土壤沉积剖面中风成沙、灰色古土壤和黑色古土壤L^*a^*b^*均值的非参数检验(Cruskal-Wallis秩和检验)显示：风成沙、灰色古土壤和黑色古土壤的L^*、b^*值存在显著差异(P＜0.05)；但风成沙和灰色古土壤的a^*值差异不明显(P>0.05)，而它们与黑色古土壤的a^*值有显著性差异(P＜0.05)，见图 3。

图 3 风成沙—古土壤沉积序列中风成沙、灰色古土壤和黑色古土壤L*a*b*均值的非参数检验 Fig. 3 Nonparametric tests of mean L*a*b* value on the paleosol-sand stratigraphic profiles across Otindag sandy land

图选项

浑善达克沙地3个典型风成沙—古土壤沉积序列的颜色在垂直方向上显示出规律性的变化(图 4)：即随着剖面深度、风成沙和古土壤层的变化，L^*、a^*、b^*值都呈现相似的变化曲线；它们在风成沙层位处都为高值，在灰色、黑色古土壤层位处会突然降低，但在黑色古土壤层位处会降低得更小，整体上风成沙的L^*值>灰色古土壤的L^*值>黑色古土壤的L^*值，a^*值和b^*值亦是如此规律，但L^*变化值>b^*变化值>a^*变化值；沙含量(>63 μm，%)在垂直方向上也表现出与L^*、a^*、b^*值相似的变化趋势和特征(图 4)。然而，粉沙含量(63~4 μm，%)、黏土含量(＜4 μm，%)、有机碳(SOC)、总氮(TN)在垂直方向上却表现出与L^*、a^*、b^*、沙含量变化趋势和特征相反的曲线：它们的值在风成沙层位处都为低值，在灰色、黑色古土壤层位处会突然升高，但在黑色古土壤层位处会升高得更大，整体上黑色古土壤的SOC >灰色古土壤的SOC >风成沙的SOC，总氮、粉沙和黏土含量亦是如此规律(图 4，剖面见图 2)。研究区其余4个风成沙—古土壤沉积序列的L^*、a^*、b^*、粒度(沙、粉沙、黏土含量)、SOC、TN在垂直方向上也遵循上述的变化趋势和特征。

图 4 沙丘剖面的颜色参数与粒度、有机质、总氮的垂直变化 Fig. 4 Vertical distribution of color parameters, grain size, SOC and TN along the typical paleosol-sand profiles

图选项

浑善达克沙地风成沙—古土壤沉积序列颜色参数的相关分析显示，它们之间存在显著的正相关(表 1)，其中L^*与b^*的相关性最显著(r=0.949，n=36，P＜0.01)，而L^*与a^*的相关性最低(图 5)。L^*、a^*、b^*与沙含量(>63 μm，%)存在显著的正相关(n=36，P＜0.05)，与粉沙(63~4 μm，%)、黏土含量(＜4 μm，%)、有机质(SOC)、总氮(TN)存在显著的负相关(n=36，P＜0.05)，但L^*与沙、粉沙、黏土含量、SOC、TN的相关性和显著性更高(n=36，P＜0.01，表 1)。进一步的研究显示：风成沙—古土壤沉积序列的L^*值随有机碳、总氮、粉沙和黏土含量对数值的增加呈线性的下降趋势，随沙含量对数值的增加呈线性增加趋势(图 6)。

图 5 风成沙—古土壤沉积序列颜色参数间的相关关系 Fig. 5 Relationship between L* value and a*,b* value of paleosol-sand profiles

图选项

图 6 风成沙—古土壤沉积序列L*值与粒度、有机质、总氮的关系 Fig. 6 The functional relationship between L* value and grain size, organic carbon, total nitrogen of paleosol-sand profiles

图选项

浑善达克沙地广泛分布的风成沙—古土壤沉积序列，在野外露头剖面上表现为风成沙层与浅灰—灰黑色的砂质土壤层交互出现(图 2)，这是它们最直观的野外现象。定量颜色测定和数据分析表明所有研究剖面的L^*(亮度)和b^*(黄度)在风成沙、灰色砂质土壤、黑色砂质土壤层位处存在显著差异，但a^*(红度)在风成沙和灰色砂质土壤层位处的差异不明显(图 3,4)。这表明，CIE/L^*a^* b^*色系的3个颜色参数中只有L^*和b^*可以指示风成沙—古土壤沉积剖面的层位变化，而a^*对这类沉积剖面的岩性变化指示不敏感(图 4)，这与a^*(红度)可指示陆相沉积中的黄土—古土壤序列^{[14, 15, 16, 17, 18, 19]}、三趾马红土^[41]、南方红土^{[3, 25]}、白垩纪地层^[26]岩性变化的结论并不一致。土壤学研究表明，赤铁矿和针铁矿等铁氧化物是土壤主要的致色矿物^{[1, 42]}，赤铁矿呈红色，使土壤呈现红色，针铁矿呈黄色—棕色，使土壤呈现明亮的黄色^[43]。前人研究显示，只需1.7%的赤铁矿就可使土壤变红色，赤铁矿含量与土壤红度指数呈显著的正相关^{[44, 45]}；何柳等^[17]对黄土高原东南缘西峰和渭南黄土剖面的颜色研究也表明黄土—古土壤序列的a^*主要受控于铁氧化物的含量和种类。土壤中的赤铁矿和针铁矿主要有两种来源，一种来自对母质矿物的继承和保留，另一种来自土壤发育过程中含Fe硅酸盐矿物和其他Fe氧化物的转化^[46]。末次冰盛期(LGM)以来研究区的物源并未发生变化^{[47, 48]}，因而成壤作用应是影响这两种矿物在风成沙—古土壤沉积剖面上变化的主要因素。中国现代表土^[1]和西北干旱区现代地表沉积物^[27]颜色与气候的定量分析都表明在干冷的干旱—半干旱区，赤铁矿的形成会受到强烈抑制，红度a^*在干旱—半干旱区是不敏感的。地处中国北方干旱半干旱区的浑善达克沙地，LGM以来一直维系着草原或荒漠的景观^{[38, 49, 50]}，在这种气候环境条件下，成壤作用较弱，所生成的赤铁矿也较少，因此使得a^*值在风成沙—古土壤沉积序列上变化不明显。

浑善达克沙地风成沙—古土壤沉积序列颜色参数之间存在正相关关系(表 1)，其中L^*和a^*值的相关性最低(R²=0.26)，而L^*和b^*值的相关性最高(R²=0.9)，它们在每个沙丘剖面的垂直方向上也表现出相似的变化趋势和特征(图 4)，这充分说明L^*和b^*值可能受控于相同的气候因子，具有较为一致的致色物质。由于b^*在风成沙—古土壤沉积剖面上的变化值低于L^*(图 4)，因此讨论颜色参数L^*的影响因素将更有实际意义。回归分析显示：研究区风成沙—古土壤沉积序列的L^*与SOC、TN、黏土和粉沙含量存在显著的负相关关系，其相关系数R²分别为0.89、0.85、0.67和0.63 (图 6)；与沙含量存在明显的正相关关系(R²=0.55)，这说明风成沙—古土壤沉积序列的L^*会受有机质含量和土壤机械组成的影响，它能指示此类沉积剖面有机质和成壤作用的变化，这与深海沉积^{[8, 9, 10, 11]}以及陆相沉积中的黄土—古土壤序列^{[17, 18, 19]}、湖泊沉积^[24]的L^*能指示有机质含量变化的结论是一致的，杨胜利等^[1]对中国表土颜色的定量研究也支持该结论。对于风成沙—古土壤沉积序列的L^*与CaCO₃含量是否存在关联，需进一步的研究。

综上所述，影响浑善达克沙地风成沙—古土壤沉积序列颜色变化的主要因素是有机质含量和成壤作用，它们与颜色参数中的L^*关系最密切。

浑善达克沙地位于中国北方干旱半干旱的温带草原带内，在这种水分限制的草地生态系统中，降水量是维系植物生长和衰亡的决定性生态因子^{[51, 52]}，在干旱的年份，降水量的减少会导致植物有效水分利用的降低，诱发大量的多年生植物死亡^[52]。因此，全球和区域降水格局的变化会影响区域荒漠化程度的改变，这也是地球上主要沙漠收缩和扩张的主要驱动力^{[53, 54]}。研究区地处现代东亚夏季风的北部边缘^[55]，末次冰盛期(LGM)以来东亚夏季风降水曾在该区域经历多次的增强与衰弱^{[29, 30, 31, 32, 33, 34, 37, 38, 56, 57]}，这表明浑善达克沙地也经历了多次的收缩和扩张。当东亚夏季风增强时，伴随而来的雨带推进给研究区带来更多的降水，提高了区域植被的盖度和生产力，促进地表有机质的积累，在风成沙层上直接发育砂质土壤，有机质含量的增加使砂质土壤颜色变深^{[1, 12, 19]}，因而L^*(亮度)就表现的小；当东亚夏季风减弱时，雨带后撤导致降水量的减少会降低区域植被的盖度和生产力，诱发研究区发生荒漠化，以前被固定的沙丘又重新活化，沙地向外扩张，形成风成砂层，有机质含量的减少和沙含量的增加使其颜色变浅^{[1, 19, 24]}，因此L^*就变大；东亚夏季风降水雨带如此进退往复就形成了研究区现今的风成沙—古土壤沉积序列。因此，浑善达克沙地风成沙—古土壤沉积序列L^*的变化受控于东亚夏季风降水强弱的变化，但更深层次的驱动机制，如北半球冰盖变化、地球轨道参数诱发太阳辐射变化和大气CO₂浓度变化等假说^[30]，还有待进一步深入研究。

浑善达克沙地风成沙—古土壤沉积序列的L^*(亮度)可以指示其岩性的变化(图 3，4)，它与SOC、TN、黏土、粉沙含量存在良好的相关关系(表 1和图 6)，能指示此类沉积序列有机质含量的变化，它们都受东亚夏季风降水的驱动；这意味着L^*具有明确的气候意义，能反映区域降水量的变化。这与陆相沉积中黄土—古土壤序列^{[16, 17, 18, 19]}、湖泊沉积^{[20, 21, 22]}和西北干旱区现代地表沉积物^[27]的L^*可以反映区域降水变化的结论是一致的。需要说明的是，风成沙—古土壤沉积序列的L^*与区域降水量呈负相关关系，L^*的低值指示区域降水量增加，可能是东亚夏季风降水增强的信号。

浑善达克沙地风成沙—古土壤沉积序列的颜色参数(L^*a^*b^*)之间存在正相关关系，其中L^*和b^*的相关性最显著(R²=0.90)，而L^*和a^*的相关性最低(R²=0.26)；L^*a^*b^*在垂直方向上具有相同的变化趋势和特征，但L^*值变化较大，a^*和b^*值变化较小；L^*(亮度)可指示风成沙—古土壤沉积序列的风成沙层和古土壤层的变化，而a^*(红度)在这类沉积序列上变化不明显，其主要原因可能是受沉积环境和区域条件的影响。影响研究区风成沙—古土壤沉积序列L^*变化的主要因素与有机质含量和成壤作用有关，它们可能受东亚夏季风降水雨带进退的影响。L^*与有机碳、总氮、黏土和粉沙含量呈显著的负相关，其相关系数R²分别为0.89、0.85、0.67和0.63，低的L^*值与高的有机质、黏土含量有良好的一致性，说明L^*的低值可指示研究区有机质含量增加和成壤作用的加强，是东亚夏季风增强的信号，反之亦然。风成沙—古土壤沉积序列的L^*(亮度)可指示其岩性变化、有机质和黏土含量以及区域降水量的变化，它具有明确的古气候意义。中国沙漠东部的呼伦贝尔沙地、科尔沁沙地、浑善达克沙地和毛乌素沙地同处现代东亚夏季风降水的北部边缘地带，末次冰盛期以来东亚夏季风降水雨带曾在这些区域内经历多次的进退，造就了如今广泛分布的风成沙—古土壤沉积序列，在重建这些区域高精度的古气候时，颜色参数L^*是一种便捷、可快速获取的气候代用指标。