Preparation and Physical Properties of LLDPE Grafted with Novel Nonionic Surfactants

Copolymers of linear low-density polyethylene (LLDPE) grafted with two novel nonionic surfactants, acrylic glycerol monostearate ester (AGMS) and acrylic polyoxyethylenesorbitan monooleate ester (ATW-EEN80), containing hydrophilic and hydrophobic groups and 1-olefin double bond were prepared by using a plasticorder at 190 degrees C. To evaluate the grafting degree, two different approaches based on H-1-NMR data were proposed, and FTIR calibration was showed to validate these methods. The rheological response of the molten polymers, determined under dynamic shear flow at small-amplitude oscillations, indicated that crosslinking formation of the chains could be decreased with increasing the monomer concentration. Their thermal behavior was studied by DSC and polarization microscope (PLM): The crystallization temperature (T-C) of grafted LLDPE shifted to higher temperature compared with neat LLDPE because the grafted chains acted as nucleating agents. Water and glycerol were used to calculate the surface free energy of grafted LLDPE films.

Robustness and Coherence of a Three-Protein Circadian Oscillator: Landscape and Flux Perspectives

Three-protein circadian oscillations in cyanobacteria sustain for weeks. To understand how cellular oscillations function robustly in stochastic fluctuating environments, we used a stochastic model to uncover two natures of circadian oscillation: the potential landscape related to steady-state probability distribution of protein concentrations; and the corresponding flux related to speed of concentration changes which drive the oscillations. The barrier height of escaping from the oscillation attractor on the landscape provides a quantitative measure of the robustness and coherence for oscillations against intrinsic and external fluctuations. The difference between the locations of the zero total driving force and the extremal of the potential provides a possible experimental probe and quantification of the force from curl flux. These results, correlated with experiments, can help in the design of robust oscillatory networks.

Single molecule dynamics and statistical fluctuations of gene regulatory networks: A repressilator

The authors developed a time dependent method to study the single molecule dynamics of a simple gene regulatory network: a repressilator with three genes mutually repressing each other. They quantitatively characterize the time evolution dynamics of the repressilator. Furthermore, they study purely dynamical issues such as statistical fluctuations and noise evolution. They illustrated some important features of the biological network such as monostability, spirals, and limit cycle oscillation. Explicit time dependent Fano factors which describe noise evolution and show statistical fluctuations out of equilibrium can be significant and far from the Poisson distribution. They explore the phase space and the interrelationships among fluctuations, order, amplitude, and period of oscillations of the repressilators. The authors found that repressilators follow ordered limit cycle orbits and are more likely to appear in the lower fluctuating regions. The amplitude of the repressilators increases as the suppressing of the genes decreases and production of proteins increases. The oscillation period of the repressilators decreases as the suppressing of the genes decreases and production of proteins increases.

PHYLOGEOGRAPHIC PATTERNS INDICATE TRANSATLANTIC MIGRATION FROM EUROPE TO NORTH AMERICA IN THE RED SEAWEED CHONDRUS CRISPUS (GIGARTINALES, RHODOPHYTA)1

Inferring how the Pleistocene climate oscillations have repopulated the extant population structure of Chondrus crispus Stackh. in the North Atlantic Ocean is important both for our understanding of the glacial episode promoting diversification and for the conservation and development of marine organisms. C. crispus is an ecologically and commercially important red seaweed with broad distributions in the North Atlantic. Here, we employed both partial mtDNA Cox1 and nrDNA internal transcribed spacer region 2 (ITS2) sequences to explore the genetic structure of 17 C. crispus populations from this area. Twenty-eight and 30 haplotypes were inferred from these two markers, respectively. Analysis of molecular variance (AMOVA) and of the population statistic Theta(ST) not only revealed significant genetic structure within C. crispus populations but also detected significant levels of genetic subdivision among and within populations in the North Atlantic. On the basis of high haplotype diversity and the presence of endemic haplotypes, we postulate that C. crispus had survived in Pleistocene glacial refugia in the northeast Atlantic, such as the English Channel and the northwestern Iberian Peninsula. We also hypothesize that C. crispus from the English Channel refugium repopulated most of northeastern Europe and recolonized northeastern North America in the Late Pleistocene. The observed phylogeographic pattern of C. crispus populations is in agreement with a scenario in which severe Quaternary glaciations influenced the genetic structure of North Atlantic marine organisms with contiguous population expansion and locally restricted gene flow coupled with a transatlantic dispersal in the Late Pleistocene.

A Holocene paleotemperature record based on radiolaria from the northern Okinawa Trough (East China Sea)

Using a radiolarian-based transfer function, mean annual sea surface temperature (SST) and seasonal temperature range are reconstructed through the last 10,500yrs in the northern Okinawa Trough. Down-core SST estimates reveal that throughout the Holocene the changes of mean annual SST display a three-step trend: (i) an early Holocene continuous warming between 10,500 and 8500yr BP which ends up with a abrupt cooling at about 8200yr BP; (ii) a relatively stable middle Holocene with high SST that lasted until 3200yr BP; and (iii) a late-Holocene distinct SST decline between 3200 and 500yr BP. This pattern is in agreement with the ice core and terrestrial paleoclimatic records in the Chinese continent and other regions of the world. Five cooling events with abrupt mean annual SST drops, which occur at similar to 300-600, 1400, 3100, 4600-5100 and 8200yr BP, are recognized during the last 10,500yrs. Comparison of our results with records of GISP2 ice core and marine sediment in North Atlantic region suggests these cooling events are strongly coupled, which implies a possible significant climatic correlation between high- and low-latitude areas. (C) 2007 Elsevier Ltd and INQUA. All rights reserved.

Paleoenvironmental change in the middle Okinawa Trough since the last deglaciation: Evidence from the sedimentation rate and planktonic foraminiferal record

Well-dated, high-resolution records of planktonic foraminifera and oxygen isotopes from two sediment cores, A7 and E017, in the middle Okinawa Trough reveal strong and rapid millennial-scale climate changes since similar to 18 to 17 thousand years before present (kyr B.P.). Sedimentation rate shows a sudden drop at similar to 11.2 cal. kyr B.P. due to a rapid rise of sea level after the Younger Dryas (YD) and consequently submergence of the large continental shelf on the East China Sea (ECS) and the retreat of the estuary providing sediment to the basin. During the last deglaciation, the relative abundance of warm and cold species of planktonic foraminifera fluctuates strongly, consistent with the timing of sea surface temperature (SST) variations determined from Mg/Ca measurements of planktonic foraminifera from one of the two cores. These fluctuations are coeval with climate variation recorded in the Greenland ice cores and North Atlantic sediments, namely Heinrich event 1 (H1), Bolling-Allerod (B/A) and YD events. At about 9.4 kyr B.P., a sudden change in the relative abundance of shallow to deep planktonic species probably indicates a sudden strengthening of the Kuroshio Current in the Okinawa Trough, which was synchronous with a rapid sea-level rise at 9.5-9.2 kyr B.P. in the ECS, Yellow Sea (YS) and South China Sea (SCS). The abundance of planktonic foraminiferal species, together with Mg/Ca based SST, exhibits millennial-scale oscillations during the Holocene, with 7 cold events (at about 1.7, 2.3-4.6, 6.2, 7.3, 8.2, 9.6, 10.6 cal. kyr BP) superimposed on a Holocene warming trend. This Holocene trend, together with centennial-scale SST variations superimposed on the last deglacial trend, suggests that both high and low latitude influences affected the climatology of the Okinawa Trough. (c) 2006 Elsevier B.V. All rights reserved.

Dynamics of intraseasonal sea level and thermocline variability in the equatorial Atlantic during 2002-03

Satellite and in situ observations in the equatorial Atlantic Ocean during 2002-03 show dominant spectral peaks at 40-60 days and secondary peaks at 10-40 days in sea level and thermocline within the intraseasonal period band (10-80 days). A detailed investigation of the dynamics of the intraseasonal variations is carried out using an ocean general circulation model, namely, the Hybrid Coordinate Ocean Model (HYCOM). Two parallel experiments are performed in the tropical Atlantic Ocean basin for the period 2000-03: one is forced by daily scatterometer winds from the Quick Scatterometer (QuikSCAT) satellite together with other forcing fields, and the other is forced by the low-passed 80-day version of the above fields. To help in understanding the role played by the wind-driven equatorial waves, a linear continuously stratified ocean model is also used. Within 3 degrees S-3 degrees N of the equatorial region, the strong 40-60-day sea surface height anomaly (SSHA) and thermocline variability result mainly from the first and second baroclinic modes equatorial Kelvin waves that are forced by intraseasonal zonal winds, with the second baroclinic mode playing a more important role. Sharp 40-50-day peaks of zonal and meridional winds appear in both the QuikSCAT and Pilot Research Moored Array in the Tropical Atlantic (PIRATA) data for the period 2002-03, and they are especially strong in 2002. Zonal wind anomaly in the central-western equatorial basin for the period 2000-06 is significantly correlated with SSHA across the equatorial basin, with simultaneous/ lag correlation ranging from-0.62 to 0.74 above 95% significance. Away from the equator (3 degrees-5 degrees N), however, sea level and thermocline variations in the 40-60-day band are caused largely by tropical instability waves (TIWs). On 10-40-day time scales and west of 10 degrees W, the spectral power of sea level and thermocline appears to be dominated by TIWs within 5 degrees S-5 degrees N of the equatorial region. The wind-driven circulation, however, also provides a significant contribution. Interestingly, east of 10 W, SSHA and thermocline variations at 10 40- day periods result almost entirely from wind-driven equatorial waves. During the boreal spring of 2002 when TIWs are weak, Kelvin waves dominate the SSHA across the equatorial basin (2 degrees S-2 degrees N). The observed quasi-biweekly Yanai waves are excited mainly by the quasi-biweekly meridional winds, and they contribute significantly to the SSHA and thermocline variations in 1 degrees-5 degrees N and 1 degrees-5 degrees S regions.

Self-organized Criticality Model for Ocean Internal Waves

In this paper, we present a simple spring-block model for ocean internal waves based on the self-organized criticality (SOC). The oscillations of the water blocks in the model display power-law behavior with an exponent of -2 in the frequency domain, which is similar to the current and sea water temperature spectra in the actual ocean and the universal Garrett and Munk deep ocean internal wave model [Geophysical Fluid Dynamics 2(1972) 225; J. Geophys. REs. 80 (1975) 291]. The influence of the ratio of the driving force to the spring coefficient to SOC behaviors in the model is also discussed.

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.

Roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean

An ocean general circulation model (OGCM) is used to study the roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean. The western boundary reflection is defined as the total Kelvin waves leaving the western boundary, which include the reflection of the equatorial Rossby waves as well as the effects of alongshore winds, off-equatorial Rossby waves, and nonlinear processes near the western boundary. The evaluation of the reflection is based on a wave decomposition of the OGCM results and experiments with linear models. It is found that the alongshore winds along the east coast of Africa and the Rossby waves in the off-equatorial areas contribute significantly to the annual harmonics of the equatorial Kelvin waves at the western boundary. The semiannual harmonics of the Kelvin waves, on the other hand, originate primarily from a linear reflection of the equatorial Rossby waves. The dynamics of a dominant annual oscillation of sea level coexisting with the dominant semiannual oscillations of surface zonal currents in the central equatorial Indian Ocean are investigated. These sea level and zonal current patterns are found to be closely related to the linear reflections of the semiannual harmonics at the meridional boundaries. Because of the reflections, the second baroclinic mode resonates with the semiannual wind forcing; that is, the semiannual zonal currents carried by the reflected waves enhance the wind-forced currents at the central basin. Because of the different behavior of the zonal current and sea level during the reflections, the semiannual sea levels of the directly forced and reflected waves cancel each other significantly at the central basin. In the meantime, the annual harmonic of the sea level remains large, producing a dominant annual oscillation of sea level in the central equatorial Indian Ocean. The linear reflection causes the semiannual harmonics of the incoming and reflected sea levels to enhance each other at the meridional boundaries. In addition, the weak annual harmonics of sea level in the western basin, resulting from a combined effect of the western boundary reflection and the equatorial zonal wind forcing, facilitate the dominance by the semiannual harmonics near the western boundary despite the strong local wind forcing at the annual period. The Rossby waves are found to have a much larger contribution to the observed equatorial semiannual oscillations of surface zonal currents than the Kelvin waves. The westward progressive reversal of seasonal surface zonal currents along the equator in the observations is primarily due to the Rossby wave propagation.

大气季节内振荡对东亚夏季风活动的影响研究

利用ERA40逐日再分析资料、NCEP/NCAR2逐日再分析资料、中国740个测站日降水资料、上海台风研究所提供的西太平洋热带气旋资料、Kaplan等重建的月平均SSTA资料、NOAA逐日长波辐射（OLR）等资料，应用离散功率谱分析、带通滤波、EOF分析等统计方法，研究了东亚夏季风(EASM）的移动特征、东亚地区季节内振荡(ISO)的基本特征、季节内振荡对东亚夏季风活动的影响、季节内振荡对东亚夏季风异常活动的影响机理。主要结论如下： （1）综合动力和热力因素定义了可动态描述东亚夏季风移动和强度的指数，并利用该指数研究了东亚夏季风的爆发和移动的季节内变化及其年际和年代际变化特征。研究发现，气候平均东亚夏季风前沿分别在28候、33候、36候、38候、40候、44候出现了明显的跳跃。东亚夏季风活动具有显著的年际变率，主要由于季风前沿在某些区域异常停滞和突然跨越北跳或南撤引起，造成中国东部旱涝灾害频繁发生。东亚夏季风的活动具有明显的年代际变化，在1965年、1980年、1994年发生了突变，造成中国东部降水由“南旱北涝”向“南涝北旱”的转变。 （2）东亚季风区季节内变化具有10~25d和30~60d两个波段的季节内振荡周期，以30-60d为主。存在三个主要低频模态，第一模态主要表征了EASM在长江中下游和华北地区活动期间的低频形势；第二模态印度洋-菲律宾由低频气旋式环流控制，主要表现了ISO在EASM爆发期间的低频形势；第三模态主要出现在EASM在华南和淮河活动期间的低频形势。第一模态和第三模态是代表东亚夏季风活动异常的主要低频形势。 （3）热带和副热带地区ISO总是沿垂直切变风的垂直方向传播。因此，在南海-菲律宾东北风垂直切变和副热带西太平洋北风垂直切变下，大气热源激发菲律宾附近交替出现的低频气旋和低频反气旋不断向西北传播，副热带西太平洋ISO以向西传播为主。中高纬度地区，乌拉尔山附近ISO以向东、向南移动或局地振荡为主；北太平洋中部ISO在某些情况下向南、向西传播。 （4）季风爆发期，伴随着热带东印度洋到菲律宾一系列低频气旋和低频反气旋， 冷空气向南输送，10~25天和30~60天季节内振荡低频气旋同时传入南海加快了南海夏季风的爆发。在气候态下，ISO活动表现的欧亚- 太平洋（EAP）以及太平洋-北美（PNA）低频波列分布特征(本文提出的EAP和PNA低频波列与传统意义上的二维定点相关得到的波列不同)。这种低频分布形式使得欧亚和太平洋中高纬度的槽、脊及太平洋副热带高压稳定、加强，东亚地区的低频波列则成为热带和中高纬度ISO相互作用影响东亚夏季风活动的纽带。不同的阶段表现不同的低频模态，30~60d低频模态的转变加快了EASM推进过程中跳跃性；30-60d低频模态的维持使得EASM前沿相对停滞。 （5）30-60d滤波场，菲律宾海域交替出现的低频气旋和低频反气旋不断向西北传播到南海-西太平洋一带。当南海-西太平洋地区低频气旋活跃时，季风槽加强、东伸，季风槽内热带气旋（TC）频数增加；当南海-西太平洋低频反气旋活跃时，季风槽减弱、西退，TC处于间歇期，生成位置不集中。 （6）在El Nino态下，大气季节内振荡偏弱，北传特征不明显，但ISO由中高纬度北太平洋中部向南和副热带西太平洋向西的传播特征显著，东亚地区ISO活动以第三模态为主，EASM集中停滞在华南和淮河流域，常伴随着持续性区域暴雨的出现，易造成华南和江淮流域洪涝灾害，长江和华北持续干旱。在La Nina态下，大气季节内振荡活跃，且具有明显的向北传播特征，PNA低频波列显著，东亚地区ISO活动以第一模态单峰为主；EASM主要停滞在长江中下游和华北地区，这些地区出现异常持续强降水，华南和淮河流域多干旱；在El Nino态向La Nina态转换期，ISO活动以第一模态双峰为主，长江中下游常常出现二度梅。

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

温度跃层是反映海洋温度场的重要物理特性指标，对水下通讯、潜艇活动及渔业养殖、捕捞等有重要影响。本文利用中国科学院海洋研究所“中国海洋科学数据库”在中国近海及西北太平洋（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%。因此说模型重构的温度廓线可以达到对我国陆架海域、深水区温跃层的较好估算。 基于对历史水文温度廓线观测资料的分析及自主温跃层统计预报模型，研制了实时可利用微机简单、快捷地进行温跃层估算及查询的可视化系统，这是迄今进行大范围海域温跃层统计与实时预报研究的较系统成果。

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed-dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3-4) and paternally inherited, pollen- and seed-dispersed chloroplast DNA (cpDNA) (rpl16 and trnS-trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G(ST) = 0.738, N-ST = 0.771) indicating low levels of seed-based gene flow and significant phylogeographical structure (N-ST > G(ST), P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north-central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north-central part of the species' range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north-central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among-population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.

Diploid and tetraploid distribution of Allium przewalskianum Regel. (Liliaceae) in the Qinghai-Tibetan Plateau and adjacent regions

Polyploidization plays an important role in generating the current high diversity of plants. Studies of distributional patterns of diploids and derivative autopolyploids have provided important insights into evolutionary processes and cryptic speciation of polyploidization within species defined on the basis of their morphology. However, few studies have been designed to examine distributions of infrageneric diploids and polyploids on the Asian Qinghai-Tibetan Plateau (QTP). Allium przewalskianum occurs widely on the QTP and in adjacent regions, at altitudes ranging from 2000m to 4500m. We collected a total of 844 individuals from 62 populations and determined their cytotypes over the entire distribution range of this species. Tetraploids tend to occur at high altitudes; however, the positive relationship between the ploidy and altitude was only marginally significant (P < 0.05). Contact zones between diploids and tetraploids were recorded on the eastern QTP from north to south. Four populations were found to harbor both cytotypes, but no triploid individuals. The wider distribution of tetraploids may be mainly due to their greater colonization ability in the new niches created by the Quaternary climatic oscillations in the QTP region. Our results offer a fundamental framework for studying evolutionary origins, adaptations and cryptic divergences of polyploids within this species complex based on molecular and/or ecological examinations in the future.

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N-ST = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N-ST > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N-ST = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.