Processes of accelerated pluvial erosion on desert hillslopes modified by vehicular traffic

Accelerated pluvial erosion on hillslopes modified by off-road vehicles (ORVs) is analysed using results from 50 rainfall simulation experiments conducted in the Mojave Desert, California. Sediment yield from 1 m² hillslope plots subjected to intense, 20-minute rainfalls is typically increased 10 to 20-fold following ORV use. Salient effects of vehicle traffic, which reduce infiltration, increase runoff sediment transport efficiency, and enhance gully formation, are further studied by combining simple theoretical relations with experimental data. This analysis helps identify factors controlling erosion on natural desert hillslopes, as well as those used by ORVs. Erosion of natural or vehicle-used desert surfaces is heavily influenced by runoff hydraulics. Calculated Darcy-Weisbach friction factors decrease by an average of 13-fold following vehicular slope modification, whereas runoff Reynolds numbers increase by an average of 5 1/2-fold. The capacity of overland flow to transport sediment is related to runoff power and its degree of localization, which usually increase considerably following ORV activity; however, the ability of overland flow to move large grains (competency) is related to a combination of factors not always systematically influenced by ORV use. Kinematic runoff routing, which is used to extrapolate experimental results to longer slope lengths, leads to the suggestion that the hydraulic roughness of desert hillslopes strongly influences their erosional behaviour.     

Ground motion scenarios consistent with probabilistic seismic hazard disaggregation analysis. Application to Mainland Portugal

Probabilistic seismic hazard for Mainland Portugal was re-evaluated in order to perform its disaggregation. Seismic hazard was disaggregated considering different spaces of random variables, namely, univariate conditional hazard distributions of M (magnitude), R (source-to-site distance) and ε (deviation of ground motion to the median value predicted by an attenuation model), bivariate conditional hazard distributions of M–R and X–Y (seismic source latitude and longitude) or multivariate conditional hazard distributions of M–R–ε and M–(X–Y)–ε. The main objective of the present work was achieved, as it was possible, based on the modal values of the above mentioned distributions, to characterize the scenarios that dominate some seismic hazard levels of the 278 Mainland Portuguese counties. In addition, results of 4D disaggregation analysis, in M–(X–Y)–ε, pointed out the existence of one geographic location shared by the dominant scenario of most analyzed counties, especially for hazard levels correspondent to high return periods. Those dominant scenarios are located offshore at a distance of approximately 70 km WSW of S. Vicente cape. On the other hand, the lower the return period the higher is the number of modal scenarios in the neighbourhood of the analyzed site. One may conclude that modal scenarios reproduce hazard target values in each site with great accuracy enabling the applications derived from those scenarios (e.g. loss evaluation) to be associated to a hazard level exceedance probability.     

Malleable disc base for direct infiltration measurements using the tension infiltrometry technique

The correct use of the tension disc infiltrometer requires the membrane of the disc base to be completely in contact with the soil surface. To achieve this contact, a thick layer of sand is commonly placed between the soil surface and the disc base. This paper presents an alternative disc (M_DB), which, by incorporating a malleable membrane, allows direct infiltration measurements without using a contact sand layer. Infiltration curves obtained with this new design in a soil under three different tillage management treatments were compared with the corresponding curves obtained with a conventional disc (C_DB) that uses a contact sand layer. The cumulative infiltration curves measured with C_DB were analysed by the differentiated linearization (DL) method, and the corresponding curves obtained with M_DB were analysed using both the DL and the cumulative linearization (CL) models. The values of hydraulic conductivity (K₀) and sorptivity (S₀) estimated with C_DB were also compared with those obtained with M_DB. Finally, the cumulative infiltration curves measured with C_DB and M_DB were compared with the corresponding modelled function for the respective K₀ and S₀ values calculated with the CL and DL models. The results show that, compared with C_DB without a contact sand layer, M_DB allows complete soil surface wetting even when non‐smoothed soil surfaces are used. The C_DB, which yielded average K₀ values 18% lower than those estimated with M_DB, gave the highest values of standard error for the hydraulic parameters calculated. Furthermore, the subjective method employed in the C_DB‐DL technique, which requires the first points of the differential infiltration line corresponding to the sand layer to be manually removed, introduces additional uncertainties in estimating S₀ and K₀. Comparison between the modelled and measured infiltration curves demonstrates that the DL or CL methods applied to M_DB gave excellent estimates of S₀ and K₀. Copyright © 2012 John Wiley & Sons, Ltd.     

M TSU : Recovering Seismic Moments from Tsunameter Records

We define a new magnitude scale, M_TSU, allowing the quantification of the seismic moment M₀ of an earthquake based on recordings of its tsunami in the far field by ocean-bottom pressure sensors (``tsunameters') deployed in ocean basins, far from continental or island shores which are known to affect profoundly and in a nonlinear fashion the amplitude of the tsunami wave. The formula for M_TSU, M_TSU = log₁₀ M₀ − 20 = log₁₀ X (ω) + C_D^TSU + C_S^TSU + C₀, where X (ω) is the spectral amplitude of the tsunami, C_D^TSU a distance correction and C_S^TSU a source correction, is directly adapted from the mantle magnitude M_m introduced for seismic surface waves by Okal and Talandier. Like M_m, its corrections are fully justified theoretically based on the representation of a tsunami wave as a branch of the Earth's normal modes. Even the locking constant C₀, which may depend on the nature of the recording (surface amplitude of the tsunami or overpressure at the ocean floor) and its units, is predicted theoretically. M_TSU combines the power of a theoretically developed algorithm, with the robustness of a magnitude measurement that does not take into account such parameters as focal geometry and exact depth, which may not be available under operational conditions in the framework of tsunami warning. We verify the performance of the concept on simulations of the great 1946 Aleutian tsunami at two virtual gauges, and then apply the algorithm to 24 records of 7 tsunamis at DART tsunameters during the years 1994–2003. We find that M_TSU generally recovers the seismic moment M₀ within 0.2 logarithmic units, even under unfavorable conditions such as excessive focal depth and refraction of the tsunami wave around continental masses. Finally, we apply the algorithm to the JASON satellite trace obtained over the Bay of Bengal during the 2004 Sumatra tsunami, after transforming the trace into a time series through a simple ad hoc procedure. Results are surprisingly good, with most estimates of the moment being over 10²⁹ dyn-cm, and thus identifying the source as an exceptionally large earthquake.     

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system,central California

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes ofM _L =5.3 andM _L =5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average.The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquakeretardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979M _L =5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes ofM _L =5.2 andM _L =6.2, respectively.Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes.Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip.     

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.     

Threshold relations for the transport of sediment by overland flow on desert hillslopes

Although the Shields relation was developed for rivers, it has been applied to sediment transport by overland flow. According to the Shields relation, where the critical boundary Reynolds number Re_*c exceeds 40, the critical Shields number F_*c is independent of both Re_*c and the ratio of the critical flow depth to particle diameter d_c/D. Analyses of data collected from runoff plots in southern Arizona reveal that F_*c is positively correlated with both Re_*c and d_c/D. Thus the Shields relation does not apply to overland flow on debris-covered desert hillslopes. Multiple regression analysis is employed to develop alternative threshold relations in which critical boundary shear stress τ_c is related to D and d_c/D (R² = 0.782) and to D and S_c (critical gradient) (R² = 0.625). The computed R² values derive in large part from the spurious correlations of d_c/D and S_c with τ_c. Nevertheless, the relations may be utilized to predict τ_c. In this regard, the latter relation is likely to prove more useful than the former because S_c is generally known, whereas d_c is not. An investigation of the functional relation between τ_c and D reveals that τ_c varies approximately with D² for overland flow on the desert hillslopes under study, whereas the Shields relation predicts a linear relation (i.e. a D exponent of 1). This result is consistent with Cheng's data which show that F_*c varies with (d_c/D)^?1 where 0.4 < d_c/D < 2 and may be explained in terms of increased energy dissipation both in separation zones downslope of particles and in distortion of the water surface as d_c/D decreases. Consequently, larger values of τ_c, and hence F_*c, are required to initiate the transport of particles of a given size D as d_c decreases.     

青藏高原西北部区域地壳形变、构造地貌与孕震构造模型研究——以2008年与2014年新疆于田7.3级地震为例

为了清晰认识发生于青藏高原西北部2008年与2014年的两次于田MS7.3地震发震构造环境与构造地貌特征,本文利用DEM(数字高程模型)数据分析"喀喇昆仑—西昆仑—康西瓦地区"的地形地貌特征,结合区域活动断裂研究资料、相对于塔里木盆地的两期GPS速度场资料和区域运动学特征等讨论两次MS7.3地震所处的青藏高原西北部区域构造环境和地壳运动学特征,分析喀喇昆仑断裂、阿尔金断裂康西瓦段、龙木错-邦达错断裂及贡嘎错断裂所围限的西昆仑地块的地质构造背景、阿尔金断裂西南端发震断裂活动性及孕震环境等发震构造基本条件;进而利用"地形剖面"方法及断裂分布特征分析震源区的地形地貌特征,给出晚第四纪以来的地貌形态与发震构造的关系,从区域构造地貌学和GPS地壳运动学的角度探讨中上地壳变形特征及孕震过程;最后讨论区域孕震构造、克尔牙张性裂谷演化过程和地球动力学背景等。通过地形剖面及区域地貌综合分析新疆于田2008年MS7.3拉张型发震构造和2014年MS7.3走滑拉张型地震的发震构造特点的区别,认为2014年发生的地震可能与2008年MS7.3地震同震库伦应力变化、触发过程及震后变形过程密切相关,并且青藏高原西北部地区存在明显的东西向拉张性构造单元,可能与青藏高原10~15 Ma以来的地壳减薄过程有关。     

2022年门源MS6.9和2016年门源MS6.4地震序列比较分析

2022年1月8日青海省海北州门源县发生M_S6.9地震,震中距离2016年1月21日门源M_S6.4地震震中约33km,两次门源地震均发生在冷龙岭断裂附近,但在震源机制、主发震断层破裂过程及地震序列余震活动等方面显著不同。针对两次门源地震序列的比较分析,对研究冷龙岭断裂及其附近区域强震序列和余震衰减特征等具有重要研究意义。通过对比分析2022年门源M_S6.9地震和2016年门源M_S6.4地震余震的时空演化特征,发现二者在震源过程和断层破裂尺度上存在明显差异,前者发震断层破裂充分,震后能量释放充分,余震丰富且震级偏高;而后者发震断层未破裂至地表,余震震级水平偏低。综合分析两次门源地震序列表现出来的差异性,认为其可能与地震发震断层的破裂过程密切相关,且同时受到区域构造环境的影响。     

Inversion of long-period body-wave data for the source process of the Gonghe, Qinghai, China earthquake

Long period body waves data recorded by the China Digital Seismograph Network (CDSN) are inverted for the seismic moment tensors of the April 26, 1990, Gonghe, QinghaiM _S=6.9 earthquake and itsM _S=5.0 after-shock occurred on May 7, 1990. In the inversion, the generalized reflection-transmission coefficient matrix method is used to generate Green’s function. From the inversion it is obtained that the rupture process of theM _S=5.0 aftershock is relatively simple, and that of the main shock is rather complex. There are at least two events during main shock rupture process with an interval about 35 seconds. The focal mechanisms of two events are roughly the same as that of the aftershock, all of them were mainly reverse dip-slipping faulting with minor left-lateral strike-slip motion. These results indicate that the Gonghe earthquake was the result of the farther extension of one NWW-SEE striking buried fault on the southern margin of Gonghe basin from shallower depth to deeper depth and from NW to SE under the action of a nearly horizontal NE direction compressive stress. Contribution No. 95A0111, Institute of Geophysics, SSB, China.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Moment tensor inversion and source rupture process of the September 27, 2003 M S=7.9 earthquake occurred in the border area of China, Russia and Mongolia

We conducted moment tensor inversion and studied source rupture process for M _S=7.9 earthquake occurred in the border area of China, Russia and Mongolia on September 27 2003, by using digital teleseismic P-wave seismograms recorded by long-period seismograph stations of the global seismic network. Considering the aftershock distribution and the tectonic settings around the epicentral area, we propose that the M _S=7.9 earthquake occurred on a fault plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of M _S=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M ₀=0.97×10²⁰ N·m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the M _S=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of aftershocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.     

Hydraulic Tomography Offers Improved Imaging of Heterogeneity in Fractured Rocks

Fractured rocks have presented formidable challenges for accurately predicting groundwater flow and contaminant transport. This is mainly due to our difficulty in mapping the fracture‐rock matrix system, their hydraulic properties and connectivity at resolutions that are meaningful for groundwater modeling. Over the last several decades, considerable effort has gone into creating maps of subsurface heterogeneity in hydraulic conductivity (K) and specific storage (S_s) of fractured rocks. Developed methods include kriging, stochastic simulation, stochastic inverse modeling, and hydraulic tomography. In this article, I review the evolution of various heterogeneity mapping approaches and contend that hydraulic tomography, a recently developed aquifer characterization technique for unconsolidated deposits, is also a promising approach in yielding robust maps (or tomograms) of K and S_s heterogeneity for fractured rocks. While hydraulic tomography has recently been shown to be a robust technique, the resolution of the K and S_s tomograms mainly depends on the density of pumping and monitoring locations and the quality of data. The resolution will be improved through the development of new devices for higher density monitoring of pressure responses at discrete intervals in boreholes and potentially through the integration of other data from single‐hole tests, borehole flowmeter profiling, and tracer tests. Other data from temperature and geophysical surveys as well as geological investigations may improve the accuracy of the maps, but more research is needed. Technological advances will undoubtedly lead to more accurate maps. However, more effort should go into evaluating these maps so that one can gain more confidence in their reliability.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.     

A Preliminary Analysis of the Seismogenic Structure of the Akto MS6.7 Earthquake Sequence on November 25, 2016

In this study,data from the Xinjiang regional network and IRIS shared global stations are used to relocate the Akto M_S6. 7 earthquake sequence on November 25,2016 by using double difference location method. Three earthquakes of M_S4. 8,M_S6. 7 and M_S5. 0 are inverted by using the g CAP method,and the focal mechanism solutions are obtained.According to the results of relocating,the location of the main shock is 39. 22°N,73. 98°E,the distribution of the earthquake sequence is about 70 km in length,and the focal depth is mainly within the range of 5-20 km. The plane and depth profiles of the earthquake sequence show that aftershocks extended in SEE direction after the main shock and the dip angle of fault plane is steep. Focal mechanism results show that the three earthquakes are characterized by strike-slip movement. Based on the results of field geological investigation,it is inferred that the seismogenic fault of the Akto earthquake is Muji fault,which is located at the northernmost end of the Kongur extensional system.The possible cause of this earthquake is that the Indian Plate continues to push northward,and during this compression process,the Indian Plate is affected by the clockwise rotation of the Tarim basin,which causes the accumulation of right-lateral action of the Muji fault,resulting in this earthquake.     

Study of the Change in Wave Velocity Ratio before and after Two Strong Earthquakes Using Airgun Source Data

Using the signals excited by the large-volume airgun source at the Binchuan transmitting seismic station from January to June, 2016, arrival-time data was acquired at four stations near the epicenter of the Eryuan M_S4.5 and M_S4.0 earthquakes on February 8, 2016, as well as the epicenter of the Yunlong M_S5.0 and Eryuan M_S4.6 earthquakes on May 18, 2016 through the waveform cross-correlation technique. The wave velocity ratio of the four stations was calculated using the single-station method. At the same time, the b-value and the focal mechanism consistency parameters of the study area were also calculated. The results show that:(1) the wave velocity ratio of each station experienced a process of decline-recovery-fast rise before the two strong earthquakes, and a significant quasi-synchronous rapid rise occurred within 3-12 days before the earthquake; (2) the timing of the rapid rise of the wave velocity ratio of the four stations before the Yunlong M_S5.0 and Eryuan M_S4.6 earthquakes were related to the epicentral distance. The station which observed the earliest increase in rapid rise is the farthest one from the epicenter, and the station where the rapid rise appeared in the latest is closest to the epicenter; (3) the form of change of the wave velocity ratio before the earthquake was different between stations located at different directions in the epicenter area; (4) the b-value and the focal mechanism consistency parameter which is commonly used to characterize the stress level both showed a downward trend before the two strong earthquakes, and were consistent with the change in the wave velocity ratio. According to the preliminary analysis, the wave velocity ratio obtained by using airgun source can better reflect the change in the stress state of the underground medium.     

新10井水位对九寨沟M_S7.0、精河M_S6.6地震同震响应

基于新10井水位对九寨沟M_S7.0、精河M_S6.6地震的同震响应特征,分别利用水震波和潮汐分析法,反演得到了新10井水位在地震波作用期间和地震波作用前后含水层参数特征,探讨了2次地震引起的新10井同震响应变化机理。结果表明,新10井水位对九寨沟地震和精河地震具有不同的同震响应形态,利用水平流模型反演所得的新10井含水层渗透系数,在九寨沟地震时为61m/d,而在精河地震时为147m/d,表明地震波作用导致水平向的渗透性增强;利用垂向流模型反演得到的新10井含水层渗透系数,在2次地震之前约为49×10~(-5)m/d,而在地震之后约为18×10~(-5)m/d,表明地震波作用导致垂直向的渗透性减弱。因此,新10井水位的同震响应变化与地震波引起的含水层渗透率的改变有关。     

2016年安徽无为MS3.0震群序列发震构造与震后趋势研究

首先通过模板匹配方法检测无为震群活动期间目录遗漏的地震事件,共识别出5次遗漏地震事件,震级为M_L0.5～1.2,得到了更为完整的地震目录;然后基于波形互相关震相检测技术标定震相到时,进而采用双差定位方法进行精定位,精定位后震群分布更加集中,未见明显的优势方位分布;采用Snoke方法计算震级较大地震的震源机制,结果表明,此次震群为NEE向的水平挤压与NNE向的水平拉张应力场作用下具逆冲分量的走滑型地震活动,严家桥-枫沙湖断裂可能为其发震构造;最后,计算了震群序列的视应力,结果显示,视应力和扣除震级影响后的差视应力随着震群序列的衰减逐渐恢复,因此,分析认为,随着无为震群序列的衰减,震源区发生更大地震的可能性不大。