Energy balance of simple elastodynamic sources

Complete relations are derived for energy and energy flux of elastic waves generated by an isotropic and double-couple source in a perfectly elastic, homogeneous, isotropic, and unbounded medium. In the energy balance of elastodynamic sources near-field waves play an essential role, transforming static energy into wave energy, andvice versa. For explosive and dislocation sources, the sources surface radiates a positive wave energy that is partially distributed to the medium transforming into static energy. For implosive and antidislocation sources, the source surface generates elastic waves, but it does not necessarily imply that it also radiates a positive wave energy. The energy transported by waves can originate in gradual transformation of the static-to-wave energy during propagation of waves through a stressed medium.On leave from Geophysical Institute, Czech Academy of Sciences, Boní II/1401, 41 31, Praha 4 Czech Republic     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.     

Cross-formational hydrocarbon fluid flows in the Tertiary deltaic system of the Beaufort–Mackenzie Basin

Molecular and isotopic compositions of crude oils in the Beaufort–Mackenzie Basin confirm three genetic end-member oil groups and suggest extensive cross-formational hydrocarbon fluid flows in the Tertiary deltaic system. Inter- and intra-fractional variations in the geochemistry of the Tertiary-reservoired oils indicate that the oil source/maturity signatures were substantially masked by biomarkers that were picked up along migration pathways. Thus, many of the previously recognized “immature non-marine oils” are in fact thermally mature, probably derived from unpenetrated deeper marine source rocks. Although the effective source rock volumes have not been evaluated and their exact stratigraphic levels remain unknown, the relative timing of oil generation versus trap formation, rather than poor source quality, may be the cause of under-filled traps in the offshore area.     

Biostratigraphy of Jurassic-Cretaceous boundary sediments in the Eastern Crimea

The Dvuyakornaya Formation section in the eastern Crimea is described and subdivided into biostratigraphic units based on ammonites, foraminifers, and ostracodes. The lower part of the formation contains first discovered ammonites of the upper Kimmeridgian (Lingulaticears cf. procurvum (Ziegler), Pseudowaagenia gemmellariana Olóriz, Euvirgalithacoceras cf. tantalus (Herbich), Subplanites sp.) and Tithonian (?(Lingulaticeras efimovi (Rogov), Phylloceras consaguineum Gemmellaro, Oloriziceras cf. schneidi Tavera, and Paraulacosphinctes cf. transitorius (Oppel)). Based on the assemblage of characteristic ammonite species, the upper part of the formation is attributed to the Berriasian Jacobi Zone. Five biostratigraphic units (zones and beds with fauna) distinguished based on foraminifers are the Epistomina ventriosa-Melathrokerion eospirialis Beds and Anchispirocyclina lusitanica-Melathrokerion spirialis Zone in the upper Kimmeridgian-Tithonian, the Protopeneroplis ultragranulatus-Siphoninella antiqua, Frondicularia cuspidiata-Saracenaria inflanta zones, and Textularia crimica Beds in the Berriasian. The Cyrherelloidea tortuosa-Palaeocytheridea grossi Beds of the Upper Jurassic and Raymoorea peculiaris-Eucytherura ardescae-Protocythere revili Beds of the Berriasian are defined based on ostracodes. A new biostratigraphic scale is proposed for the upper Kimmeridgian-Berriasian of the eastern Crimea. The Dvyyakornaya Formation sediments are considered as deepwater facies accumulated on the continental slope.     

Correlations in aftershock and seismicity patterns

Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

城市化过程中土地征用风险表征及对策

从土地征用基本概念和基本特征着手，描述了土地征用的主要风险，其表现为土地征用权力滥用、侵害失地农民的权益、诱发集体土地黑市、诱发社会矛盾等方面。分析了产生风险的“公共利益需要”缺乏明确界定、土地征用的补偿标准偏低或滞后、土地征用及相关制度不健全或存在缺陷等主要原因。提出了控制与规避风险的对策与措施：严格限定公益性用地范围，提高征地补偿标准，扩大补偿范围，完善相关制度，严格土地征用程序。     

Total gaseous mercury emissions from mercury-enriched soil in Guizhou, China

Guizhou is located in the Circum-Pacific Global Mercuriferous Belt, and mercury concentrations in soil in this area are enriched. In-situ total gaseous mercury （TGM） exchange fluxes between air and soil surface were intensively measured at four sampling sites in Guiyang from 21 May to 16 June, 2003, and five sites in the Lanmuchang mercury mining area in December 2002 and May 2003, respectively. The in-situ Hg flux measurement was conducted with a dynamic flux chamber （DFC） of quartz. Overall, net emissions were obtained from all sampling sites. Soil mercury concentration and solar radiation have been proved to be the two most important parameters to control mercury emissions from soil. Meanwhile, rain events can enhance mercury emission rate significantly.     

Origin and three-dimensional fluorescence spectrum of the chromophoric dissolved organic matter in Yangtze estuary, East China

Dissolved organic matter （DOM） is an important chemical component in natural water. Chromophoric dissolved organic matter （CDOM）, a fraction of optical properties, plays art important role in the biogeochemical cycle of nutrients in aquatic environment. People realized that DOM cycle is crucial in the global carbon and nitrogen flux, and also is inherently related to nutrients and trace metal elements. Therefore, CDOM was concerned by scientists in global oceanography and limnology fields. Water samples were collected from three sections （North Channel, South Channel and Zhuyuan） of the Yangtze （Changjiang River） estuary in March 2006 Three-dimensional excitation emission matrix （3-DEEM） fluorescence spectra were analyzed for those filtrates through Whatman GF/F filters. Dissolved organic carbon （DOC） was also measured by TOC analyzer. The tidal variety was also taken into account. The 3-D EEM fluorescence scans suggested the fluorescence characteristics of humic acid （Ex=332-344 nm, Em=439-451 nm） and fulvic acid （Ex=250-254 nm, Em=472-478 nm） were obvious, and the fluorescence group of protein-like and tyrosine （Ex=230 nm, Em=283 nm） was also found. They are mainly composed of CDOM in the Yangtze estuary. Further data analysis, especially the fluorescence index （f 450/500）, showed that terrestrial signal was rather strong （1.41-1.65） in the surface water, however, some terrestrial CDOM signals of bottom water showed excursions （1.28-1.39）. On the other hand, anthropogenic sign was impressed in the waters of Zhuyuan, which is one of the main drain outlets of Shanghai Metropolis. DOC concentrations ranged from 2.2 mg/L to 3.4 mg/L in Zhuyuan and South Channel, and from 2.0 mg/L to 2.4 mg/L in North Channel. The tide effect played a role in the composition of the CDOM measured by 3-D fluorescence scan technology.