Characteristics of relativistic solar cosmic rays during the event of December 13, 2006

The characteristics of relativistic solar protons have been obtained using the methods of optimization based on the data of ground detectors of cosmic rays during the event of December 13, 2006, which occurred under the conditions of solar activity minimum. The dynamics of relativistic solar protons during the event has been studied. It has been indicated that two populations (components) of particles exist: prompt and delayed (slow). The prompt component with a hard energy spectrum and strong anisotropy manifested itself as a pulse-shaped enhancement at Apatity and Oulu stations, which received particles with small pitch-angles. The delayed component had a wider pitch-angle distribution, as a result of which an enhancement was moderate at Barentsburg station and at most neutron monitors of the worldwide network. The energy spectra obtained from the ground-based observations are in good agreement with the direct measurements of solar protons on balloons and spacecraft.     

盈江MS6.1地震序列近场横波分裂研究

使用横波分裂系统分析方法(SAM), 对2014年5月30日盈江M_S6.1地震震区内多个近场流动台站记录到的大量波形数据进行横波分裂研究． 研究结果表明, 盈江M_S6.1地震序列的快S波偏振方向为近NS向, 与区域主压应力方向一致． 主震发生后, 由于震源区应力状态的调整, 卡场台(KAC)快S波偏振方向发生逆时针偏转, 勐弄台(MNO)快S波偏振方向离散度减小, 并且由于受到研究区内断裂的影响, MNO台偏振方向较KAC台偏振方向更加离散． KAC台和MNO台的慢S波时间延迟均表现出主震发生前短时间内突然减小, 震后逐渐增大的变化特征, 这意味着临震前震源区地壳应力的释放和震后地壳应力的增强, 预示了后续余震的持续发生． 地震序列时间延迟平均滑动曲线起伏振荡, 表明了余震的发生伴随着震源区地壳应力的不断调整．      

青海地区S波分裂研究

本文利用国家地震台网及中国地震局“十五”期间在青海布设的30个宽频带地震台站记录到的远震波形数据,分别采用最小能量法和旋转相关法对SKS、SKKS和PKS波震相进行了偏振分析,计算了台站下方介质的各向异性分裂参数：快波的偏振方向（φ）和慢波延迟时间（δt）.本文研究结果表明,研究区多数台站下方的地震各向异性参数都表现出随方位角变化而变化的特征,可以用双层各向异性模型来解释.其中上层各向异性的快波偏振方向位于N65°E～N95°E之间,可能与中下地壳物质的流动有关;而下层各向异性的快波偏振方向位于N105°E～N135°E之间,可能为祁连块体NEE向的推移及导致的岩石圈缩短有关.此外,我们还发现,与周边的台站下方各向异性分裂参数相比,昆仑断裂附近两个台站GOM和DAW的各向异性特征急剧变化,其快波方向都与该断层近乎平行,这很可能暗示昆仑断裂已经切穿整个岩石圈;阿尔金断裂附近两个台站（LEH和HTG）无效分裂事件的方位分布与阿尔金断裂走向缺乏相关性,我们推测研究区内阿尔金断裂可能为地壳尺度的断裂.     

汶川MS8.0级大震前天基与陆基电场资料联合应用研究

利用边际谱方法对2008年5月12日四川汶川地震前天基和陆基电场变化情况进行了联合分析研究,结果表明:(1)空间电场方面,由重复轨道和连续轨道电场边际谱演化看出,电离层电场从2008年4月27日左右开始,其边际谱出现了明显的增强现象(增大1~2个量级),这种现象一直持续到汶川地震发生.(2)陆基电场方面,中法合作的松山...     

An automated cross‐correlation based event detection technique and its application to a surface passive data set

In studies on heavy oil, shale reservoirs, tight gas and enhanced geothermal systems, the use of surface passive seismic data to monitor induced microseismicity due to the fluid flow in the subsurface is becoming more common. However, in most studies passive seismic records contain days and months of data and manually analysing the data can be expensive and inaccurate. Moreover, in the presence of noise, detecting the arrival of weak microseismic events becomes challenging. Hence, the use of an automated, accurate and computationally fast technique for event detection in passive seismic data is essential. The conventional automatic event identification algorithm computes a running‐window energy ratio of the short‐term average to the long‐term average of the passive seismic data for each trace. We show that for the common case of a low signal‐to‐noise ratio in surface passive records, the conventional method is not sufficiently effective at event identification. Here, we extend the conventional algorithm by introducing a technique that is based on the cross‐correlation of the energy ratios computed by the conventional method. With our technique we can measure the similarities amongst the computed energy ratios at different traces. Our approach is successful at improving the detectability of events with a low signal‐to‐noise ratio that are not detectable with the conventional algorithm. Also, our algorithm has the advantage to identify if an event is common to all stations (a regional event) or to a limited number of stations (a local event). We provide examples of applying our technique to synthetic data and a field surface passive data set recorded at a geothermal site.     

2008年汶川MS8.0地震前定点观测电磁异常回溯性分析

对2008年汶川M_S8.0地震周围定点台站观测的电磁异常的相关研究进行了简要的回顾分析,以期加深对汶川地震孕震过程中电磁异常的解读.震中周围8个地电阻率台站震前出现不同形态的异常变化,结合震后地电阻率变化形态分析,仅近邻地震破裂带的成都和江油台异常变化符合已有的地电阻率孕震异常机理.在4月24日和5月9日大致沿南北地震带出现两条南北走向的低点位移线,成都台地磁转换函数、谐波振幅比、帕金森矢量和垂直极化强度等主要反映了地下介质电性在震前出现的异常变化;断裂带附近的地电场和电磁扰动在震前出现波形畸变和能量增强,距离较远的西昌台阵和天祝—松山台阵内地电场也出现功率谱能量增加和裂隙渗透方位角扰动等短期异常现象,甚至沿龙门山断裂带NE方向1300km外的河北电磁扰动台网震前数月也出现自观测以来最大幅度的异常变化.     

利用宽频带数字地震记录测定谱震级

谱震级是表征地震强弱的量度,是地震学研究中的一个重要的地震参数。在地震预测和工程地震研究中具有重要的实际意义。谱震级测定还可得到地震辐射的能量,再通过和标量地震矩的比较,可以估计地震震源处的视应力,而视应力也是震源物理研究中一个重要的物理量。考虑到我国地震台网分布不均匀的情况,本文尝试利用单台的宽频带波形资料进行谱震级的测定。选用CDSN昆明地震台(KMI)记录的近震宽频带波形资料,分析地震图垂直分量上的P波,采用 Duda等人提出的谱震级的测定方法,改进并提出了一种能够利用宽频带波形资料,直接用速度谱值测定地震谱震级的方法。所得谱震级结果与ISC给出的Ms震级对比后,发现相差不大;所得的能量结果与NEIC用远震体波测得的能量结果具有较好的一致性。用同样的方法,将来还可对不同地区更小的地震进行谱震级的测定,这无疑对地震学的进一步研究很有帮助。     

Solar proton event in December 2006

The variations in the rigidity spectrum and anisotropy of cosmic rays in December 2006 have been studied based on the surface measurements of the cosmic ray intensity at the global network of stations, using the method of global spectrographic survey. It has been indicated that the highest degree of anisotropy (to ～50%) with the maximal intensity of particles with a rigidity of 4 GV in the direction from the Sun (an asymptotic direction of about ?25° and 160°) was observed at 0400 UT on December 13. The parameters of the cosmic ray rigidity spectrum, which reflect the electromagnetic characteristics of the heliospheric fields during the studied period, have been determined when the surface and satellite measurements of protons in the energy range from several megaelectronvolts to several tens of gigaelectronvolts were jointly analyzed. The observed anisotropy and variations in cosmic rays in a wide energy range have been explained based on an analysis of the results.     

中天山及邻区S波分裂研究及其动力学意义

本文利用天山及其邻区布设的37个宽频带地震台站记录到的远震波形数据,分别采用最小能量法和旋转相关法对SKS和SKKS波震相进行了偏振分析,计算出了台站下方介质的S波分裂参数：快波的偏振方向（φ）和慢波延迟时间（δt）.本文研究结果表明：中天山内部大多数台站的各向异性快波方向呈NEE-SWW向,与天山走向平行,慢波时间延迟为0.4～1.7 s,这是塔里木、哈萨克斯坦的南北双向俯冲及其导致的天山地区岩石圈地幔南北向缩短变形的直接反映.本文研究发现中天山北部部分台站下方地震各向异性快波方向与慢波延时随方位角呈现规律性的变化,可能暗示该区上地幔各向异性不能仅用单层水平各向异性这一简单模式来解释.75°E以西的天山地区台站下方S波快波方向和延时具有强烈的横向变化,可能与研究区下方存在的小规模对流有关.中天山不同地段地震台站下方各向异性明显不同,进一步证实了天山地区构造变形的复杂性.     

南北构造带南段上地幔各向异性特征

对布设在南北构造带南段的中国地震科学探测台阵项目一期350个宽频带流动台站和中国地震台网90个宽频带固定台站记录的远震XKS(SKS、SKKS和PKS)波形资料作偏振分析,采用最小切向能量的网格搜索法和"叠加"分析方法求得每一个台站的XKS波的快波偏振方向和快、慢波的时间延迟,获得了南北构造带南段上地幔各向异性图像.结果显示研究区的各向异性具有明显的南北分区特征,北部的快波方向为近N-S方向,而南部主要表现为近E-W方向,且北部的平均时间延迟小于南部.分析表明,具有厚岩石圈的北部的各向异性主要由岩石圈变形引起,是一种垂直连贯变形模式;具有薄岩石圈的南部的各向异性主要由软流圈地幔流引起,缅甸和巽达板片的后撤/回转作用产生了指向西南的软流圈地幔流,在岩石圈底部和软流圈之间产生了一个水平差异运动,产生了一个与简单剪切一致的软流圈变形结构,从而产生了南部观测的各向异性.     

基于陆态网络GPS数据的电离层空间天气监测与研究

中国大陆构造环境监测网络(简称陆态网络)是以全球卫星导航定位系统(GNSS)为主,辅以多种空间观测技术,实时动态监测大陆构造环境变化,探求其对资源、环境和灾害的影响的地球科学综合观测网络.基于陆态网络约200个基准站的GPS观测数据,本文探讨了其在电离层空间天气监测与研究方面的应用.包括磁暴期间电离层暴扰动形态,大尺度电离层行进式扰动,太阳耀斑引起的电离层骚扰和低纬电离层不规则体结构等.研究结果表明:陆态网络布局合理,观测数据质量良好,完全可用于中国及周边地区电离层空间天气监测与研究,为进一步开展我国电离层空间天气预警和预报奠定了观测基础.     

Ionospheric response to the geomagnetic storm on 13–17 April 2006 in the West Pacific region

This paper presents an investigation of geomagnetic storm effects in the equatorial and middle-low latitude F-region in the West Pacific sector during the intense geomagnetic storm on 13–17 April, 2006. The event, preceded by a minor storm, started at 2130 UT on April 13 while interplanetary magnetic field (IMF) B_z component was ready to turn southward. From 14–17 the ionosphere was characterized by a large scale enhancement in critical frequency, foF2 (4～6 MHz) and total electron content (TEC) (～30TECU, 1TECU=1×10¹⁶el/m²) followed by a long-duration negative phase observed through the simultaneous ionospheric sounding measurements from 14 stations and GPS network along the meridian 120°E. A periodic wave structure, known as traveling ionospheric disturbances (TIDs) was observed in the morning sector during the initial phase of the storm which should be associated with the impulsive magnetospheric energy injection to the auroral. In the afternoon and nighttime, the positive phase should be caused by the combination of equatorward winds and disturbed electric fields verified through the equatorial F-layer peak height variation and modeled upward drift of Fejer and Scherliess [1997. Empirical models of storm time equatorial electric fields. Journal of Geophysical Research 102, 24,047–24,056]. It is shown that the large positive storm effect was more pronounced in the Southern Hemisphere during the morning-noon sector on April 15 and negative phase reached to lower magnetic latitudes in the Northern Hemisphere which may be related to the asymmetry of the thermospheric condition during the storm.     

Multipoint observations of Pi2 pulsations and correlation with dynamic processes in the near-Earth magnetotail on March 18, 2009

Simultaneous measurements from THEMIS spacecraft,GOES-11 and ground stations(Canadian Array for Realtime Investigations of Magnetic Activity or CARISMA,and 210°magnetic meridian or MM)on March 18,2009 allow the study of dynamic processes in the near-Earth magnetotail and corresponding Pi2 pulsations on the ground in great detail.Fast earthward flows along with traveling Alfvén waves and fast mode waves in the Pi2 band were observed by three Time History of Events and Macroscale Interactions during Substorms(THEMIS)probes(P3,P4 and P5)in the near-Earth plasmasheet.At the mid-to high-latitude nightside,the CARISMA stations located near the foot points of the three probes recorded Pi2s with two periods,about 80 s after the earthward fast flows observed by the P4 probe.The long-period Pi2(140–150 s)belongs to the transient response Pi2(TR Pi2),since the travel time of the Alfvén waves between the plasma sheet and CARISMA stations is very close to half the period of the long-period Pi2.The short-period Pi2(60–80 s)has the same period band as the perpendicular velocity of the fast flows,which indicates that it may relate to the inertial current caused by periodic braking of the earthward fast flows.The 210°MM stations located at the low-latitude duskside also observed Pi2s with the same start time,waveform and frequency,about~120 s after the earthward fast flows.Strong poloidal oscillations are shown by GOES-11(~23 MLT)and the compressional component(Bb)is highly correlated with H components of the 210°MM stations,whereas the other two components(Br and Be)are not.These results confirm that the low-latitude Pi2s are generated by cavity mode resonance,which is driven by an impulsive broadband source in the near-Earth magnetotail.     

基于希尔伯特-黄变换的汶川地震记录分析

利用HHT方法对6个台站的汶川8.0级地震加速度记录进行EMD分解,进而对地震波输入结构的总能量进行分解,得到地震动能量在IMF分量之间的分布关系。同时通过对各分量边际能量谱的分析,可以从另一方面得到地震波的频谱特性,边际能量谱峰值频率与地震波傅立叶谱的卓越频率相同,边际能量谱也可以作为一种估算震动卓越周期的方法。     

Free‐field rotations during earthquakes: Relevance on buildings

Soil rotations around horizontal axes, during an earthquake, are studied through records collected by closely spaced arrays of strong motion accelerometers. The cross power spectrum of accelerations at nearby stations has been generally utilized to describe the spatial distribution of the motion. A number of cross spectra have been obtained during the training of these arrays. To take profit of these elaborations, a mathematical relation is established between the cross power spectrum and the power spectrum of the rotation. Rotation data presented by Liu et al, concerning 52 earthquake records collected at a single station in Taiwan, are compared with rotation data computed according to our procedure. The two series of data are suitably normalized to the peak horizontal acceleration. The data are shown in function of the distance from epicentre. The same ratio, computed according to our procedure, is in good agreement with the average value of these data. Direct measurements and the present approach have lead to evaluations of rotation higher than those predicted by mathematical investigations on the basis of the wave propagation theory, for comparable circumstances. The relevance of this input motion for relatively tall structures is examined, with reference to the structural effects that the horizontal motion concurrently provides. Meaningful will be ranked those effects of the order of magnitude of 20% or higher than those implied by the horizontal excitation. For understanding the relevance on building structures, the procedure has two areas of concern: 1) the coherence implicit in the cross power spectra, which depends on the interpolation process of the original records, collected in the arrays of instruments, and 2) the relative importance of the vertical to the horizontal input motion. As to the second item, the relevance of the rotation component on structures largely depends on the relative importance of the vertical to the horizontal input motion. When the records in an area a few km from the epicenter are considered, the response spectrum of vertical motion is comparable and in some records even higher, than that of horizontal motion, over the entire range of frequencies. This has been observed as well for the 2009 earthquake event of L'Aquila, Italy, and that at the Christchurch (New Zealand) 2011. When the response spectrum of vertical motion is comparable to that of horizontal motion, the effects of rotational motions on most engineering structures can be meaningful. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydromagnetic spectroscopy of the magnetosphere with Pc3 geomagnetic pulsations along the 210° meridian

Analysis of Pc3 observational data along the 210° magnetic meridian showed a complicated frequency-latitude structure at middle latitudes. The observed period-latitude distributions vary between events with a “noisy source”: the D component has a colored-noise spectrum, while the spectrum of H component exhibits regular peaks that vary with latitude, and events with a “band-limited source”: the spectral power density of the D component is enhanced at certain frequencies throughout the network. For most ULF events a local gap of the H component amplitude has been exhibited at both conjugate stations at L ≃ 2.1. A quantitative interpretation has been given assuming that band-limited MHD emission from an extra-magnetospheric source is distorted by local field line resonances. Resonant frequencies had been singled out with the use of the asymmetry between spectra of H and D components. Additionally, a local resonant frequency at L ≃ 1.6 was determined by the quasi-gradient method using the data from nearly conjugate stations. The experimentally determined local resonance frequencies agree satisfactorily with those obtained from a numerical model of the Alfven resonator with the equatorial plasma density taken by extrapolation of Carpenter-Anderson model. We demonstrate how simple methods of hydromagnetic spectroscopy enable us to monitor simultaneously both the magnitude of the IMF and the magnetospheric plasma density from ULF data.     

Backtracking of Noble Gas Measurements Taken in the Aftermath of the Announced October 2006 Event in North Korea by Means of PTS Methods in Nuclear Source Estimation and Reconstruction

The announced October 2006 nuclear test explosion in the Democratic People’s Republic of Korea (DPRK) has been the first real test regarding the technical capabilities of the verification system built up by the Vienna-based Provisional Technical Secretariat (PTS) of the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) to detect and locate a nuclear test event. This paper enhances the resolution of the DPRK events’ xenon source reconstruction published by Saey et al. (2007, “A long distance measurement of radioxenon in Yellowknife, Canada, in late October 2006”, GRL, Vol. 34, L20802) that was based solely on radio-xenon measurements taken at the remote radionuclide station in Yellowknife, Canada by involving additional measurements taken by a mobile noble gas system deployed quite close to the event location in the Republic of Korea (ROK). Moreover the horizontal resolution of the forward and backward atmospheric transport modelling methods applied for the source scenario reconstruction has been enhanced appropriately to reflect the considerably shorter source-receptor distances examined in comparison to the previously published source reconstruction. It is shown that the ¹³³Xe measurements in Yellowknife could register ¹³³Xe traces from the nuclear explosion during the first 3 days after the event, while the mobile measurements were rather sensitive to releases during days 2–4 after the explosion. According to the analysis, the most likely source scenario would consist of an initial (possibly up to 21 h delayed) venting of 1 × 10^?15 Bq ¹³³Xe during the first 24 h, followed by a two orders of magnitude weaker seepage during the following 3 days. Both measurements corroborate the scenario of a rather rapid venting and soil diffusion of the ¹³³Xe yielded during the explosion. While the Swedish mobile measurements were crucial to enhancement of the reconstruction of the source scenario, given the installation status of the IMS xenon network at the time of the event, a sensitivity analysis revealed that the fully developed network would have been able to detect ¹³³Xe traces from the Korean explosion at a number of stations and allowed for an even better constraint on the release function. The station Ussuriysk, Russia, being in operation in 2006, would have registered ¹³³Xe within 1 day and with a three orders of magnitudes stronger signal compared to the detection at Yellowknife.     

Coherence between the fields of low-frequency seismic noise in Japan and California

The relationship is considered between the statistics of the field of low-frequency seismic noise which was synchronously recorded by two broadband seismic networks in Japan (78 stations) and California (81 stations). The analysis is based on the data for seven years of observations (2008–2014). For each network, the daily time series of the median values are constructed for five parameters of seismic noise: kurtosis (excess), minimal normalized entropy of the distribution of the squared wavelet coefficients, generalized Hurst exponent, support width of the singularity spectrum, and index of linear predictability. The median values for each parameter were calculated on a daily basis over all the stations of the networks and resulted in a time series containing 2557 data points of the integral characteristics of the noise with a daily time step. The use of the median values of the noise parameters avoids considering the effects of the gaps in recording by individual stations and provides the continuous time series as the integral characteristic of the whole network. Next, for each network, the aggregate signals were calculated for the obtained five-variate time series. By construction, the aggregate signal is a scalar signal which maximally accumulates the most general variations that are simultaneously present in all the analyzed signals and, at the same time, rejects the components that are only characteristic of a single process. The final step of the analysis consists in estimating the evolution of the quadratic spectrum in the moving time window with a length of one year. It is shown that during the considered interval of the observations, the coherence is characterized by the increasing linear trend, which independently supports the previous conclusion about the enhancement of the synchronization between the parameters of the global seismic noise.     

Preliminary analysis of crustal shear-wave splitting in the Sanjiang lateral collision zone of the southeast margin of the Tibetan Plateau and its tectonic implications

The collision of the Indian and Eurasian plates, to the east of the eastern Himalayan syntaxes, forms the Sanjiang lateral collision zone in the southeast margin of the Tibetan Plateau, where there are intense crustal deformation, active faults, earthquakes, as well as a metallogenic belt. Given the lack of adequate seismic data, shear-wave splitting in this area has not been studied. With seismic data from a temporary seismic linear array, as well as permanent seismic stations, this paper adopts the identification on microseismic event to pick more events and obtains shear-wave splitting parameters from local earthquakes. From the west to the east, the study area can be divided into three subzones. The “fast” polarization (i.e. the polarization of the fast shear wave) varies gradually from NNW to NS to NNE in these three subzones. The time delay of the slow shear wave (i.e. the time difference between the two split shear waves) also increases in the same direction, indicating the presence of seismic anisotropy above 25 km in the crust. Both shear-wave splitting parameters are closely related to stress, faults and tectonics. The scatter and the “dual” (i.e. two) dominant orientations of the fast polarizations at several stations indicate strong distortions caused by nearby faults or deep tectonics. The anisotropic parameters are found to be related to some degree to the metallogenic belt. It is worth to further analyse the link between the anisotropic pattern and the metallogenic area, which suggests that shear-wave splitting could be applied to study metallogeny. This paper demonstrates that the identification on microseismic event is a useful tool in detecting shear-wave splitting details and exploring its tectonic implications.     

Localized structure in the cusp and high-latitude ionosphere: a modelling study

The ionospheric signature of a flux transfer event (FTE) seen in EISCAT radar data has been used as the basis for a modelling study using a new numerical model of the high-latitude ionosphere developed at the University of Sheffield, UK. The evolution of structure in the high-latitude ionosphere is investigated and examined with respect to the current views of polar patch formation and development. A localized velocity enhancement, of the type associated with FTEs, is added to the plasma as it passes through the cusp. This is found to produce a region of greatly enhanced ion temperature. The new model can provide greater detail during this event as it includes anisotropic temperature calculations for the O⁺ ions. This illustrates the uneven partitioning of the energy during an event of this type. O⁺ ion temperatures are found to become increasingly anisotropic, with the perpendicular temperature being substantially larger than the parallel component during the velocity enhancement. The enhanced temperatures lead to an increase in the recombination rate, which results in an alteration of the ion concentrations. A region of decreased O⁺ and increased molecular ion concentration develops in the cusp. The electron temperature is less enhanced than the ions. As the new model has an upper boundary of 10 000 km the topside can also be studied in great detail. Large upward fluxes are seen to transport plasma to higher altitudes, contributing to the alteration of the ion densities. Plasma is stored in the topside ionosphere and released several hours after the FTE has finished as the flux tube convects across the polar cap. This mechanism illustrates how concentration patches can be created on the dayside and be maintained into the nightside polar cap.