青藏高原东缘青杨（Populus cathayana Rehd.）遗传多样性研究

青杨（Populus cathayana Rehd.）是青杨派杨树的主要树种之一，为我国特有乡土树种，其主要分布区之一是我国的青藏高原，集中分布地带在甘肃省中部及青海省东部，四川省西北部岷江上游和松潘等地区。本研究以青藏高原东缘青杨天然分布区的6个群体143个个体为材料，用AFLP、SSR和叶绿体SSR分子标记分析青杨天然群体的遗传多样性，分析其遗传结构和分化，比较6个群体间遗传多样性的高低和群体间的遗传关系。旨在为青杨基因资源评价、保护与保存、遗传改良策略制定等提供科学理论依据。通过以上研究，得出如下主要研究结果： 1 AFLP分子标记研究结果 采用4对选择性引物对6个青杨天然群体143个个体进行分析，扩增谱带分析共检测到175个位点，其中173个位点表现为多态，多态位点百分率高达98.9%。从整体上表现出较高的遗传多样性，Nei’s基因多样度(h)水平为0.306。从青杨天然群体位点分布来看，有高达20%的位点（32位点）为群体所特有，仅有9.14%的位点（16位点）在所有群体中存在。群体间的遗传分化极大，所有遗传变异中，有48.9%的遗传变异存在于群体间。在个体群丛（Individuals cluster）和主坐标（PCO analysis）分析中，青杨各群体未呈现任何地理模式，Mantel检测也显示各群体间遗传距离与地理距离无明显相关。研究认为，由于地理和空间上大尺度的隔离和地形地貌复杂使得群体间无法进行基因交流，导致群体间遗传分化极大，另外各群体在不同的选择压力下，经历各自独立的进化历程，这些都可能导致群体间遗传距离与地理距离的不相关。 2 SSR分子标记研究结果 在SSR分析中，7个位点在6个青杨天然群体143个个体中共检测到79个等位基因，每位点检测到的等位基因数在5-16之间，平均11.3个，总体上多态位点百分率达100%。平均观察杂合度和期望杂合度分别为0.792和0.802。Hardy-Weinberg平衡检验表明青杨大部分群体都处于非平衡状态，群体大部分位点都是偏离哈迪-温伯格平衡（76.3%），只有23.7%的测验满足哈迪-温伯格平衡。分析青杨天然群体内和群体间的遗传变异，基因分化系数（GST）为0.373，即有62.7%的遗传变异存在群体内，37.3%的遗传变异存在群体间。群体内的遗传变异高于群体间水平。根据各群体遗传距离UPGMA聚类分析，有来自相临分布区、近似气候类型的群体聚在一起的趋势，但Mantel检测反映遗传距离与地理距离间并无明显相关性。 3 cpSSR分子标记研究结果 分析来自青藏高原东缘6个青杨天然群体，所用cpSSR引物中有5对cpSSR引物（CCMP2、CCMP5、SCUO01、SCU03、SCU07）都表现较高的多态性，单个引物检测的片段数都在4以上。5对cpSSR引物共检测片段数26个，组成了12种叶绿体DNA单倍型。各群体的单倍型分布和频率有较大差异，群体单倍型多样性范围为0-0.4926，TS、JZ、PW和SHY群体单倍型多样性高于QHY和LED群体水平。本研究发现，分布在青藏高原东缘的青杨天然群体，群体间不存在共享的单倍型，各群体间存在极大的遗传分化（GST=0.9223）。从青藏高原东缘地区经历的地质历史事件来看，第四纪的冰期气候变迁可能是造成青杨现今遗传结构模式的主要因素之一。根据单倍型在各群体的分布情况，进行青杨群体聚类分析结果，各群体无明显的分组现象，青杨各群体也未呈现任何清晰地理模式。 由于不同分子标记在对群体遗传多样性检测能力与效率上存在差异，所以三种标记检测的青藏高原东缘青杨天然群体遗传多性水平也不尽一致，但在与用同种方法检测其它物种或同一物种不同种源群体比较，三种分子标记方法都揭示了青藏高原东缘青杨天然群体具有中等偏上的遗传多样性水平。结果分析表明，群体间遗传分化极大，这是由于青杨天然群体分布于青藏高原东缘，既有高原又有高山峡谷，由于地理和空间上大尺度的隔离和地形地貌复杂导致了基因流物理上的阻隔。三种分子标记研究结果经Mantel分析检测，遗传距离与地理距离之间都无明显相关性。较为一致的解释是，青杨分布区域地理和空间上大尺度的隔离和和地形地貌复杂导致群体之间不存在均匀扩散现象，另外各群体在不同的选择压力下，经历各自独立的进化历程，这些都可能导致群体间遗传距离与地理距离的不相关。 The wide geographical and climatic distribution of P. cathayana Rehd. indicates that there is a large amount of genetic diversity available, which can be exploited for conservation, breeding programs and afforestation schemes. The results are as follows: 1 Research results of AFLP genetic diversity In present study, genetic diversity was evaluated in the natural populations of P. cathayana originating from southern and eastern edge of the Qinghai-Tibetan Plateau of China by means of AFLP markers. For four primer combinations, a total of 175 bands were obtained, of which 173 (98.9%) were polymorphic. Six natural populations of P. cathayana possessed different levels of genetic diversity, high level of genetic differentiation existed among populations (GST=0.489) of P. cathayana. Individuals cluster and PCO analysis based on Jaccard’s similarity coefficient also showed evident population genetic structure with high level population genetic differentiation. The long evolutionary process coupled with genetic drift within populations, rather than contemporary gene flow, are the major forces shaping genetic structure of P. cathayana populations. Moreover, there is no correspondence between geographical and genetic distances in the populations of P. cathayana, seldom gene exchange among populations and different selection pressures may be the causes. Our finding of different levels of genetic diversity within population and high level of genetic differentiation among populations provided promising condition for further breeding or conservation programs. 2 Research results of SSR genetic diversity In this study, the genetic diversity of P. cathayana was investigated using microsatellite markers. In a total of 150 individuals collected from six natural populations in the southeastern part of the Qinghai-Tibetan Plateau in China, a high level of microsatellite polymorphism was detected. At the seven investigated microsatellite loci, the number of alleles per locus ranged from 5 to 16, with a mean of 11.3, the observed heterozygosities across populations ranged from 0.408 to 0.986, with a mean of 0.792, and the expected heterozygosities across populations ranged from 0.511 to 0.891, with a mean of 0.802. The proportion of genetic differentiation among populations accounted for 37.3% of the whole genetic diversity. The presence of such a high level of genetic diversity could be attributed to the features of the species and the habitats where the sampled populations occur: The southeastern part of the Qinghai-Tibetan Plateau is regarded as the natural distribution and variation center of the genus Populus in China. Variation in environmental conditions and selection pressures in different populations, and topographic dispersal barriers could be factors associated with the high level of genetic differentiation found among populations. The populations possessed significant heterozygosity excesses, which may be due to extensive population mixing at the local scale. The cluster analysis showed that the populations are not strictly grouped according to their geographic distances but the habitat characteristics also influence the divergence pattern. In addition, we suggest that population SHY should be regarded as an ecologically divergent species of P. cathayana. 3 Research results of cpSSR genetic diversity Genetic diversity of six natural populations of P. cathayana originating from the southeastern part of the Qinghai-Tibetan Plateau in China was studied by use of cpSSR markers. Based on 5 pairs of polymorphic primers screened from 12 pairs of primers, twenty-six different length fragments and twelve different kinds of haplotypes were reduced in 143 samples. There were significant variant haplotypes among the populations．There were no shared haplotypes found among populations, analysis of molecular variance indicated that a high proportion of the total genetic variance was attributable to variations among populations (92.23%). The pattern of genetic structure which is associated with spatial separation, variation in environmental conditions and selection pressures in different populations, is also the result of geological historical factor. A molecular phylogenetic tree based on the 12 haplotypes showed that the populations are not strictly grouped according to their geographic distances.

太阳活动变化对中低纬电离层的影响

With the variations of solar activity, solar EUV and X-ray radiations change over different timescales (e.g., from solar cycle variation to solar flare burst). Since solar EUV and X-ray radiations are the primary energy sources for the ionosphere, theirs variations undoubtedly produce significant and complicated effects on the ionosphere. So the variations of solar activity significantly affect the ionosphere. It is essential for both ionospheric theory and applications to study solar activity effects on the ionosphere. The study about solar activity variations of the ionosphere is an important part of the ionospheric climatology. It can enhance the understanding for the basic processes in the ionosphere, ionospheric structure and its change, ionosphere/thermosphere coupling, and so on. As for applications, people need sufficient knowledges about solar activity variations of the ionosphere in order to improve ionospheric models so that more accurate forecast for the ionospheric environments can be made. Presently, the whole image about the modalities of ionospheric solar activity variations is still unknown, and related mechanisms still cannot be well understood. This paper is about the effects of the 11-year change in solar activity to the low- and mid-latitude ionosphere. We use multi-type ionospheric observations and model to investigate solar activity effects on the electron density and ionospheric spatial structure, and we focus on discussing some related mechanisms. The main works are as follows: Firstly, solar activity variations of ionospheric peak electron density (NmF2) around 1400 LT were investigated using ionosonde observations in the 120°E sector. The result shows that the variation trend of NmF2 with F107 depends on latitudes and seasons. There is obvious saturation trend in low latitudes in all seasons; while in middle latitudes, NmF2 increases linearly with F107 in winter but saturates with F107 at higher solar activity levels in the other seasons. We calculated the photochemical equilibrium electron density to discuss the effects induced by the changes of neutral atmosphere and dynamics processes on the solar activity variations of NmF2. We found that: (1) Seasonal variation of neutral atmosphere plays an important role in the seasonal difference of the solar activity variations of NmF2 in middle latitudes. (2) Less [O]/[N2] and higher neutral temperature are important for the saturation effect in summer, and the increase of vibrational excited N2 is also important for the saturation effect. (3) Dynamics processes can significantly weaken the increase of NmF2 when solar activity enhances, which is also a necessary factor for the saturation effect. Secondly, solar activity variations of nighttime NmF2 were investigated using ionosonde observations in the 120°E sector. The result shows that the variation trends of NmF2 with F107 in nighttime are different from that in daytime in some cases, and the nighttime variation trends depend on seasons. There is linear increase trend in equinox nighttime, and saturation trend in summer nighttime, while the increase rate of NmF2 with F107 increases when solar activity enhances in winter nighttime (we term it with “amplification trend”). We discussed the possible mechanisms which affect the solar activity variations of nighttime NmF2. The primary conclusions are as follows: (1) In the equatorial ionization anomaly (EIA) crest region, the plasma influx induced by the pre-reversal enhancement (PRE) results in the change of the variation trend between NmF2 and F107 from “saturation” to “linear” after sunset in equinoxes and winter; while the recombination process at the F2-peak is the primary factor that affects the variation trend of NmF2 with F107 in middle latitudes. (2) The recombination coefficient at the F2-peak height reaches its maximum at moderate solar activity level in winter nighttime, which induces NmF2 attenuates more quickly at moderate solar activity level. This is the main reason for the amplification trend. (3) The change of the recombination process at the F2-peak with solar activity depends on the increases of neutral parameters (temperature, density et al.) and the F2-peak height (hmF2). The seasonal differences in the changes of neutral atmosphere and hmF2 with solar activity are the primary reasons for the seasonal difference in the variation trend of nighttime NmF2 with F107. Finally, we investigated the solar activity dependence of the topside ionosphere in low latitudes using ROCSAT-1 satellite (at 600 km altitude) observations. The primary results and conclusions are as follows: (1) Latitudinal distribution of the plasma density is local time, seasonal, and solar activity dependent. In daytime, there is a plasma density peak at the dip equator. The peak is obviously enhanced at high solar activity level, and the strength of the peak strongly depends on seasons. While at sunset, two profound plasma density peaks (double-peak structure) are found in solar maximum equinox months. (2) Local time dependence of the latitudinal distribution is due to the local time variation of the equatorial dynamics processes. Double-peak structure is attributed to the fountain effect induced by strong PRE. Daytime peak enhances with solar activity since the plasma density increases with solar activity more strongly at the dip equator due to the equatorial vertical drift, and its seasonal dependence is mainly due to the seasonal variations of neutral density and the equatorial vertical drift. In the sunset sector, seasonal and solar activity dependences of the latitudinal distribution are related to the seasonal and solar activity variations of PRE. (3) The variation trend of the plasma density with solar activity shows local time, seasonal, and latitudinal differences. That is different from the changeless amplification trend at the DMSP altitude (840 km). Profound saturation effect is found in the dip equator region at equinox sunset. This saturation effect in the topside ionosphere is realated to the increase of PRE with solar activity. Solar activity variation trend of the topside plasma density was discussed quantitatively by Chapman-α function. The result shows that the effect induced by the change of the scale height is dominant at high altitudes; while the variation trend of ROCSAT-1 plasma density with solar activity is suggested to be related to the changes of the peak height, the scale height, and the peak electron density with solar activity.

A generalized wave action conservative equation for the dissipative dynamical system in the nearshore region

To describe the various complex mechanisms of the dissipative dynamical system between waves, currents, and bottoms in the nearshore region that induce typically the wave motion on large-scale variation of ambient currents, a generalized wave action equation for the dissipative dynamical system in the nearshore region is developed by using the mean-flow equations based on the Navier-Stokes equations of viscous fluid, thus raising two new concepts: the vertical velocity wave action and the dissipative wave action, extending the classical concept, wave action, from the ideal averaged flow conservative system into the real averaged flow dissipative system (that is, the generalized conservative system). It will have more applications.

Non-relativistic conformal symmetries in fluid mechanics

The symmetries of a free incompressible fluid span the Galilei group, augmented with independent dilations of space and time. When the fluid is compressible, the symmetry is enlarged to the expanded Schrodinger group, which also involves, in addition, Schrodinger expansions. While incompressible fluid dynamics can be derived as an appropriate non-relativistic limit of a conformally invariant relativistic theory, the recently discussed conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry. This is explained by the subtleties of the non-relativistic limit.

COASTAL SEA-LEVEL AND THE LARGE-SCALE CLIMATE STATE - A DOWNSCALING EXERCISE FOR THE JAPANESE ISLANDS

A major problem which is envisaged in the course of man-made climate change is sea-level rise. The global aspect of the thermal expansion of the sea water likely is reasonably well simulated by present day climate models; the variation of sea level, due to variations of the regional atmospheric forcing and of the large-scale oceanic circulation, is not adequately simulated by a global climate model because of insufficient spatial resolution. A method to infer the coastal aspects of sea level change is to use a statistical ''downscaling'' strategy: a linear statistical model is built upon a multi-year data set of local sea level data and of large-scale oceanic and/or atmospheric data such as sea-surface temperature or sea-level air-pressure. We apply this idea to sea level along the Japanese coast. The sea level is related to regional and North Pacific sea-surface temperature and sea-level air pressure. Two relevant processes are identified. One process is the local wind set-up of water due to regional low-frequency wind anomalies; the other is a planetary scale atmosphere-ocean interaction which takes place in the eastern North Pacific.

Small-scale species richness and its spatial variation in an alpine meadow on the Qinghai-Tibet Plateau

We investigated how the high small-scale species richness of an alpine meadow on the Qinghai-Tibet Plateau, China, is maintained. This area is characterized by strong wind and severe cold during long winters. In winter, most livestock is grazed on dead leaves in small pastures near farmers' residences, whereas in the short summer, livestock is grazed in mountainous areas far from farmers' residences. The number of plant species and the aboveground biomass were surveyed for three adjacent pastures differing in grazing management: a late-winter grazing pasture grazed moderately from 1 February to 30 April, an early-winter grazing pasture grazed lightly from 20 September to late October, and a whole-year grazing pasture grazed intensively throughout the entire year. In each pasture, we harvested the aboveground biomass from 80 or 100 quadrats of 0.01 m(2) along a transect and classified the contents by species. We observed 15.5-19.7 species per 0.01 m(2), which is high richness per 0.01 m(2) on a worldwide scale. The species richness in the two winter grazing pastures was higher than that in the whole-year grazing pasture. The spatial variation in species richness and species composition in the two winter grazing pastures in which species richness was high was greater than that in the whole-year grazing pasture in which species richness was lower. Most of the leaves that are preserved on the winter grazing pastures during summer are blown away by strong winds during winter, and the remaining leaves are completely exhausted in winter by livestock grazing. A pasture with a high richess is accompanied by a high spatial variation in species richness and species composition. There is a high possibility that the characteristic of spatial variation is also caused by traditional grazing practices in this area.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Influences of strain rate and temperature variation on material intrinsic length of plasticity strain gradient theory

Cowper-Symonds and Johnson-Cook dynamic constitutive relations are used to study the influence of both strain rate effect and temperature variation on the material intrinsic length scale in strain gradient plasticity. The material intrinsic length scale decreases with increasing strain rates, and this length scale increases with temperature.

Microdamage Evolution, Energy Dissipation and Their Trans-Scale Effects on Macroscopic Failure

The process of damage evolution concerns various scales, from micro- to macroscopic. How to characterize the trans-scale nature of the process is on the challenging frontiers of solid mechanics. In this paper, a closed trans-scale formulation of damage evolution based on statistical microdamage mechanics is presented. As a case study, the damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that the following dimensionless numbers: reduced Mach number M, damage number S, stress wave Fourier number P, intrinsic Deborah number D*, and the imposed Deborah number De*, govern the whole process of deformation and damage evolution. The evaluation of P and the estimation of temperature increase show that the energy equation can be ignored as the first approximation in the case of spallation. Hence, apart from the two conventional macroscopic parameters: the reduced Mach number M and damage number S, the damage evolution in spallation is mainly governed by two microdamage-relevant parameters: the Deborah numbers D* and De*. Higher nucleation and growth rates of microdamage accelerate damage evolution, and result in higher damage in the target plate. In addition, the mere variation in nucleation rate does not change the spatial distribution of damage or form localized rupture, while the increase of microdamage growth rate localizes the damage distribution in the target plate, which can be characterized by the imposed Deborah number De*.

Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.