Relative paleointensity of the geomagnetic field during the past 200 ka from the West Philippine Sea and its chronological significance

Our rock magnetic analysis of core Ph05 from the West Philippine Sea demonstrates that the core preserves a strong, stable remanent magnetization and meets the magnetic mineral criteria for relative paleointensity (RPI) analyses. The magnetic minerals in the sequence are dominated by pseudosingle-domain magnetite, and the concentration of magnetic minerals is at the same scale. Both the conventional normalizing method and the pseudo-Thellier method were used in conjunction with the examination of the rock magnetic properties and natural remanent magnetization. Susceptibility (chi), anhysteretic remnant magnetization (ARM) and saturation isothermal remnant magnetization (SIRM) were used as the natural remanent magnetization normalizer. However, coherence analysis indicated that only ARM is more suitable for paleointensity reconstruction. The age model of core is established based on oxygen isotope data and AMS(14)C data, which is consistent with the age model estimated from RPI records. The relative paleointensity data provide a continuous record of the intensity variation during the last 200 ka, which correlates well with the global references RPI stacks. Several prominent low paleointensity values are identified and are correlated to the main RPI minima in the SINT-200 record, suggesting that the sediments have recorded the real changes of geomagnetic field.

P wave velocity of the pumice from the Okinawa Trough at high temperature and high pressure.

P wave velocity of the pumice sample from the middle Okinawa Trough and andesite sample from vicinity Yingdao volcanic island, Kyushu Japan were measured at temperature (from room temperature to 1500 C) and pressure (from room pressure to 2.4GPa) using a multi-anvil pressure apparatus called the YJ-3000 press. The measured data shows that at low temperature and low pressure (<1GPa, <800degreesC), the P wave velocity of pumice is lower than that of andesite, while at high temperature and high pressure (>1GPa, >800degreesC) the P wave velocity of pumice and andesite. becomes consistent (5.9km/s). The paper points out that 1GPa/800degreesC is the point of thermodynamic phase transformation Okinawa Trough pumice and vicinity andesite, and the point is deeper than 18km.

Impact of spatial resolution and spatial difference accuracy on the performance of Arakawa A-D grids

This paper alms at illustrating the impact of spatial difference scheme and spatial resolution on the performance of Arakawa A-D grids in physical space. Linear shallow water equations are discretized and forecasted on Arakawa A-D grids for 120-minute using the ordinary second-order (M and fourth-order (C4) finite difference schemes with the grid spacing being 100 km, 10 km and I km, respectively. Then the forecasted results are compared with the exact solution, the result indicates that when the grid spacing is I kin, the inertial gravity wave can be simulated on any grid with the same results from C2 scheme or C4 scheme, namely the impact of variable configuration is neglectable; while the inertial gravity wave is simulated with lengthened grid spacing, the effects of different variable configurations are different. However, whether for C2 scheme or for C4 scheme, the RMS is minimal (maximal) on C (D) grid. At the same time it is also shown that when the difference accuracy increases from C2 scheme to C4 scheme, the resulted forecasts do not uniformly decrease, which is validated by the change of the group A velocity relative error from C2 scheme to C4 scheme. Therefore, the impact of the grid spacing is more important than that of the difference accuracy on the performance of Arakawa A-D grid.

On the western boundary currents in the Philippine Sea

In this study we describe the velocity structure and transport of the North Equatorial Current (NEC), the Kuroshio, and the Mindanao Current (MC) using repeated hydrographic sections near the Philippine coast. A most striking feature of the current system in the region is the undercurrent structure below the surface flow. Both the Luzon Undercurrent and the Mindanao Undercurrent appear to be permanent phenomena. The present data set also provides an estimate of the mean circulation diagram (relative to 1500 dbar) that involves a NEC transport of 41 Sverdrups (Sv), a Kuroshio transport of 14 Sv, and a MC transport of 27 Sv, inducing a mass balance better than 1 Sv within the region enclosed by stations. The circulation diagram is insensitive to vertical displacements of the reference level within the depth range between 1500 and 2500 dbar. Transport fluctuations are, in general, consistent with earlier observations; that is, the NEC and the Kuroshio vary in the same phase with a seasonal signal superimposed with interannual variations, and the transport of the MC is dominated by a quasi-biennial oscillation. Dynamic height distributions are also examined to explore the dynamics of the current system.

A subsurface countercurrent along the east coast of Luzon

A hydrographic section in the region east of Luzon was repeated 14 times during the period from 1986 to 1991. The data revealed the existence of a subsurface countercurrent located on the shoreward side of the Kuroshio with its upper boundary at about 500 m. The countercurrent, which should be called the Luzon Undercurrent (LUG), was only about 50 km wide, which is comparable to the baroclinic radius of deformation. Despite considerable variabilities both in velocity profile and intensity, the LUC appears to be a permanent feature. Over the period of observations, the maximum speed in the LUC calculated from the mean temperature and salinity by assuming geostrophy (relative to 2500 db) was 7 cm s(-1) at about 700 m and its mean geostrophic volume transport was 3.6 Sv (1 Sv = 10(6) m(3) s(-1)). About 28% of this transport was composed of the low-salinity North Pacific Intermediate Water (NPIW) advected to the south along the coast of Luzon. (C) 1997 Elsevier Science Ltd.

南海北部陆缘深水区构造演化及其资源效应

南海北部陆缘深水区（水深>300m）蕴藏着丰富的资源，我国对深水区的地质研究刚刚起步，但相关领域已成为科研热点。深水油气盆地的构造演化是油气勘探中最重要的基础性研究之一，因此针对我国南海北部陆缘深水区开展构造演化及其资源效应的研究具有重要的理论意义和实际意义。 本文利用钻井和地震资料并结合区域地质资料，重点研究了珠江口盆地深水区的结构和构造演化，取得如下创新性成果：1）首次利用半地堑分析方法系统解剖了研究区的结构、各构造单元发育特征，在此基础上指出五个有利油气运聚带；2）采用回剥法并利用最新资料进行校正，得到了研究区更为可靠的构造沉降曲线，重新划分了裂陷期和裂后期的分界，认为32Ma南海海底扩张之后裂陷作用仍在持续，直到23Ma左右才开始大规模裂后热沉降，并进一步解释了裂陷期延迟的形成机制；3）应用非连续拉张模型计算拉张系数的方程计算了研究区的壳幔拉张系数，指出了深水区地幔相对于地壳的优势伸展作用；首次运用平衡剖面技术重建了研究区的构造发育史，计算了各构造期的拉张率和沉积速率，指出研究区新生代整体呈现持续拉张，拉张系数在1.1-1.24之间；4）精细刻画了水合物钻采区的地质构造特征，建立了该区天然气水合物成藏的概念模式；建立了一套根据地震叠加速度计算流体势的方法，为水合物成藏规律的研究提供了新的思路。

Numerical study on the velocity structure around tidal fronts in the Yellow Sea

The velocity components across tidal fronts are examined using the Blumberg and Mellor 3-D nonlinear numerical coastal circulation model incorporated with the Mellor and Yamada level 2.5 turbulent closure model based on the reasonable model output of the M-2 tide and density residual currents. In the numerical experiments, upwelling motion appears around all the fronts with different velocity structures, accounting for surface cold water around the fronts. The experiments also suggest that the location and formation of fronts are closely related to topography and tidal mixing, as is the velocity structure around the front.

Control mechanisms of diel vertical migration: Theoretical assumptions

We explore control mechanisms underlying the vertical migration of zooplankton in the water column under the predator-avoidance hypothesis. Two groups of assumptions in which the organisms are assumed to migrate vertically in order to minimize realized or effective predation pressure (type-I) and to minimize changes in realized or effective predation pressure (type-II), respectively, are investigated. Realized predation pressure is defined as the product of light intensity and relative predation abundance and the part of realized predation pressure that really affects organisms is termed as effective predation pressure. Although both types of assumptions can lead to the migration of zooplankton to avoid the mortality from predators, only the mechanisms based on type-II assumptions permit zooplankton to undergo a normal diel vertical migration (morning descent and evening ascent). The assumption of minimizing changes in realized predation pressure is based on consideration of DVM induction only by light intensity and predators. The assumption of minimizing changes in effective predation pressure takes into account, apart from light and predators also the effects of food and temperature. The latter assumption results in the same expression of migration velocity as the former one when both food and temperature are constant over water depth. A significant characteristic of the two type-II assumptions is that the relative change in light intensity plays a primary role in determining the migration velocity. The photoresponse is modified by other environmental variables: predation pressure, food and temperature. Both light and predation pressure are necessary for organisms to undertake DVM. We analyse the effect of each single variable. The modification of the phototaxis of migratory organisms depends on the vertical distribution of these variables. (C) 2001 Academic Press.

Embryology of Swertia (Gentianaceae) relative to taxonomy

The embryological features of three species of Swertia (s.l.) - S. erythrosticta, S. franchetiana, and S. tetraptera were characterized, and the observations were used, together with previously gathered data on other species, to evaluate a recently proposed polyphyly, based on molecular data, of Swertia s.l. Comparisons of species within the genus showed that they have diversified embryologically, and there are significant between-species differences. Notable features that vary between species include the number of cell layers that form the anther locule wall, the construction of the wall of the mature anther, tapetum origin, the cell number in mature pollen grains, the structure of the fused margins of the two carpels, the ovule numbers in placental cross-sections, the shape of the mature embryo sac, the degree of ovule curvature, antipodal variation and the presence of a hypostase, and seed appendages. They share characters that are widely distributed in the tribe Gentianeae, such as a dicotyledonous type of anther wall formation, a glandular tapetum with uninucleate cells, simultaneous cytokinesis following the meiosis of the microsporocytes, tetrahedral microspore tetrads, superior, bicarpellary and unilocular ovaries, unitegmic and tenuinucellar ovules, Polygonum-type megagametophytes, progamous fertilization, nuclear endosperm, and Solanad-type embryogeny. The presence of variation in embryological characters amongst the species of Swertia s.l. strongly supports the view that Swertia s.l. is not a monophyletic group. Frasera is better separated from Swertia s.l. as an independent genus, and is only distantly related to Swertia s. s. judging from the numerous differences in embryology. Swertia tetraptera is very closely related to Halenia, as they show identical embryology. (C) 2007 The Linnean Society of London, Botanical Journal of the Linnean Society, 2007, 155, 383-400.

Global Translation Preferred Active Contour

具有全局平移优先属性的主动轮廓更适于目标跟踪。演化轮廓具有的全局平移优先性可以理解为沿轮廓的速度场具有相等的倾向。根据此思想,通过定义在曲线扰动集合上的新内积空间导出了一种简单,具有平移优先的梯度流。新的内积空间由于是通过向H0主动轮廓对应的內积空间引入曲线扰动的方差获得,所以此主动轮廓称为方差主动轮廓。方差主动轮廓是将H0主动轮廓与其对应的平均梯度流通过加权求和获得,而H1主动轮廓则是通过H0主动轮廓与特定类型的核函数进行卷积得到。因此方差主动轮廓实现时更简单和快速。最后给出了H0,H1和方差主动轮廓在频域与时域的分析。

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

二维自组织介质中地震波的传播特性

The topic of this study is simulation in the two dimensional self-organized media. The study in complexity of the earth plays an important role in structures, sources and energy seismic detection. And it mainly focuses on vertical or horizontal heterogeneous, anisotropic and linear media. Based on 2D self-organized velocity model and four-order finite-difference method, we simulate different types self-organized media and the same type mode with various parameters such as horizontal relative length, vertical relative length, variations, and velocity background gradient. Also we analyze the seismograms with complexity methods with instant information including amplitude, energy and frequency. The results can be summarized as the fallows: (1) The waveforms fluctuate with the velocity variations; (2) Different type self-organized media bring different effects on the amplitudes, energy and waveforms; (3) Different parameters also produce various influences to seismograms. (4) The layer contains their self-organized features, from which we can investigate the quality of the earth.

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

微动勘查方法在煤矿采区勘探中的应用研究

The conventional microtremor survey is based on the single point of exploration, which includes collecting field data，estimating the phase velocity, investing the dispersion curve and obtaining the S–wave velocity structure. In the case of large-scale exploration, and when making the two-dimensional velocity section, the inversion is quite time-consuming, laborious and its precision depends on the subjective interpretation, which makes the results differently from person to person. In fact，we do not need the S-wave velocity values but only need the relative variation of velocity. For these reasons, this paper is desired to calculate the apparent S-wave velocity (Vx) to replace the S-wave velocity inversion and to obtain the relative variation of the S-wave velocity. Using this method, we can decrease the analysts’ effect, shorten the data processing time and improve work efficiency. The apparent S-wave velocity is a variable of the surface wave property, which can clearly reflect the downcast columns, mined-out areas and other unusual geological bodies. In this paper, Matlab is used to establish the three-dimensional data volume of the apparent S-wave velocity, from which we can get any apparent S-wave velocity section we need. Through the application case, the designed method is proved to be reliable and effective. The downcast columns, mined-out areas and other unusual geological bodies can be clearly showed in the apparent S-wave velocity section. And from the contour of the apparent S-wave velocity, the interface shape of the major target layers can be controlled basically.

人工源宽角地震资料中的广角效应及其对地壳速度成像的影响

The seismic wide-angle reflection/refraction method is the one of the most effective method for probing the crustal and upper mantle structure. It mainly uses the wide-angle reflection information from the boundary in the crust and the top boundary of the upper mantle to rebuild the crust and upper mantle structure. Through analyzing the reflection and transmission coefficients of various incident waves on the interface, we think relative to the pre-critical angle reflection information the post critical angle reflection information that received by wide-angle seismic data exists a time-shift effect with the offset variation, and then it must cause the error for velocity analysis and structure image. The feature of the wide-angle seismic wave field of the fourteen representative crust columns tell us that the wide-angle effects in the different representative tectonic units for the interface depth and the interval velocity in crust. We studied the features of the wide-angle seismic wave field through building the crust model and inverse its travel time by GA method to know the wide-angle influence on crustal velocity image. At last we finished the data processing of the Tunxi-Wenzhou wide-angle seismic profile. The results are as following: (1) Through building crust model, we labeled the travel time for all the phases by ray tracing method and remove wide-angle effects method, it revealed the wide-angle effect exists in the seismic data. (2) The travel time inversion by GA method can tell us that the depth by traditional ray tracing method is shallower than the result by remove wide-angle effects method, the latter can recover the crust structure model in effect. (3) We applied the two method mentioned before to the fourteen representative crust columns in China. It indicates that the removed wide-angle effect method in travel time inversion is reasonable and effective. (4) The real data processing from Tunxi-Wenzhou wide-angle seismic profile give us the basic structure through the two ways. The main influence exhibits in the difference of the interval velocity of the curst, and the wide-angle effects in shallow interface are stronger than the deep interface.