Grain-size records at ODP site 1146 from the northern South China sea: Implications on the east Asian monsoon evolution since 20 Ma

273 samples from Ocean Drilling Program (ODP) Site 1146 in the northern South China Sea (SCS) were analyzed for grain-size distributions using grain-size class vs. standard deviation method and end-member modeling algorithm (EMMA) in order to investigate the evolution of the East Asian mon-soon since about 20 Ma. 10-19 mu m/1.3-2.4 mu m, the ratio of two grain-size populations with the highest variability through time was used to indicate East Asian winter monsoon intensity relative to summer monsoon. The mass accumulation rate of the coarsest end member EM1 (eolian), resulting from EMMA, can be used as a proxy of winter monsoon strength and Asian inland aridity, and the ratio of EM1/(EM2+EM3) as a proxy of winter monsoon intensity relative to summer monsoon. The combined proxies show that a profound enhancement of East Asian winter monsoon strength and winter monsoon intensity relative to summer monsoon occurred at about 8 Ma, and it is possible that the summer monsoon simultaneously intensified with winter monsoon at 3 Ma. Our results are well consistent with the previous studies in loess, eolian deposion in the Pacifc, radiolarians and planktonic foraminifera in the SCS. The phased uplift of the Himalaya-Tibetan Plateau may have played a significant role in strengthening the Asian monsoon at 8 Ma and 3 Ma.

Generation and evolution of magma beneath the East Pacific Rise: Constraints from U-series disequilibrium and plagioclase-hosted melt inclusions

Mid-ocean ridge basalt (MORB) samples from the East Pacific Rise (EPR 12 degrees 50'N) were analyzed for U-series isotopes and compositions of plagioclase-hosted melt inclusions. The Ra-226 and Th-230 excesses are negatively correlated; the Ra-226 excess is positively correlated with Mg# and Sm/Nd, and is negatively correlated with La/Sm and Fe-8; the Th-230 excess is positively correlated with Fe-8 and La/Sm and is negatively correlated with Mg# and Sm/Nd. Interpretation of these correlations is critical for understanding the magmatic process. There are two models (the dynamic model and the "two-porosity" model) for interpreting these correlations, however, some crucial parameters used in these models are not ascertained. We propose instead a model to explain the U-series isotopic compositions based on the control of melt density variation. For melting either peridotite or the "marble-cake" mantle, the FeOt content, Th-230 excess and La/Sm ratio increases and Sm/Nd decreases with increasing pressure. A deep melt will evolve to a higher density and lower Mg# than a shallow melt, the former corresponds to a long residence time, which lowers the Ra-226 excess significantly. This model is supported by the existence of low Ra-226 excesses and high Th-230 excesses in MORBs having a high Fe-8 content and high density. The positive correlation of Ra-226 excess and magma liquidus temperature implies that the shallow melt is cooled less than the deep melt due to its low density and short residence time. The correlations among Fe-8, Ti-8 and Ca-8/Al-8 in plagioclase-hosted melt inclusions further prove that MORBs are formed from melts having a negative correlation in melting depths and degrees. The negative correlation of Ra-226 excess vs. chemical diversity index (standard deviation of Fe-8, Ti-8 and Ca-8/Al-8) of the melt inclusions is in accordance with the influence of a density-controlled magma residence time. We conclude that the magma density variation exerts significant control on residence time and U-series isotopic compositions. (c) 2010 Elsevier B.V. All rights reserved.

Equivalent-barotropic definition of tropospheric mean temperature

An equivalent-barotropic (EB) description of the tropospheric temperature field is derived from the geostrophic empirical mode (GEM) in the form of a scalar function Gamma(p, phi), where p is pressure and phi is 300-850-mb thickness. Baroclinic parameter phi plays the role of latitude at each longitudinal section. Compared with traditional Eulerian-mean methods, GEM defines a mean field in baroclinic streamfunction space with a time scale much longer than synoptic variability. It prompts an EB concept that is only based on a baroclinic field. Monthly GEM fields are diagnosed from NCEP-NCAR reanalysis data and account for more than 90% of the tropospheric thermal variance. The circumglobal composite of GEM fields exhibits seasonal, zonal, and hemispheric asymmetries, with larger rms errors occurring in winter and in the Northern Hemisphere (NH). Zonally asymmetric features and planetary deviation from EB are seen in the NH winter GEM. Reconstruction of synoptic sections and correlation analysis reveal that the tropospheric temperature field is EB at the leading order and has a 1-day phase lag behind barotropic variations in extratropical regions.

Mean sea surface heights of the South and East China Seas from ocean circulation model and geodetic leveling

The mean sea surface heights (sea surface topography) of the South China, East China, Yellow and Bohai Seas are derived from an ocean general circulation model and surface air pressure. The circulation model covers the global oceans, with fine grid (1/6degrees) covering the East Asian marginal seas and coarse grid (31) covering the rest part of the global oceans. The result shows that the China 1985 National Altitude Datum is 24.7 cm above the me-an sea surface height of the world oceans. The mean sea surface in the coastal ocean adjacent to China is higher in the south than in the north. Intercomparison of the model results with the geodetic leveling measurements at 28 coastal tidal stations shows a standard deviation of 4.8 cm and a fitting coefficient of 95.3%. After correction through linear regression, the standard deviation is reduced to 4.5 cm. This indicates that the accuracy of model results is sufficient for practical application. Based on the model results, the mean sea surface heights for the study area with a resolution of 1/6 degree are given. This result also links the mean sea levels at islands with those on the mainland coast and gives the mean sea surface heights at tidal stations in the Taiwan Island, the Dongsha Islands, the Yisha Islands and the Nansha Islands relative to the China 1985 National Altitude Datum.

Intraseasonal variability of Indian Ocean sea surface temperature during boreal winter: Madden-Julian Oscillation versus submonthly forcing and processes

[ 1] Intraseasonal variability of Indian Ocean sea surface temperature (SST) during boreal winter is investigated by analyzing available data and a suite of solutions to an ocean general circulation model for 1998 - 2004. This period covers the QuikSCAT and Tropical Rainfall Measuring Mission (TRMM) observations. Impacts of the 30 - 90 day and 10 - 30 day atmospheric intraseasonal oscillations (ISOs) are examined separately, with the former dominated by the Madden-Julian Oscillation (MJO) and the latter dominated by convectively coupled Rossby and Kelvin waves. The maximum variation of intraseasonal SST occurs at 10 degrees S - 2 degrees S in the wintertime Intertropical Convergence Zone (ITCZ), where the mixed layer is thin and intraseasonal wind speed reaches its maximum. The observed maximum warming ( cooling) averaged over ( 60 degrees E - 85 degrees E, 10 degrees S - 3 degrees S) is 1.13 degrees C ( - 0.97 degrees C) for the period of interest, with a standard deviation of 0.39 degrees C in winter. This SST change is forced predominantly by the MJO. While the MJO causes a basin-wide cooling ( warming) in the ITCZ region, submonthly ISOs cause a more complex SST structure that propagates southwestward in the western-central basin and southeastward in the eastern ocean. On both the MJO and submonthly timescales, winds are the deterministic factor for the SST variability. Short-wave radiation generally plays a secondary role, and effects of precipitation are negligible. The dominant role of winds results roughly equally from wind speed and stress forcing. Wind speed affects SST by altering turbulent heat fluxes and entrainment cooling. Wind stress affects SST via several local and remote oceanic processes.

Statistical distribution of nonlinear random wave height

A statistical model of random wave is developed using Stokes wave theory of water wave dynamics. A new nonlinear probability distribution function of wave height is presented. The results indicate that wave steepness not only could be a parameter of the distribution function of wave height but also could reflect the degree of wave height distribution deviation from the Rayleigh distribution. The new wave height distribution overcomes the problem of Rayleigh distribution that the prediction of big wave is overestimated and the general wave is underestimated. The prediction of small probability wave height value of new distribution is also smaller than that of Rayleigh distribution. Wave height data taken from East China Normal University are used to verify the new distribution. The results indicate that the new distribution fits the measurements much better than the Rayleigh distribution.

潮波模型的优化开边界方法研究

数值模式是潮波研究的一种有利手段，但在研究中会面临各种具体问题，包括开边界条件的确定、底摩擦系数和耗散系数的选取等。数据同化是解决这些问题的一种途径，即利用有限数量的潮汐观测资料对潮波进行最优估计，其根本目的是迫使模型预报值逼近观测值，使模式不要偏离实际情况太远。本文采用了一种优化开边界方法，沿着数值模型的开边界优化潮汐水位信息，目的是设法使数值解在动力约束的意义下接近观测值，获得研究区域的潮汐结果。边界值由指定优化问题的解来定，以提高模拟区域的潮汐精度，最优问题的解是基于通过开边界的能量通量的变化，处理开边界处的观测值与计算值之差的最小化。这里提供了辐射型边界条件，由Reid 和Bodine（本文简称为RB）推导，我们将采用的优化后的RB方法（称为ORB）是优化开边界的特殊情况。 本文对理想矩形海域（ E- E， N- N， 分辨率 ）进行了潮波模拟，有东部开边界，模式采用ECOM3D模式。对数据结果的误差分析采用，振幅平均偏差，平均绝对偏差，平均相对误差和均方根偏差四个值来衡量模拟结果的好坏程度。 需要优化入开边界的解析潮汐值本文采用的解析解由方国洪《海湾的潮汐与潮流》（1966年）方法提供，为验证本文所做的解析解和方文的一致，本文做了其第一个例子的关键值a,b,z,结果与其结果吻合的相当好。但略有差别，分析的可能原因是两法在具体迭代方案和计算机保留小数上有区别造成微小误差。另外，我们取m=20，得到更精确的数值，我们发现对前十项的各项参数值，取m=10,m=20各项参数略有改进。当然我们可以获得m更大的各项参数值。 同时为了检验解析解的正确性讨论m和l变化对边界值的影响，结果指出，增大m，m=20时，u的模最大在本身u1或u2的模的6%；m=100时，u的模最大在本身u1或u2的模的4%；m再增大，m=1000时，u的模最大在本身u1或u2的模的4%，改变不大。当l<1时， =0处u的模最大为2。当l=1时， =0处u的模最大为0.1，当l>1时，l越大，u的模越小，当l=10时，u的模最大为0.001，可以认为为0。 为检验该优化方法的应用情况，我们对理想矩形区域进行模拟，首先将本文所采用的优化开边界方法应用于30m的情况，在开边界优化入开边界得出模式解，所得模拟结果与解析解吻合得相当好，该模式解和解析解在整个区域上，振幅平均绝对偏差为9.9cm,相位平均绝对偏差只有4.0 ，均方根偏差只有13.3cm，说明该优化方法在潮波模型中有效。 为验证该优化方法在各种条件下的模拟结果情况，在下面我们做了三类敏感性试验： 第一类试验：为证明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，我们来比较ORB条件与RB条件的优劣，我们模拟用了两个不同的摩擦系数，k分别为：0,0.00006。 结果显示，针对不同摩擦系数，显示在开边界上使用ORB条件的解比使用RB条件的解无论是振幅还是相位都有显著改善，两个试验均方根偏差优化程度分别为84.3%，83.7%。说明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，大大提高了模拟水平。上述的两个试验得出， k=0.00006优化结果比k=0的好。 第二类试验，使用ORB条件确定优化开边界情况下，在东西边界加入出入流的情况，流考虑线性和非线性情况，结果显示，加入流的情况，潮汐模拟的效果降低不少，流为1Sv的情况要比5Sv的情况均方根偏差相差20cm,而不加流的情况只有0.2cm。线性流和非线性流情况两者模式解相差不大，振幅，相位各项指数都相近, 说明流的线性与否对结果影响不大。 第三类试验，不仅在开边界使用ORB条件，在模式内部也使用ORB条件，比较了内部优化和不优化情况与解析解的偏差。结果显示，选用不同的k,振幅都能得到很好的模拟，而相位相对较差。另外，在内部优化的情况下，考虑不同的k的模式解， 我们选用了与解析解相近的6个模式解的k,结果显示，不同的k,振幅都能得到很好的模拟，而相位较差。 总之，在开边界使用ORB条件比使用RB条件好，振幅相位都有大幅度改进，在加入出入流情况下，流的大小对模拟结果有影响，但线形流和非线性流差别不大。内部优化的结果显示，模式采用不同的k都能很好模拟解析解的振幅。

Determination of tetrodotoxin by liquid chromatography coupled with electrospray ion-trap mass spectrometry

Two methods for tetrodotoxin analysis using liquid chromatography coupled with electrospray iontrap mass spectrometry (LC-ESI-MS) have been established with C,, reversed phase column and hydrophilic interaction liquid chromatography (HILIC) column, respectively. Sensitivity and reproducibility of the methods were compared. The method using C-18 column in selected ion monitoring (SIM) mode had a detection limit (S/N = 3) of 120 pg, and a good linearity of the calibration curve was obtained for tetrodotoxin (r = 0. 9992). High reproducibility of the method was observed, with a relative standard deviation (RSD) below 10%. The method using HILIC column in SIM mode and selected reaction monitoring (SRM) mode had detection limits (S/N = 3) of 15 and 3.75 pg, respectively. Good linearity of the calibration curves was obtained for tetrodotoxin (r = 0. 9996 and 0. 9998 in SIM and SRM mode, respectively). T he reproducibility was high in SIM mode but relatively poor in SRM mode. Based on the results, the method using HILIC column in SIM mode was suggested for the analysis of tetrodotoxin with LC-MS system.

Yellow perch strain evaluation I: Genetic variance of six broodstock populations

As a prelude to strain selection for domestication and future marker assisted selection, genetic variation revealed by microsatellite DNA was evaluated in yellow perch, Perca flavescens, from four wild North American populations collected in 2003-2004 (Maine, New York, North Carolina, and Pennsylvania,), and two captive populations (Michigan and Ohio). For the loci examined, levels of heterozygosity ranged from H-e=0.04 to 0.88, genetic differentiation was highly significant among all population pairs, and effective migration ranged from low (N(e)m=0.3) to high (N(e)m=4.5). Deviation from Hardy-Weinberg equilibrium was regularly observed indicating significant departures from random mating. Instantaneous measures of inbreeding within these populations ranged from near zero to moderate (median F=0.16) and overall inbreeding levels averaged F-IS=0.18. Estimates of genetic diversity, Phi(ST), and genetic distance were highest between Michigan and all other broodstock groups and lowest between New York and Ohio. Genetic differentiation among groups did not correlate with geographic distance. Overall, the patterns of variation exhibited by the captive (Michigan and Ohio) populations were similar to patterns exhibited by the other wild populations, indicating that spawning and management practices to date have not significantly reduced levels of genetic variation. (C) 2007 Elsevier B.V. All rights reserved.

Detection of carbohydrates using new labeling reagent 1-(2-naphthyl)-3-methyl-5-pyrazolone by capillary zone electrophoresis with absorbance (UV)

A novel labeling reagent 1-(2-naphthyl)-3-methyl-5-pyrazolone (NMP) coupled with capillary electrophoresis (CE) with DAD detection for the determination of carbohydrates has been developed. The chromophore in the 1-phenyl-3-methyl-5-pyrazolone (PMP) reagent is replaced by naphthyl functional group, which results in a reagent with very high molar absorptivity (epsilon(251nm) = 5.58 x 10(4) L mol(-1) cm(-1)). This pen-nits NMP-labeled carbohydrates to be detected with UV absorbance in standard 50-mu m-i.d. fused silica capillaries by zone electrophoresis. in this mode, nanomolar concentrations of detection limits are obtained. The method for the derivatization. of carbohydrates with NMP is simplified. The derivatization reaction is rapid and mild in the presence of ammonia catalyst without further transfer steps. Nine monosaccharide derivatives such as mannose, galacturonic acid, glucuronic acid, rhamnose, glucose, galactose, xylose, arabinose and fucose can successfully be detected in CE mode. Good reproducibility can be obtained with relative standard deviation (R.S.D.) values of the migration times and peak area, respectively, from 0.44 to 0.48 and from 3.2 to 4.8. Furthermore, the developed method has been successfully applied to the analysis of carbohydrates in the hydrolyzed rape bee pollen samples. (C) 2008 Published by Elsevier B.V.

轮式移动机器人轨迹跟踪的最优控制

深入分析了轮式移动机器人的运动状态,建立了WMR路径偏差系统的非线性数学模型。应用小偏差线性化理论,将该多输入多输出非线性系统简化成一个单输入单输出线性系统。然后基于线性二次型调节器理论进行了系统最优控制器的设计,并针对该理论中加权矩阵Q与R难以确定的问题,从控制效果出发,采用自适应遗传算法对其进行了优化。实现了移动机器人对预定轨迹的满意鲁棒跟踪,同时满足了实时性要求。实验结果证明了该方法的正确性与实用性。

AUV自主导航航位推算算法的分析研究

地球形状的不规则性,各种导航传感器本身的误差,以及仪器的安装偏差等,使得AUV(自治水下机器人)在进行远距离自主航行时,自主导航的精度大大下降。针对以上问题及实际工程需要,论文对AUV自主导航的航位推算算法做了进一步研究并加以改进,以提高其自主导航精度。最后,利用2004年中国科学院沈阳自动化所水下机器人研究中心进行AUV湖试所获得的数据,对文中提出的算法进行了验证。结果表明,AUV的自主导航精度得到大大提高,可以用于修正原来的自主导航算法。

二维流形三角网格模型顶点法矢量估计

提出一种以顶点的一邻域中三角形在该顶点处的顶角与对应三角形的面积比值加权三角面法矢量估计二维流形三角网格模型顶点法矢量的方法.回顾了现有的五种顶点法矢量估计方法,然后给出了新的方法.设计了利用理论法矢量与估计法矢量的夹角作为误差评价标准的实验,应用球体和椭球体模型分析了所涉及的6种估计方法的性能。

利用MC-ICP-MS测定硼同位素及其在西藏地热和盐湖中的初步应用

Based on previous studies, boron can be separated from aqueous samples with Amberlite IRA-743 resin. Experiments on the elute temperature, elute volume and the dynamic resin exchange capacity have been performed in this study. Results show that the dynamic exchange capacity of the resin is 4.2mg B/g and at room temperature, boron fixed on the resin within this capacity level can be extracted quantitatively by using 5ml 2%HNO3. A new procedure has been developed for the measurement of boron isotope ratios in water samples using a Neptune MC-ICP-MS, after resolving the memory effect, which is a key problem, and investigating the impacts of mass bias and Si matrix effect. Using this method, it usually takes 20 min to perform one measurement on 0.1ppm boron solution with a precision of 0.23‰ (SD). If the relative deviation between a sample and the standard is large, the washout time needs to be doubled to achieve a higher precision. δ11B values of water samples from Yangbajing geothermal field vary from -10.53 to -9.13‰. Owing to the large difference B concentration and the small B isotope difference between deep geothermal water and surface water, B isotope ratios of the shallow geothermal fluids are dominated by the deep end member rather than the shallower one in the mixing process. As a consequence, δ11B-B relation is indicative basically of a dilution process. Vapor-liquid separation and calcite scaling also greatly influence B isotope fractionation. δ11B values of water samples from Dagejia geothermal field are from -15.98‰ to -11.67‰. Boron in Changma River near the field has two sources, freshwater lakes (Dajiamang Lake and Canke Lake) and geothermal waters. Finally, a preliminary discussion is included on boron geochemical characteristics of the salt lakes in Shuanghu area and other geothermal fields, to provide information for future studies on boron isotope geochemistry of geothermal systems and salt lakes in Tibet.

高温高压下H2O-CH4体系PVTx性质的分子动力学模拟研究

Geological fluids are important components in the earth system. To study thephysical chemistry properties and the evolution of fluid system turns out to be one of the most challenging issues in geosciences. Besides the conventional experimental approaches and theoretical or semi-theoretical modeling, molecular level computer simulation(MLCS) emerges as an alternative tool to quantificationally study the physico-chemical properties of fluid under extreme conditions in order to find out the characteristics and interaction of geological fluids in and around earth. Based on our previous study of the intermolecular potential for pure H2O and thestrict evaluation of the competitive potential models for pure CH4 and the ab initio fitting potential surface across H2O-CH4 molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH4 and the H2O-CH4 mixtures. Comparison of 1941 simulations with experimental PVT data for pure CH4 shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the H2O-CH4 mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase,indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region. After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the H2O-CH4 system covering 673 to 2573 K and 0.01 to 10 GPa. Isochores for compositions < 4 mol% CH4 up to 773 K and 600 MPa are also determined in this thesis.