Thermal and compositional stratification of the inner core

The improvements of the knowledge of the seismic structure of the inner core and the complexities thereby revealed ask for a dynamical origin. Sub-solidus convection was one of the early suggestions to explain the seismic anisotropy, but it requires an unstable density gradient either from thermal or compositional origin, or from both. Temperature and composition profiles in the inner core are computed using a unidimensional model of core evolution including diffusion in the inner core and fractional crystallisation at the inner core boundary (ICB). The thermal conductivity of the core has been recently revised upwardly and, moreover, found to increase with depth. Values of the heat flow across the core mantle boundary (CMB) sufficient to maintain convection in the whole outer core are not sufficient to make the temperature in the inner core super-isentropic and therefore prone to thermal instability. An unreasonably high CMB heat flow is necessary to this end. The compositional stratification results from a competition of the increase of the concentration of light elements in the outer core with inner core growth, which makes the inner core concentration also increase, and of the decrease of the liquidus, which makes the partition coefficient decrease as well as the concentration of light elements in the solid. While the latter (destabilizing) effect dominates at small inner core sizes, the former takes over for a large inner core. The turnover point is encountered for an inner core about half its current size in the case of S, but much larger for the case of O. The combined thermal and compositional buoyancy is stabilizing and solid-state convection in the inner core appears unlikely, unless an early double-diffusive instability can set in.     

基于包络检波和STFT谱分析的探地雷达土壤分层信息识别

土壤分层信息,特别是表土层结构,对土地生产力具有重要影响,是评价土壤质量的一个重要指标。为了快速、准确地获取土壤分层信息,本文利用探地雷达对分层土壤进行了回波信号采集,并分别在时域和频域分析土壤层位置和层厚信息。首先在信号预处理的基础上,借助包络检波方法确定在土壤分层界面在时域上的位置;然后获取电磁波速度,得到土壤分层厚度。考虑到土壤介电常数与电磁波在土壤中传播速度的相关性,采用短时傅里叶变换方法（Short-time Fourier Transform,STFT）获取各土壤层时频域特征值,并利用回归分析建立特征值与介电常数之间的数学关系,实现对各土壤层的介电常数估算,从而计算出电磁波传播速度,进而确定土壤各层厚度。为验证算法的有效性,分别对理想模拟实验环境和农田环境进行了探地雷达实验,结果表明利用包络检波对探地雷达回波信息进行分析,土壤层检出率达到94.5%,借助STFT谱分析进行探地雷达回波速度估计,对于70 cm深度以上土层厚度计算误差大都保持在10%以下,但随着土壤深度的增加,误差变大。总体来说,本方法能有效识别浅层土壤的分层信息,可应用于实际生产中耕层厚度的估测。     

ANALYSIS OF THE STRATIFICATION FEATURES OF THE YELLOW SEA AND THE LOCATION OF ITS GENERAL SHELF FRONT IN THE WARM HALF-YEAR

A stratification parameter ,defined as theamount of mechanical energy required to bring about vertical mixing, has been calculated for theYellow Sea using available data over the past ten years.T he monthly distributions of Log are obtained to explain the features of the Yellow Sea stratification.Fronts of the shallow shelf sea are often inseparably related with its stratifications. The front of the Yellow Sea in the warm half-year is generated in May and disappears in November. The shelf front moves shoreward and becomes strong in the heating season, but becomes weak in the cooling season upon return.     

北京周口店雾迷山组碳酸盐岩风暴沉积序列研究*

北京周口店地区中元古界雾迷山组以含硅质条带的白云岩为主,发育不规则硅质砾屑。这些硅质砾屑形态各异,具有强烈扰动的迹象,目前对其成因解释一直存在争议,如泄水构造、冲刷构造、震积岩、风暴岩等。文中选择北京周口店地区黄山店村恒顺厂剖面为研究对象,通过野外精细的沉积学分析,系统描述了岩性及沉积构造特征,确定了该套沉积地层发育典型的碳酸盐岩风暴沉积序列。该风暴沉积序列自下往上可分为5段: A段,为风暴前正常沉积的中厚层状泥晶白云岩;B段,具口袋构造的侵蚀冲刷面及硅质砾屑,为风暴高潮期产物;C段,灰色硅质条带白云岩,普遍发育平行层理和波状、丘状、洼状交错层理,是风暴衰减期产物;D段,灰白色薄层状泥晶白云岩,代表了风暴间歇期的正常沉积;E段,含硅质团块白云岩,为正常天气下海水局部扰动形成。该剖面上递变层理几乎不发育,结合侵蚀冲刷面特征、砾屑高含量以及杂乱堆积方式,认为该套沉积为典型的近原地风暴沉积。结合古地理资料,推测研究区在雾迷山组沉积时期可能处于热带海洋环境,热带气旋引发的频繁风暴潮导致了近原地碳酸盐岩风暴沉积的形成。上述研究成果不仅补充了燕山地区雾迷山组的风暴沉积记录,而且为该时期华北板块的古地理、古纬度和古气候研究提供了重要的沉积学证据。     

缓倾角塔柱状危岩压裂损伤-突变失稳预测

陡高边坡灰岩地区缓倾角塔柱状危岩常具有底部压裂诱发整体失稳破坏模式,其损伤-突变失稳机制属于山地灾害关键科学问题之一。以重庆南川区甑子岩W12#危岩崩塌为例,建立了荷载、水致弱化效应耦合下的损伤突变地质力学模型,获得了基于应变等价原理的损伤本构方程及总损伤度演化方程,并将水致弱化函数改进为软化系数的一元三次函数。将地质力学模型概化为等效弹簧模型,通过能量平衡原理建立损伤-折叠突变模型,获得了塔柱状危岩系统压裂失稳判据及临界突变位移特征值表达式。计算表明：W12#危岩压裂失稳时,该系统的控制变量为–0.003 251,小于0,表示该系统进入不稳定状态;计算理论突变位移起始值148.70 mm比实测位移第一拐点值154.34 mm偏小,相对误差约3.65%,偏于安全;理论损伤本构曲线、理论损伤演化曲线与现有文献数值模拟结果趋于一致,所构建的理论预测模型具有较好的适用性。其研究成果可用于预判塔柱状危岩压裂失稳损伤演化过程及突发失稳特征位移值,为灰岩区陡高边坡危岩崩塌监测预警和防灾减灾提供理论依据。     

气候变暖对湖泊热力及溶解氧分层影响研究进展

气候变暖对湖泊物理、化学、生物和生态系统有着复杂而深刻的直接和间接影响,而具体影响随研究区域和水体表现不尽相同。气候变暖通过改变湖泊热力和溶解氧分层进而影响湖泊生物过程和生态系统结构与功能。从全球湖泊变暖趋势、长期缓慢气温上升、极端高温事件以及气候情景模拟等方面详细综述了气候变暖对湖泊热力及溶解氧分层影响的研究进展。研究表明,全球不同区域湖泊均存在不同程度的变暖趋势;长期缓慢气温上升和短期极端高温均会造成湖泊热力分层提前,分层结束推迟,分层时间延长,混合层和温跃层深度下降,以及热稳定性增加;相伴随的是溶解氧扩散深度和氧跃层深度明显下降,加剧了湖泊底部好氧和厌氧环境。除了这种直接影响外,气候变暖引起的流域降水、入湖物质的变化以及风速的变化也会对湖泊热力和溶解氧分层产生许多间接的影响,因此未来仍然需要更多的实验证据、经验和数值模型来验证和预测气候变暖对湖泊热力及溶解氧分层的影响。     

两次秸秆焚烧导致辽宁重污染过程的气象条件对比分析

利用环境监测站大气污染物数据、地面自动气象站观测资料、L波段加密探空资料和0.125°×0.125°的EC再分析资料,结合MODIS遥感火点监测和HYSPLIT4后向轨迹模拟结果,对比分析了2015年11月8日和2016年11月5日的两次由于东北地区秸秆焚烧导致辽宁重污染天气过程的大气边界层特征、气象扩散条件和大气污染物输送来源等。结果表明:两次过程地面PM_(2.5)浓度均出现快速上升和下降,其中2015年11月8日重污染过程的污染强度较2016年11月5日强,且持续时间更长。2015年11月8日重污染过程的混合层高度较低,其上层的中性层结转变为逆温层结,抑制混合层高度的发展。同时低层冷平流不断侵入到暖平流下方,使得大气层结稳定性增强,维持时间较2016年11月5日重污染过程更长,低层下沉运动和黑龙江西南部、吉林西部污染物的远距离输送增强使得辽宁地面污染物浓度快速累积。而2016年11月5日重污染天气过程主要受深厚冷空气影响,东北地区西部污染物的区域输送和地面风场辐合是地面污染物浓度快速上升的主要原因。     

杭州市典型雨转雪天气成因及预报模型

利用2008—2018年的NCEP(1°×1°)再分析资料、常规气象观测资料和降雪加密观测资料,选出杭州地区10次典型的雨转雪天气过程,从大尺度环流背景和动力、水汽以及热力因子等物理量场结构方面展开研究,最终得出杭州冬季典型雨转雪天气的预报模型:①大尺度环流配置需满足能为雨转雪天气的形成提供有利的水汽、动力抬升以及中低层上暖下冷的逆温或等温层结条件;②水汽和动力因子等物理量须满足产生纯雪的特定条件;③杭州温度层结须为T_(2m)≤1.5℃、T_(925)≤-4.0℃、T_(850)≤0℃、T_(700)≤-1.0℃和T_(500)≤-10.0℃。此外,进一步补充了杭州可能产生大雪甚至暴雪量级降雪的特定条件。最终选取2019年初的2次典型降水过程进行预报回报检验。     

“尤特”强降水过程大气层结分析

1311号强台风“尤特”登陆后给广东带来持续性大范围强降水,对流降水特征显著.分析了“尤特”影响期间大尺度环流背景,重点讨论了此次持续性强降水过程中大气层结问题.发现低空急流向广东输送强的暖平流,是广东大气层结不稳定得以持续维持的根本原因.进一步分析发现,低空急流本身并不是“暖”的,当“尤特”趋向陆地时,陆地上的暖气团在“尤特”环流强迫下向南传播扩散,低空急流穿越这一暖区时温度升高才具备“暖”的特性.这一事实在以前并未被关注到.通过个例反查,在许多登陆后造成连续强降水的台风过程中均发现了这一特征.因此,台风登陆引起环境温度场的演变以及与低空急流的配置需引起重视.     

The general solution of the HENON–HEILES problem

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.