利用经验格林函数谱比法估计2013年芦山地震序列的应力降

本文以芦山地震余震为例,分析了基于强震动观测记录的经验格林函数（EGF）谱比法估计地震的拐角频率和应力降的可行性。首先给出能够可靠估计地震拐角频率和应力降的经验格林函数谱比曲线的质量标准;然后利用其估算出17次M_L3.8—5.4芦山强余震的拐角频率;最后,参考其它研究给出的地震矩震级,计算出地震应力降。结果显示,芦山强余震的拐角频率主要介于1.0—2.0 Hz之间,应力降平均值为9.98 MPa,且地震应力降表现出明显的震级相关性。      

姚安地震序列与大姚地震序列震源参数的对比研究

2000年1月15日在云南省姚安县发生了M_S6.5地震, 2003年7月21日和10月16日在云南省大姚县发生了M_S6.2、 M_S6.1地震。 这两个地震序列均发生在滇中块体, 震中位置相距42.5 km, 构造、 应力背景相似, 但序列类型不同, 时间间隔短。 为深入研究姚安地震和大姚地震序列的特征, 根据昆明区域数字台网记录的数字地震波资料, 利用遗传算法计算了2000年姚安地震和2003年大姚地震序列的震源参数。 研究结果表明, 姚安地震序列和大姚地震序列的地震矩和近震震级有很好的线性关系, 它们的地震矩M₀均在10¹²~10¹⁴ N·m, 地震矩与震源半径之间呈线性关系。 应力降随震级变化不明显。 拐角频率和地震矩之间有明显的依赖关系, 用最小二乘法拟合拐角频率和地震矩的关系式, 利用此关系式可计算出给定地震矩的拐角频率估计值f_a, 进而可分析实际计算的拐角频率与估计值之差(f_c-f_a)的时间变化曲线。 姚安地震序列和大姚地震序列实际计算的拐角频率与估计值之差(f_c-f_a)及应力降随时间的变化表现出了不同的特征: 姚安地震序列的应力降和拐角频率差值在4.3级和4.6级强余震前均有一个上升-下降过程, 反映余震区应力场有一个增高-下降过程; 大姚地震序列应力降和拐角频率差值在M_S≥4.6强余震发生之前均会出现高值异常, 但在M_S4.6～4.7地震之前出现的高值异常持续时间短, 在M_S6.1地震前出现的高值异常持续较长。 文中结合两个地震的发震构造, 序列类型等特征进行了分析讨论。     

2017年8月8日九寨沟MS7.0地震震源谱及震中区域品质因子

2017年8月8日四川省九寨沟县发生M_S7.0地震,国家强震动台网中心正式发布了此次地震中收集的66组强震动记录,采用其中震源距不超过150 km且台站场地反应已确定的8组记录,基于S波傅氏谱反演方法确定了此次地震的震源谱及震中区域的地壳介质品质因子.结果为：品质因子Q（f）=84.9f^0.71,说明该区域的S波非弹性衰减较为强烈;根据ω^-2震源谱模型,确定该地震的地震矩为9.42×10¹⁸Nm,其对数标准差为0.12,相应的矩震级M_W=6.616±0.079,拐角频率为0.131±0.011 Hz,应力降为3.854 MPa,该值明显低于全球板内地震的应力降平均水平（4.89 MPa）,但高于2013年芦山M_W6.6地震的应力降.     

2013年岷漳地震和2014年景谷地震的能矩比对比研究

2013年岷漳地震和2014年景谷地震两个地震的矩震级均为M_W6.1,但相同等级烈度区域的面积差异显著,前者均明显大于后者.地面震动的强烈程度与震源破裂过程中辐射的地震波高频成分密切相关,作者系统地测定与研究了这两个地震的震源谱和反映其震源特征的能矩比等震源参量.结果表明,虽然岷漳地震和景谷地震释放的地震矩基本相同,分别为1.6×10¹⁸Nm和1.8×10¹⁸Nm,而且前者还比后者略小些,但其释放的地震波辐射能却差异显著,分别为1.3×10¹⁴J和0.6×10¹⁴J;因而前者的能矩比（8.1×10^-5）是后者能矩比（3.3×10^-5）的大约2.5倍.相应地,岷漳地震的视应力和应力降也都明显大于景谷地震的视应力和应力降.此外,这两个地震震源谱结果均显示出其高频趋势均按角频率ω的-2幂次衰减,符合"ω²模式";采用理论震源谱公式计算结果分别约占实际结果的54%和50%,两个地震的地震波辐射能解析解之间的差异依然约2.3倍.震源破裂特征的显著差异导致岷漳地震单位地震矩释放的地震波辐射能更多,即岷漳地震能矩比更大,是导致这两个矩震级均为M_W6.1地震的地震烈度差异显著的一个重要原因.     

2014年于田M_S7.3地震序列的频谱特征分析及其前震识别

利用新疆区域数字地震台网的波形资料,对2014年2月12日于田MS7.3地震前震和余震序列的频谱特征进行了分析。基于Brune模型,对观测记录进行仪器响应、传播路径和场地响应的校正,使用遗传算法对于田地震序列103次ML≥3.0地震的震源谱进行了反演计算。结果表明:1前震序列的拐角频率相对偏低,扣除地震矩影响后,余震的拐角频率有明显增高的趋势;2距离震中最近的YUT台资料分析显示,于田地震几次前震活动的台站观测谱相对较为独立,与主震以及彼此之间的相关系数均不高;3序列事件和主震的台站观测谱相关系数的变化反映了序列发展的动态破裂过程,相关系数较低的事件基本上偏离了主震破裂区而发生在分支断裂上;4前震的震源谱中存在异常的高频成份。     

1995年陡河地震前后小震震源参数及其相互关系

在震源谱的多台多震尾波解法的基础上,应用Brune震源模型,求解1995年10月6日陡河ML5.9地震前后在其周围发生的48次小地震的震源因子和震源波谱,进而估算了地震矩、拐角频率和应力降,并讨论了它们之间的相互关系和在陡河地震前后随时间的变化．结果表明:小震震源因子表现出相当大的稳定性,震源因子的峰值频率和随频率的变化形态都很相似．发现震前一年左右曾多次出现应力降较高的事件．需要指出的是,考虑到应力降由地震矩和拐角频率计算得到以及地震矩和拐角频率的相关关系,实际上,这里的高应力降指的是拐角频率高于正常值的情形．      

一种矩张量反演新方法及其对2008年汶川MS8.0地震序列的应用

本文利用全球地震台网（GSN）的宽频带与长周期地震波形资料,采用矩张量反演的新方法,反演得到了2008年5月12日四川省汶川县MS8.0地震及其7个较大余震（MS5.0-6.0）的矩张量解与震源时间函数等震源参数.文章首先简要叙述矩张量反演新方法的理论背景和技术途径,并以汶川大地震的一个余震为例阐述了具体的实现过程;然后给出包括主震在内的8次地震的矩张量解和震源时间函数;最后分析探讨这些结果的构造意义.本文提出的矩张量反演新方法,不但与全球矩心矩张量（GCMT）一样,可以给出点源矩张量解,而且还可以给出点源的震源时间函数. 反演得到的汶川大地震的7个较大余震的震源时间函数表明,即使是中等强度的地震也可能有复杂的震源过程;汶川大地震的多数余震发生在以逆冲为主、兼具小量走滑分量的龙门山断裂带的主断裂上,但很可能有些余震则发生在主断裂附近的次级走滑断裂上.      

近场全波形反演：芦山7.0级地震及余震矩张量解

2013年4月20日四川芦山发生了7.0级地震.该地震是自2008年汶川8.0级地震以来发生在龙门山断裂带上的最强地震,截至5月5日,四川省数字遥测台网记录到了52个（M_L≥4.0）余震.我们搜集了信噪比相对较高的宽频带数据,根据已有的三个分别适合于四川盆地和龙门山断裂东南地区、西南山区、西北山区速度模型,利用TDMT_INVC方法计算矩张量,获得了芦山主震及其16个余震的高质量矩张量解;根据P轴、B轴和T轴的倾角对主震及余震的断层类型进行了分类,结果表明：芦山地震的震级为M_w6.46（矩震级）,地震矩为0.61×10¹⁹Nm,震源深度为13 km,震源机制类型为逆断层.芦山余震类型简单,主要为逆冲型;芦山地震序列震中分布于龙门山断裂带西南端的彭县—灌县断裂上,本研究获得的主震和余震的震源机制主要为逆冲型,这表明了龙门山断裂西南端主要为逆冲特征.芦山地震震源性质与汶川地震相近,发震间隔时间短,且均发生于同一构造单元内,表明了该地震与汶川地震可能有密切关系.     

唐山地震前后京津唐张地区的应力场

作者从断裂力学观点研究了圆盘形纯剪切破裂产生的远场位移,导出了远场位移的解析表达式.用此表达式计算了远场位移脉冲波形及其谱.求得了P波和S波拐角频率F_cα、f_cβ与震源破裂半径的关系式,t_2α和t_2β分别是P波与S波的震源脉冲宽度在震源球上的平均值.V_r是破裂速度,取0.775β,α、β是P波与S波速度.介质吸收的影响,假定在略去路径滞后因子e-（iωr）/c后是最小相移滤波器,应用希尔伯特（Hilbert）变换导出了介质的频率响应公式 式中, r为震源距,c为波速,Q为介质品质因数,ω_m为高频截止频率.它与富特曼（Futterman）公式是不一样的.在考虑了介质吸收影响与仪器频率特性影响后,在频率域中求出了震源尺度2a和地震矩M₀.在时间域中,利用快速傅立叶反变换求出介质和仪器的脉冲响应后与震源脉冲折积,求得了合成地震图,用此图得到了地震图初动半周期T_2P与震源尺度2a或震源脉冲宽度t_2a之间的关系量板.用此量板在时间域中也求得震源尺度2a,与频率域中求得的结果基本一致. 利用从断裂力学观点研究震源破裂过程导出的关系式求出了京津唐张地区的区域剪切应力τ₀,式中η为地震辐射效率,取0.05;v为泊松比,取0.252;μ为刚性模量,取3.3×10¹¹达因/厘米².结果表明,1976年7月28日唐山大震前,唐山-丰南一带应力值最高，达870巴，离开唐山仅100公里左右的昌黎地区的应力就只有200巴左右.唐山大震后，地震破裂带中部的唐山一丰南一带，t₀降至100巴左右，而破裂带两头应力值仍很高，尤其是西南头的宁河一带，达400-500巴.直到11月15日此处发生6.9级强余震后应力值才降到200巴左右.     

汶川地震科学钻探3号井孔附近微震震源参数确定及其意义

基于微震监测仪器频带及波形记录的时频特征,我们筛选出汶川地震科学钻探3号井孔周围微震台阵2012年记录的218个t_s-t_p1 s的微震.通过盖戈法与和达法相结合确定微震震源的几何参数,发现这些微震分布大体呈NE-SW展布,与龙门山断裂带的走向基本一致.无论是采用《地震台站观测规范》中的量规函数(量规函数GF)还是李学政等(2003)的量规函数(量规函数LXZ),近震震级均与矩震级呈现出较好的线性关系,向震级小的一端延伸时都表现为M_LM_W,但采用李学政等的量规函数时该趋势更加明显.同时采用Brune和Boatwright震源谱衰减模型对观测震源谱的拟合表明拐角频率具有模型依赖性:基于Brune模型拟合得到的拐角频率大于基于Boatwright模型得到的拐角频率.无论基于哪一种模型,矩震级与拐角频率、破裂半径的对数线性关系均较弱,与应力降、视应力的对数则存在较好的线性关系.这些关系不支持应力降、折合能量、视应力为常数的观点,表明微震与大地震的震源物理过程存在差异.视应力与应力降成比例,比例系数小于0.5,表明破裂动力学模式符合SavageWood模式.近震震级与矩震级拟合关系Mw=a+bM_L中b的大小与应力降和地震矩的关系有关,△σ∝M_0~γ,则b=1/(1+γ),因此从b的大小可以粗略地判断应力降与地震矩的关系.本文对应于γ≈1的情况,与基于量规函数LXZ得到的近震震级与矩震级关系中系数b=0.53吻合,这说明仅从辐射能量的角度考虑量规函数LXZ较量规函数GF准确.     

Contribution to the knowledge of the magmatic evolution by the study of the variation of the coefficient of viscosity

The determination of the coefficient of viscosity of eruptive products gives useful elements to the knowledge of possible variations of composition and physical conditions of the magmas.     

Pictorial Series of the Manifestations of the Dynamics of the Earth

Studia Geophysica et Geodaetica -     

Evolution of the Source of the Active Region of the Sun

Geomagnetism and Aeronomy - The magnetic evolution of the active region of the Sun at the level of the photosphere has been studied. Magnetic synoptic maps of the 2007 Kitt Peak Observatory were...     

Magnetic stratigraphy of the Ordovician in the lower reach of the Kotuy River: the age of the Bysy-Yuryakh stratum and the rate of geomagnetic reversals on the eve of the superchron

Until recently, the existing data prevented the geophysicists from accurately dating the Bysy-Yuryakh stratum, which outcrops in the middle reach of the Kotuy River, constraining the time of its formation to a wide interval from the end of the Late Cambrian to the beginning of the Silurian. The obtained paleomagnetic data unambiguously correlate the Bysy-Yuryakh stratum to the Nyaian regional stage and constrain its formation, at least a considerable part of it, by the Tremadocian. This result perfectly agrees with the data on the Bysy-Yuryakh conodonts studied in this work and yields a spectacular example of the successful application of paleomagnetic studies in solving important tasks of stratigraphy and, correspondingly, petroleum geology. Within the Bysy-Yuryakh stratum, we revealed a large normal-polarity interval corresponding to the long (>1 Ma) period when the geomagnetic reversals were absent. This result, in combination with the data for the Tremadocian and Middle–Upper Cambrian sequences of the other regions, indicates that (1) the rate of occurrence of the geomagnetic reversals on the eve of the Ordovician Moyero superchron of reversed polarity was at most one reversal per Ma; (2) the superchron does not switch on instantaneously but is preceded by a certain gradual change in the operation conditions of the dynamo mechanism which, inter alia, manifests itself by the reduction of the frequency of geomagnetic reversals with the approach of the superchron. This finding supports the views according to which a process preparing the establishment of the superchrons takes place at the core–mantle boundary.     

On the determination of the non-linearity of the gravimeter scale

Summary The results of two different methods of the determination of the non-linearity of the gravimeter scale are compared; the economy is discussed as well.     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.     

The effect of the pollution of the ground layer of the atmosphere on the albedo

Summary A simple model is proposed suitable for studying the effect of the ground layer of the atmosphere, polluted by aerosol, on the albedo. This model is founded on solving the equation of transfer of radiative energy. The numerical results are discussed, particular attention being paid to the analysis of errors due to neglecting the multiple reflection of solar radiation on the aerosol particles. A method which would also include the multiple reflection is proposed, and the conditions under which the presence of the aerosol is responsible for an increase or decrease of the solar radiation balance on the Earth's surface, are analysed.     

On the solution of the reverse problem of geodesics with the help of the chord

Резюме Дается вывод формул для вычислений длины хорды, длины и азимутов нормаляных сечений и длины геодезической линии эллипсоида. Применяются тригонометрические функции данных величин и постоянные принятого ргфгргнц-эллипсоида, что является целесообразным при использовании вычислителQjных машин, особенно автоматов. Рабочие формулы рекомендуются в форме (8–11).

---

---

Dedicated to Professor František Fiala on His 85th Birthday

---

Address: Veveři 95, Brno.     

On the maintenance of the axisymmetric part of the flow in the atmosphere

Summary The maintenance of the axisymmetric component of the flow in the atmosphere is investigated by means of a steady-state, quasi-geostrophic formulation of the meteorological equations. It is shown that the meridional variations in the time-averaged axisymmetric variables can be expressed as the sum of three contributions, one being due to the eddy heat transport, another to the eddy momentum transport, and a third to the convective-radiative equilibrium temperature which enters the problem through the specification of a Newtonian form of diabatic heating. The contributions by the large scale eddies are evaluated through the use of observed values for the eddy heat and momentum transports.The contributions from each of the three forcing mechanisms to the temperature and zonal wind fields are invstigated individually and found to be of about equal importance. The sum of the three contributions are also presented for the temperature, the zonal wind, the stream function associated with the mean meridional circulation and the corresponding vertical motion. Although the results fail to reproduce the main observed features of the lower stratosphere, they are found to be in good agreement with observations in the middle latitude troposphere. At any pressure level, for example, the computed mean zonal wind has a jet-like profile and the axis of the jet is found to slope to the south with height, as observed in the atmosphere.Based in part on a thesis submitted by the first author as partial fulfillment of the requirements for the Ph.D. degree at the University of Michigan. — Publication No. 194 from the Department of Meteorology and Oceanography, The University of Michigan.