“12·14”山东暴雪过程降水相态的多源观测分析

基于毫米波云雷达、降水天气现象仪、风廓线雷达、双偏振雷达、自动气象站等多源观测资料及欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析（ECMWF Reanalysis v5,ERA5）资料,对2023年12月14日山东暴雪过程的环流背景和降水粒子微物理参数进行分析,探讨新型观测资料在降水相态监测与预报中的应用。结果表明：（1）此次过程受高空槽、低空西南急流和江淮气旋影响,伴随地面气温下降和中层暖层消退,出现雨、雨夹雪、冰粒和雪等相态。（2）降水天气现象仪探测发现,雪和雨的下落末速度均较小,雪粒子直径超过8 mm,雨粒子直径大多在4 mm以下。（3）毫米波云雷达观测到反射率因子、径向速度、谱宽和垂直液态水含量降低时,雨夹雪转为雪。（4）风廓线雷达显示雨夹雪和冰粒阶段对应强的低空西南急流和最大垂直速度（4~5 m·s^-1）,转雪时3 km以下垂直速度降低至2 m·s^-1左右。（5）相关系数（C_c）、差分反射率（Z_DR）和水平极化反射率因子（Z_h）等双偏振参量可判断融化层亮带,亮带的下降和消失可指示此次过程雨向雪的转换。     

2021年2月黄河中下游两次暴雪的相态转换特征及成因

利用常规高空资料、地面加密自动站、双偏振多普勒天气雷达、微波辐射计与ERA5再分析数据等多源资料,分析了2021年2月下旬黄河中下游两次暴雪过程的相态演变及形成机理。结果表明:两次过程的大尺度影响系统基本一致,只是影响系统的强度和位置不同导致两次过程存在些许差异。两次过程均存在相态转换,过程Ⅰ中存在相态逆转(由雨转雪再转雨),而过程Ⅱ中只存在雨转雪的转换。在太行山以西的山区,当地面2 m气温低于0.5 ℃时,降水相态以雪为主,在0.5~1 ℃之间时,多为雪或雨夹雪并存;在平原地区,当2 m气温为1~2 ℃时,降水相态为雨或雨夹雪,在0~1 ℃之间时,则为雪与雨夹雪并存,低于0 ℃时,降水相态为雪。在降雨阶段,双偏振雷达产品相关系数(CC)值约在0.98以下,差分反射率(Z_DR)在0.6 dB以上,差分传播相移率(K_DP)值约在0.2 (°)·km^-1以上;在降雪阶段,CC值在0.98~0.99之间,Z_DR值在0.2~0.8 dB之间,K_DP 值约在0.2 (°)·km^-1以下;但在降水相态由雨转雨夹雪时,Z_DR、CC与K_DP 值没有明显变化。     

山东一次典型持续冰雹过程的闪电活动和雷达特征分析

2020年6月1日下午至夜间,山东部分地区出现降雹天气,利用卫星和雷达以及闪电探测资料详细分析此次降雹过程的闪电和雷达参量特征。结果表明：（1）鲁北和鲁中雹暴闪电活动的主要差异为：鲁中雹暴的对流强度强于鲁北雹暴,其正地闪比例和正地闪平均强度明显高于鲁北雹暴。其相同特征为：在开始阶段正地闪比例和云闪占总闪电的比例较高;在地面降雹开始之前总闪电频数均出现跃增,其峰值提前降雹6~18 min,同时云闪表现更加活跃;地面出现降雹之后,云闪频数快速下降,地闪占总闪电比增加,甚至超过50%。（2）雹暴的闪电活动均发生在云顶温度小于-50 ℃的云区内,且密集出现在云顶温度小于-60 ℃的范围内。闪电主要发生在30 dBZ以上回波区域,密集的云闪对应强对流区,表明云内垂直运动剧烈。负地闪与差分相移率K_DP＞0.5 （°）·km^-1和差分反射率Z_DR＞2.0 dB的区域对应关系非常一致,表明负地闪与强降水区密切相关。雹暴回波穹隆区的Z_DR值较高、零滞后相关系数（C_c）值较低,很好地指示了上升气流区。较高的水平反射率因子（Z_H）和较低的C_c以及较低的Z_DR区域对应冰雹粒子区域。闪电很好地对应于云内霰、湿雪、冰雹、干雪等大的冰相粒子区。（3）0 ℃以上Z_DR≥1.0 dB的体积、Z_H≥30 dBZ的体积、Z_H≥30 dBZ的冰水含量与总闪电频数的时间演变趋势基本一致,其与总闪电频数的相关系数分别为0.756、0.780和0.710,进一步证实大的冰相粒子在起电过程中发挥着主导作用。     

2023年和2005年山东两次极端海效应暴雪的对比

利用常规气象观测资料、降水天气现象仪资料、美国国家海洋和大气管理局（National Oceanic and Atmospheric Administration,NOAA）资料、欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析（ECMWF Reanalysis v5,ERA5）资料等对2023年12月山东半岛的海效应特大暴雪和2005年12月持续性海效应强降雪过程的高空形势、海气温差、低层切变线、大气水凝物等进行了综合分析。结果表明：（1）亚欧大陆强大的经向高压脊和脊前冷涡环流是发生海效应暴雪的大尺度环流背景。2005年高压脊宽广、稳定少动是海效应强降雪持续时间长的主要原因,2023年高压脊、冷涡的经向特征更加显著,冷空气爆发力强。（2）2005、2023年渤海海面温度较常年偏高;2023年偏高2.5 ℃的范围更广,12月20—21日海气温差超过30 ℃。（3）出现暴雪时,山东半岛北部存在偏西风和北—东北风之间的风场辐合;荣成站出现暴雪时,低层需要更强劲的引导气流。（4）2023年12月20—21日,云体主要由冰晶和雪晶构成,较多雪晶位于上升气流上方,与冰粒子分布区域重叠,说明除了水汽凝华外,冰雪晶粒子之间的聚合作用对雪粒子的增长有很大帮助,聚合产生的枝状雪粒子有利于积雪深度增大,文登站雨滴谱也表明,21—22日雪粒子直径偏大的特征更明显。     

“23·12”山东半岛特大海效应暴雪特征及成因

采用地面气象观测站、多普勒天气雷达、闪电、积雪深度人工加密观测资料、常规观测及欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析（ECMWF Reanalysis v5,ERA5）资料,对2023年12月15—22日山东半岛特大海效应暴雪过程的降雪特征及极端性成因进行了分析。结果表明：（1）此次过程有4站积雪深度突破本站历史极值,有1站2 d的日降雪量为山东半岛海效应降雪有气象记录以来的最大值,文登积雪深度达74 cm,超过山东所有国家级地面气象观测站纪录,是一次极端海效应暴雪事件。（2）欧亚中高纬度阻塞形势下两次异常强冷空气持续影响渤海和山东半岛地区,850 hPa温度最低降至-21~-20 ℃,冷空气强度明显强于往年12月海效应暴雪过程,造成降雪持续时间长、累计降雪量大。异常强冷空气是此次极端暴雪过程产生的关键因素,渤海海面温度（简称“海温”）异常偏高是有利的海温背景。（3）冷空气强、海温偏高造成海气温差偏大,700 hPa以下产生对流不稳定,使得降雪强度大;强降雪发生在海气温差快速增大阶段。（4）925 hPa以下存在来自渤海的北—东北风与内陆地区的西北风构成的切变线,产生强上升运动,切变线长时间维持形成“列车效应”。（5）主要降雪时段强垂直上升运动、高相对湿度层的温度为-20~-12 ℃,适宜树枝状冰晶形成和维持,有利于产生大的积雪和降雪含水比;2 m气温持续低于-5 ℃,0 cm地温在降雪开始时即降至0 ℃以下,且两次强降雪过程仅间隔1 d,均有利于降雪累积产生极端积雪深度。     

“12·14”山东暴雪过程的极端性特征及成因

利用国家级地面气象观测站、风廓线雷达、X波段双偏振相控阵雷达等多源观测资料和欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析（ECMWF Reanalysis v5,ERA5）资料,总结了2023年12月13—14日山东大范围暴雪、局地大暴雪（简称“12·14”暴雪）极端性的特征和成因,并与2021年11月7日山东极端暴雪过程（简称“11·7”暴雪）对比分析了降雪量和雪水比差异的原因。结果表明：（1）典型的暖平流型天气形势是产生极端暴雪有利的环流背景条件,低层切变线和风速辐合区在鲁西北叠加,形成强烈而持久的上升运动。（2）低空急流异常偏强,降水强度不仅与低空急流的强度有关,而且与其厚度有关。当3.0 km高度保持低空急流的强度时,10 m·s^-1风速到达的高度越低,降雪强度越大。（3）700 hPa比湿超过4 g·kg^-1、850 hPa比湿超过3 g·kg^-1的持续时间均长达10 h以上,为极端暴雪过程提供了充足的水汽。850 hPa比湿和比湿平流远远高于“12·14”暴雪过程,是“11·7”过程最大累计降雪量大于“12·14”过程的原因之一。（4）“12·14”暴雪过程垂直运动旺盛,最大上升速度位于不稳定层顶的前沿,处于-20~-10 ℃层,有利于树枝状冰晶生长,降雪效率高。对流层整层温度低于0 ℃,雪花下落过程中没有融化,雪水比高,积雪深度大。“11·7”暴雪过程初期最大上升运动中心高度低,形成更多柱状冰晶,经过暖层时融化,雪水比低,积雪深度小于“12·14”暴雪过程。     

“12·14”山东暴雪过程积雪深度特征及成因

利用常规观测、积雪深度逐时加密观测资料和欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析（ECMWF Reanalysis v5,ERA5）资料,对2023年12月13—15日山东一次极端暴雪天气过程积雪特征及其成因进行分析,得到以下结论：（1）此次过程是一次江淮气旋暴雪天气过程,具有持续时间长、降水相态复杂、基础温度低、降温幅度大和积雪深度厚等特征。（2）最大小时新增积雪深度可达8 cm;过程平均雪水比为0.7 cm·mm^-1,呈“西大东小” 的分布特征。（3）有积雪的站近地面温度从开始降雪到地面产生积雪,气温和雪面温度均呈下降趋势,0 cm地温在降雪前期降温明显,积雪形成后地温不再明显变化。无积雪的站在整个降雪时段内近地面温度可分为4种情况。（4）雪水比随气温变化最明显;积雪形成之后地温对雪水比大小影响不大;当雪水比小于0.75 cm·mm^-1时,雪水比随雪面温度降低而增大,当雪水比大于0.76 cm·mm^-1后,雪面温度不再有明显变化。     

中国黄渤海真江蓠(Agarophyton vermiculophyllum)群体遗传多样性研究

真江蓠是原产于西北太平洋的重要经济红藻。我们利用10对微卫星引物检测中国黄渤海地区真江蓠的群体遗传多样性和结构。10个微卫星位点在12个群体中共检测到65个等位基因,每个位点的等位基因（N_a）为1~28,有效等位基因（N_e）为1.0~9.6。每个群体的平均等位基因（N_a）、平均有效等位基因（N_e）、平均香浓指数（I）、平均观察杂合度（H_o）和平均预期杂合度（H_e）分别为2.4、1.6、0.419、0.133和0.227,显示较低的群体遗传多样性。中国黄渤海12个真江蓠群体间遗传分化较大（F_st=0.398 7）,基因流有限（N_m=0.377 1）,近交系数为正（F_is=0.391 3,F_it=0.634 0）,表明可能存在近交和杂合子缺失现象。Structure和UPGMA系统进化分析一致将12个群体分为两个遗传组,并在黑石礁群体（HS）和石岛群体（SD）中发现明显的遗传混杂现象。AMOVA分析显示遗传变异主要来自于群体内（73.27%）。该研究可为黄渤海地区真江蓠自然资源保护和管理提供科学依据。     

济宁S波段双线偏振多普勒雷达数据质量分析

以济宁S波段双线偏振多普勒雷达观测的层状云降水资料为例,分析了信噪比(signal-to-noise ratio,SNR)、相关系数(ρhv)对偏振参量的影响以及近距离双偏振参数的特征。利用雷达最高扫描仰角的层状云降水资料,对小雨和干雪两种自然目标物差分反射率因子(Z_DR)的订正效果进行分析。研究表明:(1)SNR小于21 dB时,Z_DR、ρ_hv、差分传播相位偏差(ΔΦ_DP)的平均值和标准差随SNR的变化波动明显,Z_DR、ρ_hv、差分传播相位(Φ_DP)不稳定,数据质量差。SNR大于21 dB时,Z_DR、ρ_hv、Φ_DP才具有可信度,且随着SNR的增大,Z_DR、ρ_hv、Φ_DP的稳定性越来越好。(2)SNR大于21 dB、ρ_hv＜0.965时,ΔΦ_DP的平均值和标准差随ρhv的变化而波动剧烈;SNR大于21 dB、ρhv＞0.965时,ρhv越大,Φ_DP的值越稳定。(3)受地面噪声干扰、近距离地物以及天线旁瓣的影响,0~15 km的近距离偏振参数不稳定,数据质量差,业务应用中需特别注意。(4)利用小雨和干雪两种自然目标物对ZDR系统误差进行标定,通过对比发现,小雨粒子的ZDR系统误差的平均值和标准差比干雪粒子小很多,且小雨粒子ZDR系统误差订正曲线比干雪粒子稳定;小雨粒子相比于干雪粒子,其ZDR系统误差订正效果更稳定;经过系统误差订正后的降水粒子的ZDR值较订正前更加合理。     

“12·14雷打雪”过程的闪电特征及成因

2023年12月14日在山西和河北南部、河南北部、山东西部出现一次大到暴雪过程,并伴随频繁的闪电活动。利用欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析（ECMWF Reanalysis v5,ERA5）资料以及闪电、大气电场仪、双偏振多普勒天气雷达等观测资料对此次“雷打雪”天气过程进行详细分析,得到以下研究结果：（1）此次“雷打雪”是一次明显的高架对流天气过程,低层一直存在逆温层,850~600 hPa为显著的暖湿平流,925 hPa以下近地面东北风形成的“冷垫”楔入到条件性对称不稳定层结下方,抬升暖湿空气,触发对流天气,释放不稳定能量,产生较强的上升运动。（2）降雪过程的雷电活动比较活跃,闪电频数高达99次·（10 min）^-1,云闪和地闪比例为1.7∶1,正地闪的占比为19.4%。此次“雷打雪”天气过程的闪电空间分布和走向与强降雪落区表现出良好的一致性,闪电密集区域对应强降雪区域,74%的暴雪站周围30 km范围内出现闪电,100%的大暴雪站周围30 km内发生了闪电。（3）高空正K_DP区域与地面降水率的增大密切相关。总闪电与6 km等高层雷达回波的空间一致性很好,基本分布在大于20 dBZ的云区内,闪电落区的回波顶高几乎都在5 km以上。雷达发现在5~8 km的高度范围出现Z_DR和K_DP的正值区,揭示云内存在明显的垂直上升运动和过冷却液态水,存在良好云内起电条件,结合地面电场观测资料和闪电资料,推测此次降雪云系的电荷结构为上正下负的电荷结构。     

渤海海效应暴雪微物理过程的数值模拟

采用RAMS4.4中尺度数值模拟结果,分析了2008年12月4～6日1次渤海海效应暴雪的微物理过程。结果表明:(1)在降雪初期,云中水物质包含云水、雨水、霰、冰晶和雪晶,及地水物质为雨水和霰。随着温度的降低,中后期仅存冰晶和雪晶,产生降雪。由于整个过程以降雪为主,降雨时间短暂,通常忽略降雨,称为降雪过程。(2)本次海效应暴雪的微物理过程表现在两方面:一是"播撒-反馈"机制,二是合适的冰相过程。这两种过程均有利于降雪增幅。(3)西风槽前产生的环境云和冷空气流经渤海暖海面时形成的海效应云之间在合并时发生"播撒-反馈"作用,前者是中云,后者是低云,前者从上层播撒冰晶和雪晶到下层,使得降雪增强。(4)微物理过程另一个有利因素是环境温度,本次强冷空气使得降水云中的温度在-10～-15℃之间,有利于树枝状冰晶的增长,从而产生强降雪。强降雪发生在强上升运动、高相对湿度和适宜的温度的叠置区域。     

山东深秋一次暴雪过程雪水比影响因子分析

利用地面气象观测站资料、加密地面观测资料和欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts,ECMWF）第五代大气再分析数据（ECMWF Reanalysis v5,ERA5;分辨率为0.25°×0.25°）逐小时资料,对山东2021年11月6—8日极端暴雪过程雪水比影响因子进行研究。结果显示：此次暴雪过程平均雪水比分布总体呈“北大南小、西大东小”的分布特征,降雪初期产生的雪水比小,降雪中后期产生的雪水比大;温度偏高、云内液态水含量较高的地区雪水比较小,温度偏低、云内液态水含量较低的地区雪水比较大;雪水比与地面气温、地表温度呈负相关,地面气温与雪水比的相关性最大,积雪产生之后地表温度与雪水比变化无明显相关。     

Evaluation of hydraulic conductivity through particle shape and packing density characteristics of sand–silt mixtures

Abstract

This study aims to evaluate the relationship between saturated hydraulic conductivity with particle shape and packing density characteristics of silty sand soils. The article presents a series of hydraulics tests performed on three kinds of sand with different particles shapes (Chlef rounded sand, Fontainebleau sub-rounded sand and Hostun sub-angular sand) mixed with low plastic rounded Chlef silt in the range of 0–30% fines content. The sand–silt mixture samples were tested in the constant-head permeability device at a loose relative density (D_r = 18%) and a constant room temperature (T?=?20?°C). The obtained results indicate that the measured saturated hydraulic conductivity (K_s) correlates very well with the fines content (F_c), packing density in terms of [maximum void ratio “e_max,” minimum void ratio “e_min,” predicted maximum void ratio “e_max”_pr and predicted minimum void ratio “e_min”_pr] and particle shape characteristics ratios in terms of roundness ratio (R_r = R_hs/R_mixture) and sphericity ratio (S_r = S_hs/S_mixture) of the silty sand materials under consideration. Moreover, the analysis of the available data show a noticeable success in exploring the prediction of the saturated hydraulic conductivity (K_s) based on the particle shape and packing density characteristics (R_r, S_r, e_max, and e_min) of the studied sand–silt mixture samples.     

Experiment Study of the Evolution of Coral Sand Particle Clouds in Water

The motion of particle clouds (i.e., sediment clouds) usually can be found in engineering applications such as wastewater discharge, land reclamation, and marine bed capping. In this paper, a series of laboratory tests are conducted on coral sand to investigate the shape feature of the single particle and the mixing processes of the coral sand particle clouds. The shape of coral sand particle is measured and quantified. The experimental results demonstrate that the shape of coral sand particles tends to be spherical as the particle size decreases, and empirical equations were established to explain the variation of D₅₀ and f_S,50 of coral sand. Compared with the silica sand, the evolution of the coral sand particle cloud still experiences three stages, but the threshold for the Reynolds number of particle clouds entering the next stage changes. Further, the normalized axial distance of the coral sand particle clouds is 58% smaller. The frontal velocity exhibits similar varying tendency for the coral sand particle cloud. Considering the difference in shape between coral sand particles and silica sand particles, a semi-empirical formula was proposed based on the original silica sand prediction formula by adding the shape factor and the experimental data of 122 µm⩽D₅₀⩽842 µm. It can predict the frontal velocity of the coral sand particle clouds.

     

位涡在冷流暴雪短时预报中的应用

利用NCAR/NCEP逐日6 h间隔的再分析资料,对2005 年冬季发生在山东半岛的冷流暴雪过程中的位涡演变特征及其对降雪的作用机制进行了分析.结果表明:(1) 中高纬度阻塞环流形势经历了转型,其发展演变过程有利于冷空气频繁南下从不同的路径入侵渤海和山东半岛地区,造成持续性冷流暴雪;(2) 高位涡从高层下传,降雪强度与下伸高度有关;(3) 高位涡的移动可以很好地示踪冷空气的源地和路径,此次持续性强降雪过程的冷空气来源于4种路径;(4) 对流层中层的高位涡区的强度和影响时间可作为冷流暴雪预报的有益指标,500 hPa等压面的高位涡信号明显,与强降雪的出现有很好的同位相对应关系,700 hPa的位涡场有时候也表现出强劣信号,业务应用中可同时比较分析二者的高位涡区;(5) 高层的高位涡冷空气下沉时产生位势不稳定和对称不稳定层结,对降雪起到增强作用,500 hPa出现明显槽时可形成"不稳定-稳定-不稳定"的大气层结稳定垂直结构,在短时间内产生强降雪.     

Estuarine phytoplankton dynamics and shift of limiting factors: A study in the Changjiang (Yangtze River) Estuary and adjacent area

Environmental factors in estuaries are highly variable in terms of both spatial and temporal dimensions and hence phytoplankton biomass, as well as community structure, is dynamic. Two cruises were carried out in the Changjiang (Yangtze River) Estuary and adjacent area in spring and summer. The result of CHEMTAX calculation suggests that in spring diatoms and chlorophytes contribute equally to phytoplankton biomass, while phytoplankton community structure is mainly composed of diatoms in summer. We encountered blooms in summer with chlorophyll a (CHLa) over 10 μg l⁻¹ off the Changjiang Estuary and they were mainly caused by diatoms (>90%). Based on the HPLC analysis of samples collected, phytoplankton pigments mainly concentrated beyond the front between 122.5°E and 123°E where nutrients and turbidity were best balanced. Euphotic depth (Z_eu, calculated from Secchi disk depth) to surface mixed layer depth (Z_mix) ratio (i.e. Z_eu/Z_mix) were comparable in spring (average value 1.2) and the ratio increased to 5.2 in summer. Variation of the ratio indicates an apparent shift of light and physical conditions from spring to summer. Correspondingly, CHLa was positively related to Z_eu/Z_mix ratio (r² = 0.83) in spring, indicating the light limitation over the whole investigation area. On the other hand, the relationship of CHLa and Z_eu/Z_mix ratio became unclear when Z_eu/Z_mix ratio >3 in summer. This is probably due to the combination of both light limitation before the front and nutrient limitation beyond the front. In addition, evidence was found that light condition can impact the diagnostic pigments in the Changjiang Estuary.     

Wave friction factor rediscovered

The wave friction factor is commonly expressed as a function of the horizontal water particle semi-excursion (A _wb) at the top of the boundary layer. A _wb, in turn, is normally derived from linear wave theory by \fracU_\textwbT_\textw2p \frac{{{U_{\text{wb}}}{T_{\text{w}}}}}{{2\pi }} , where U _wb is the maximum water particle velocity measured at the top of the boundary layer and T _w is the wave period. However, it is shown here that A _wb determined in this way deviates drastically from its real value under both linear and non-linear waves. Three equations for smooth, transitional and rough boundary conditions, respectively, are proposed to solve this problem, all three being a function of U _wb, T _w, and δ, the thickness of the boundary layer. Because these variables can be determined theoretically for any bottom slope and water depth using the deepwater wave conditions, there is no need to physically measure them. Although differing substantially from many modern attempts to define the wave friction factor, the results coincide with equations proposed in the 1960s for either smooth or rough boundary conditions. The findings also confirm that the long-held notion of circular water particle motion down to the bottom in deepwater conditions is erroneous, the motion in fact being circular at the surface and elliptical at depth in both deep and shallow water conditions, with only horizontal motion at the top of the boundary layer. The new equations are incorporated in an updated version (WAVECALC II) of the Excel program published earlier in this journal by Le Roux et al. Geo-Mar Lett 30(5): 549–560, (2010).     

Experimental study on interaction between compressibility and permeability of a crushed sandstone–mudstone particle mixture

The present study focuses on the compressibility and permeability characteristics of a crushed sandstone–mudstone particle mixture (SMPM). Two type of laboratory tests, which are compressibility–permeability test (CPT) and compressibility test only (CTO), are performed. Based on the test data, the effects of the seepage action on the compressibility and ones of the void ratio (e) on the permeability are analyzed. The rate of consolidation of the crushed SMPM is also discussed. The values of compressibility index (C_c) obtained from the CPT are greater about 1.32–4.81% than ones obtained from the CTO, but the values of preconsolidation stress (σ_p) obtained from the CPT are smaller about 2.34–9.83% than ones obtained from the CTO. The slope of fitting line of e～logK (where K is the coefficient of permeability, and log is the logarithm to base 10), defined as the permeability index (K_c), ranges from 0.146 to 0.337 with an average of 0.226. The value of C_c/K_c, used to evaluate the rate of consolidation, ranges from 0.905 to 1.250 with an average of 1.031. The rate of consolidation of the crushed SMPM may be analyzed by Terzaghi’s theory due to the average value of C_c/K_c very close to 1.0.