接收函数复谱比的最大或然性估计及非线性反演

利用Ｓｈａｍｗａｙ和Ｄｅｒ的多道最大或然性反褶积原理,讨论了从单台三分量远震Ｐ波波形中分离接收函数径向与垂向分量复谱比的方法．根据Ｔａｒａｎｔｏｌａ的波形反演理论发展了接收函数复谱比的非线性反演方法．本文的反演方法除了拟合接收函数的复谱比之外,还需拟合时间域中接收函数垂向与径向分量的初至振幅比．合成波形数据反演的结果表明本文方法的反演结果不依赖初始模型．利用该方法和ＣＤＳＮ台网兰州台记录的宽频带远震Ｐ波波形数据,研究了该台站的接收函数及其随方位角的变化．接收函数复谱比非线性反演给出了该台站下方１４０ｋｍ深度岩石层的Ｓ波速度结构,得到兰州台下方地壳厚度为５４ｋｓ,岩石层厚度为９４ｋｍ．     

The Kapoeta howardite: Implications for the regolith evolution of the howardite-eucrite-diogenite parent body

Abstract— A large hand sample and numerous polished thin sections, made from the hand sample, of the Kapoeta howardite and its many diverse lithic clasts were studied in detail by optical microscopy and electron microprobe techniques in an attempt to understand the surface processes that operated on the howardite-eucrite-diogenite (HED) parent body (most likely the asteroid 4 Vesta). Four unique, unusually large clasts, designated A (mafic breccia), B (granoblastic eucrite), D (howardite) and H (melt-coated breccia), were selected for detailed study (modal analysis, mineral microprobe analysis, and noble gas measurements). Petrographic studies reveal that Kapoeta consists of a fine-grained matrix made mostly of minute pyroxene and plagioclase fragments, into which are embedded numerous different lithic and mineral clasts of highly variable sizes. The lithic clasts include pyroxene-plagioclase (eucrite), orthopyroxenite (diogenite), howardite, impact-melt, metal-sulfide-rich, and carbonaceous chondrite clasts. The howardite clasts include examples of lithic clasts that constitute breccias-within-breccias, suggesting that at least two regolith generations are represented in the Kapoeta sample we studied. The clast assemblage suggests that repeated shock lithification was an important process during regolith evolution. Noble gas analyses of clast samples fall into two populations: (a) solar-gas-rich clasts H (rim only) and D and (b) clasts A and B, which are essentially free of solar gases. The concentrations of solar noble gases in the two matrix samples differ by a factor of ～40. It appears that clast D is a true regolith breccia within the Kapoeta howardite (breccia-within-breccia), while clast H is a regolith breccia that has been significantly impact reworked. Our data indicate that the Kapoeta howardite is an extraordinarily heterogeneous rock in modal mineral and lithic clast abundances, grain size distributions, solar-wind noble gas concentrations and cosmic-ray exposure ages. These results illustrate the repetitive nature of impact comminution and lithification in the regolith of the HED parent body.     

Mid-crustal Metasomatic Reaction Veins ina Spinel Peridotite

In Central Dronning Maud Land, East Antarctica, rare metre-sizedlenses of spinel peridotite are enclosed in high-grade metamorphicrocks. The rocks experienced a medium-P granulite-facies metamorphismat 575 Ma and a low-P amphibolite-facies overprint at 530 Ma.The latter is probably related to extensive granitoid magmatismbetween 530 and 500 Ma, which produced large volumes (abouthalf of the outcrops today) of granitic to syenitic rocks aswell as abundant K-feldspar–quartz pegmatites. One ofthe spinel peridotite lenses in the Schirmacher Oasis of CentralDronning Maud Land is crosscut by several small (up to 10 cmwide) veins with a characteristic zoned sequence of mineralassemblages, which was formed by reaction of a hydrous, SiO₂-saturatedfluid or pegmatitic melt with the peridotite. The zoned sequenceconsists of the following mineral assemblages (from the centreof the vein towards the outer margin): zone 0, plagioclase +quartz; zone 1, green biotite intergrown with zircon + clinoamphibole;zone 2, cummingtonite + dark brown biotite intergrown with rutile+ clinoamphibole; zone 3, cummingtonite + light brown biotite+ spinel; zone 4, olivine + orthopyroxene + spinel ±clinopyroxene (unaltered peridotite). This sequence was investigatedwith respect to its conditions of formation, modal mineralogy,mineral chemistry, fluid inclusions, and oxygen and hydrogenisotope compositions of selected minerals. Based on the stabilityof cummingtonite and on equilibrium calculations in the MgO–SiO₂–H₂Osystem and on quartz–biotite oxygen isotope thermometry,the reaction vein formed at 650°C, which is in accord withtypical pegmatite crystallization temperatures. The pegmatiteof zone 0 is interpreted to have formed in an open fissure whereas,on textural grounds, zone 3 replaces former peridotite. On thebasis of mass balance constraints, the boundary between zones1 and 2 is interpreted to approximately represent the formerboundary between peridotite and the open fissure before reaction.Oxygen isotope systematics show that the infiltrating fluidhad an isotopic composition of 9–10 SMOW. All mineralsof the reaction vein with the exception of the inherited spineland olivine in the adjacent peridotite are in equilibrium withsuch a fluid. Spinel in the peridotite is depleted in ¹⁸O comparedwith coexisting olivine, which suggests isotopic disequilibrium.Spinel in zone 3 has a distinctly different isotopic compositioncompared with that in the peridotite, apparently approachingbut not reaching equilibrium. The combination of mineral chemistryand mass balance constraints of the modal mineralogy constrainsthe volume change during metasomatism and the direction of elementaldiffusion. It is indicated that Mg, Cr and Ni always diffusedtowards the vein, whereas Si, Al, K, Na, H₂O and possibly Fediffused into the peridotite. KEY WORDS: peridotite; metasomatism; pegmatite; diffusion; reaction     

THE DETERMINATION OF THE SEISMIC QUALITY FACTOR Q FROM VSP DATA: A COMPARISON OF DIFFERENT COMPUTATIONAL METHODS1

Ten methods for the computation of attenuation have been investigated, namely: amplitude decay, analytical signal, wavelet modelling, phase modelling, frequency modelling, rise-time, pulse amplitude, matching technique, spectral modelling and spectral ratio. In particular, we have studied the reliability of each of these methods in estimating correct values of Q using three synthetic VSP seismograms for plane P-waves with different noise contents. The investigations proved that no single method is generally superior. Rather, some methods are more suitable than others in specific situations depending on recording, noise or geology. The analytical signal method has been demonstrated to be superior if true amplitude recordings are available. Otherwise spectral modelling or, in the ‘ noise-free’ case the spectral ratio method, is optimal. Finally, two field VSPs in sediments are investigated. Only in the case of the highest quality VSP can significant information be deduced from the computed attenuation.     

中国数字地震台网的接收函数及其非线性反演

利用中国数字台网（CDSN）记录的85个远震事件的宽频带P波波形数据和分离接收函数的最大或然性反褶方法,获得了CDSN台网10个台站不同方位的岩石层接收函数.利用这些台站不同方位的平均接收函数和非线性接收函数反演方法,获得了上述各台站下方100km深度范围内的岩石层S波速度结构.结果表明CDSN台网各台站下方的地壳厚度和岩石层速度结构存在明显的差别.其中,拉萨台下方的地壳厚度为66kin,壳幔界面较模糊,而余山台下方的地壳厚度为34km,壳慢界面两侧速度反差明显.刮用本文方法估计的地壳厚度与已有的结果基本一致,由于CDSN台网覆盖了中国大陆的各主要构造单元,本文的结果为研究中国大陆的岩石层S波速度结构及其横向变化提供了新的证据.