读刊有感

《云南林业》一奇葩， 知识园地职工家。 内容丰富图文茂， 健笔常开鲜艳花。

Using multispectral landsat and sentinel-2 satellite data to investigate vegetation change at Mount St. Helens since the great volcanic eruption in 1980

Long-term analyses of vegetation succession after catastrophic events are of high interest for an improved understanding of succession dynamics. However, in many studies such analyses were restricted to plot-based measurements. Contrarily, spatially continuous observations of succession dynamics over extended areas and timeperiods are sparse. Here, we applied a change vector analysis(CVA) to investigate vegetation succession dynamics at Mount St. Helens after the great volcanic eruption in 1980 using Landsat. We additionally applied a supervised random forest classification using Sentinel-2 data to map the currently prevailing vegetation types. Change vector analysis was performed with the normalized difference vegetation index(NDVI) and the urban index(UI) for three subsequent decades after the eruption as well as for the whole observation time between 1984 and 2016. The influence of topography on the current vegetation distribution was examined by comparing altitude, slope angles and aspect values of vegetation classes derived by the random forest classification. WilcoxRank-Sum test was applied to test for significant differences between topographic properties of the vegetation classes inside and outside of the areas affected by the eruption. For the full time period, a total area of 516 km2 was identified as re-vegetated, whereas the area and magnitude of re-growing vegetation decreased during the three decades and migrated closer to the volcanic crater. Vegetation losses were mainly observed in regions unaffected by the eruption and related mostly to timber harvesting. The vegetation type classification reached a high overall accuracy of approximately 90%. 36 years after the eruption, coniferous and deciduous trees have established at formerly devastated areas dominating with a proportion of 66%, whereas shrubs are more abundant in riparian zones. Sparse vegetation dominates at regions very close to the crater. Elevation was found to have a great influence on the reestablishment and distribution of the vegetation classes within the devastated areas showing in almost all cases significant differences in altitude distribution. Slope was less important for the different classes-only representing significantly higher values for meadows, whereas aspect seems to have no notable influence on the reestablishment of vegetation at Mount St. Helens. We conclude that major vegetation succession dynamics after catastrophic events can be assessed and characterized over large areas from freely available remote sensing data and hence contribute to an improved understanding of succession dynamics.

Vertical Distribution and Pollution Assessment of Heavy Metals in Mangrove Sediments of Dongzhai Harbor, Hainan

This study was conducted to reveal the vertical distribution and concentration variation of heavy metals（Cr, Co, Ni, Cu, Zn, As, Cd and Pb） in the mangrove sediments of Dongzhai Harbor, Hainan by using inductively-coupled plasma mass spectrometry（ICP-MS）. The surface sediment samples were collected from 15 sites, and a geoaccumulation index（Igeo） was estimated to quantify the pollution status of the mangrove sediments. The results showed that the concentrations of heavy metals in sediments followed the order of Cr〉 Zn〉 Ni〉 Cu 〉Pb 〉Co〉 Cd. A common pattern was found in some heavy metals（Cu, Zn, As and Pb） as their concentrations were significantly and negatively correlated with depth. Cr concentration in upper layer was slightly lower than that in middle and lower layers,while the concentrations of Ni and Cd in middle layer were slightly higher than in upper and lower layers. The geoaccumulation index indicated that the surface sediments at sites QT-1-01, QT-1-02, QT-2, GN-1, GN-2 and GN-3 were strongly polluted by Cd, and unpolluted to moderately polluted or moderately polluted by As.

Effects of slope aspects on forest compositions, community structures and soil properties in natural temperate forests of Garhwal Himalaya

The present study was undertaken in seven natural forest types of temperate zone （1500 to 3100 m, a.s.l.） of Garhwal Himalaya to understand the effects of slope aspects viz., north-east （NE）, north-west （NW）, south-east （SE） and south-west （SW）, on the forest structure, composition and soil characteristics of each selected forest type. The sample plots in each forest type were laid out by using stratified random approach. The indices i.e., the Importance Value Index, Shannon-Wiener diversity index, Simpson’s concentration of Dominance, Simpson diver-sity index, Pielou equitability and Margalef species richness index were calculated statistically using standard softwares to elucidate the differ-ences in forest structure and composition of forest types on different slope aspects of the sites. The composite soil samples were taken from each forest stand and the physico-chemical properties of the soil i.e., moisture content （MC）, water holding capacity （WHC）, pH, organic carbon （OC）, phosphorus （P）, potassium （K） and available nitrogen （N） were analyzed. The results show that the higher values of total basal cover （74.4 m2·ha-1 in Quercus semecarpifolia forest）, Concentration of dominance （0.85 in Pinus roxburghii forest） and Tree diversity （1.81 in Quercus floribunda forest） in the forests were recorded in the northern aspects. MC （40.8% in Quercus leucotrichophora forest）, WHC （48.9% in Cupressus torulosa forest）, OC （3.8% in Cedrus deodara forest）, P （31.9 kg·ha-1 in Quercus leucotrichophora forest） and N （1.0% in Pinus roxburghii forest） had also higher values in the soils of northern aspects. Consequently the higher productivity of the forests was also noticed on the northern aspects.

Changes in spatial variations of sap flow in Korean pine trees due to environmental factors and their effects on estimates of stand transpiration

It is difficult to scale up measurements of the sap flux density(J_S) for the characterization of tree or stand transpiration(E) due to spatial variations in J_S and their temporal changes.To assess spatial variations in the sap flux density of Korean pine(Pinus koraiensis) and their effects on E estimates,we measured the J_S using Granier-type sensors.Within trees,the J_S decreased exponentially with the radial depth,and the J_S of the east aspects were higher than those of the west aspects.Among trees,there was a positive relationship between J_S and the tree diameter at breast height,and this positive relationship became stronger as the transpiration demand increased.The spatial variations that caused large errors in E estimates(i.e.,up to 110.8 % when radial variation was ignored) had varied systematically with environmental factors systematic characteristics in relation to environmental factors.However,changes in these variations did not generate substantial errors in the E estimates.For our study periods,the differences in the daily E(E_D) calculated by ignoring radial,azimuthal and tree-to-tree variations and the measured E_D were fairly constant,especially when the daily vapor pressure deficit(D_D)was higher than 0.6 k Pa.These results imply that the effect of spatial variations changes on sap flow can be a minor source of error compared with spatial variations(radial,azimuthal and tree-to-tree variations) when considering E estimates.

Biomass production of herb species in broad leaf forests in Kumaun Himalaya,India

The study focuses on the dynamics and biomass production （g·m-2） of understory （herbaceous） plant community under broad-leaf forests consisting of Quercus leucotrichophora （Banj-oak）, Quercus floribunda （Tilonj-oak） and Quercus semicarpifolia （Kharsu-oak）, respectively in central Himalaya, India. With increasing altitudes, the den-sity and biomass decreased significantly across the three types of forests. Banj-oak forest harboured the maximum density and biomass among the other sampling sites. The mean density of herb species in two contrasting orientation differed significantly （p0.05）, showing relatively higher density on slope orientation （west）. Across the sites, total production declined significantly with increasing altitude of the sites （p0.05）, and Banj-oak forest presented the highest production.

Estimation of wear performance of AZ91 alloy under dry sliding conditions using machine learning methods

The wear behavior of AZ91 alloy was investigated by considering different parameters,such as load(10−50 N),sliding speed(160−220 mm/s)and sliding distance(250−1000 m).It was found that wear volume loss increased as load increased for all sliding distances and some sliding speeds.For sliding speed of 220 mm/s and sliding distance of 1000 m,the wear volume losses under loads of 10,20,30,40 and 50 N were calculated to be 15.0,19.0,24.3,33.9 and 37.4 mm3,respectively.Worn surfaces show that abrasion and oxidation were present at a load of 10 N,which changes into delamination at a load of 50 N.ANOVA results show that the contributions of load,sliding distance and sliding speed were 12.99%,83.04%and 3.97%,respectively.The artificial neural networks(ANN),support vector regressor(SVR)and random forest(RF)methods were applied for the prediction of wear volume loss of AZ91 alloy.The correlation coefficient(R2)values of SVR,RF and ANN for the test were 0.9245,0.9800 and 0.9845,respectively.Thus,the ANN model has promising results for the prediction of wear performance of AZ91 alloy.

Geochronological Significance of the Post-Orogenic Mafic-Ultramafic Rocks in the Hongqiling Area of Jilin Province, Northeast China

Mafic-uitramafic complexes are widespread in Hongqiling in central Jiulin Province, NE China. The Hongqiling complex comprises pyroxenite, olivine websterite, lberzolite, gabbro and leucogabbro. Based on the latest geochronological results, these intrusions yield a zircon U-Pb age of about 216 Ma, implying that they emplaced in the late Triassic period and that the mafic-uitramafic complexes are post-orogenic in origin. These ages are coeval with the emplacement of A-type granites in the area, but slightly younger than syn-orogenic granitic magmatism. The composition of the parent magma during the equilibration of the accumulated olivine crystallizing is basaltic with high MgO, and it comes from depleted or juvenile lithospheric mantle, according to the results of trace elements and Sr-Nd isotopes. Amount of crustal material contaminated the magmatic source, implying that the mafic-uitramafic rocks originate from the mixing product of crust- and mantle-derived magma. Therefore, the magmatic source was contaminated by an amount of crustal material, and the subsequent crystal fractionation resulted in the Cu-Ni mineralization. Dynamic analyses indicate that, after the collision of the North China Craton and Jiamusi Massif, the lithospheric delamination during post-orogeny resulted in the upwelling of asthenosphere and the intruding of the mass and underplating of the mantle-derived magma, which led to the formation of the primary mafic-uitramafic magma due to the mantle-derived magma and partial melting of the lower crust. This result suggests that the mafic-uitramafic complexes belong to a typically important magmatic suite that evolved during post-orogenic processes, and they were exposed as maficultramafic dyke swarms that existed in the post-orogenic extension environment. It also implies that the orogenic process finished ultimately in this regional geological setting. The widespread occurrence of mafic-uitramafic complexes in the Xing＇an-Mongolia Orogenic Belt in NE China and in the Altay-Tianshan-Junggar Orogenic Belt in Northern Xinjiang indicates that mafic intrusions are an important magmatic suite that evolved during post-orogenic processes. Portions of this mafic magma could have underplated the lower crust, and served as a heat source for associated late-stage granitic magmas.

Formulation of Green's Function in Solving Schroedinger Equation forBipartite System with Kinetic Coupling Gained by Virtue of Entangled States

Using the entangled state representation we present a formulation of Green'sfunction in solving Schrodinger equation for bipartite system with kinetic coupling.

Effects of windbreak width in wind direction on wind velocity reduction

The variations of drag force acting on the windbreak and the bulk drag coefficients for different windbreak widths were studied experimentally in the Eiffel-type non-circulating wind tunnel at the Hydraulic Engineering Laboratory, Saitama University, Japan, to elucidate the effects of windbreak width in the wind direction on wind velocity reduction behind a windbreak. The variations of flow field for different windbreak widths were studied numerically by using the two-dimensional Reynolds-averaged Navier-Stokes （RANS） equation with a k-c turbulence closure model. Results show that the total drag force to wind increased with increasing windbreak width, but the bulk drag coefficient decreased slightly. The relationship between the bulk drag coefficient Cd and the windbreak width W and height H can be presented by the equation of Ca= kd （W/H）^-b （kd, b： constants）. The result of the numerical simulation shows that the windbreak width greatly affects the location and the value of the minimum wind velocity. The wind velocity decreased by 15%-22% as the windbreak width increased.

Prostate Cancer Risk Prediction and Online Calculation Based on Machine Learning Algorithm

Objective To build a prostate cancer(PCa) risk prediction model based on common clinical indicators to provide a theoretical basis for the diagnosis and treatment of PCa and to evaluate the value of artificial intelligence(AI) technology under healthcare data platforms.Methods After preprocessing of the data from Population Health Data Archive,smuothly clipped absolute deviation(SCAD) was used to select features.Random forest(RF),support vector machine(SVM),back propagation neural network(BP),and convolutional neural network(CNN) were used to predict the risk of PCa,among which BP and CNN were used on the enhanced data by SMOTE.The performances of models were compared using area under the curve(AUC) of the receiving operating characteristic curve.After the optimal model was selected,we used the Shiny to develop an online calculator for PCa risk prediction based on predictive indicators.Results Inorganic phosphorus,triglycerides,and calcium were closely related to PCa in addition to the volume of fragmented tissue and free prostate-specific antigen(PSA).Among the four models,RF had the best performance in predicting PCa(accuracy:96.80%;AUC:0.975,95% CI:0.964-0.986).Followed by BP(accuracy:85.36%;AUC:0.892,95% CI:0.849-0.934) and SVM(accuracy:82.67%;AUC:0.824,95% CI:0.805-0.844).CNN performed worse(accuracy:72.37%;AUC:0.724,95% CI:0.670-0.779).An online platform for PCa risk prediction was developed based on the RF model and the predictive indicators.Conclusions This study revealed the application value of traditional machine learning and deep learning models in disease risk prediction under healthcare data platform,proposed new ideas for PCa risk prediction in patients suspected for PCa and had undergone core needle biopsy.Besides,the online calculation may enhance the practicability of AI prediction technology and facilitate medical diagnosis.

Asymmetric Green's functions for exponentially graded transversely isotropic substrate–coating system

By virtue of a complete set of two displacement potentials,an analytical derivation of the elastostatic Green’s functions of an exponentially graded transversely isotropic substrate–coating system is presented.Three-dimensional point–load and patch–load Green’s functions for stresses and displacements are given in line-integral representations.The formulation includes a complete set of transformed stress–potential and displacement–potential relations,with utilizing Fourier series and Hankel transforms.As illustrations,the present Green’s functions are degenerated to the special cases such as an exponentially graded half-space and a homogeneous two-layered half-space Green’s functions.Because of complicated integrand functions,the integrals are evaluated numerically and for numerical computation of the integrals,a robust and effective methodology is laid out which gives the necessary account of the presence of singularities of integration.Comparisons of the existing numerical solutions for homogeneous two-layered isotropic and transversely isotropic half-spaces are made to confirm the accuracy of the present solutions.Some typical numerical examples are also given to show the general features of the exponentially graded two-layered half-space Green’s functions that the effect of degree of variation of material properties will be recognized.

Magnetic Properties of Spin-1／2 Quantum Heisenberg Ferromagnet with Spatially Anisotropy

The spin-1/2 quantum Heisenberg ferromagnet with anisotropic spin exchange interactions in three dimensions are investigated by means of the Green's function method. In the Tyablikov approximation, the correlation function s, the magnetization, and the susceptibilities are computed. The magnetic properties of this model are found to be dependent of anisotropy.

Multi-domain Boundary Element Method with Dissipation

The wave diffraction and radiation around a floating body is considered within the framework of the linear potential theory in a fairly perfect fluid. The fluid domain extended infinitely in the horizontal directions but is limited by the sea bed, the body hull, and the part of the free surface excluding the body waterplane, and is subdivided into two subdomains according to the body geometry. The two subdomains are connected by a control surface in fluid. In each subdomain, the velocity potential is described by using the usual boundary integral representation involving Green functions. The boundary integral equations are then established by satisfying the boundary conditions and the continuous condition of the potential and the normal derivation across the control surface. This multi-domain boundary element method （MDBEM） is particularly interesting for bodies with a hull form including moonpools to which the usual BEM presents singularities and slow convergence of numerical results. The application of the MDBEM to study the resonant motion of a water column in moonpools shows that the MDBEM provides an efficient and reliable prediction method.

Omni-Directional Distribution Patterns of Montane Coniferous Forest in the Helan Mountains of China

Slope aspect has significant effect on the development and distribution of montane forest, especially in arid and semiarid regions. This paper, using SPOT5 images and 1：5o,ooo DEM, digitally extracts and analyzes the spatial information of montane coniferous forest （mainly Qinghai spruce）, and thereby explores how the upper and lower limits and the altitudinal range of coniferous forest vary and how the area of coniferous forest is related with annual insolation with all aspects in the Helan Mountains. The results show that： 1） In the eastern flank, the lower limit of coniferous forest is between 1,6oo m and a,ooo m a.s.l., and the upper limit between 2,800 m and 3,ooo m a.s.1.; in the western flank, the lower limit of coniferous forest is between a,ooo m and 2,300 m a.s.l., and the upper limit between 9,80o m and 3,1oo m a.s.1.2） The altitudinal ranges of coniferous forest are 8o6-1,435 m, 438-1,14o m for eastern flank and western flank, respectively. 3） The area of coniferous forest takes on a normal distribution with aspect, and it has a close relationship with annual insolation. This distribution model developed in this paper quantitatively reveals the significant effect of slope aspect on the distribution of coniferous forest in arid and semi-arid land.

Delineation of Haigou gold metallogenic province and its metallogenic prospection

Based on the accumulated data for the gold deposits in the central Jilin Province in recent years and our understanding of the gold metallogenic province,the Haigou gold metallogenic province is delineated and the denudation degree of gold deposits in this province is discussed. The potential and the ore-searching direction of the province are also considered. The Haigou gold metallogenic province occurs as an independent province with low denudation degree and high ore-producing potential. Regional fault belts and small basic intrusions are two ore-constrains and could serve as the ore-searching indictors in the province.

Study of Phase Transition Properties in Epitaxial Ferroelectric Film

Based on the transverse Ising model and using decoupling approximation to the Fermi-type Green＇s function, we study the phase transition properties of the epitaxial ferroeleetric film with one substrate. A general recursive equation of the ferroelectric thin film with two n-layer materials is obtained, which enables us to study the phase transition properties for any number layers for epitaxial ferroelectric thin film. With the help of this equation, we analyze the effect of the exchange interaction and the transverse field in the phase diagram, the influence to the polarizations and Curie temperature numerically. The results show that epitaxial ferroelectric film are able to induce a strong increase or decrease of Curie temperature to different exchange interactions and transverse fields within the epitaxial film layers. The theoretical results are in reasonable accordance with experimental data of different ferroelectric thin film.

Spatial Differentiation Characteristics and Driving Forces of Forest Transition:A Case Study of Zunyi City,Guizhou

Spatial differentiation in forest transition was measured in terms of space transition and function transition using the exploratory spatial data analysis method（ESDA） and data from 2004—2014 for Zunyi city,Guizhou province,China.The validity of factors affecting forest transition was analyzed by constructing radial basis function neural networks（RBFNN） based on the data processing system（DPS）.Our results will provide references for scientific understanding of the potential mechanism underlying forest transition in mountainous areas.We found that Global Moran＇s I of space transition and function transition of forest land was 0.0336 and 0.2323,respectively.This suggests a significant positive correlation in spatial distribution of space transition and function transition of forest land,and significant spatial aggregation.The Global Moran＇s I of function transition was higher than that of space transition,and the spatial aggregation characteristics of function transition were more significant than for space transition.The Global Moran＇s I at each time period tended to increase,and the spatial aggregation degree of the function transition and space transition was further enhanced.Hot and cold spots of space transition of forest land stably evolved,suggesting a minor spatial difference in forest land among different administrative units at the county level.The number of hot spots at the county level in function transition increased.Hot spots were intensively distributed at the western edge and continuously distributed in the northeast.The space transition and function transition of forest land were both greatly influenced by urbanization rate and second and third industries.The development of urbanization and industrialization was the main factor driving forest transition,suggesting a positive role of economic growth on forest transition in mountainous areas.The development of urbanization and industrialization is an effective approach to forest transition in mountainous areas.

Man-machine verification of mouse trajectory based on the random forest model

Identifying code has been widely used in man-machine verification to maintain network security.The challenge in engaging man-machine verification involves the correct classification of man and machine tracks.In this study,we propose a random forest(RF)model for man-machine verification based on the mouse movement trajectory dataset.We also compare the RF model with the baseline models(logistic regression and support vector machine)based on performance metrics such as precision,recall,false positive rates,false negative rates,F-measure,and weighted accuracy.The performance metrics of the RF model exceed those of the baseline models.

Dynamics of the Alpine Treeline Ecotone under Global Warming:A Review

The alpine treeline ecotone is defined as a forest-grassland or forest-tundra transition boundary either between subalpine forest and treeless grassland,or between subalpine forest and treeless tundra.The alpine treeline ecotone serves irreplaceable ecological functions and provides various ecosystem services.There are three lines associated with the alpine treeline ecotone,the tree species line(i.e.,the highest elevational limit of individual tree establishment and growth),the treeline(i.e.,the transition line between tree islands and isolated individual trees)and the timber line(i.e.,the upper boundary of the closed subalpine forest).The alpine treeline ecotone is the belt region between the tree species line and the timber line of the closed forest.The treeline is very sensitive to climate change and is often used as an indicator for the response of vegetation to global warming.However,there is currently no comprehensive review in the field of alpine treeline advance under global warming.Therefore,this review summarizes the literature and discusses the theoretical bases and challenges in the study of alpine treeline dynamics from the following four aspects:(1)Ecological functions and issues of treeline dynamics;(2)Methodology for monitoring treeline dynamics;(3)Treeline shifts in different climate zones;(4)Driving factors for treeline upward shifting.