A deep-learning-based approach for seismic surface-wave dispersion inversion (SfNet) with application to the Chinese mainland

Surface-wave tomography is an important and widely used method for imaging the crust and upper mantle velocity structure of the Earth. In this study, we proposed a deep learning (DL) method based on convolutional neural network (CNN), named SfNet, to derive the v_S model from the Rayleigh wave phase and group velocity dispersion curves. Training a network model usually requires large amount of training datasets, which is labor-intensive and expensive to acquire. Here we relied on synthetics generated automatically from various spline-based v_S models instead of directly using the existing v_S models of an area to build the training dataset, which enhances the generalization of the DL method. In addition, we used a random sampling strategy of the dispersion periods in the training dataset, which alleviates the problem that the real data used must be sampled strictly according to the periods of training dataset. Tests using synthetic data demonstrate that the proposed method is much faster, and the results for the v_S model are more accurate and robust than those of conventional methods. We applied our method to a dataset for the Chinese mainland and obtained a new reference velocity model of the Chinese continent (ChinaVs-DL1.0), which has smaller dispersion misfits than those from the traditional method. The high accuracy and efficiency of our DL approach makes it an important method for v_S model inversions from large amounts of surface-wave dispersion data.     

Variation of Q value before and after the 1999 Xiuyan, Liaoning Province, M =5.4 earthquake on the basis of analysis on attenuation dispersion of P waves

A method for determining medium quality factor is developed on the basis of analyzing the attenuation dispersion of the arrived first period P wave. In order to enhance signal to noise ratio, improve the resolution in measurement and reduce systematic error we applied the data resampling technique. The group velocity delay of P wave was derived by using an improved multi-filtering method. Based on a linear viscoelastic relaxation model we deduced the medium quality factor Q _m, and associated error with 95% confidence level. Applying the method to the seismic record of the Xiuyan M=5.4 earthquake sequences we obtained the following result: (1) High Q _m started to appear from Nov. 9, 1999. The events giving the deduced high Q _m value clustered in a region with their epicenter distances being between 32 and 46 km to the Yingkou station. This Q _m versus distance observation obviously deviates from the normal trend of Q _m linearly increasing with distance. (2) The average Q _m before the 29 Dec. 1999 M=5.4 earthquake is 460, while the average Q _m between the M=5.4 event and the 12 Jan. 2000 M=5.1 earthquake is 391, and the average Q _m after the M=5.1 event is 204. Foundation item: State Key Project of Science and Technology during the Tenth Five-year Plan (2004BA601B01-03-01).     

基于洱海水生态特征的流域最大日负荷总量控制

基于洱海水生态历史数据及现状资料,采用概率密度分布曲线法及水生生物基准相结合的方式计算洱海总氮(TN)、总磷(TP)、高锰酸盐指数(COD_(Mn))和氨氮(NH_3-N)的控制目标,目标值分别为0.36、0.026、4和0.28 mg/L.再根据该水质目标,得到洱海的最大日负荷(TMDL)总量.TMDL总量采用线性规划法计算,其中污染物响应系数矩阵通过MIKE 21二维水质模型计算所得.安全容余(MOS)则通过一阶误差分析法确定.经过一系列的计算,最终确定洱海的TMDL总量控制计划.计算结果表明,洱海TN、TP、COD_(Mn)和NH_3-N的TMDL总量分别为2005.989、149.671、19258.844和1348.119 kg/d,其MOS所占比例分别为6.152%、5.570%、4.380%和5.021%,表明洱海农业面源污染为该流域主要的污染形式,其允许最大排放量约占全部允许排放量的90%左右.     

Estimation of codaQ using the maximum likelihood method

A seismic coda wave from local earthquakes is an indication of heterogeneity in the crust and upper mantle, and codaQ (Q _c ) is the parameter defining the temporal decay of the coda wave. ButQ _c estimates obtained with the conventional least-square method are not based on any solid statistical background, assuming the Gaussian distributiona priori. In this study, we propose a statistically reliable estimation method of estimatingQ _c using the maximum likelihood method, and show its validity and usefulness with the data from the 1986 Joint Seismological Research in the western Nagano Prefecture. We found first that theQ _c estimation with the maximum likelihood method is statistically valid and its reliability can be checked with the -square test. Next,Q _c around Ohtaki village, within the studied area, is estimated using the maximum likelihood method. The averageQ _c value at low frequencies (up to 2 Hz) is one of the lowest in the world, which is in harmony with the geological setting of the site condition directly above an active fault.     

WEIGHTED VERTICAL STACKING IN CRUSTAL SEISMIC REFLECTION STUDIES ON THE CANADIAN SHIELD*

Seismic reflection methods are being developed at the University of Manitoba to aid in determining fine crustal structure in the Precambrian of Manitoba and northwestern Ontario. Present-day environmental concern as well as mineshaft conditions necessitate the detonation of several smaller charges repeated, say, I times and followed by ‘vertical’ stacking. To obtain the familiar √I improvement in signal-to-noise (S:N) amplitude ratio applying the straight-sum (SS) method, one assumes, among other things, that both S:N ratio and signal variance are the same on all traces. Dropping these assumptions, as we must for our data, it becomes necessary to apply weighting coefficients to optimize the S:N ratio of the stacked trace. We still assume the signal shapes to be the same for repeated shots, so for the jth trace on the record of the ith shot we model the time series as: t_ij=a_i (s_j+n_ij); where a_i is a scaling factor. The proper weights w_i are then shown to be proportional to σ_si/σ²_ni where σ² is variance, or to γ_i/a_i where γ_i is S:N power ratio. Applying the weighted-stack (WS) method gives S:N amplitude ratios which are, on average, 55% of the optimal ratios expected from WS theory compared with only 24% for the SS method. The 45% shortfall in WS performance is ascribed mainly to trace-alignment (or time-delay) errors. Varying noise levels on individual traces, slight dissimilarity of signal shape, and correlated noise may also contribute to a lesser extent (in decreasing order of significance). This WS method appears to strike a good practical balance between S:N improvement and processing efficiency.     

Intrinsic absorption and scattering attenuation in the southern part of the Netherlands

We present the first systematic study of attenuation derived from the S-wave coda in the frequency range 1-32 Hz for the southern part of the Netherlands and its surroundings. For this we used two methods, the codaQ (Q _c) method and the Multiple Lapse Time Window (MLTW) method. In the interpretation of the results both single and multiple scattering in a half space are considered. Our aim is to validate these interpretations in our region and to try to identify theeffects of attenuation due to intrinsic absoprtion (Q _i)and scattering attenuation (Q _s). For this we analyzedmore than 100 3-component high-quality digital seismograms from 43 crustalevents and 23 different stations in the Netherlands, Germany and Belgium.Coda Q results show smaller Q _c (=Q ₀fⁿ) values for epicentral distances shorter than 25 km (Q ₀=90) compared to larger epicentral distances (Q ₀=190), but similar frequency dependence (f^-0.9). Interpretation of MLTW results provided a seismic albedo smaller then 0.5, suggesting that the intrinsic absorption dominates over scattering in this region. Both Q _i and Q _s show similar frequency dependences as Q _c. These results are comparable to those obtained in other areas, but we also show that more sophisticated models are required to remove ambiguities in the interpretation. For short lapse times and shortevent-station distances we find for the simple half space model a correspondinginterpretation of both methodologies, where Q _c correspondsto Q _t, suggesting that a model with single scattering in ahalf space is appropriate. For long lapse times and long event station distances, however, we find that the S-wave coda is, most probably, too much influenced by crust-mantel heterogenities and more sophisticated Qinversion models using larger data sets are required for more reliable attenuation estimates.     

The Conversion Relationship between the Local Earthquake Magnitude and Surface Wave Magnitude in the Inner Mongolia Digital Seismic Network

Using 116 earthquakes over M_L3.8 in the Inner Mongolia region from 2008 to 2015, the local earthquake magnitude M_L and surface wave magnitude M_S are remeasured. Based on norm linear regression(SR1 and SR2) and norm(OR) orthogonal regression method, we established the conversion relationship between M_L and M_S. The results were tested with Gaussian disturbance. The result shows that the orthogonal regression method(OR) result has the best fitting curve, and the conversion relation is M_S=0.96 M_L-0.10. The difference between our result and Guo Lücan's(M_S=1.13 M_L-1.08) may be caused by regional tectonic characteristics. M_(S Inner Mongolia) value is significantly higher than the M_(S empirical) value, with an average difference of 0.23, the difference distribution of empirical relation and the rectified relation is in the range of 0.2-0.3.     

A Type-Curve Method for the Analysis of Pumping Tests with Piecewise-Linear Pumping Rates

Aquifer hydraulic parameters are commonly inferred from constant-rate pumping tests, while variable pumping rates are frequently encountered in actual field conditions. In this study, we propose a generally applicable dimensionless form of the analytical solution for variable-rate pumping tests in confined aquifers. In particular, we adopt a piecewise-linear fitting of variable pumping rates and propose a new type-curve method for estimating the hydraulic conductivity (K ) and specific storage (S_s ) of the investigated confined aquifer. For each test, a series of type curves, which depend on the variable pumping rates, the location of observation wells and the introduced first dimensionless inflection time, need to be provided for matching the observed drawdown data on a log-log graph. We first demonstrate the applicability and robustness of this method through a synthetic pumping test. Subsequently, we apply this method to analyze drawdown data from four pumping tests conducted within a multilayered aquifer/aquitard system in Wuxi city, Jiangsu Province, China. The parameter estimates are then compared with those reported by PEST. The K and S_s values estimated by the new type-curve method are found to be quite close to PEST-based estimates. Parameter estimation results demonstrate the difference in K and S_s values between observation wells. The difference could be attributed to the spatial heterogeneity in K and S_s . A future research topic may focus on the characterization of K and S_s heterogeneity with the currently available drawdown data from variable-rate pumping tests.     

IMPROVED PHASE-ELLIPSE METHOD FOR IN SITU GEOPHONE CALIBRATION*

For amplitude and phase response calibration of moving-coil electromagnetic geophones two parameters are needed, namely, the geophone natural frequency, f₀, and the geophone upper resonance frequency f_u. The phase-ellipse method is commonly used for the in situ determination of these parameters. For a given signal-to-noise ratio, the precision of the measurement off₀ andf_u depends on the phase sensitivity, f(δφ/δf) For some commercial geophones f(δφ/δf) atf_u can be an order of magnitude less than the sensitivity atf₀. In this paper we present an improved phase-ellipse method with increased precision. Compared to measurements made with the existing phase-ellipse methods, our method shows a 6- and 3-fold improvement in the precision, respectively, on measurements of f₀ andf_u on a commercial geophone.     

Hyperspectral feature analysis of chlorophyll a and suspended solids using field measurements from Taihu Lake,eastern China

Abstract

We report results of three field campaigns conducted at 39 stations. At each station, we measured reflectance spectra in situ and collected water samples for measuring chlorophyll a (CHL) and suspended solids (SS) concentrations in the laboratory. To identify the indicative bands and develop suitable estimation models for CHL (C _CHL) and SS (C _SS) concentrations in Taihu Lake, a spectral-feature method and a derivative method were applied. The following conclusions were drawn: (a) the critical C _CHL and C _SS probably causing their spectral variation are, respectively: 0, 10, 50 and 75 μg L^?1, and 0, 10, 50 and 100 mg L^?1; (b) the derivative method is better than the spectral-feature method for estimating C _CHL and C _SS; (c) the optimal variable for CHL is a reflectance second-order derivative at 501 nm or a reflectance first-order derivative at 698 nm; the optimal variable for SS can change when its concentration is low and the range is narrow; otherwise, the optimal variable is a reflectance first-order derivative at 878 nm; and (d) the CHL and SS have an effect on one another's retrieval. The C _CHL estimation accuracy would benefit from narrowing the C _SS range. With C _CHL increasing and its range broadening, the corresponding C _SS estimation accuracy decreases gradually.     

Application of scaling-rule theory in crustal rock fracture to studying characteristics of seismological precursors associated with M=6.1 Shandan-Minle earthquake

In the paper, we introduce Allegre＇s scaling-rule theory of rock fracture and the probability to develop a method for predicting earthquake occurrence time on its basis. As an example, we study the characteristics of seismological precursors （seismic spatial correlation length and coda Qc） associated with the earthquake （M=6.1） occurred in Shandan-Minle, Gansu Province. The results show an increasing trend of seismic spatial correlation length and coda Qc before the earthquake. And a power exponent relation is used to fit the increasing variation form of these two parameters. The study has provided a basis for creating a method and finding indexes to predict the earthquake occurrence time by using the monitored seismic spatial correlation length and coda Qc.     

Study on the processing method of nighttime CO2 eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the “true value” of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO₂ flux data) appearing in post-processing CO₂ flux data. The ratio of effective and reliable nighttime EC CO₂ flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO₂ flux data and its effect analysis on estimating CO₂ flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO₂ flux using EC method, and introduce the general theory and method for processing nighttime CO₂ flux data. By analyzing the relationship between nighttime CO₂ flux and air fraction velocity u*, we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u*_c) for screening the effective nighttime CO₂ flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO₂ flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.

     

Decentralized static output‐feedback H∞ controller design for buildings under seismic excitation

In this work, we present a new method in designing static output‐feedback H_∞ controllers suitable for vibrational control of buildings under seismic excitation. The method produces a Linear Matrix Inequality (LMI) formulation that allows obtaining static output‐feedback controllers with different information structure constraints by imposing a convenient zero–nonzero structure on the LMI variables. The application of the proposed methodology is illustrated by designing centralized and decentralized velocity‐feedback H_∞ controllers to mitigate the seismic response of a five‐story building. Copyright © 2011 John Wiley & Sons, Ltd.     

Constraining the crustal structure under the central and western Tian Shan based on teleseismic receiver functions and gravity anomalies

The Tian Shan is a vast range that spans several countries in Asia. Understanding its evolutionary history may provide valuable insights into intracontinental orogenic dynamics. In this study, we explored the crustal characteristics of the Tian Shan and their relationships to the tectonic evolution of the region. A new H-stacking method that combines the P receiver function and gravity anomalies was used to estimate the thickness and ratio of P- to S-wave velocities (v_P/v_S) for 91 broadband seismic stations in the central and western Tian Shan. Our results revealed significant lateral variations in crustal thickness and v_P/v_S. A ~45-km-thick crust and an intermediate-high v_P/v_S (~1.74–1.84) were found in the Kazakh Shield and Tarim Basin, which we interpreted to indicate a mafic crystalline basement and lower crust. The central Tian Shan varied greatly in crustal thickness (40–64 km) and v_P/v_S ratio (1.65–2.00), which may be due to crustal shortening, mafic underplating, and crustal melting. In contrast, we observed a relatively thin crust (42–50 km) with an intermediate v_P/v_S ratio (~1.78) in the western Tian Shan. The differences in the crustal structures between the western and central Tian Shan imply that the Talas-Fergana Fault may be trans-lithospheric.     

Qβ tomography under the crust and upper mantle in eastern China

On the basis of data of long period Rayleigh surface wave, we select 43 two-station paths which cover the eastern China thoroughly. By using the improved method of multi-filtration, we obtain the group velocity and amplitude spectrum, and then get attenuation factor for each paths. We employ Talentola inversion method to get local attenuation factor, and further invert the three-dimension Q _β image under the crust and upper mantle in the eastern Chinese continent. The Q _β image shows the following basic characters. There is correlation between the seismic activity and Q _β structure under the crust and upper mantle in North China region. The Yangtze block begins to collide with and subduct to the North China block from the southern border of the Qinling in the southern Shaanxi. In the large part of Yangtze quasi-platform appear an obvious high Q _β area at 88 km deep. In the east of Sichuan depression platform, the juncture of Sichun and Guizhou, and the Jiangnan block near the juncture of Guizhou and Hunan, the lateral variation of Q _β in the crust is little, and there is a high-Q _β layer no thinner than 40 km in the top mantle. In the Dian-Qian fold and fracture region between Yunnan and Guizhou, the vertical variation of Q _β at the region of the crust and upper mantle is little, there is a low-Q _β layer in the top mantle, about 40 km thick, low-Q _β layer of the upper mantle begins to appear at about 95 km deep. In the east of Yangtze quasi-platform and the central and eastern part of the South China fold system, the Moho is smooth, the lateral variation of Q _β in the crust is also little, low-Q _β layer of the upper mantle begins to appear at about 85 km deep.     

Generating ground motion by two new techniques of adding harmonic wave in the time domain and approximating to response spectrum as a whole

ＧｅｎｅｒａｔｉｎｇｇｒｏｕｎｄｍｏｔｉｏｎｂｙｔｗｏｎｅｗｔｅｃｈｎｉｑｕｅｓｏｆａｄｄｉｎｇｈａｒｍｏｎｉｃｗａｖｅｉｎｔｈｅｔｉｍｅｄｏｍａｉｎａｎｄａｐｐｒｏｘｉｍａｔｉｎｇｔｏｒｅｓｐｏｎｓｅｓｐｅｃｔｒｕｍａｓａｗｈｏｌｅＣＨＡＮＧ...     

2020年和林格尔M_L4.5地震微震匹配定位及发震构造探讨

内蒙古和林格尔地处鄂尔多斯块体北缘阴山地震带内,历史上6级以上强震频发。2020年3月30日和林格尔发生M_L4.5地震,打破了自2005年以来阴山地震带M_L4.0以上地震的长期平静。研究此次地震序列的发震构造对区域应力状态和地震危险性分析有重要作用,然而内蒙古地震台网台站较为稀疏,相对于华北其他地区地震监测能力较低,对和林格尔地震序列的分析造成不利影响。本文采用匹配定位检测方法(ML)检测区域台网目录遗漏的微震,并对检测到的地震事件进行精定位。匹配定位方法共检测到序列中可定位的地震事件61个,约为台网目录的1.3倍,可定位地震约为台网目录的2.9倍。地震序列重定位结果中余震整体呈现NE向的分布,与使用CAP方法得到主震震源机制NE向节面走向基本一致,破裂为正断走滑型,显示和林格尔M_L4.5地震序列发震构造为NE走向、NW倾向的拉张性断层。结合计算结果和区域地震地质构造特征,认为M_L4.5地震序列发震构造为NE走向的岱海-黄旗海盆地南缘断裂。     

基于速率-状态依从摩擦定律的前郭震群余震活动率及模型参数相关性研究

基于速率-状态依从摩擦定律的地震活动率时空预测模型,以同震库伦应力变化作为模型初始应力扰动,模拟了2013年吉林前郭MS5.8震群的余震活动率变化。考虑模型参数相关性,在模拟中采用2种不同的拟合方案,一是余震持续时间t_a不固定条件下的拟合,二是余震持续时间t_a固定条件下的拟合。结果显示,t_a不固定条件下的拟合方式可获得较好的AIC评价,适用于震后早期的趋势判定; t_a固定条件下的拟合计算耗时更短,拟合误差更小,理论模拟结果与前郭震群实际地震时序特征更为吻合。采用该方案对截至2016年10月24日的余震活动率变化进行了回溯性预测检验,结果显示模型预期的余震日频次与实际记录呈较好的正相关关系。研究还发现,主震破裂面附近的同震应力影区导致震后早期模型预测值相对于实际偏低,说明前郭序列余震活动可能还存在其他触发机制。     

2014年金寨M_L3.9震群序列遗漏地震检测及发震构造分析

针对2014年8月—2015年1月安徽金寨发生的M_L3.9震群,利用匹配滤波技术补充台网目录遗漏的地震事件,再利用波形互相关震相检测技术标定P波和S波到时,进一步采用双差定位方法对震群进行重定位,结合震源机制解等分析此次震群活动可能的发震构造。计算结果显示,通过互相关扫描检测到1376个地震台网常规分析遗漏的地震,数量约为台网目录给出的585个事件的2.35倍。检测到的遗漏地震震级估算为M_L0～2.3,通过震级-频次统计分析,加入遗漏地震后地震目录的完整性在M_L0～1.5范围内有较明显的改善。重定位后地震走时残差更小,水平位置更集中,沿NNE向断裂F和NW向青山-晓天断裂呈现近直立的条带状分布。结合地质构造、震源机制解和水库因素,推测2014年金寨M_L3.9震群可能是由周边水库水下渗引起NW向青山-晓天断裂与NNE向断裂F慢滑动而触发的。     

Statistical Detection and Characterization of a Deviation from the Gutenberg-Richter Distribution above Magnitude 8

— We present a quantitative statistical test for the presence of a crossover c₀ in the Gutenberg-Richter distribution of earthquake seismic moments, separating the usual power-law regime for seismic moments less than c₀ from another faster decaying regime beyond c₀. Our method is based on the transformation of the ordered sample of seismic moments into a series with uniform distribution under condition of no crossover. A simulation method allows us to estimate the statistical significance of the null hypothesis H₀ of an absence of crossover (c₀=infinity). When H₀ is rejected, we estimate the crossover c₀ using two different competing models for the second regime beyond c₀ and the simulation method. For the catalog obtained by aggregating 14 subduction zones of the Circum-Pacific Seismic Belt, our estimate of the crossover point is log(c₀)=28.14 ± 0.40 (c₀ in dyne-cm), corresponding to a crossover magnitude m_W=8.1 ± 0.3. For separate subduction zones, the corresponding estimates are substantially more uncertain, so that the null hypothesis of an identical crossover for all subduction zones cannot be rejected. Such a large value of the crossover magnitude makes it difficult to associate it directly with a seismogenic thickness as proposed by many different authors. Our measure of c₀ may substantiate the concept that the localization of strong shear deformation could propagate significantly in the lower crust and upper mantle, thus increasing the effective size beyond which one should expect a change of regime.