黄土粉尘搬运过程的动力学物理模型.pdf

书书书第 卷第 期 年 月第四纪研究 ，：文章编号 （）黄土粉尘搬运过程的动力学物理模型秦小光蔡炳贵穆燕宁波殷志强（中国科学院地质与地球物理研究所新生代地质与环境重点实验室，北京 ）摘要通过对黄土粉尘重力沉降过程的动力学分析，给出了黄土粉尘粒度分布的数学表达，讨论了粉尘沉积通量随搬运距离、粒径变化的物理过程，首次确定了估算粉尘搬运距离和风力强度的计算方法，为区分粉尘搬运距离和风力强度对粒度的影响及其它们在冰期间冰期中的差异提供了物理学的判别依据。分析结果表明：）在重力沉降作用下，粉尘沉降通量随搬运距离的变化服从几何分布，具有沉降通量随搬运距离的增加迅速减小，越粗的粉尘颗粒其沉降通量初值越大，同时下降速度也越迅速的特点；）如果用携粉尘气流的水平通量作代表风力强度，则粉尘搬运距离与粒度分布曲线上重力沉降部分最高点的粉尘沉降通量成反比，风力强度与该点对应粒径的平方和粉尘搬运距离成正比。因此，根据该点的粉尘沉降通量和对应粒径，可以估算粉尘的搬运距离和风力强度。根据上述理论对渭南阳郭中学 黄土 古土壤剖面进行了粉尘搬运距离和风力强度的估算：首先，从粒度分布中提取出 个对数正态分布的独立组份；然后利用粗粒组份的参数计算粉尘搬运距离和风力强度。分析结果表明粉尘搬运距离具有冰期近、间冰期远的特点，风力强度的变化则具有冰期弱、间冰期强的特点，时段的风力强度比 阶段大，但小于全新世适宜期，而 时期粉尘搬运距离并未明显减小，因此，可能黄土粒度的变化并非反映了冬季风的变化，而是反映了夏季风的变化，夏季风是通过影响粉尘源区来影响粉尘粒度的变化。主题词黄土粉尘搬运距离风力强度物理模型中图分类号 ，文献标识码第一作者简介：秦小光男 岁副研究员 博士第四纪地质与数学地质专业 ：国家自然科学基金项目（批准号：和 ）和中国科学院知识创新工程重要方向项目（批准号：）资助 收稿，收修改稿 前言在全球变化研究中，黄土 古土壤一直是人们研究东亚地区古气候变化的重点对象之一。在各种古气候替代 性指标 中黄土 粒度 被 人 们 所 广 泛 采用 。黄土的粒度变化最早被解释为受控于风力强度的 改 变，因 而 认 为 反 映 了 东 亚 冬 季 风 的 变化 。近年来人们开始认识到除风力强度外，源区距离、风化程度对黄土粒度也会产生影响，黄土粒度受到冬季风强度和沙漠进退的共同影响，认为黄土粒度应该作为季风 沙漠系统变化的替代性指标，而不是单一的风力强度指标 ，。最近人们开始尝试利用黄土不同粒径粒子的含量反映沙漠 黄土过渡带的进退变化，在假设风力稳定情况下拟合粉尘含量与源区距离的关系 。但是有什么因素影响源区距离以及它们又如何影响源区距离等问题仍然困扰着人们，由于缺乏对粉尘沉降动力学过程的了解，制约了人们对源区距离研究的深入。近年来 人 们 开 始 重 视 粒 度 构 成 的 古 气 候 意义 。我们也利用大气粉尘沉降动力学原理，对黄土中的 种悬浮组份进行了研究分析，对其成因和古气候意义进行了合理解释 ，为估算源区距离奠定了基础。本文正是以黄土粗、中、细三组份为基础，推导了粉尘搬运距离的物理模型，然后利用黄土粗粒组份的数字特征，对渭南地区 时期黄土粉尘沉积的源区距离、风力强度进行估算，并探讨了源区距离、风力强度与冬夏季风的关系。样品与预处理研究 剖 面 位 于 陕 西 省 渭 南 市 阳 郭 镇 中 学（，），属黄土高原南部。剖面高 ，以 间距、黄土 以 间距（部分 ）连续采集了 个样品。样品自然风干后用英国产 磁化率仪测量低频磁化率。对比研究剖面与渭南标准剖面的磁化率时间序列 ，获得主要的地层界限年龄控制点，然后根据年龄控制点深度与对应年龄拟合得到深度与年龄之间的非线性方程，最后插值得到期秦小光等：黄土粉尘搬运过程的动力学物理模型研究剖面的时间标尺。粒 度 分 析 的 前 处 理 采 用 比 较 彻 底 前 处 理 方法 ，。然后用英国 公司的 型激 光 粒 度 仪 完 成，仪 器 测 量 范 围 为 ，重复测试误差 。结果为等效体积百分比，等效于重量百分比。黄土粒度分布混合组份筛分采用了累积概率坐标的对数正态组份筛分方法 ，。黄土粉尘搬运过程的物理模型黄土与古土壤的粒度分布中普遍存在着粗、中、细 个组份，其动力学成因分析表明粗粒组份是重力沉降速度远大于大气湍流强度的粉尘组份，细粒组份是重力沉降速度远小于湍流强度的粉尘组份，而中粒组份是重力沉降速度与湍流强度相当的颗粒组份 。粉尘在水平风力带动下运移，并在重力影响下沿途沉降，由于沉降速度的差异，不同粒径粉尘的百分含量沿途不断改变，因此各组份的含量变化包含了粉尘搬运距离和风力强度的信息。而粉尘源区距离、风力强度作为古风场的重要指标对了解古大气环境特征、揭示古大气环流空间格局有着重要意义。粉尘沉降通量的运移模型设源区粉尘扬起后粉尘到达高度为 米，在单位面积上方含粉尘空气柱的体积为 ，有 ，空气柱中粉尘的初始总量为，初始粉尘浓度 。为了简化分析，假设为：）当底部近地表的大气粉尘沉降后，体积 上部的粉尘可对底部粉尘亏损迅速弥补，使粉尘浓度均一化。显然粉尘搬运距离越远，这个高度越大；）粉尘在体积 中均匀分布，同时没有侧向的粉尘扩散；）粉尘扬起后沿途只有沉降而没有新粉尘的加入；）从源区到沉降区平均风速不变，大气湍流强度的变化也可忽略不计；）一次沙尘暴过境等同于无数这样的空气柱过境，而所有空气柱的特征相同；）一个黄土 古土壤样品是多年内无数尘暴的粉尘堆积而成，而所有尘暴的特征都相同，样品反映了多年内大气粉尘的平均状态。另外，设粉尘沉降速度为 ，每个时刻的粉尘沉降时间长度为 ，水平平均风速为珔，这样单位面积上粉尘的沉降时间 单位长度 珔 珔。如果粉尘从源区到达沉降点的搬运距离为 ，历时 ，为时刻数，也即经过的单位面积数，有 珔和 。根据流体力学沉积通量原理，在开始沉降的第一个时刻，粉尘沉降量（）为：（）这时空气柱中剩余的粉尘量（）为：（）（）空气柱中剩余粉尘的浓度（）为：（）（）在时刻 ，粉尘沉降量（）为：（）（）空气柱中剩下的粉尘量（）为：（）（）（）（）空气柱中剩下粉尘的浓度（）为：（）（）在时刻 ，粉尘沉降量（）为：（）（）空气柱中剩下的粉尘量（）为：（）（）空气柱中剩下粉尘的浓度（）为：（）（）令 （珔 ），显然 ，公式（）可写成：（）（）令 （），则 ，是 的 函 数。表达了粉尘的搬运分选过程或沉积贡献率。如果忽略大气湍流沉降的影响，粉尘颗粒的沉降速度 等于由 公式确定的重力沉降速度 ，：()（）其中，是重力加速度，是空气粘滞系数，是颗粒 密 度，是 空 气 密 度，是 颗 粒 直 径，()。因此，（珔）（）这里 珔，（珔）。表示单位面积上方含粉尘空气柱（对地面粉尘沉降有贡献部分）在单位时间内的通过量，表示了含粉尘大气的水平通量，因此可以视为风力强度，单位是 。公式（）显示 与颗粒粒径的平方正相关，与风速和粉尘到达高度负相关。当 值很大时，于是公式（）变为：（）（）（）公式（）刻画了粉尘沉降通量与粉尘粒径、风第四纪研究 年速、粉尘到达高度和原始粉尘浓度的关系。粉尘沉降通量与粒径、搬运距离的关系我们分别考察沉降贡献率 随搬运距离 和粉尘粒径的变化。设各个粒径的粉尘原始浓度 均为相同的常数（一般等于 ），并取平均风速珔 ，。（）不同粒径粉尘沉降通量 随搬运距离 的变化（书书书图 ）对于给定粒径的粉尘，随搬运距离的增加，粉尘沉降通量呈指数形式迅速减小。越粗的粉尘颗粒其沉降通量初值越大，同时沉降通量下降速度也越迅速。如果将搬运距离取对数以观察近距离部分的细节，可以看到沉降通量的变化呈 字形，以内，粉尘沉降通量变化较小，之内沉降通量迅速降低，尤其是粗颗粒沉降通量下降迅速，而颗粒越细，其沉降通量的降低速率越小，搬运距离也越大。局部放大书书书图 ，显示粗颗粒沉降通量的降低速率快于细颗粒的，大约 以后，粗颗粒的沉降通量开始小于细颗粒的。可见 是黄土粉尘粒径具有明显变化的搬运距离，因此也是估算粉尘源区的最有效距离范围。的粉尘，搬运距离的增加而沉降通量变化很小，这说明 细粒粉尘受重力沉降的影响很小，是长期悬浮的 粉尘成 分。而 的粉 尘相比 的粉尘，在 以内其 值的减小幅度很小，近于稳定，说明在这个距离范围内 的粉尘含量可以作为衡量粗粒粉尘沉降量的标准。不同的风速对给定粒径颗粒沉降贡献率的影响，实际上类似于粒径的影响，公式（）表 明风速增加的影响等同于粒径减小的影响（书书书图 ），风速增大，粉尘沉降通量将降低，并以较低的水平维持更远的距离；相反，风速较小时，粉尘会在较短距离内以较高的粉尘沉降通量迅速沉降。同样粉尘可到达高度的影响也与风速相同，高度越大，粉尘沉降通量越小，能被携带的距离越大；高度越小，粉尘沉降通量越大，粉尘搬运距离越小。如果以时刻数 （或搬运距离 ）为变量，则由 的定义，可知 服从几何分布，其数学期望（即均值）为：（）（珔 ）（）它表示粒径为 的粉尘在给定风速和高度条件下能够被搬运的平均距离。该公式表明粉尘的平均搬运距离与风力强度成正比，与粉尘粒径的平方成图 不同条件下粉尘沉降通量随搬运距离的变化（）不同粒径粉尘沉降通量随搬运距离的变化（）对数坐标下不同粒径粉尘沉降通量随搬运距离的变化（）不同风速下，粉尘沉降通量随搬运距离的变化 （），（），（），反比。这个关系与我们的一般认识完全一致，可见该公式在理论上刻画了粉尘平均搬运距离、风力强度和粉尘粒径之间的物理学关系。（）不同搬运距离下粉尘沉降通量随粒径的变化（书书书图 ）期秦小光等：黄土粉尘搬运过程的动力学物理模型图 给定搬运距离下粉尘沉降通量的粒径分布（）近距离搬运下粉尘沉降通量的粒径分布（）长距离搬运下粉尘沉降通量的粒径分布 （），（）分析给定搬运距离上不同粒径粉尘的沉降通量变化，我们发现在粒度对数坐标图上，沉降通量曲线具有负偏的单峰分布特征，与黄土的粒度分布曲线总体形态极为相似，虽然没有出现 报道的黄土粒度分布中的细粒（中值粒径 ）和中粒（中值粒径在 之间）组份，但沉降通量曲线很好解释了黄土粒度分布基本形态，证明其粗粒部分含量变化形态是粉尘在风力搬运过程中受重力作用影响而自然分选形成的。另一方面也说明黄土粉尘中、细粒组份的主要形成原因还与大气湍流或其他物理、化学过程有关。另一个特点是随搬运距离的增加，粉尘沉降通量曲线不仅峰值相应减小，而且峰值对应的粒径也随之减小。这表明通过比较粒度分布曲线的峰值位置可以相应估算粉尘的搬运距离。（）粉尘搬运距离与风力强度模型因此如果对公式（）取导数，则书书书图 中沉降通量曲线峰值点的导数等于零，令该峰值点处的粒径为 和原始粉尘含量为，于是有如下公式：（）（）（）（）（）（）因此，结合 的定义，有：（珔）（）（）（）最后得到：（珔 ）（）（）由公式（）可见，风力强度 与粉尘搬运距离和粉尘粒径的平方成正比。分析风力强度与粉尘粒度和搬运距离的关系，显然，风力能够搬运的粉尘越粗，则风力强度越大，粉尘能搬运得越远，风力强度也越大，因此，公式（）客观地刻画了这种关系，是风力强度、搬运距离和粉尘粒径之间合理的物理学表达。该式实际上从公式（）中也可推出，这表明粒径为 的粉尘达到平均搬运距离时，相比其他粒径的粉尘，其沉降通量最大。将公式（）代 入（）式，得到下式：（）（）（）（）（）如果取粒径 ，对应粒径粉尘的沉降通量为，则有：（）（）公式（）表明，粉尘搬运距离与粒度分布曲线上最高点的粉尘沉降通量成反比。这样我们就可以根据粒度分布曲线上最高点的粉尘沉降通量，利用公式（）估 算 出 粉 尘 的 搬 运 距 离，再 利 用 公 式（）估算出风力强度。（）黄土粉尘参数估计由于实测的黄土粒度分布曲线包含了粗、中、细个组份，只有粗粒组份完全受控于重力沉降作用，因此必须先提取出粗粒组份，然后再利用粗粒组份的有关参数来估算搬运距离。公式（）显示沉降通量 随 粒 径 的 分 布 类 似 分 布，但 并 非 分布，而是服从一个以 为变量的伽马 分布。可见不管用对数正态分布还是 分布来拟合粉尘粒度曲线都只是伽马分布的一种近似。而伽马函数存在均值和方差，根据中心极限定理，其样本均值服从正态分布。公式（）显示粉尘沉降通量随粒径的分布主要受控于顶点粒径，搬运距离 主要影响沉降通量的幅度，没有其他影响沉降通量方差的变量。而实际观测表明，粉尘源区越近（即搬运距离越小），黄土粉尘的粒度分布越粗、峰值越高，如邙山的黄土主要来自北侧紧邻的黄河滩涂 ，因此其粒度粗、第四纪研究 年峰值高、峰区宽；相反，粉尘源区越远，黄土粉尘的粒度分布越细、峰值越低、峰区越宽，如位于黄土高原南端的渭南黄土 。实际上黄土粉尘必然来自一个相当面积的区域，不会是一个点，而且伴随粉尘的搬运，沿途会不断有新的粉尘加入，包括其他源区的粉尘混入，因此，粉尘搬运距离越大，其物源区范围必然越大。由于不同搬运距离的粉尘具有不同的峰值粒径，大的源区范围必然导致峰值粒径在较大区间范围内波动。根据中心极限定理，其均值将服从正态分布，结果是造成峰值粒径方差的增加，但粒度分布的实际峰值粒径仍然可视为代表了粉尘的平均搬运距离。由于公式（）只是描述了给定原始浓度的粉尘在搬运过程中的分选效应，并未刻划所有来源粉尘的综合效应，因此，一般来说直接用公式（）拟合实测粒度曲线，并不能保证获得最佳的拟合效果。目前用来提取粒度组份的分布函数主要有对数正态分布和 分布 ，如前所述这两种分布都是一种近似，理论上根据中心极限定理，采用对数正态分布能够更好地提取峰值粒径，因此，我们采用对数正态分布提取粒度的组份特征可以得到细、中、粗 个组份的含量、均方差 和均值（，）等参数，而粗粒组份是估算粉尘搬运距离和风力强度的依据。由于粒度测量中得到的数据是各粒级粉尘颗粒的体积百分含量，并非绝对沉降量，因此，需要进行换算才能利用公式（）和（）来估算粉尘搬运距离和风力强度。黄土剖面中第 个样品点的粉尘沉降速率可用公式 计算，这里 是采样间距，是样品年龄间隔。如果假设每年的沙尘暴次数相同，样品粗粒组份最高点的绝对沉降量则可用下式表达：槡 （）将 作为，代入公式（），即可计算出给定初始粉尘浓度 情况下的粉尘搬运时间 （或距离）。计算中需要先根据现代粉尘源区的距离确定合理的 值，然后再应用于整个剖面的计算中。结果与讨论图 是渭南剖面 万年以来 个组份含量变化曲线，其中书书书图 中粉尘搬运距离曲线显示末次冰期中 大约相当于 的时期源区距离最小，时期粉尘源区距离略大于 时期。自 以来粉尘源区距离持图 渭南黄土剖面 以来粉尘的三组份特征与搬运距离、风力强度（）个组份的百分含量（）个组份的中值粒径（）粉尘搬运距离与风力强度 （），（），（）续增加，的末次冰盛期并非源区最近时期。进入冰消期后，粉尘源区已远于深海氧同位素第三阶段时期（）。全新世适宜期源区距离最大，其中有几次快速短暂的源区距离急剧减小事件。适宜期以后，大约 以来，源区距离波动中持续减小，并延续到现在。这种源区距离的变化特征基本符合冰期间冰期沙漠进退的基本事实，同时又反映出一些全新的信息，即冰盛期 并非粉尘源区最近的时期。风力强度的变化则展现了一种出乎意料的特征，即冰期风力强度弱，而全新世中、早期强。这与过去人们所想象的“冰期风力强，间冰期风力弱”相差甚远。对比源区距离和风力强度两条曲线，可以发现它们具有很大的相似性。事实上我们很容易理解，因为只有风力强时才能将粉尘带到更远的地方，期秦小光等：黄土粉尘搬运过程的动力学物理模型公式（）反映了这种关系。因此，显然冰期粉尘粗的现象实际上是由于粉尘源区近所造成。分析风力强度曲线的细节，可以发现 时期，风力强度确实有明显增强，出现了一个峰值，说明 时期虽然源区距离没有显著缩小，但风力强度显著加强，导致了这个时期的黄土快速堆积和粒度变粗。另一个细节是全新世晚期，不仅源区距离在缩小，而且风力强度也有明显加强，形成了一个小的峰值。源区距离和风力强度的上述特点表明，源区距离的变化是冰期间冰期旋回中最明显的变化，是导致粉尘粒度变化的主要原因。而风力强度并非冰期强、间冰期弱，相反具有冰期弱、间冰期强的特点，因此，风力强度不是粉尘粒度冰期间冰期变化的主要影响因素，但对一些短尺度的气候事件有显著影响。而造成粉尘源区距离变化的实际原因是源区植被的变化，当原来的粉尘源区（如毛乌素沙漠）开始生长大量的植被后，这里就不再是粉尘的主要源区，粉尘源区必然向环境仍然恶劣的西北方向迁移，造成粉尘源区距离增加；而原来的草原区如果植被退化，裸露出下伏沙土，就会成为新的粉尘源区。因此，源区距离的变化实际反映了植被的进退变化，而植被的进退在中国北方受控于夏季风降雨量的变化，降雨量增加必然导致植被的茂盛，因此，源区距离的变化应该是夏季风强弱的反映。这似乎告诉我们黄土粒度的变化并非反映了冬季风的变化，而是反映了夏季风的变化，夏季风是通过影响粉尘源区来影响粉尘粒度的变化。这个结论颠覆了过去黄土粒度变化反映了冬季风强度变化的观点，但具有合理的大气动力学理论基础，也完全能够解释各种已有的观测结果，我们认为是可以接受的，但还需要其他资料的进一步证实。上述粉尘搬运距离模型是一个基于大气动力学的理想模型，虽然较好解释了黄土粒度分布的基本特征，并给出了搬运距离和风力强度的量化估算公式，但仍然存在一些不足。首先，模型没有考虑粉尘混合问题，粉尘不仅是相当范围内源区物质的高度混合物，而且在搬运过程中还会不断有新的粉尘加入，因此公式中的粉尘浓度除了沉降去除外，还应该有起尘物质的加入，但由于无法获取新加入粉尘的信息，我们忽略了这部分粉尘的影响，这势必造成模型中的伽马 分布不能完全反映粉尘粒度的实际分布特 征；其 次，不 管 采 用 对 数 正 态 分 布，还 是 分布来拟合粉尘的粒度分布都只是实际粒度分布的一种近似，因此，拟合的结果本身也会带来误差；第三，原始粉尘浓度 对于不同粒径的粉尘是不一样的，我们的敏感性试验结果显示，不同粒径原始粉尘浓度的差异似乎对沉降通量分布形式的影响不大，因此，计算中我们假设沉积分异作用是导致黄土粒度变化的主要原因，由此设定 是一个固定值。但这个假设还是会对结果产生一定的影响。实际上如果 是粒径的函数，并服从某种分布，则公式（）中的 也是粒径的函数；第四，我们假设一年中的沙尘暴次数不变，也可能导致估算结果产生偏差。尽管如此，我们的分析仍然第一次从理论上明确了粉尘搬运距离与风力强度的关系，提供了各种地质观测和分析的物理学依据，对于正确理解冰期间冰期中东亚季风的变化具有有益的启示作用。结论（）在重力沉降作用下，粉尘沉降通量随搬运距离的变化服从几何分布，具有沉降通量随搬运距离的增加迅速减小，越粗的粉尘颗粒其沉降通量初值越大，同时下降速度也越迅速的特点。（）黄土粒度分布的基本型式主要是粉尘重力沉降导致的沉积分异结果。随搬运距离的增加，粉尘沉降通量的粒度分布曲线不仅峰值相应减小，而且峰值对应的粒径也随之减小，因此通过比较粒度分布曲线的峰值位置及其粉尘含量可以估算相应的粉尘搬运距离。（）粉尘搬运距离与粒度分布曲线上重力沉降部分最高点的粉尘沉降通量成反比，风力强度与该点对应粒径的平方和粉尘搬运距离成正比。因此根据该点的粉尘沉降通量和对应粒径，可以估算粉尘的搬运距离和风力强度。（）粉尘搬运距离具有冰期小、间冰期大的特点，自 以来粉尘源区距离持续增加，冰盛期 并非粉尘源区最近的时期，全新世适宜期源区距离最大，大约 以来源区距离持续减小，并延续到现在。（）风力强度的变化具有冰期弱、间冰期强的特点。时期，粉尘源区距离没有明显减小，但风力强度明显大于深海氧同位素 阶段时期。全新世中、早期风力强度大，晚期风力强度有明显小幅度加强。（）源区距离的变化是冰期间冰期旋回中最明显的变化，而风力强度冰期弱、间冰期强的特点表明它不是粉尘粒度冰期间冰期变化的主要影响因素，但对一些短尺度气候事件有明显关系。因此可能黄第四纪研究 年土粒度的变化并非反映了冬季风的变化，而是反映了夏季风的变化，夏季风是通过影响粉尘源区来影响粉尘粒度的变化。致谢本文在撰写过程中得到了旺罗的有益提示，肖举乐提供了粒度实验室，刘秀平、曹瑜参加了野外工作，特予致谢。参考文献（）刘东生等黄土的物质成分和结构北京：科学出版社，：，刘东生等著黄土与环境北京：科学出版社，：，：丁仲礼，余志伟，刘东生中国黄土研究新进展（三）时间标尺第四纪研究，（）：，（）：，（）：丁仲礼，杨石岭，孙继敏等 前后大气环流重构的黄土红粘土沉积证据第四纪研究，（）：，（）：孙东怀，鹿化煜晚新生代黄土高原风尘序列的粒度和沉积速率与中国北方大气环流演变 第四纪研究，（）：，（）：熊尚发，朱园健，周茹等 白水黄土 红粘土化学风化强度的剖面特征与粒度效应第四纪研究，（）：，（）：，（）：，：，（）：，：，（）：，：，：孙东怀，鹿化煜，等 中国黄土粒度的双峰分布及其古气候意义沉积学报，（）：，（）：张小曳亚洲粉尘的源区分布、释放、输送、沉降与黄土堆积第四纪研究，（）：，（）：张小曳有关中国黄土高原黄土物质的源区及其输送方式的再评述第四纪研究，（）：，（）：强明瑞，陈发虎，周爱锋等苏干湖沉积物粒度组成记录尘暴事件的初步研究第四纪研究，（）：，（）：张崧，勒春胜等 年 月 日北京降尘的粒度分布与磁学特征第四纪研究，（）：，（）：秦小光，蔡炳贵，刘东生深海氧同位素第 阶段时期东亚季风区大气湍流特征第四纪研究，（）：，（）：秦小光，吴金水，蔡炳贵等全新世时期北京 张家口地区与黄土高原地区风成系统的差异第四纪研究，（）：，（）：刘嘉麒，陈铁梅，聂高众等 渭南黄土剖面的年龄测定及十五万年来高分辨时间序列的建立第四纪研究，（）：期秦小光等：黄土粉尘搬运过程的动力学物理模型 ，（）：鹿化煜，安芷生前处理方法对黄土沉积物粒度测量影响的实验研究科学通报，（）：，（）：（），（）：赵鹏大，胡旺亮，李紫金矿床统计预测北京：地质出版社，：，：，：，：，殷志强，秦小光，李玉梅等源区距离对黄土粒度多组份分布特征的影响中国地质，（）：，（）：（，），