Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

邹平火山岩盆地铜矿地质成因地球物理特征探讨

邹平铜矿处于齐河-广饶深大断裂带南部的邹平火山岩盆地中,形成于破火山口火山通道充填的石英正长闪长岩岩颈中央上部,包括伟晶岩型铜矿和细脉浸染状斑铜矿床两种类型。前者矿体较小,但品位高;后者品位较低,但规模中等。含矿石英正长闪长岩等密度小、磁性弱,故在火山岩系中呈现高背景重力场上的重力低和杂乱高磁场背景中的低负异常,即“重磁同低”,且高极化。重磁同低异常区和高极化率异常带,是本区寻找铜矿的有利部位。     

The Glueckstadt Graben of the North-German Basin: new insights into the structure from 3D and 2D gravity analyses

The structure of the Glueckstadt Graben has been investigated by use of 3D gravity backstripping technique and by 2D gravity and magnetic modelling. Subtracting the gravity effects of the Meso-Cenozoic sediments together with Permian salt reveals a positive residual anomaly within the Glueckstadt Graben. This anomaly includes two local maxima over the Westholstein and Eastholstein Troughs. The 2D gravity models point to the presence of a high-density body within the lower crust of the Glueckstadt Graben. In addition, the results of 2D magnetic modelling indicate that the central part of the high-density body is overlain by an area with high susceptibility. Most probable, the formation of this high-density body is a result of complex poly-phase tectonic history of the study area. Finally, the results of gravity modelling indicate that Permian salt is not homogeneous. 3D gravity analysis and, especially, 2D gravity modelling have distinguished the differences in degree of salt saturation in salt-rich bodies, and elucidate the proportion of Rotliegend salt.     

New perspectives on Mars’ crustal magnetic field

The planet Mars lacks, today, a planetary, dynamic magnetic field, but strong, intense, localized magnetic fields of lithospheric origin, one to two orders of magnitude larger than the terrestrial lithospheric field, are present. This lithospheric magnetic field is the result of magnetization processes in the presence of a magnetic dynamo and of demagnetization processes after the dynamo shutdown, such as impact or volcanoes. This crude scenario can be more accurately specified by interpreting global and local models of the current magnetic field of Mars. Some specific areas are studied, including the intensely magnetized Terra Sirenum, as well as the magnetic anomaly associated with Apollinaris Patera. Magnetic minerals could be of primary and/or secondary origin; this latter would imply an early hydration of a basaltic crust. A scenario, in which Mars experienced a major polar wander due to the Tharsis bulge, prior to the cessation of its dynamo, is proposed and discussed.     

地面核磁共振模型约束反演含水层参数

导电模型的地面核磁共振感应电动势是含水量的非线性函数.引入模型约束的迭代反演方法求解该非线性问题的反演.在反演过程中,根据理论公式计算迭代过程中灵敏度矩阵,并采用平坦模型和光滑模型两种约束.对均匀半空间、层状导电模型和实际数据进行了反演模拟,结果表明,模型约束迭代反演方法能从地面核磁共振感应电动势获得含水层较为合理的含水量及分布,且该结果可以从作为初始模型的均匀含水量分布反演得到.对无噪音数据,平坦模型和光滑模型约束对反演的含水量分布影响不大;但当数据存在一定的噪音时,平坦模型约束将比光滑模型约束获得更为精确的含水层参数.     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

勘查油气田化探新方法——土壤相态烃测量法

烃气测量是油气化探最主要的方法。目前应用的烃气测量方法有两类:一类是测量土壤中游离烃气如K—V指纹法、吸附丝法、物上气法及直接取土壤气测量法等。此类方法由于需特殊设备及壤气测量方法本身固有的缺陷,目前还未得到普遍应用,另一类为土壤吸附烃气法即国内外目前普遍使用的酸解烃气法、此法在多数油气田上取得好效果,但也有相当多的地区如黄土、荒漠及水网发育等特殊景观区效果不好或不显著。本文介绍一种新的烃气测量法──土壤相态烃测量法。这是用一种特殊的技术,只选择性提取由油气藏渗漏过程中保存在运移通道围岩及近地表土壤各种次生矿物结构中以"穴居"态存在的轻烃组分。"穴居"态烃为累积性的指标,不受外界温度、修水等的影响,保存较牢固,测定结果稳定。相态烃法避免了酸解烃法存在的生物成因烃、矿物同生烃及碳酸盐的干扰影响,提高了测量指标的信噪比,突出了应有的异常。本方法已在黄土厚覆盖区两个已知油气田和一个风成沙覆盖的断裂发育油田上试验获得成功,结果表明相态烃法优于酸解烃法。     

Strain analysis in Jurassic argillites of the Monte Sirino area (Lagonegro Zone,southern Apennines,Italy) and implications for deformation paths in pelitic rocks

Analysis of strain in Jurassic argillites forming part of the folded and thrusted sedimentary succession of the Lagonegro basin (southern Italian Apennines) has been carried out using ellipsoid-shaped reduction spots as strain markers. Most of the determined finite strain ellipsoids are of oblate type and show a peculiar distribution of the maximum extension direction (X), with maxima either subparallel or subperpendicular to the local fold axes. Using the strain matrix method, two different deformation histories have been considered to assist the interpretation of the observed finite strain pattern. A first deformation history involved vertical compaction followed by horizontal shortening (occurring by a combination of true tectonic strain and volume loss), whereby all strain is coaxial and there is no change in the intermediate axis of the strain ellipsoid. By this type of deformation sequence, which produces a deformation path where total strain moves from the oblate to the prolate strain field and back to the oblate field, prolate strain ellipsoids can be generated and may be recorded where tectonic deformation has not been large enough to reverse pretectonic compaction. This type of deformation history may be of local importance within the study area (i.e. it may characterize some fold hinge regions) and, more generally, is probably of limited occurrence in deformed pelitic rocks. A second deformation sequence considered the superposition of pre-tectonic compaction and tectonic strain consisting of initial layer-parallel shortening followed by layer-parallel shear (related to flexural folding). Also in this instance, volume change during tectonic deformation and tectonic plane strain have been assumed. For geologically reasonable amounts of volume loss due to compaction and of initial layer-parallel shortening, this type of deformation history is capable of producing a deformation path entirely lying within the oblate strain field, but still characterized by a changeover, during deformation, of the maximum extension axis (X) from a position parallel to the fold axis to one perpendicular to it. This type of deformation sequence may explain the main strain features observed in the study area, where most of the measured finite strain ellipsoids, determined from the limb regions of flexural folds, display an oblate shape, irrespective of the orientation of their maximum extension direction (X) with respect to the local structural trends. More generally, this type of deformation history provides a mechanism to account for the predominance of oblate strains in deformed pelitic rocks.