Plio-Pleistocene climatic change drives allopatric speciation and population divergence within the Scrophularia incisa complex (Scrophulariaceae) of desert and steppe subshrubs in Northwest China.

Numerous temperate plants and animals on the Qinghai-Tibet Plateau (QTP) are hypothesized to have differentiated due to vicariant allopatric speciation associated with the geologic uplifts. However, this hypothesis has rarely been tested through a phylogeographic study of relative species in a broader geographic context, including the QTP, Tianshan Mountains, Mongolian Plateau, and surrounding regions. To understand the speciation and diversification process of plants across this wide area, phylogeographic analysis were examined from Scrophularia incisa and two other closely relative species comprising S. kiriloviana and S. dentata. Thirty-two populations of the three close relatives were genotyped using chloroplast DNA fragments and nuclear microsatellite loci to assess population structure and diversity, supplemented by phylogenetic dating, ancestral area reconstructions and species distribution modelings, as well as niche identity tests. Our chloroplast DNA (cpDNA) phylogeny showed that this monophyletic group of desert and steppe semi-shrub is derived from a Middle Pliocene ancestor of the Central Asia. Lineages in Central Asia vs. China diverged through climate/tectonic-induced vicariance during Middle Pliocene. Genetic and ENM data in conjunction with niche differentiation analyses support that the divergence of S. incisa, S. dentata and S. kiriloviana in China lineage proceeded through allopatric speciation, might triggered by early Pleistocene climate change of increase of aridification and enlargement of deserts, while subsequent climate-induced cycles of range contractions/expansions enhanced the geographical isolation and habit fragmentation of these taxa. These findings highlight the importance of the Plio-Pleistocene climate change in shaping genetic diversity and driving speciation in temperate steppes and deserts of Northwestern China.

Development and Application of Genomic Resources in an Endangered Palaeoendemic Tree, Parrotia subaequalis (Hamamelidaceae) From Eastern China.

Parrotia subaequalis is an endangered palaeoendemic tree from disjunct montane sites in eastern China. Due to the lack of effective genomic resources, the genetic diversity and population structure of this endangered species are not clearly understood. In this study, we conducted paired-end shotgun sequencing (2 × 125 bp) of genomic DNA for two individuals of P. subaequalis on the Illumina HiSeq platform. Based on the resulting sequences, we have successfully assembled the complete chloroplast genome of P. subaequalis, as well as identified the polymorphic chloroplast microsatellites (cpSSRs), nuclear microsatellites (nSSRs) and mutational hotspots of chloroplast. Ten polymorphic cpSSR loci and 12 polymorphic nSSR loci were used to genotype 96 individuals of P. subaequalis from six populations to estimate genetic diversity and population structure. Our results revealed that P. subaequalis exhibited abundant genetic diversity (e.g., cpSSRs: Hcp = 0.862; nSSRs: HT = 0.559) and high genetic differentiation (e.g., cpSSRs: RST = 0.652; nSSRs: RST = 0.331), and characterized by a low pollen-to-seed migration ratio (r ≈ 1.78). These genetic patterns are attributable to its long evolutionary histories and low levels of contemporary inter-population gene flow by pollen and seed. In addition, lack of isolation-by-distance pattern and strong population genetic structuring in both marker systems, suggests that long-term isolation and/or habitat fragmentation as well as genetic drift may have also contributed to the geographic differentiation of P. subaequalis. Therefore, long-term habitat protection is the most important methods to prevent further loss of genetic variation and a decrease in effective population size. Furthermore, both cpSSRs and nSSRs revealed that P. subaequalis populations consisted of three genetic clusters, which should be considered as separated conservation units.

Use of portable X-ray fluorescence spectroscopy and geostatistics for health risk assessment.

Laboratory analysis of trace metals using inductively coupled plasma (ICP) spectroscopy is not cost effective, and the complex spatial distribution of soil trace metals makes their spatial analysis and prediction problematic. Thus, for the health risk assessment of exposure to trace metals in soils, portable X-ray fluorescence (PXRF) spectroscopy was used to replace ICP spectroscopy for metal analysis, and robust geostatistical methods were used to identify spatial outliers in trace metal concentrations and to map trace metal distributions. A case study was carried out around an industrial area in Nanjing, China. The results showed that PXRF spectroscopy provided results for trace metal (Cu, Ni, Pb and Zn) levels comparable to ICP spectroscopy. The results of the health risk assessment showed that Ni posed a higher non-carcinogenic risk than Cu, Pb and Zn, indicating a higher priority of concern than the other elements. Sampling locations associated with adverse health effects were identified as 'hotspots', and high-risk areas were delineated from risk maps. These 'hotspots' and high-risk areas were in close proximity to and downwind from petrochemical plants, indicating the dominant role of industrial activities as the major sources of trace metals in soils. The approach used in this study could be adopted as a cost-effective methodology for screening 'hotspots' and priority areas of concern for cost-efficient health risk management.

[Magnetic Response of Heavy Metal Pollution in Playground Dust of an Industrial Area].

The use of mineral magnetic measurements as a potential pollution proxy is explored as an alternative means of monitoring heavy metal pollution in playground dust in industrial areas. Dust samples were collected from playgrounds in an industrial area in Nanjing. The magnetic properties of the samples were analyzed, and the total and acid-extractable concentrations of heavy metals were measured.Mineral magnetic concentration parameters[mean magnetic susceptibility (χlf) 939.31×10-8 m3·kg-1 and saturation isothermal remanent magnetization (SIRM) 16618.74×10-5 A·m2·kg-1] indicate that the dust samples contain high concentrations of magnetic minerals when compared to the non-industrial zone of the city. High and stable S-ratio values (mean 0.97, standard deviation 0.01) indicate the prominence of ferromagnetic minerals. Low values of anhysteretic remanent magnetization (χARM)/χlf(2.35) and χARM/SIRM (0.12×10-3 m·A-1) indicate predominantly coarse pseudo single domain (PSD) and multi-domain (MD) magnetic grain size assemblages. The playground dust is considered moderately to highly polluted, and the enrichment factor values for Cd, Zn, Cu, and Pb indicate extremely high contamination levels in the playground dust. Mineral magnetic parameters[χlf, SIRM, hard isothermal remanent magnetization (HIRM), χARM, and χARM/SIRM] reveal significant correlations with total and acid-extractable concentrations of the trace metals, and the correlations with total concentrations are stronger. The most significant correlations were found with mineral magnetic concentrations (χARM and SIRM) and total concentrations of the elements Mn, Ni, and Cr and the Tomlinson pollution load index (PLI) (0.69 ≤ r ≤ 0.86, P<0.01). Consistent spatial characteristics between them were also found, showing that the total concentrations of heavy metals and magnetic parameter values are higher in the playgrounds near and downwind of the power plant. Results of a principal component analysis indicate that anthropogenic activities (mainly from coal consumption by power plants) are the main sources for both heavy metals and magnetic minerals. In summary, significant correlations and consistent spatial characteristics between heavy metal contents and magnetic parameters and the same anthropogenic source for heavy metals and magnetic minerals confirm that magnetic parameters could be considered as efficient indicators for heavy metal contamination in the playground dust in industrial areas.

[Structural characteristics of Sophora moorcroftiana community on wind-sandy land in middle reaches of Yaluzangbu River].

Based on quadrat investigation and hierarchical cluster analysis, this paper studied the structural characteristics of Sophora moorcroftiana community at its different restoration succession stages on the wind-sandy land in middle reaches of Yaluzangbu River. The plant community on the wind-sandy land could be divided into 6 types, i.e., S. moorcroftiana community, Orinus thoroldii community, S. moorcroftiana < O. thoroldii community, S. moorcroftiana-O. thoroldii community, Artemisia weiibyii+S. moorcroftiana community, and S. moorcroftiana+A. weiibyii community. The restoration succession began with the sand-fixing S. moorcroftiana distributed on moving dune, followed by O. thoroldii invasion, and then, S. moorcroftiana-O. thoroldii community formed gradually. Another restoration succession approach began with O. thoroldii community, S. moorcroftiana occurred then, and finally, S. moorcroftiana-O, thoroldii community formed. On semi-fixed dune, A. weiibyi could invade the community, and S. moorcroftiana+A. weiibyi community formed. The climax community on the wind-sandy land was the shrub community dominated by S. moorcroftiana. With increasing slope grade, the coverage and height of S. moorcroftiana community tended to decrease. Biological crust had definite improvement effect on the soil physical and chemical properties. At the places with biological crust, surface soil water content was significantly higher, compared with that at the places with other restoration succession stages of S. moorcroftiana community. The crown diameter, plant height, and the sprouting branch number of S. moorcroftiana community were the highest on moving dune, followed by on fixed dune, and on semi-fixed dune. The dead branch number at different restoration succession stages indicated that S. moorcroftiana population had a declining trend after the sand was fixed.

[Dynamics and influence width of dry valley landscape boundary in upper reaches of Minjiang River].

Based on RS and GIS techniques, the distribution extent of the dry valley in upper reaches of Minjiang River was determined, with the dynamics and influence width of its landscape boundary analyzed. The results indicated that the dry valley had a gradually expanded distribution, and the total length of its boundary was increased. The boundary complexity increased from 1974 to 1995 but decreased from 1995 to 2000, and the upper boundary of the dry valley ascended rapidly along mountain slope. The uppermost elevation of the dry valley boundary was 3128 m in 1974, 3167 m in 1995, and 3181 m in 2000. In the last 26 years, the upper boundary of the dry valley ascended 53 m in total, being about 2 m per year. The influence width of the dry valley boundary was 800 m, which could be regarded as a buffer zone to disturbances. The area between 800-1200 m from the dry valley boundary to exterior was a transition zone intensely influenced by the dry valley and its surrounding landscape, which should be addressed in the restoration of the dry valley.

[Effects of sub-watershed landscape patterns at the upper reaches of Minjiang River on soil erosion].

Based on GIS, the spatial distribution of soil loss and sediment yield in Heishui and Zhenjiangguan sub-watersheds at the upper reaches of Minjiang River was simulated by using sediment delivery-distribution (SEDD) model, and the effects of land use/cover types on soil erosion and sediment yield were discussed, based on the simulated results and related land use maps. A landscape index named location-weighted landscape contrast index (LCI) was calculated to evaluate the effects of landscape components' spatial distribution, weight, and structure of land use/cover on soil erosion. The results showed the soil erosion modulus varied with land use pattern, and decreased in the order of bare rock > urban/village > rangeland > farmland > shrub > forest. There were no significant differences in sediment yield modules among different land use/covers. In the two sub-watersheds, the spatial distribution of land use/covers on slope tended to decrease the final sediment load at watershed outlet, hut as related to relative elevation, relative distance, and flow length, the spatial distribution tended to increase sediment yield. The two sub-watersheds had different advantages as related to landscape components' spatial distribution, but, when the land use/cover weight was considered, the advantages of Zhenjiangguan sub-watershed increased. If the land use/cover structure was considered in addition, the landscape pattern of Zhenjiangguan subwatershed was better. Therefore, only the three elements, i.e., landscape components' spatial distribution, land use/cover weight, and land use/cover structure, were considered comprehensively, can we get an overall evaluation on the effects of landscape pattern on soil erosion. The calculation of LCI related to slope suggested that this index couldn' t accurately reflect the effects of land use/cover weight and structure on soil erosion, and thus, needed to be modified.