甘孜-玉树断裂带晚第四纪走滑速率与滑动分解作用

滑动速率是研究断裂运动学特征、地震活动性和区域应变分配的重要参数和依据。前人关于甘孜-玉树断裂带滑动速率的研究结果存在较大差异,因此,其晚第四纪滑动速率有待进一步调查研究。本文基于卫星影像解译和野外实地考察,对甘孜-玉树断裂带西段（玉树断裂）上典型断错地貌点进行测量分析,得到玉树断裂晚第四纪走滑速率为6.6±0.1-7.4±1.2mm/a。通过与前人对甘孜-玉树断裂带东段（甘孜断裂）滑动速率的研究结果进行对比,发现甘孜-玉树断裂带东、西段滑动速率不一致,其原因是甘孜断裂的左旋滑移在向西传递的过程中,一部分应变被分配到了巴塘盆地南缘断裂上。巴塘盆地南缘断裂的存在很好地解释了玉树断裂的走滑速率比甘孜断裂偏低的原因。但是,从区域变形来看,巴塘盆地南缘断裂分配的滑动速率恰好说明了甘孜-玉树断裂带东、西段及鲜水河断裂带的水平构造变形是协调一致的。     

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

Effect of compressive strength on the performance of the NEMO-LIM model in Arctic Sea ice simulation

Satellite records show that the extent and thickness of sea ice in the Arctic Ocean have significantly decreased since the early 1970s. The prediction of sea ice is highly important, but accurate simulation of sea ice variations remains highly challenging. For improving model performance, sensitivity experiments were conducted using the coupled ocean and sea ice model (NEMO-LIM), and the simulation results were compared against satellite observations. Moreover, the contribution ratios of dynamic and thermodynamic processes to sea ice variations were analyzed. The results show that the performance of the model in reconstructing the spatial distribution of Arctic sea ice is highly sensitive to ice strength decay constant (C^rhg). By reducing the C^rhg constant, the sea ice compressive strength increases, leading to improved simulated sea ice states. The contribution of thermodynamic processes to sea ice melting was reduced due to less deformation and fracture of sea ice with increased compressive strength. Meanwhile, dynamic processes constrained more sea ice to the central Arctic Ocean and contributed to the increases in ice concentration, reducing the simulation bias in the central Arctic Ocean in summer. The root mean square error (RMSE) between modeled and the CryoSat-2/SMOS satellite observed ice thickness was reduced in the compressive strength-enhanced model solution. The ice thickness, especially of multiyear thick ice, was also reduced and matched with the satellite observation better in the freezing season. These provide an essential foundation on exploring the response of the marine ecosystem and biogeochemical cycling to sea ice changes.

     

Numerical Analysis of Groundwater Ridging Processes Considering Water‐Air Flow in a Hillslope

In this study, a water‐air two‐phase flow model was employed to investigate the formation, extension, and dissipation of groundwater ridging induced by recharge events in a hypothetical hillslope‐riparian zone, considering interactions between the liquid and gas phases in soil voids. The simulation results show that, after a rain begins, the groundwater table near the stream is elevated instantaneously and significantly, thereby generating a pressure gradient driving water toward both the stream (the discharge of groundwater to the stream) and upslope (the extension of groundwater ridging into upslope). Meanwhile, the airflow upslope triggered by the advancing wetting front moves downward gradually. Therefore, the extension of groundwater ridging into upslope and the downward airflow interact within a certain region. After the rain stops, groundwater ridging near the stream declines quickly while the airflow in the lower part of upslope is still moving into the hillslope. Thus, the airflow upslope mitigates the dissipation of groundwater ridging. Additionally, the development of groundwater ridging under different conditions, including rain intensity, intrinsic permeability, capillary fringe height, and initial groundwater table, was analyzed. Changes in intrinsic permeability affect the magnitude of groundwater ridging near the stream, as well as the downward speed of airflow, thereby generating highly complex responses. The capillary fringe is not a controlling factor but an influence factor on the formation of groundwater ridging, which is mainly related to the antecedent moisture. It was demonstrated that groundwater ridging also occurs where an unsaturated zone occurs above the capillary fringe with a subsurface lateral flow.     

强震临震微波动现象初步研究(二)

为了验证和检验强震临震微波动现象及其主要特征,利用甘、青、川、滇、藏等区域台网200余个宽频带数字地震台站的波形资料,建立了实时监控技术系统。通过对2012~2014年间青藏块体发生的24次5级以上地震及部分无震区域台站观测的全程实时跟踪、动态监控,检验和验证了临震微波动现象的重现性和客观性,验证了临震微波动现象的频谱、时间、空间、方向性等特征,进一步得出了临震微波动现象可能与地震强度有关,且震级越大,震前出现持续临震微波动的可能性越大的结论。     

唐山马家沟矿井水位异常分析

地下水动态受水文因素影响较大,对地震和构造活动具有较灵敏的响应。判别并排除各种水文干扰,确认地下水在地震前的异常变化,对提高地震分析预报能力,具有重要作用。马家沟矿井水位动态观测层与地下水开采层为同一含水层,井水位于2010年出现破年变异常,加速持续上升,截至2015年,最大上升幅度约30 m。依据该井水文地质环境特征,根据唐山市区2001-2015年地下水位、降雨量、地下水开采量实测资料,建立合理多元回归模型和三维地下水流动模型,发现地下水开采量减少应为影响马家沟矿井水位动态的可能因素。文中采取的异常识别与分析方法,可为其他类似井孔的地下水动态异常识别及判定提供一定借鉴。     

镇江地震台网理论监测能力评估

基于近震震级公式,对镇江地震台网理论监测能力进行评估,绘制地震监测能力图,以确定监测区不同地点发生地震时能被台网有效分析和定位的地震最小震级。镇江市句容县在该台网理论监测震级最小,为ML1.2—1.3;在镇江及相邻城市,发生M_L 1.8分析认为,以上地震均能被有效监测。     

以颗粒物理原理认识汶川MW 7.9地震发生前临震预滑和震颤现象

发生在孕震区周边地块上的临震预滑和震颤现象,对破坏性地震预测有一定前兆意义,是值得地震学界关注的问题。选取2008年5月12日汶川M_W 7.9地震发生前,临夏和湟源地震台分量应变仪记录与临夏、恩施和西安地震台数字地震仪记录以及临夏和周至地震台深井水位仪记录,分析发现,在临震前数天至数小时,上述各地震台不同学科观测仪器均记录到一些"跃变"和"震颤"震相。文中试图以颗粒物理原理,来认识不同距离、不同台站、不同学科的观测仪器在临震前相近时间段内记录的低频和高频震相,可能是不同地块在临震前发生预滑错动后激发的预滑震相Xp和地下气体在裂隙内流动激发的震颤震相Tp。观测结果表明：2008年5月8日03时至主震发生,各地震台所处地块在相近时段内逐次发生次数不等的预滑错动,其中1-2次较大错动可在噪声背景中被识别;各地震台预滑错动方向指向或背向主震震中。据此认为：汶川M_W 7.9地震前,上述各地震台所处地块在不同大小、不同方向的力链驱动下,发生指向或背向主震震中的临震预滑现象。     

大尺度雁列断层粘滑运动数值模拟

通过ANSYS软件建立二维雁列断层模型,采用粘弹性力学参数,建立摩擦系数接触单元,模拟10万年时间大尺度雁列走滑断层活动（地震相关活动）,分析挤压条件下的断层活动特征,研究结果表明：①断层时空区域不同,所受应力也不同,接近贯通区域重合地区所受应力较大,粘滑运动剧烈,由贯通区域向断层两侧逐渐延伸,所受应力依次减小;②摩擦系数越小,地震发生频率越高,地震运动周期越短,则位错量相对较小,易发生震级较小地震;摩擦系数越大,地震发生频率越低,地震运动周期越长,则位错量相对较大,易发生震级较大地震;③给定的边界挤压速率越大,地震发生频率越高,地震运动周期越短,则位错量增大,易发生震级较大地震。     

Source-independent wave-equation based microseismic source location using traveltime inversion

A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.