The cDNA cloning and mRNA expression of heat shock protein 70 gene in the haemocytes of bay scallop (Argopecten irradians, Lamarck 1819) responding to bacteria challenge and naphthalin stress

Heat shock protein 70 (HSP70) is an important member of the heat shock protein superfamily, and it plays a key role in the process of protecting cells, facilitating the folding of nascent peptides and responding to stress. The cDNA of bay scallop Argopecten irradians HSP70 (designated AIHSP70) was cloned by the techniques of homological cloning and rapid amplification of cDNA end (RACE). The full length of AIHSP70 cDNA was 2651 bp in length, having a 5' untranslated region (UTR) of 96 bp, a 3' UTR of 575 bp, and an open reading frame (ORF) of 1980 bp encoding a polypeptide of 659 amino acids with an estimated molecular mass of 71.80 kDa and an estimated isoelectric point of 5.26. BLAST analysis revealed that the AIHSP70 gene shared high identity with other known HSP70 genes. Three classical HSP signature motifs were detected in AIHSP70 by InterPro, analysis. 3-D structural prediction of AIHSP70 showed that its N terminal ATPase activity domain and,C terminal substrate-binding domain shared high similarity with that in human heat shock protein 70. The results indicated that the AIHSP70 was a member of the heat shock protein 70 family. A semi-quantitive RT-PCR method was used to analyse the expression of AIHSP70 gene after the treatment of naphthalin which is one kind of polycyclic aromatic hydrocarbon (PAH) and the challenge of bacteria. mRNA expression of AIHSP70 in scallop was up-regulated significantly after the stimulation of naphthalin and increased with increasing naphthalin concentration. A clearly time-dependent expression pattern of AIHSP70 was observed after the scallops were infected by Vibrio anguillarum, and the mRNA expression reached a maximum level at 8 h and lasted to 16 h, and then dropped progressively. The results indicated that AIHSP70 could play an important role in mediating the environmental stress and immune response in scallop. (c) 2006 Elsevier Ltd. All rights reserved.

The first molluscan TCTP in Venerupis philippinarum: Molecular cloning and expression analysis

Translationally controlled tumor protein (TCTP) is one of the abundant and ubiquitously expressed proteins in metazoans In the present study, the first molluscan TCTP (denoted as VpTCTP) was identified from Venerupis philippinarum haemocytes by EST and RACE approaches The full-length cDNA of VpTCTP consisted of 1148 nucleotides with an open-reading frame of 555 bp encoding 184 amino acids The deduced amino acid sequence of VpTCTP shared high similarity with TCTPs from other species, indicating that VpTCTP should be a new member of TCTP family Several highly conserved motifs, including 5'terminal ologopyrimidine (5'TOP) starting sequence and rich AU and AUUT elements in 3'UTR, were also identified in VpTCTP The tissue and temporal expression of VpTCTP after Vi boo anguillarum challenge was recorded by quantitative real-time RT-PCR. VpTCTP transcript could be detected in all examined tissues with the highest expression level in haemocytes and the lowest in hepatopancreas Concerning the time-course expression in haemocytes, the relative expression of VpTCTP mRNA was down-regulated sharply from 6 h to 12 h post-infection. Then, the expression level was obviously up-regulated and reached 3.4-fold to that in the control group at 48 h post challenge As time progressed, the expression of VpTCTP recovered to the original level at 96 h. All these results indicated that VpTCTP was an acute-phase protein involved in the Immune response of V philippinarum (C) 2010 Elsevier Ltd. All rights reserved.

Molecular cloning and expression analysis of Crustin-like gene from Chinese shrimp Fenneropenaeus chinensis

A new member of antimicrobial protein genes of the Crustin family was cloned from haemocytes of the Chinese shrimp Fenneropenaeus chinensis by 3' and 5' RACE. The full-length cDNA of Crustin-like gene contains a 390 bp open reading frame, encoding 130 amino acids. The deduced peptide contains a putative signal peptide of 17 amino acids and mature peptide of 113 amino acids. The molecular mass of the deduced mature peptide is 12.3 ku. It is highly cationic with a theoretical isoelectric point of 8.5. The deduced amino acids sequence of this Crustin showed high homology with those of Penaeus (Litopenaeas) setferus. Northern blotting showed that the cloned Crustin gene was mainly expressed in haemocytes, gill, intestine, and RNA in situ hybridization indicated that the Crustin gene was constitutively expressed exclusively in haemocytes of these tissues. Capillary electrophoresis RT-PCR analysis showed that Crustin was up-regulated dramatically from 12 to 48 h after a brief decrease of mRNA during first 6 h in response to microbe infection. The level of Crustin mRNA began to restore at 72 h post-challenge. This indicated that Crustin gene might play an important role when shrimps are infected by bacterial pathogen.

Tide-influenced acidic hydrothermal system offshore NE Taiwan

Elemental sulfur and hydrogen sulfide emitted offshore of northeastern Taiwan known to local fishermen for generations, but never studied until recently, are found to have originated from a cluster of shallow (< 30 m depth) hydrothermal vents. Among the mounds is a massive 6 m high chimney with a diameter of 4 m at the base composed of almost pure sulfur and discharging hydrothermal fluid containing sulfur particles. The sulfur in the chimney has a delta(34)S= 1.1 parts per thousand that is isotopically lighter than seawater. A yellow smoker at shallow depths with such characteristics has never been reported on anywhere else in the world. Gas discharges from these vents are dominated by CO2 (> 92%) with small amounts of H2S. Helium isotopic ratios 7.5 times that of air indicate that these gases originate from the mantle. High temperature hydrothermal fluids have measured temperatures of 78-116 degrees C and pH (25 degrees C) values as low as 1.52, likely the lowest to be found in world records. Low temperature vents (30-65 degrees C) have higher pH values. Continuous temperature records from one vent show a close correlation with diumal tides, suggesting rapid circulation of the hydrothermal fluids. (c) 2005 Elsevier B.V. All rights reserved.

Reconstruction of vertical thermal structure from several subsurface temperatures in the China Seas and adjacent waters

Empirical Orthogonal Function (EOF) analysis is used in this study to generate main eigenvector fields of historical temperature for the China Seas (here referring to Chinese marine territories) and adjacent waters from 1930 to 2002 (510 143 profiles). A good temperature profile is reconstructed based on several subsurface in situ temperature observations and the thermocline was estimated using the model. The results show that: 1) For the study area, the former four principal components can explain 95% of the overall variance, and the vertical distribution of temperature is most stable using the in situ temperature observations near the surface. 2) The model verifications based on the observed CTD data from the East China Sea (ECS), South China Sea (SCS) and the areas around Taiwan Island show that the reconstructed profiles have high correlation with the observed ones with the confidence level > 95%, especially to describe the characteristics of the thermocline well. The average errors between the reconstructed and observed profiles in these three areas are 0.69A degrees C, 0.52A degrees C and 1.18A degrees C respectively. It also shows the model RMS error is less than or close to the climatological error. The statistical model can be used to well estimate the temperature profile vertical structure. 3) Comparing the thermocline characteristics between the reconstructed and observed profiles, the results in the ECS show that the average absolute errors are 1.5m, 1.4 m and 0.17A degrees C/m, and the average relative errors are 24.7%, 8.9% and 22.6% for the upper, lower thermocline boundaries and the gradient, respectively. Although the relative errors are obvious, the absolute error is small. In the SCS, the average absolute errors are 4.1 m, 27.7 m and 0.007A degrees C/m, and the average relative errors are 16.1%, 16.8% and 9.5% for the upper, lower thermocline boundaries and the gradient, respectively. The average relative errors are all < 20%. Although the average absolute error of the lower thermocline boundary is considerable, but contrast to the spatial scale of average depth of the lower thermocline boundary (165 m), the average relative error is small (16.8%). Therefore the model can be used to well estimate the thermocline.

趋磁细菌AMB-1培养条件优化及磁小体变化过程研究

趋磁细菌(Magnetotactic bacteria)的研究是国际微生物学研究热点之一。趋磁细菌体内含有纳米单磁畴的氧化铁/硫化铁(Fe3O4或Fe3S4)晶体，称为磁小体。由于趋磁细菌营养条件要求苛刻，在环境中需要微好氧条件，且营养类型属于化能自养，使得培养趋磁细菌时常遇到问题。 本研究首先通过正交试验优化趋磁细菌AMB-1菌株培养条件，在培养条件铁源为奎尼酸铁0.02 mmol/L，装瓶量75% ，pH值6.7，温度25 ℃时，AMB-1 OD600达到0.440(1.166×109 cells/ml)。同时运用磁收集传代法，使带有磁小体的AMB-1细胞比例占95％以上(Cmag值稳定在1.9-2.0)。 在AMB-1具有较好的生物量，同时又具有较好的含磁小体细胞比例后，研究磁小体的变化过程。通过透射电镜观察磁小体变化过程，发现培养24 h细菌体内已有较小晶体形成(平均27 nm，n=188)且沿长轴分布；48 h晶体长大(平均43 nm，n=203)且形成分段链沿长轴排列；72 h晶体进一步成熟(平均50 nm，n=191)仍以分段链沿长轴排列；随后细菌逐渐衰亡磁小体变小，168 h可见部分自溶细菌中仍有磁小体链(平均37 nm，n=186)；192 h细菌自溶磁小体链(平均33 nm，n=184)分散到环境中。 通过透射电镜在细胞水平上研究趋磁细菌细胞分裂时发现，磁小体在细菌分裂时采用两种分离方式：一种为磁小体分配到两个子细胞；另一种为磁小体只分配到一个子细胞。无磁小体的子细胞，在随后的生长过程又分为两种情况：一种为细胞逐渐产生磁小体，另一种为不再产生磁小体。这种现象的发现，解释了随着传代次数的增多，细菌磁性有所下降的原因(Cmag值降低)。 在对趋磁细菌磁小体合成机制的研究中，常使用基因敲除的办法获得缺陷型，并与野生型对比进行研究。但是，利用基因敲除获得缺陷型不仅操作繁琐并且所得缺陷型不稳定。本研究利用特殊的磁富集传代法，先将带有磁小体的菌体收集并连续传代，筛选获得了高磁菌株；利用这种方法，收集不含磁小体的菌体并连续传代，筛选获得了无磁菌株。 趋磁细菌磁小体在医疗、环保等领域具有广阔应用价值，但是目前由于趋磁细菌难以大规模培养，并且磁小体纯化存在成本高等原因，将磁小体真正实际应用尚有一段距离。通过研究磁小体在趋磁细菌中的变化过程发现，AMB-1菌株在培养192 h后自溶，并且磁小体随着细胞的破碎释放到环境中去。

中国近海和西北太平洋温跃层时空变化分析、模拟及预报

温度跃层是反映海洋温度场的重要物理特性指标，对水下通讯、潜艇活动及渔业养殖、捕捞等有重要影响。本文利用中国科学院海洋研究所“中国海洋科学数据库”在中国近海及西北太平洋（110ºE-140ºE,10ºN-40ºN）的多年历史资料（1930-2002年，510143站次），基于一种改进的温跃层判定方法，分析了该海域温跃层特征量的时空分布状况。同时利用Princeton Ocean Model(POM)，对中国近海，特别是东南沿海的水文结构进行了模拟，研究了海洋水文环境对逆温跃层的影响。最后根据历史海温观测资料，利用EOF分解统计技术，提出了一种适于我国近海及毗邻海域，基于现场有限层实测海温数据，快速重构海洋水温垂直结构的统计预报方法，以达到对现场温跃层的快速估计。 历史资料分析结果表明，受太阳辐射和风应力的影响，20°N以北研究海域，温跃层季节变化明显，夏季温跃层最浅、最强，冬季相反，温跃层厚度的相位明显滞后于其他变量，其在春季最薄、秋季最厚。12月份到翌年3月份，渤、黄及东海西岸，呈无跃层结构，西北太平洋部分海域从1月到3月份，也基本无跃层结构。在黄海西和东岸以及台湾海峡附近的浅滩海域，由于风力搅拌和潮混合作用，温跃层出现概率常年较低。夏季，海水层化现象在近海陆架海域得到了加强，陆架海域温跃层强度季节性变化幅度(0.31°C/m)明显大于深水区(约0.05°C/m)，而前者温跃层深度和厚度的季节性变化幅度小于后者。20°N以南研究海域，温跃层季节变化不明显。逆温跃层主要出现在冬、春季节（10月-翌年5月）。受长江冲淡水和台湾暖流的影响，东南沿海区域逆温跃层持续时间最长，出现概率最大，而在山东半岛北及东沿岸、朝鲜半岛西及北岸，逆温跃层消长过程似乎和黄海暖流有关。多温跃层结构常年出现于北赤道流及对马暖流区。在黑潮入侵黄、东、南海的区域，多温跃层呈现明显不同的季节变化。在黄海中部，春季多温跃层发生概率高于夏季和秋季，在东海西部，多跃层主要出现在夏季，在南海北部，冬季和春季多温跃层发生概率大于夏季和秋季。这些变化可能主要受海表面温度变化和风力驱动的表层流的影响。 利用Princeton Ocean Model(POM)，对中国东南沿海逆温跃层结构进行了模拟，模拟结果显示，长江冲淡水的季节性变化以及夏季转向与实际结果符合较好，基本再现了渤、黄、东海海域主要的环流、温盐场以及逆温跃层的分布特征和季节变化。通过数值实验发现，若无长江、黄河淡水输入，则在整个研究海域基本无逆温跃层出现，因此陆源淡水可能是河口附近逆温跃层出现的基本因素之一。长江以及暖流（黑潮和台湾暖流）流量的增加，均可在不同程度上使逆温跃层出现概率及强度、深度和厚度增加，且暖流的影响更加明显。长江对东南沿海逆温跃层的出现，特别是秋季到冬季初期，有明显的影响，使长江口海域逆温跃层位置偏向东南。暖流对于中国东南沿海的逆温跃层结构，特别是初春时期，有较大影响，使长江口海域的逆温跃层位置向东北偏移。 通过对温跃层长期变化分析得出，黄海冷水团区域，夏季温跃层强度存在3.8年左右的年际变化及18.9年左右的年代际变化，此变化可能主要表现为对当年夏季和前冬东亚地区大气气温的热力响应。东海冷涡区域，夏季温跃层强度存在3.7年的年际变化，在El Nino年为正的强度异常，其可能主要受局地气旋式大气环流变异所影响。谱分析同时表明，该海域夏季温跃层强度还存在33.2年的年代际变化，上世纪70年代中期，温跃层强度由弱转强，而此变化可能与黑潮流量的年代际变化有关。 海洋水温垂直结构的统计预报结果显示，EOF分解的前四个主分量即能够解释原空间点温度距平总方差的95%以上，以海洋表层附近观测资料求解的特征系数推断温度垂直结构分布的结果最稳定。利用东海陆架区、南海深水区和台湾周边海域三个不同区域的实测CTD样本廓线资料，对重构模型的检验结果表明，重构与实测廓线的相关程度超过95%的置信水平。三个区重构与实测温度廓线值的平均误差分别为0.69℃，0.52℃，1.18℃，平均重构廓线误差小于平均气候偏差，统计模式可以很好的估算温度廓线垂直结构。东海陆架海区温度垂直重构廓线与CTD观测廓线获得的温跃层结果对比表明，重构温跃层上界、下界深度和强度的平均绝对误差分别为1.51m、1.36m和0.17℃/m，它们的平均相对误差分别为24.7%、8.9%和22.6%，虽然温跃层深度和强度的平均相对误差较大，但其绝对误差量值较小。而在南海海区，模型重构温跃层上界、下界和强度的平均绝对预报误差分别为4.1m、27.7m和0.007℃/m，它们的平均相对误差分别为16.1%、16.8%和9.5%，重构温跃层各特征值的平均相对误差都在20%以内。虽然南海区温跃层下界深度平均绝对预报误差较大，但相对于温跃层下界深度的空间尺度变化而言（平均温跃层下界深度为168m），平均相对误差仅为16.8%。因此说模型重构的温度廓线可以达到对我国陆架海域、深水区温跃层的较好估算。 基于对历史水文温度廓线观测资料的分析及自主温跃层统计预报模型，研制了实时可利用微机简单、快捷地进行温跃层估算及查询的可视化系统，这是迄今进行大范围海域温跃层统计与实时预报研究的较系统成果。

The unique skeleton of siliceous sponges (Porifera; Hexactinellida and Demospongiae) that evolved first from the Urmetazoa during the Proterozoic: a review

Sponges (phylum Porifera) had been considered as an enigmatic phylum, prior to the analysis of their genetic repertoire/tool kit. Already with the isolation of the first adhesion molecule, galectin, it became clear that the sequences of sponge cell surface receptors and of molecules forming the intracellular signal transduction pathways triggered by them, share high similarity with those identified in other metazoan phyla. These studies demonstrated that all metazoan phyla, including Porifera, originate from one common ancestor, the Urmetazoa. The sponges evolved prior to the Ediacaran-Cambrian boundary (542 million years ago [myr]) during two major "snowball earth events", the Sturtian glaciation (710 to 680 myr) and the Varanger-Marinoan ice ages (605 to 585 myr). During this period the ocean was richer in silica due to the silicate weathering. The oldest sponge fossils (Hexactinellida) have been described from Australia, China and Mongolia and are thought to have existed coeval with the diverse Ediacara fauna. Only little younger are the fossils discovered in the Sansha section in Hunan (Early Cambrian; China). It has been proposed that only the sponges possessed the genetic repertoire to cope with the adverse conditions, e.g. temperature-protection molecules or proteins protecting them against ultraviolet radiation. The skeletal elements of the Hexactinellida (model organisms Monorhaphis chuni and Monorhaphis intermedia or Hyalonema sieboldi) and Demospongiae (models Suberites domuncula and Geodia cydonium), the spicules, are formed enzymatically by the anabolic enzyme silicatein and the catabolic enzyme silicase. Both, the spicules of Hexactinellida and of Demospongiae, comprise a central axial canal and an axial filament which harbors the silicatein. After intracellular formation of the first lamella around the channel and the subsequent extracellular apposition of further lamellae the spicules are completed in a net formed of collagen fibers. The data summarized here substantiate that with the finding of silicatein a new aera in the field of bio/inorganic chemistry started. For the first time strategies could be formulated and experimentally proven that allow the formation/synthesis of inorganic structures by organic molecules. These findings are not only of importance for the further understanding of basic pathways in the body plan formation of sponges but also of eminent importance for applied/commercial processes in a sustainable use of biomolecules for novel bio/inorganic materials.

The manganese superoxide dismutase gene in bay scallop Argopecten irradians: Cloning, 3D modelling and mRNA expression

A novel manganese superoxide dismutase (MnSOD) was cloned from bay scallop Argopecten irradians by 3' and 5' rapid amplification of cDNA ends (RACE) PCR. The full-length cDNA of MnSOD was of 1207 bp with a 678 bp open reading frame encoding 226 amino acids. The deduced amino acid sequence contained a putative signal peptide of 26 amino acids. Sequence comparison showed that the MnSOD of A. irradians shared high identity with MnSOD in invertebrates and vertebrates, such as MnSOD from abalone Haliotis discus discus (ABG88843) and frog Xenopus laevis (AAQ63483). Furthermore, the 3D structure of bay scallop MnSOD was predicted by SWISS-MODEL Protein Modelling Server and compared with those of other MnSODs. The overall structure of bay scallop MnSOD was similar to those of zebrafish Danio rerio, fruit fly Drosophila melanogaster, Chinese shrimp Fenneropenaeus chinensis, human Homo sapiens, and had the highest similarity to scallop Mizuhopecten yessoensis and abalone H. discus discus. A quantitative real-time PCR (qRT-PCR) assay was developed to detect the mRNA expression of MnSOD in different tissues and the temporal expression in haemocytes following challenge with the bacterium Vibrio anguillarum. A higher-level of mRNA expression of MnSOD was detected in gill and mantle. The expression of MnSOD reached the highest level at 3 h post-injection with V. anguillarum and then slightly recovered from 6 to 48 h. The results indicated that bay scallop MnSOD was a constitutive and inducible protein and thus could play an important role in the immune responses against V anguillarum infection. (c) 2008 Elsevier Ltd. All rights reserved.

Isolation and characterization of a novel variant of HMW glutenin subunit gene from the St genome of Pseudoroegneria stipifolia

The x- and y-type high molecular weight (HMW) glutenin subunits are conserved seed storage proteins in wheat and related species. Here we describe investigations on the HMW glutenin subunits from several Pseudoroegneria accessions. The electrophoretic mobilities of the HMW glutenin subunits from Pd. stipifolia, Pd tauri and Pd strigosa were much faster than those of orthologous wheat subunits, indicating that their protein size may be smaller than that of wheat subunits. The coding sequence of the Glu-1St1 subunit (encoded by the Pseudoroegneria stipifolia accession PI325181) was isolated, and found to represent the native open reading frame (ORF) by in vitro expression. The deduced amino acid sequence of Glu-1St1 matched with that determined from the native subunit by mass spectrometric analysis. The domain organization in Glu-1St1 showed high similarity with that of typical HMW glutenin subunits. However, Glu-1St1 exhibited several distinct characteristics. First, the length of its repetitive domain was substantially smaller than that of conventional subunits, which explains its much faster electrophoretic mobility in SDS-PAGE. Second, although the N-terminal domain of Glu-1St1 resembled that of y-type subunit, its C-terminal domain was more similar to that of x-type subunit. Third, the N- and C-terminat domains of Glu-1St1 shared conserved features with those of barley D-hordein, but the repeat motifs and the organization of its repetitive domain were more similar to those of HMW glutenin subunits than to D-hordein. We conclude that Glu-1St1 is a novel variant of HMW glutenin subunits. The analysis of Glu-1St1 may provide new insight into the evolution of HMW glutenin subunits in Triticeae species. (C) 2007 Elsevier Ltd. All rights reserved.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

中低纬电离层暴的统计分析与经验建模

As a key issue of ionospheric weather study, systemic studies on ionospheric storms can not only further improve our understanding of the response of the ionosphere to solar and geomagnetic disturbances, but also help us to reveal the chemical, dynamic and electro-dynamic mechanisms during storms. Empirical modelling for regional ionospheric storm is also very useful, because it can provide us with tools and references for the forecasting and further practical application of ionospheric activity. In this thesis, we focus on describing and forecasting of ionospheric storms at middle and low latitudes. The main points of my investigations are listed as follows. (1) By using magnetic storms during the period over 50 years, the dependence of the type, onset time and time delay of the ionospheric storms on magnetic latitude, season and local time at middle and low latitudes in the East-Asian sector are studied. The results show that the occurrences of the types of ionospheric disturbances differ in latitude and season. The onset of the ionospheric storms depends on local time. At middle latitudes, most negative phase onsets are within the local time interval from night to early morning, and they rarely occurred in the local noon and afternoon sectors. At low latitudes, positive phases commence most frequently in the daytime sector as well as pre-midnight sector. The average time delays for both the positive and negative ionospheric storms increase with descending latitudes. The time delay has significant dependence on the local time of main phase onset (MPO). The time delay of positive response is shorter for daytime MPO and longer for night-time MPO, whereas the opposite applies for negative response. (2) Based on some previous researches, a primary empirical model for mid-latitude ionospheric disturbance is set up. By fitting to the observed data, we get a high accuracy with a mean RMSE of only 12-14% in summer and equinox. The model output has been compared with the output of STORM model, and the results show that, our model is much better than STORM in summer and a little better for some mid-latitude stations at equinox. Especially, for the type of two-step geomagnetic storm, our model can present twice descending of foF2 very well. In addition, our model can forecast positive ionospheric storms.

通过第一性原理计算研究内地核的结构与成分

The determination of the composition and structure of the Earth’s inner core has long been the major subject in the study of the Earth’s deep interior. It’s widely believed that the Earth’s core is formed by iron with a fraction of nickel. However, light elements must exist in the inner core because the earth core is less dense than pure iron-nickel alloy (~2-3% in the solid inner core and ~6-7% in the liquid outer core). The questions are what and how much light element is there in the iron-nickel alloy. Besides the composition, the crystal structure of the iron with or without light element is also not well known. According to the seismological observations, the sound waves propagate 3-4% faster along the spin axis than in the equatorial plane. That means the inner core is anisotropic. The densest structure of iron-nickel alloy should be h.c.p structure under the very high pressures. However, the h,c,p structure does not propagate waves anisotropic ally. Then what is the structure of the iron-nickel alloy or the iron-nickle-light element alloy. In this study, we tried to predict the composition and the structure of the inner core through ab initio calculation of the Gibbs free energy, which is a function of internal energy, density and entropy. We conclude that the h.c.p structure is more stable than the b.c.c structure under high pressure and 0 K, but with the increase of temperature, the free energy of the b.c.c structure is decreasing much faster than the h.c.p structure caused by the vibration of the atomics, so the b.c.c structure is more stable at high temperatures. With the addition of light elements (S or Si or both), the free energy of b.c.c. decreases even faster, about 3at% of Si not only explains why the inner core is about 2-3 % lighter than the iron-nickle alloy, but also reasons why the inner core is anisotropic, since the b.c.c. structure becomes more stable than the h.c.p structure at 5500-6000K and b.c.c. is anisotropic in propagating seismic waves. Therefore, we infer that the inner core of the earth is formed by b.c.c iron and a fraction of nickel plus ~3at.% Si, with a temperature higher than 5500K, which is consistent with the studies from other approaches.