海洋天然气水合物的地震识别方法研究

天然气水合物作为21世纪新的自然能源将为人类的生存发展服务。20世纪60年代证实，俄罗斯西伯利亚的麦索亚哈气田为典型的天然气水合物形成的气田，70年代又在海底发现了固体天然气水合物岩样。1971年，RStoll首先将地震剖面中的似海底反向层解释为海洋天然气水合物存在的标志，后来被深海钻探证实，从此地震方法成为大面积研究天然气水合物的重要手段。天然气水合物既是潜在能源，也是影响环境和形成灾害的因素之一，因此，研究天然气水合物是人类在21世纪的重要课题。探讨海洋天然气水合物的地震识别方法，由于这项工作刚刚起步，还没有做出具体的成果，在此只能根据我们仅有的工作和参照国外分开的出版物，以及出国访问得到的有关资料进行分析，提出我们的一些基本设想，与各位专家探讨。     

Hydrodynamic modeling of flushing time in a small estuary of North Bay, Florida, USA

Freshwater fraction method is popular for cost-effective estimations of estuarine flushing time in response to freshwater inputs. However, due to the spatial variations of salinity, it is usually expensive to directly estimate the long-term freshwater fraction in the estuary from field observations. This paper presents the application of the 3D hydrodynamic model to estimate the distributions of salinity and thus the freshwater fractions for flushing time estimation. For a case study in a small estuary of the North Bay in Florida, USA, the hydrodynamic model was calibrated and verified using available field observations. Freshwater fractions in the estuary were determined by integrating freshwater fractions in model grids for the calculation of flushing time. The flushing time in the North Bay is calculated by the volume of freshwater fraction divided by the freshwater inflow, which is about 2.2 days under averaged flow conditions. Based on model simulations for a time series of freshwater inputs over a 2-year period, a power regression equation has been derived from model simulations to correlate estuarine flushing time to freshwater inputs. For freshwater input varying from 12 m³/s to 50 m³/s, flushing time in this small estuary of North Bay changes from 3.7 days to 1.8 days. In supporting estuarine management, the model can be used to examine the effects of upstream freshwater withdraw on estuarine salinity and flushing time.     

Velocity-Interface Structure of the Southwestern Ryukyu Subduction Zone from EW9509-1 OBS/MCS Data

A wide-angle seismic survey, combining ocean-bottom seismometers (OBS) and multi-channel seismic (MCS) profiling, was implemented in the southwestern Ryukyu subduction zone during August and September 1995. In this paper, we present the data analysis of eight OBSs and the corresponding MCS line along profile EW9509-1 from this experiment. Seismic data modeling includes identification of refracted and reflected arrivals, initial model building from velocity analysis of the MCS data, and simultaneous and layer-stripping inversions of the OBS and MCS arrivals. The velocity-interface structure constructed along profile EW9509-1 shows that the northward subduction of the Philippine Sea Plate has resulted in a northward thickening of the sediments of the Ryukyu Trench and the Yaeyama accretionary wedge north of the trench. The boundary between the subducting oceanic crust and the overriding continental crust (represented by a velocity contour of 6.75 km/s) and a sudden increase of the subducting angle (from 5 degrees to 25 degrees) are well imaged below the Nanao Basin. Furthermore, velocity undulation and interface variation are found within the upper crust of the Ryukyu Arc. Therefore, the strongest compression due to subduction and a break-off of the slab may have occurred and induced the high seismicity in the forearc region. This revised version was published online in November 2006 with corrections to the Cover Date.     

广东红海湾、珠江口及电白—阳江近海第四纪沉积

根据物探和钻探资料，研究区内第四纪沉积物可分为两部分，下部地震反射结构复杂，斜层理发育，粒度较粗，为晚更新世产物；上部的地震反射结构简单，呈水平状或近水平状结构，粒度较细，为全新世沉积。     

Southern Segment of the European Geotraverse — a wide-angle seismic refraction experiment in the Sardinia Channel

The Sardinia Channel dataset was collected as part of the European Geotraverse (EGT)—a 4000 km seismic refraction line running from Northern Norway to the Sahara, designed to investigate the structure of the lithosphere beneath Europe. Wideangle seismic data recorded by ocean bottom seismometers deployed in the Sardinia Channel as part of the Southern Segment of the EGT, together with gravity data, were used to constrain the final crustal model. In the centre of the Channel the crust is identified as thinned continental in nature, with a crystalline thickness of 10 km overlain by 4 km of sediments and 2.5 km of water in the most extended region. High velocities in the lower crust in the central region are thought to represent an area of underplating or intrusion by igneous material caused by extension related to the opening of the Tyrrhenian Sea. The crust overlies an anomalously low velocity upper mantle.     

近58a来影响和登陆浙江的热带气旋统计特征分析

本文利用1949-2006年热带气旋资料,对58a来影响和登陆浙江热带气旋的登陆位置、登陆时间的年际、月际和日变化、登陆后的移速及路径变化与降水中心分布的关系、登陆后维持和衰减情况等进行了统计分析.统计结果表明：影响和登陆浙江的热带气旋有较明显的年际、月际和日变化特征;登陆浙江的热带气旋陆上维持时间与登陆时的强度正相关,登陆后12h内热带气旋衰减较快;登陆后移向偏北分量大的或移速加快的热带气旋,主要降水区出现在路径右侧的可能性大.     

Correspondence between the distribution of hydrodynamic time parameters and the distribution of biological and chemical variables in a semi-enclosed coral reef lagoon

Hydrodynamic modeling can be used to spatially characterize water renewal rates in coastal ecosystems. Using a hydrodynamic model implemented over the semi-enclosed Southwest coral lagoon of New Caledonia, a recent study computed the flushing lag as the minimum time required for a particle coming from outside the lagoon (open ocean) to reach a specific station [Jouon, A., Douillet, P., Ouillon, S., Fraunié, P., 2006. Calculations of hydrodynamic time parameters in a semi-opened coastal zone using a 3D hydrodynamic model. Continental Shelf Research 26, 1395–1415]. Local e-flushing time was calculated as the time requested to reach a local grid mesh concentration of 1/e from the precedent step. Here we present an attempt to connect physical forcing to biogeochemical functioning of this coastal ecosystem. An array of stations, located in the lagoonal channel as well as in several bays under anthropogenic influence, was sampled during three cruises. We then tested the statistical relationships between the distribution of flushing indices and those of biological and chemical variables. Among the variables tested, silicate, chlorophyll a and bacterial biomass production present the highest correlations with flushing indices. Correlations are higher with local e-flushing times than with flushing lags or the sum of these two indices. In the bays, these variables often deviate from the relationships determined in the main lagoon channel. In the three bays receiving significant riverine inputs, silicate is well above the regression line, whereas data from the bay receiving almost insignificant freshwater inputs generally fit the lagoon channel regressions. Moreover, in the three bays receiving important urban and industrial effluents, chlorophyll a and bacterial production of biomass generally display values exceeding the lagoon channel regression trends whereas in the bay under moderate anthropogenic influence values follow the regressions obtained in the lagoon channel. The South West lagoon of New Caledonia can hence be viewed as a coastal mesotrophic ecosystem that is flushed by oligotrophic oceanic waters which subsequently replace the lagoonal waters with water considerably impoverished in resources for microbial growth. This flushing was high enough during the periods of study to influence the distribution of phytoplankton biomass, bacterial production of biomass and silicate concentrations in the lagoon channel as well as in some of the bay areas.     

High parallelization efficiency in barotropic-mode computation of ocean models based on multi-grid boundary ghost area

In general, barotropic-mode computation requires the largest communication time in ocean models because of its iterative nature, when parallel computation is performed based on regional partitioning. In this study, reduction in the overhead time included in the communication in the parallel computation of the barotropic mode is considered to achieve a high parallelization efficiency of ocean models. We verify that the reduction in the communication frequency based on the multi-grid boundary ghost area reduces the total communication time. We find that this is because the sum of several kinds of overhead time for communication occupies a fairly large part of the total communication time. We discuss the trade-off between the decrease in communication time and the increase in computation time due to increased boundary area in such cases, leading to a determination of a suitable width of the multi-grid boundary that minimizes the total required time. We also discuss the efficiencies of the one- and two-dimensional partitioning of the model domain, when a multi-grid boundary is used. In general, two-dimensional partitioning is more efficient than one-dimensional for large partitioning numbers.     

A technique for conducting seismic refraction experiments on the ocean bed using bottom shots

A technique has been devised for firing arrays of bottom shots on the ocean bed in depths upto 4000 m or more. Ten kilogram explosive charges are dropped from the surface while the shooting ship is navigated acoustically. They are detonated at preset times by an electronic timer which initiates an electrical detonator, detonating cord and cast PETN/TNT explosive. Ranges to ocean bottom seismographs, and the shot instants, can be calculated from the arrival-time differences of the direct and surface-reflected water waves. The accuracy, which is dependent on water-depth and range, was better than 22 m in range and 14 msec in shot instant for our experiments.     

冲绳海槽天然气水合物BSR的地震研究

根据多道地震反射资料分析,在冲绳海槽南部和中部发现了拟海底反射层(BSR)现象。通过对海底异常反射层的振幅特征、速度异常和AVO属性分析,说明该BSR可能反映了天然气水合物的存在,并发现冲绳海槽断层与天然气水合物的形成有密切关系。