Tomography for Q of the Eastern Section of the Tianshan Area from High-frequency Attenuation of S-wave

Based on the waveform data of 5,076 local earthquakes recorded at 25 stations in Xinjiang during the period from 2009 to 2014 and the observation reports provided by the Xinjiang Digital Seismic Network,a data set of 19,140 attenuation factors t*is obtained by fitting the high-frequency attenuation of S-wave spectra with a genetic algorithm. The spatial distribution of Q_S is determined by inverting the t*data with seismic tomography. The results show that the average Q0 in eastern Tianshan is 520,and there is a significant correlation between the Q_S value distribution or attenuation characteristics it disclosed and the surface structure of the study area. The Q_S value is lower in the intersection area of the mountain basin which is located on the north and south sides of the Tianshan Mountains,and the high Q_S distribution is more concentrated inside the Tianshan orogenic belt. The M≥6. 0 earthquakes have been basically located in the Low-Q_S region since 1900. 24 high heat flow points in eastern Tianshan are located at the north and south of Tianshan Mountains where low Q_S exists,indicating a negative correlation. In addition,there is a positive correlation between the velocity structure and the attenuation structure in the study area,which reflects the consistency of the 2-D attenuation structure with the velocity structure and the two-dimensional density structure.     

Seismic Ground Motion Zoning Maps of the Pangxi Region

The seismotectonic environment and seismic activity in Southwest China region were studied based on new data and new results obtained during the Eighth and Ninth Five-Year Plans,the seismic areas and zones and potential seismic source zones were determined.and the relation between seismic activity parameters and ground motion attenuation was determined.Finally the seismic gound motion zoning maps of the Pangxi region was compiled by using the multi-parameter and multi-scheme method.     

S wave inelastic attenuation and site effect in the northern Tianshan area

In this paper,according to the relationship of observation spectrum with source spectrum,travel-path attenuationfactor and site effect,we use the digital seismic data of moderate-small earthquakes to invert the soft rock site ef-fect and S wave inelastic attenuation under the constraint of site effect,and determine the parameters of path andsite.The parameters obtained by this method exclude the topographic effect of seismic stations,and truly reflectthe seismic effect of free soft rock sites and the S wave inelastic attenuation in the studied region.The average softrock site effect is about 1.5 times in the frequency domain of 2～4 Hz and the S wave quality factor is obtained asQS=278 f0.346,which can be directly applied to the stochastic modeling of ground motions in the studied region.     

Macroseismic intensity attenuation in Iran

Macroseismic intensity data plays an important role in the process of seismic hazard analysis as well in developing of reliable earthquake loss models. This paper presents a physical-based model to predict macroseismic intensity attenuation based on 560 intensity data obtained in Iran in the time period 1975–2013. The geometric spreading and energy absorption of seismic waves have been considered in the proposed model. The proposed easy to implement relation describes the intensity simply as a function of moment magnitude, source to site distance and focal depth. The prediction capability of the proposed model is assessed by means of residuals analysis. Prediction results have been compared with those of other intensity prediction models for Italy, Turkey, Iran and central Asia. The results indicate the higher attenuation rate for the study area in distances less than 70 km.     

Study on the Attenuation Features of Earthquake Intensity in the Yunnan Region

This research on the attenuation features of intensity in the Yunnan region has been conducted by using intensity isoseismals of 127 earthquakes, and the following three conclusions have been put forward : ( 1 ) The average ratio value of long axis to short axis of innermost isoseismals in the Yunnan region is larger than that of eastern China and smaller than that of western China.The velocity of intensity attenuation varies in different directions; (2) The zoning feature of intensity attenuation in the Yunnan region is obvious; (3) In Northwest China and North China, the ratios of long axis to short axis of low-magnitude isoseismals are larger than those of high magnitude ones. Compared with that phenomenon, the ratios in the Yunnan region are more complicated. Finally, some relevant questions are discussed.     

Study on Ground Motion Attenuation Relation in Shanghai and Its Adjacent Region

Based on intensity data in Shanghai and its adjacent region, the intensity attenuation relation is determined. Selecting the western United States as a reference area where there are rich strong ground motion records and intensity data, and by determining ground motion attenuation relation in an area lacking in strong ground motion data, we obtain the ground motion attenuation relation in Shanghai and its adjacent region.     

Attenuation relations for horizontal peak ground acceleration and response spectrum in northeastern Tibetan Plateau region

Introduction Developing local attenuation relations of ground motion is one of the key steps in seismic hazard assessment. Because of inadequate strong ground motion records in China, the attenuation relations used in China are usually developed by using the transforming method (Hu, Zhang, 1983; HU, ZHANG, 1984). To use this method, we need to have both the attenuation relation of seismic intensity for the studied region and the attenuation relations of seismic intensity and ground mo-tion…     

Comparative study on the deviation characteristics of post-earthquake intensity evaluation based on two intensity attenuation relationship models

We collect the isoseismal data of 488 earthquakes in western China and 182 earthquakes in eastern China after 1900. The least square method is used to establish two models of the seismic intensity attenuation relationship partitions in China, based on the major and minor axis ellipse models, and the area and the major axis radius,respectively. The two models are applied to the calculation of the intensity circle of the earthquake events with a magnitude M_S≥5.0 from 2008 to 2019, and the actual intensity is compared with the model intensity value as an indicator to verify the consistency between the actual intensity and the value from the empirical statistical model.Three real earthquake results are selected to calculate the major and minor axis and area of the ellipse using the two intensity attenuation relationship models. After comparison, we summarize the deviation characteristics of the intensity value, and put forward corresponding improvement suggestions.     

The Study on the Strong Ground Motion Attenuation Relationship in the Pishan Area, Xinjiang Uygur Autonomous Region

Small earthquake data from the Pishan MS6.5 aftershocks is collected by the Xinjiang Regional Digital Seismic Observation Network.Five parameters of the focal region are obtained by micro genetic inversion:stress dropΔσof 75.95 bars,quality factor parameters Q0of 186.33 andηof 0.26,geometric attenuation parameters R1of 72.18km and R2of 139.70km.We calculate the Fourier spectrum and combine it with the random phase spectrum to get the ground motion time history,and build the strong motion acceleration attenuation relationship.The strong ground motion acceleration attenuation of the Pishan area is thus obtained.Because of the insufficiency of strong ground motion records,we added the records from the Wuqia MS6.9 earthquake on October 5,2008,the Akto MS6.2 earthquake on October 6,2008,and the Lop MS6.0 earthquake on March 9,2012 to the data.The comparison of the calculation results and the empirical attenuation relationships with strong ground motion records reveal that the strong motion data of Pishan and Lop earthquakes is higher than the empirical attenuation relationships.The Wuqia MS6.9 earthquake strong motion data is consistent with Yu Yanxiangs(2013)short axis result,and lower than the present result.     

Temporal and Spatial Images of Ambient Noise Intensity in the Chinese Mainland

Successive waveforms of the vertical component recorded by 888 broadband seismic stations in the China Seismography Network from January,2010 to June,2011 are used to investigate the temporal and spatial distribution of ambient noise intensity,and the images of ambient noise intensity at the period of 10 s in the Chinese Mainland are obtained. The temporal variation of ambient noise intensity shows some seasonal and periodic characteristics. The maximum ambient noise intensity occurred from January,2011 to March,2011. The spatial distribution images of ambient noise intensity show obvious zoning features,which doesnt correlate with surface geology,suggesting that the noise field is stronger than the site factors. The strength in southeastern coastal areas reaches its maximum and generally decreases toward to inland areas,and arrives at the minimum in the Qinghai-Tibetan Plateau. The zonal intensity distribution is probably correlated with ocean tides from the Philippine Ocean and the Pacific Ocean. It also shows that the influence from the Indian Ocean seems small. However, the ambient noise intensity increases to a certain degree in the Xinjiang area,indicating that the main source of ambient noise in the western area of the Chinese Mainland is not derived from the East and South China Sea,but rather from the deep interior of the Eurasian continent. The ambient noise intensity obtained in this study can supply reference for seismology research based on ambient noise correlation. Moreover,it can supply basic data for attenuation research based on ambient noise, and thus help achieve the object of retrieving the attenuation of Rayleigh waves from ambient noise.     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

Solute concentrations in the lower niger river and the source rock contribution

Water quality analyses for the Niger River for the 1980/81 hydrological year are presented. The samples were collected from the main river at Lokoja, and from two main tributaries, the Kaduna and the Benue Rivers. Different water types were distinguished by the concentrations of major ions. The type Ca > Na > Mg > K - HCO₃ > SO₄ > Cl was represented at all stations during at least part of the year. Chloride was found to dominate the sulphate ion in the Kaduna and Niger, while the Benue maintained a higher concentration of sulphate relative to chloride all year round. Distinct patterns of seasonal variation in the ion concentrations were observed, particularly for the samples collected at Lokoja. Low ion concentrations were prominent during periods of high discharge, while low flow periods coincided with high dissolved ion concentrations. The contribution of rainwater to the total dissolved solids in the river waters was assessed indirectly using rainwater chemistry data from the Gulf of Guinea. The estimated rainwater contribution to the Lower Niger amounts to 5.15 mg 1^?1. Geochemical weathering calculations involving reactions of the four major minerals of granitic rocks - anorthite, biotite, albite, and K-feldspar - with carbon dioxide and water, can account for the average water composition of the Lower Niger. The proportion of the ionic components was also related to the occurrence of the respective element in the minerals.     

Scenario predictions for potential near-field and far-field earthquakes affecting Hong Kong

Historical earthquakes noted in the written records of the South China region, including Hong Kong, are not well delineated along identified prominent fault sources. Despite the lack of any definitive, localised trend in the spatial distribution of seismic activity in the region, there does appear to be some major disparity in the seismic activity rates (especially for large magnitude earthquakes) between the near-field and the far-field regions of Hong Kong. Despite this observation, previous studies of the regional seismic activity and seismic ground motion hazard (the latter using a probabilistic seismic hazard assessment, PSHA) have considered very broad source zone regions, in which uniform levels of seismic activity have been assumed. The present paper further scrutinises this broad source zone (BSZ) approach by adopting a novel expanding circular disc (ECD) method to determine the rates of earthquake recurrence. Such a method is intended to counter-check previously developed models by determining earthquake scenario events in terms of magnitude–distance (M–R) pairs or combinations, having defined values of average return period. Unlike the BSZ approach, the ECD method specifically accounts for the supposed variations in the seismic activity rates between events in the near-field and the far-field of Hong Kong.The form of the developed method is particularly suited to the determination of design-level earthquake ground motions for bedrock sites, since it assumes a directionally-independent attenuation model as described in the companion paper. It is found that, whilst the BSZ approach may indicate the overall average levels of hazard that are representative of the South China region as a whole, it does not capture the large disparity in seismic activity rates between near-field and far-field events. This important feature is expected to have a significant impact upon engineering assessments of the seismic safety of structures in Hong Kong and elsewhere in the South China region. For example, it is found that for events with M≥6, the seismic activity rate (normalised by time and area) in the very far-field is around 3.5 times larger than in the near-field and medium-field of Hong Kong. The resulting design M-R combinations, covering a range of return periods from 70 to 2500 years, are limited, for very long return periods and for distant events, by the maximum credible earthquake (MCE) magnitude. Intensive research to determine this seismic hazard parameter is recommended, in order to refine further the results of the ECD analysis, which presently conservatively assumes the MCE to range between M=6 in the near-field of Hong Kong to M=8 in the very far-field, at distances greater than 280 km from Hong Kong.     

Depositional ages and characteristics of Middle–Upper Jurassic and Lower Cretaceous lacustrine deposits in southeastern Mongolia

Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Late Permian post‐ophiolitic trondhjemites from Central Iran: a mark of subduction role in growth of Paleozoic continental crust

Late Permian trondhjemites in the Anarak area occur as stocks and dykes, which cross cut the Anarak ophiolite and its overlying metasedimentary rocks, and are exposed along the northern Anarak east–west main faults. These leucocratic intrusive bodies have enclaves of all ophiolitic units and metamorphic rocks. They are composed of amphibole, plagioclase (oligoclase), quartz, zircon and muscovite. Secondary minerals are chlorite (pycnochlorite), epidote, albite, magnetite and calcite. Whole‐rock major‐ and trace‐element analyses reveal that they are characterized by high SiO ₂ (67.8–71.0 wt%), Al ₂ O ₃ (14.9–17.1 wt%) and Na ₂ O (5.3–8.6 wt%), low K ₂ O (0.1–1.5 wt%; average: 0.8 wt%), low Rb/Sr ratio (0.01–0.40; average: 0.09), low Y (3–6 ppm), negative Ti, Nb and Ta anomalies, slightly negative or positive Eu anomaly, LREE enrichment and fractionated HREE. These rocks present 2 to 40 times enrichment in inclined chondrite‐normalized REE patterns. Geochemical characteristics of the Anarak trondhjemites all reflect melting of a mafic protolith at more than 10 kbar. The field evidence and whole‐rock chemistry reveal that these rocks have been crystallized from magmas derived from melting of subducted Anarak oceanic crust. This study reveals that melting of garnet amphibolite was an important element of continent formation in the study area.