高原涡和西南涡影响的两次四川暴雨过程的对比分析

为了进一步研究高原涡、西南涡对西南地区暴雨的影响,本文用中国气象局自动站与CMORPH降水数据融合的逐时降水资料、国家卫星气象中心的逐时FY-2E卫星的云顶亮温（TBB）资料、欧洲气象资料中心（ERA-interim）的再分析资料,通过天气学诊断分析方法以及拉格朗日轨迹模式HYSPLITv4.9,对发生在四川盆地的有高原涡东移影响西南涡发展引发暴雨的两次过程进行对比分析,发现:（1）两次暴雨过程的降水强度和分布有明显区别,并且TBB活动特征显示在过程一中有MCC（Mesoscale Convective Complex）的产生和发展,过程二则没有。（2）对于过程一,500 hPa上,高原涡逐渐减弱为高原槽并伸展到四川盆地上空,850 hPa上,在鞍型场附近有MCC的产生和发展,200 hPa上,高原涡在南亚高压北部偏西风急流下方的强辐散区内,位于南亚高压东南侧急流区下方稳定少动,偏东风急流北部有辐散中心,有利于西南涡的加强。对于过程二,500 hPa高原涡东移在四川盆地上空与西南涡耦合,形成一个稳定且深厚的系统,这也是过程二的暴雨强度比过程一强的最主要原因。200 hPa上,四川盆地始终位于南亚高压东侧的西北气流中,“抽吸作用”明显。（3）在过程一中,位涡逐渐东传且位涡增加的地方对应强降水区与MCC发展区,反映了暴雨和位涡的发展基本一致。在过程二中,中层位涡高值区从高原上东移并下传至盆地上空,两涡耦合使得上下层打通,位涡值比耦合之前单独的两涡强度更强。 MCC产生的必要条件是中层大气要有强正涡度、强辐合和强上升运动,在未产生MCC前,过程一与过程二在盆地上空的动力条件甚至是相反的;从热力条件看,过程一中有明显的干冷空气入侵,增强不稳定条件,有利于MCC的产生并引发强降水;另一方面,本文也应证了二阶位涡的水平分布与暴雨落区有较好的对应关系。（4）通过拉格朗日方法的水汽轨迹追踪模式和聚类分析方法分析可得两次暴雨过程的水汽输送源地和通道也有明显区别,过程一主要有两条水汽通道,通道一来自阿拉伯海和孟加拉湾洋面的底层,通道二来自四川南部750 m以下高度;而过程二的主要水汽输送通道有三条,通道一来自西方地中海、黑海和里海上空1500～2500 m高度附近,通道二来自阿拉伯海和印度洋的底层,通道三的水汽从孟加拉湾低层绕过云贵高原直接输送到四川盆地。     

延庆-张家口地区复杂地形冬季山谷风特征分析

基于2016年12月—2017年2月和2017年12月—2018年2月两年冬季的近地面自动气象站逐时观测数据以及张家口探空数据分析延庆-张家口一带（包括张家口崇礼、赤城、海坨、小五台山区,延怀、怀涿、洋河、蔚县盆地以及北京延庆、昌平、怀柔部分平原地区）复杂地形的风场精细化时、空分布特征,揭示不同复杂地形下局地风场的时、空变化规律,加深对复杂地形动力、热力作用对近地面风场影响的认识,为冬季山区风场预报以及复杂地形数值模式改进提供参考。结果表明:晴朗小风天风持续性作为矢量平均风速和标量平均风速的比值,可以作为研究风场变化规律的重要参数。根据风持续性的日变化特征,可以将研究区域内所有站点分为10种类型,分别代表不同局地地形特征的影响,风持续与风向变化的相关也很强。研究区域主要有3种类型的地形风:斜坡风、峡谷风以及较大尺度的山区平原风。不同地形特征下的风场、风持续性存在明显不同的日变化特征,山风和谷风相互转化的时间也不同,山区最早,盆地次之,平原区最晚;山风时段持续时间较谷风时段长,风速小;晴朗小风天实测风反映了实际风场的特征,而排除环境背景风场,弱化地形动力作用后整个冬季的局地风作为理论山谷风,更能反映热力作用下的山谷风特征。      

北京山区与平原冬季近地面风的精细观测特征

利用2009—2018年冬季北京地区200多个自动气象站逐时10 m风速、风向观测数据,分典型区域（山区、山区与平原过渡区、平原区、城区）研究北京地区冬季近地面风的精细特征,并使用有完整记录的2 a（2017和2018年）冬季延庆高山区不同海拔高度10 m风逐时观测数据,多视角分析高山区不同海拔高度近地面风的特征和成因,以深刻认识北京地区复杂地形条件下冬季近地面风的特征和规律。结果表明:（1）北京地区冬季近地面平均风受西部北部地形、城市下垫面粗糙度和冷空气活动共同影响,平均风速沿地形梯度分布,山区高平原低,平原中又以城区风速最小;盛行西北风和北风,在城区东、西两侧盛行风出现扰流,在山区和过渡区一些地方还存在与局地地形环境明显关联的其他盛行风向。（2）4个典型区域冬季近地面风速日变化均表现为白天风速大于夜间,午间风速最大的“峰强谷平”单峰特征,这一特征的稳定性在城区高、山区低。（3）4个区域冬季弱风（< 1 m/s）频率为31%—42%,城区较高、山区较低;强风（> 10.8 m/s）频次则是山区多、城区少,强风风向主要表现为偏西—偏北,与冷空气活动密切关联;城区、平原区和过渡区偏南风频率均为极小,暗示北京“山区—平原”风模态在冬季是“隐式”的、不易被直接观测到。（4）近地面风的水平尺度代表范围在延庆高山区高海拔处明显大于低海拔处,海拔1500 m附近（平均的边界层顶高度）是延庆高山近地面风速日变化特征的“分水岭”,低于该海拔高度时近地面风速日变化表现为前述“峰强谷平”单峰特征,而高于该海拔高度时近地面风速日变化则呈现相反特征,即夜间大白天小、午间最小的“峰平谷深”特征,这是由边界层湍流活动的日变化及伴随的低层自由大气动量向边界层内下传所致。（5）延庆高山近地面风速大体上随观测高度而增大,高海拔站点日平均风速数倍于低海拔站点。白天—前半夜,海拔约2000 m的站点冬季盛行偏西风,风向变化不大,但风速为2—12 m/s;1000 m左右的低海拔站则风速比较稳定（< 6 m/s）,风向从午间至傍晚相对多变。      

我国中东部平原地区临界气温条件下降水相态判别分析

基于2001—2013年地面观测和探空资料,对地面气温位于0～2℃(以下称临界气温)我国降雪的时空分布及其与降雨的垂直热力特征进行了研究,引入了决策树判别方法对上述条件下雪和雨进行了判别分析,结果表明:临界气温下降雪出现频率总体高于降雨、雨夹雪出现频率,且在我国华北南部至江南北部的中东部地区分布较多,年均可达7.69～15.38站次;临界气温下,降水相态为雨或雪对应的平均温度廓线最大差异位于650 hPa附近,且地面气温较低时,平均温度差异更明显,平均湿度廓线差异则主要位于低层,且在地面气温较高时,平均湿度差异更明显;临界气温下,降水相态为雨时,地面上空存在暖层样本占比,较降水相态为雪时更高,且降雨时暖层主要位于中层,降雪时暖层则主要位于低层,降雨时其暖层强度显著大于降雪时暖层强度;在临界气温下雨雪判别分析中,地面气温能显著提升判别准确率,湿球温度能在一定程度上提升判别准确率,基于云顶温度、中层融化参数、低层湿球温度构建的决策树判别模型,判别准确率达到91.86%,能较好地解决临界气温下雨和雪的判别问题。     

气象因子对半冬性小麦灌浆速度的影响效应研究

选用黄淮海冬麦区4个半冬性小麦品种郯麦98、山农18、徐麦33、皖麦52为试验材料,通过分期播种试验,利用方差分析、相关分析、逐步回归和通径分析等方法,分析半冬性小麦籽粒灌浆速度变化趋势和气象因子对灌浆速度的影响。结果表明,正常播期冬小麦灌浆速度波动性最小、千粒重最大,迟播10 d冬小麦灌浆速度波动性最大、千粒重最小;华北区品种郯麦98灌浆速度表现最稳定、千粒重最高,而黄淮区品种皖麦52灌浆速度最大;半冬性小麦灌浆持续期为35～39 d;南北气候差异是影响各品种冬小麦灌浆速度不同的原因之一。半冬性小麦各播期灌浆速度的变化趋势一致,灌浆速度变化与相关显著气象因子的变化规律相符合;灌浆速度峰值期一般出现在开花后15～25 d,迟播冬小麦最大灌浆速度出现时间较对照处理提前,不利于提高粒重;气温条件对冬小麦灌浆速度影响显著,其中最高气温要素是影响不同播期品种灌浆速度的共有关键因子。通径分析表明,最高气温对灌浆速度的作用由自身的直接效应决定,而日照时数与最低气温对灌浆速度的作用与间接效应一致;最高气温平均值对灌浆速度的影响最重要,日照时数和最低气温平均值对灌浆速度的影响较弱;最高和最低气温平均值、日照时数均为灌浆速度的限制因子,其中最高气温平均值对灌浆速度变化的决策作用最大。     

Explicit integration of a porosity‐dependent hydro‐mechanical model for unsaturated soils

An adaptive substepping explicit integration scheme is developed for a porosity‐dependent hydro‐mechanical model for unsaturated soils. The model is referred to as the modified σ –Θ model in this paper, which features the employment of the subloading surface plasticity and the stress–saturation approach. On numerical aspects, convex/nonconvex subloading surfaces in the σ –Θ space may result in incorrect loading–unloading decisions during the integration. A new loading–unloading decision method is developed here to solve the problem and then embedded into the explicit integration scheme for the modified σ –Θ model. In addition, to enhance the accuracy of the explicit integration, local errors from both hydraulic and mechanical components are included in the error control for each substep. A drift correction method is also developed to ensure the state point lies on the subloading surface in the σ –Θ space within a set error level. The performance of the loading–unloading decision method for the modified σ –Θ model is discussed through comparing it with the conventional loading–unloading decision method. The importance of involving the hydraulic component in the error control is also demonstrated. The accuracy and efficiency of the proposed adaptive substepping explicit integration scheme for the modified p–Θ model are also studied via several numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Structural and metamorphic evolution of the Moornambool Metamorphic Complex,western Lachlan Fold Belt,southeastern Australia

The wedge‐shaped Moornambool Metamorphic Complex is bounded by the Coongee Fault to the east and the Moyston Fault to the west. This complex was juxtaposed between stable Delamerian crust to the west and the eastward migrating deformation that occurred in the western Lachlan Fold Belt during the Ordovician and Silurian. The complex comprises Cambrian turbidites and mafic volcanics and is subdivided into a lower greenschist eastern zone and a higher grade amphibolite facies western zone, with sub‐greenschist rocks occurring on either side of the complex. The boundary between the two zones is defined by steeply dipping L‐S tectonites of the Mt Ararat ductile high‐strain zone. Deformation reflects marked structural thickening that produced garnet‐bearing amphibolites followed by exhumation via ductile shearing and brittle faulting. Pressure‐temperature estimates on garnet‐bearing amphibolites in the western zone suggest metamorphic pressures of ～0.7–0.8 GPa and temperatures of ～540–590°C. Metamorphic grade variations suggest that between 15 and 20 km of vertical offset occurs across the east‐dipping Moyston Fault. Bounding fault structures show evidence for early ductile deformation followed by later brittle deformation/reactivation. Ductile deformation within the complex is initially marked by early bedding‐parallel cleavages. Later deformation produced tight to isoclinal D₂ folds and steeply dipping ductile high‐strain zones. The S₂ foliation is the dominant fabric in the complex and is shallowly west‐dipping to flat‐lying in the western zone and steeply west‐dipping in the eastern zone. Peak metamorphism is pre‐ to syn‐D₂. Later ductile deformation reoriented the S₂ foliation, produced S₃ crenulation cleavages across both zones and localised S₄ fabrics. The transition to brittle deformation is defined by the development of east‐ and west‐dipping reverse faults that produce a neutral vergence and not the predominant east‐vergent transport observed throughout the rest of the western Lachlan Fold Belt. Later north‐dipping thrusts overprint these fault structures. The majority of fault transport along ductile and brittle structures occurred prior to the intrusion of the Early Devonian Ararat Granodiorite. Late west‐ and east‐dipping faults represent the final stages of major brittle deformation: these are post plutonism.     

Vanadium in magnetite gabbros and its behaviour during lateritic weathering,Windimurra Complex,Western Australia

Shephards Discordant Zone is a 500–600 m thick interlayered sequence of deformed, altered and metamorphosed magnetite metagabbro and about 50 layers or lenses of magnetitite (> 80–90% magnetite). The sequence shows progressive magmatic fractionation upwards: Ti and Ti/Fe increase, and V, V/Ti and Cr decrease upwards in magnetite and in whole‐rock compositions. The main magnetite‐rich sequence (about 400 m thick) is deeply weathered, with 40 m of saprolite showing vertical zonation of weathering minerals due to progressive weathering. Magnetitites (average 1% V₂O₃) are resistant to weathering and show little chemical change, but magnetite gabbros (average 0.27% V₂O₃) are extensively weathered and show progressive loss of Ca, Na, Mg and S. Plagioclase, magnetite (1.37% V₂O₃), chlorite (up to 0.35% V₂O₃), actinolite, epidote and minor sulfides in unweathered rocks weather to kaolinite, hematite, goethite and minor vermiculite, ilmenite remaining largely unaffected. Vanadium is essentially immobile during weathering and is unaffected during weathering of magnetitites (1% V₂O₃), but is slightly depleted during weathering of magnetite gabbros (0.23% V₂O₃).     

Petrogenetic modelling of strongly residual metapelitic xenoliths within the southern Platreef,Bushveld Complex,South Africa

Xenoliths of quartz‐absent Fe‐rich aluminous metapelite are common within the platinum group element‐rich mafic/ultramafic magmatic rocks of the Platreef. Relative to well‐characterized protoliths, the xenoliths are strongly depleted in K₂O and H₂O, and have lost a substantial amount of melt (>50 vol.%). Mineral equilibria calculations in the NCKFMASHTO system yield results that are consistent with observations in natural samples. Lower‐grade rocks that lack staurolite constrain peak pressures to ～2.5 kbar in the southern Platreef. Smaller xenoliths and the margins of larger xenoliths comprise micro‐diatexite rich in coarse acicular corundum and spinel, which record evidence for the metastable persistence of lower‐grade hydrous phases and rapid melting consequent on a temperature overstep of several hundred degrees following their incorporation in the mafic/ultramafic magmas. In the cores of larger xenoliths, temperatures increased more slowly enabling progressive metamorphism by continuous prograde equilibration and the loss of H₂O by subsolidus dehydration; the H₂O migrated to xenolith margins where it may have promoted increased melting. According to variations in the original compositional layering, layers became aluminosilicate‐ and/or cordierite‐rich, commonly with spinel but only rarely with corundum. The differing mineralogical and microstructural evolution of the xenoliths depends on heating rates (governed by their size and, therefore, proximity to the Platreef magmas) and the pre‐intrusive metamorphic grade of the protoliths. The presence or absence of certain phases, particularly corundum, is strongly influenced by the degree of metastable retention of lower‐grade hydrates in otherwise identical protolith bulk compositions. The preservation of fine‐scale compositional layering that is inferred to be relict bedding in xenolith cores implies that melt loss by compaction was extremely efficient.