Loading...
Welcome to Chinese Journal of Applied Ecology! Today is Share:

Table of Content

    15 December 2021, Volume 32 Issue 12
    Responses of soil erosion to changes in landscape pattern and its evolution in watershed in the loess hilly region under characteristic management and development
    LUO Jia-ru, LI Bin-bin, ZHANG Feng-bao, CONG Pei-juan, WANG Hai-yan, YANG Ming-yi
    2021, 32(12):  4165-4176.  doi:10.13287/j.1001-9332.202112.006
    Asbtract ( 497 )   PDF (4308KB) ( 339 )  
    Understanding the relationship between soil erosion and the changes in landscape patterns is important for guiding the management and development of watersheds. The Nangou watershed in Ansai County, Shaanxi on the Loess Plateau, is an area with the implement of “Grain for Green”, ecological agriculture, ecological tourism and the demonstration of science and technology for landscape management. We quantified the spatial and temporal variations of landscape pattern and soil erosion from 1981 to 2018 using the GIS and the universal soil loss equation (USLE). The relationship between the soil erosion modulus and nine landscape pattern indices in three categories was analyzed using the principal component analysis at both plot and landscape levels. The results showed that, among the five landscape types, the spatial and temporal changes of cultivated land and woodland dominated the evolution of landscape patterns, which affected the concentration and distribution degree of the whole watershed. Soil erosion in the Nangou watershed decreased annually, with erosion area, erosion modulus and soil erosion intensity decreasing by 29.7%, 61.2%, and 73.4% from 1981 to 2018, respectively. The variation of cultivated land and forest land areas determined the changes of soil erosion modulus of the whole watershed. The change trend of landscape pattern index was consistent with that of soil erosion. “Grain for Green” Project was the major dri-ving force for the changes in the landscape pattern and for the reduction of soil erosion. The characteristic development and management could weaken soil erosion intensity in parts of the study area. The rational configuration of landscape types could effectively control soil erosion in a watershed. The combination of rational configuration and characteristic management could help achieve the goals for sustainable and high quality development of the watershed.
    Impacts of seepage flow and soil thaw depth on hillslope snowmelt erosion in Chinese Mollisol region
    WANG Lun, ZHENG Fen-li, SHI Hong-qiang, ZHAO Lu-you, MO Shuai-hao, QIN Qi-shan, GENG Hua-jie, ZHAO Ya-jun
    2021, 32(12):  4177-4185.  doi:10.13287/j.1001-9332.202112.020
    Asbtract ( 424 )   PDF (2051KB) ( 79 )  
    Snowmelt erosion is an important way of soil loss in Chinese Mollisol region. However, little is known about the effects of seepage flow and soil thaw depth on hillslope snowmelt runoff erosion. An indoor simulated experiment was conducted to analyze the impacts of seepage flow and soil thaw depth on hillslope snowmelt erosion. There were two snowmelt flow rates (1 and 4 L·min-1), two soil thaw depths (5 and 10 cm), and two near-surface hydrological conditions (with and without seepage flow). The results showed that hillslope runoff depth and soil erosion amount in the treatment with seepage flow were 1.1 to 1.2 times and 1.3 to 1.9 times of those in the treatment without seepage flow, respectively. Under two snowmelt flow rates, when soil thaw depth increased from 5 cm to 10 cm, hillslope runoff depth and soil erosion amount increased by 10.0% to 13.5% and 15.4% to 37.1% in the treatment without seepage flow, respectively. In the treatment with seepage flow, when soil thaw depth shifted from 5 cm to 10 cm, hillslope runoff depth increased by 6.5% to 8.5%, and soil erosion amount remained stable. Moreover, hillslope rill development was comprehensively influenced by seepage flow, soil thaw depth, and snowmelt flow rate, with rill erosion amount occupying more than 72% of hillslope snowmelt erosion amount. Compared with the treatment without seepage flow, flow velocity and shear stress under the treatment with seepage flow increased by 20.3% to 23.2% and 37.0% to 51.3%, respectively; but Darcy-Weisbach friction coefficient reduced by 9.0% to 21.4%, which caused an increase of hillslope snowmelt erosion. In addition, seepage flow enhanced rill development, which caused rill erosion amount to increase by 43.6% to 69.9% compared with the treatment without seepage flow, and it further resulted in the increase of hillslope snowmelt erosion amount. The main reason for soil thaw depth enhancing hillslope snowmelt erosion amount under the treatment without seepage flow was that both sloping runoff erosivity and erodible materials increased with increasing soil thaw depth. Furthermore, soil thaw depth had a significant impact on hillslope rill morphology development under the treatment with seepage flow. Rill widening process was dominated when soil thaw depth was 5 cm, whereas rill incision process was dominant when soil thaw depth was 10 cm. This study could improve the understanding of hillslope snowmelt erosion mechanism in Chinese Mollisol region and provide theoretical guidance for the development of water erosion model.
    Effects of biological soil crusts on the initial abstraction ratio of SCS-CN model in the Loess Plateau
    GU Kang-min, ZHAO Yun-ge, GAO Li-qian, YANG Kai, SUN Hui, GUO Ya-li
    2021, 32(12):  4186-4194.  doi:10.13287/j.1001-9332.202112.005
    Asbtract ( 215 )   PDF (1676KB) ( 39 )  
    Hydrological model is an effective tool for hydrological research. The initial abstraction ratio (λ) is a key parameter of SCS-CN model, a commonly used runoff model of great significance to simulate the hydrological process at the watershed scale. In order to examine the effects of biological soil crusts (biocrusts) on λ and improve the accuracy of the model used in the restored grasslands where biocrusts widely presented in the Loess Plateau region, we firstly determined the relationship between the amount of the potential maximum infiltration (S) and the amount of the actual infiltration (F), and then investigated the effects of biocrust coverage on λ by using the simulated rainfall experiment in the Yingwoshanjian watershed in Dingbian County, Shaanxi Province. The revised model was verified by the runoff results of the simulated rainfall experiments in the Zhifanggou watershed in Ansai County, Shaanxi Province. The results showed that the relationship between S and F on biocrust slope was described as S/F=2.5×60/T (where T was the rainfall duration). There was a negative correlation between λ and biocrust coverage (CBSC) described as λ=0.0791×e(-0.015×CBSC), R2=0.60. Compared with that using the standard value of λ, the efficiency coefficient of the model was increased by 338.7% and the qualified rate was increased by 16.1% after revising λ according to the biocrust coverage. The results provided a scientific basis for the calibration of λ on biocrust slopes in the Loess Plateau region, and were of great significance to accurately assess the hydrological effects of the implementation of the “Grain for Green” Program on the Loess Plateau.
    Responses of subalpine meadow to climatic factors and the time lag effects in Wuyi Mountains from 2000 to 2019
    AN De-shuai, XU Dan-dan, PU Yi-han, WANG Hao-bin, LIU Yan-qing, ZHU Jian-qin
    2021, 32(12):  4195-4202.  doi:10.13287/j.1001-9332.202112.001
    Asbtract ( 277 )   PDF (2601KB) ( 125 )  
    Subalpine meadow is extremely sensitive to climate change. Few studies, however, focused on the responses of subalpine meadow to climatic factors in subtropical regions. It is still a challenge to extract the biophysical parameters from optical remote sensing imagery during the growing season. Based on the normalized difference vegetation index (NDVI) time series dataset from the MOD13Q1 vegetation index products and meteorological data, we analyzed the changes of vegetation growth of subalpine meadow at Huanggang Mountain in the top of Wuyishan National Park from 2000 to 2019, its responses to climate factors and the time lag effects. The results showed that NDVI in summer increased insignificantly during 2000-2019, and that NDVI in the growing season, spring, and autumn increased significantly. The enhancement of NDVI was mainly contributed by the increasing temperature (0.026 ℃·a-1) from 2000 to 2019. The increasing temperature in spring and autumn influenced meadow growth more than that in summer and growing season. NDVI of the growing season in subalpine meadow was sensitive to precipitation,indicating that the growth of subalpine meadow was strongly affected by precipitation even in the subtropical region with sufficient precipitation. Temperature and precipitation in different growth periods had different time lag effects on the NDVI of subalpine meadowo. The time lag effects of temperature on subalpine meadow were 0-1 month, and that of precipitation were 2-3 month.
    Temporal and spatial variations of soil moisture in North China
    LI Xuan, GUO Zhi-feng, WU Men-xin, HE Yan-bo
    2021, 32(12):  4203-4211.  doi:10.13287/j.1001-9332.202112.008
    Asbtract ( 340 )   PDF (6407KB) ( 76 )  
    Soil moisture is an important hydrological parameter and a basic element for research in water cycle and climate change. Using hourly recorded soil moisture data of 374 stations from the new generation of automatic soil moisture encryption observation network constructed by China Meteo-rological Administration, we analyzed the spatial and temporal variations of soil moisture in North China from 2013 to 2019 and the relationship with precipitation and temperature. The results showed that soil moisture in 10-100 cm layer decreased fluctuatly as a whole, with the decreasing at the 100 cm being serious. The spatial distribution of soil moisture at different depths was characterized by high in the southeast and low in the northwest. About 63% of the surface layer was short of moisture. Soil moisture at different depths changed significantly with the seasons. In summer, soil moisture of each layer reached the highest and soil entropy of each layer was suitable, while it reached a low point in spring. Soil moisture was closely correlated with precipitation and temperature, but the correlation got weaker gradually with the increase of soil depth. Soil moisture was more sensitive to precipitation than to temperature.
    Growth stability of four drought resistant plant species in different regions
    YU Xiao, JI Ruo-xuan, CHANG Yuan, SHEN Chao, GUO Hui-hong, XIA Xin-li, YIN Wei-lun, LIU Chao
    2021, 32(12):  4212-4222.  doi:10.13287/j.1001-9332.202112.011
    Asbtract ( 263 )   PDF (880KB) ( 65 )  
    There are abundant germplasm resources of drought resistant trees in China. It is difficult for foresters to evaluate and screen excellent germplasm that is suitable for various drought adverse circumstances. In this study, four tree species from different provenances, namely Amygdalus davi-diana, Prunus sibirica, Salix gordejevii, and Caryopteris mongolica, were used as test materials. Four soil regions, namely Dalad Banner, Siziwang Banner, Kouhezi town and Liujiazi town of Kulun Banner in the Central and Eastern Inner Mongolia Autonomous Region were selected as multi-point experiment sites to analyze the growth and physiological status of different tree species and provenances. The additive main effects and multiplicative interaction model was used to evaluate the regional adaptability and stability of the germplasm. The growth and physiological indices of four tree species differed between provenances and locations. Soil conditions (potassium content, nitrogen content, and pH) and climate conditions (annual average temperature, precipitation, and potential evaportranspiration) in different locations all influenced the growth of different provenance species. Concerning tree species, S. gordejevii and C. mongolica are more adapted to the sandy loam and chestnut soil of Dalad Banner and Siziwang Banner. A. davidiana and P. sibirica L. are more adapted to the loess and aeolian sandy soil of Kouhezi town and Liujiazi town. Concerning tree provenances, A. davidiana of Tuzuo, P. sibirica of Ningcheng and Yuanzhou, S. gordejevii of Lanqi, and C. mongolica of Jingbian displayed higher regional stability and better growth adaptability, indicating their suitability for afforestation in similar areas.
    Spatiotemporal variations of wet-bulb temperature and its impact factors of Nanjing urban neighborhood
    LU Yang, YANG Yi-chen, WU Xin-ran, WANG Lei, WAN Zi-wen, ZHENG Zhe-jun, XU Jia-ping, CAO Chang
    2021, 32(12):  4223-4236.  doi:10.13287/j.1001-9332.202112.012
    Asbtract ( 262 )   PDF (11838KB) ( 75 )  
    Global climate change and local urban heat islands enhance urban heat stress. Studies focused at the urban neighborhood scale are limited. Wet-bulb temperature represents the combined effects of both temperature and humidity, and therefore can more accurately reflect human thermal comfort. In this study, air temperature, relative humidity and geographic information of different times, seasons, and sky conditions of the Nanjing Jiangbei New Area were obtained based on mobile measurements. The spatiotemporal variation of wet-bulb temperature at the urban neighborhood scale and the effects of sky conditions, land cover and urban morphology (sky view factor, SVF) were further analyzed. The results showed that: 1) the spatiotemporal variations of wet-bulb temperature at the Nanjing urban neighborhood scale were consistent with that of air temperature. Compared with vapor pressure, air temperature played a dominant role. The extremely high values of wet-bulb temperature in this area were mostly caused by the synergy between air temperature and vapor pressure. 2) The correlation between SVF and wet-bulb temperature was significantly positive in the daytime and negative at night. An increase in the vegetation fraction could reduce wet-bulb temperature, while impervious surfaces had the opposite effect. The wet-bulb temperature significantly decreased and its spatial distribution was much more homogeneous under overcast sky conditions. 3) The horizontal scale effect showed diurnal and seasonal differences and was more sensitive to sky conditions during nighttime than during daytime. Compared with vegetation, the horizontal effect of impervious surfaces was much larger in winter than in the other two seasons. The horizontal scale effects of vege-tation and impervious surfaces on wet-bulb temperature were similar to those of air temperature. These results could provide effective scientific support and a theoretical basis for improving and optimizing the thermal environment of urban neighborhoods, as well as alleviating urban heat stress.
    Soil properties and microbial respiration activities of riparian forest wetland in the north of permafrost zone, the Great Hing’an Mountains, Northeast China
    WANG Xian-wei, TAN Wen-wen, SONG Chang-chun, DU Yu, ZHANG Hao, CHEN Ning
    2021, 32(12):  4237-4246.  doi:10.13287/j.1001-9332.202112.017
    Asbtract ( 235 )   PDF (1851KB) ( 69 )  
    Riparian wetlands in permafrost regions are critical regions for hydrological, ecological, and biochemical processes. We studied the soils of riparian and transition wetlands and analyzed physicochemical properties, stoichiometry, and microbial respiration activities (microbial biomass carbon, basal respiration, microbial entropy, and metabolic entropy) of the humus layer and diffe-rent soil layers. The results showed that the main differentiation of soil physical and chemical pro-perties in riparian forest wetlands was below 20 cm. Compared to the wetlands of transition zone, total carbon content, total nitrogen content, C/P and N/P decreased significantly with soil depth in riparian forest wetlands. These changes in soil stoichiometry were mainly caused by soil nitrogen content. Such a result meant that the transferring of nitrogen was relatively fast and that there was nitrogen limitation. The main differentiation of Na, Mg, K and Ca in soil occurred in the 30 cm layer of the transition zone and the 20 cm layer of the riparian forest wetlands. The correlations between soil Mg content and total C, total N, total P contents were significant. It meant that the soil Mg was an important element to riparian wetlands in the Great Hing’an Mountains. Microbial respiration activities of the humus layer in riparian forest wetlands and transition zone were higher than those in the other soil layers, indicating that the content of labile carbon fractions was high. The correlations between soil microbial respiration activities and soil properties, stoichiometry, nutrient elements were different in riparian wetland and transition zone. Soil total nitrogen contents were significantly correlated with soil microbial respiration activities in riparian wetland, indicating that soil microbial respiration activities were limited by nitrogen in riparian wetland of the Great Hing’an Mountains.
    Effects of wetland reclamation on amino sugar accumulation in soils of the Sanjiang Plain
    CHEN Qi, DING Xue-li, ZHANG Bin
    2021, 32(12):  4247-4253.  doi:10.13287/j.1001-9332.202112.018
    Asbtract ( 255 )   PDF (1131KB) ( 80 )  
    Microbial residues play important roles in the formation and stabilization of soil organic matter. The responses of soil microbial residues to wetland reclamation, however, remain unclear. In this study, we collected soil samples from a native wetland dominated by Calamagrostis angustifolia as well as three agricultural lands cultivated for 5, 15, and 25 years to examine the effects of wetland reclamation on the content of microbial residues (amino sugar analysis). Results showed that soil amino sugar contents were significantly reduced after wetland reclamation, with a positive relationship between the reduction and the duration of reclamation. After 25 years of reclamation, the content of glucosamine, galactosamine, and muramic acid in agricultural soils decreased by 38.0%, 38.1%, and 35.9%, respectively, compared to the natural wetland. The reduction of muramic acid (25.8%) was stronger than glucosamine (14.9%) after 5 years of reclamation, indicating that bacterial-derived microbial residues were more sensitive to wetland reclamation than fungal-derived counterparts in the short term. Total amino sugars were decreased by 21.1%, 34.0%, and 38.0% after 5, 15, and 25 years of wetland reclamation, respectively. The proportion of total amino sugars in soil organic matter was significantly decreased from 4.8% in natural wetland to 4.4% after 25 years of reclamation, indicating that long-term wetland reclamation accelerated the depletion of microbial-derived organic components, and thus changed the composition of soil organic matter. Such changes may affect long-term stabilization of soil organic matter and ecosystem functioning.
    Effect of cropping patterns onbacterial community diversity and metabolic function in rhizosphere soil of Angelica sinensis
    JIANG Xiao-feng, GUO Feng-xia, CHEN Yuan, GUO Jian-guo, LIU Xiao-feng
    2021, 32(12):  4254-4262.  doi:10.13287/j.1001-9332.202112.031
    Asbtract ( 364 )   PDF (1368KB) ( 59 )  
    The sustainable development of Angelica sinensis industry is seriously restricted by continuous cropping obstacles. In order to explore an efficient cultivation technique for A. sinensis, an experiment with five cropping patterns [A: Pisum sativum (Ps)-A. sinensis (As)-As, control); B: Ps-Triticum aestivum (Ta)-As; C: Ps-Mongolia astragalus (Ma)-As; D: Ps-Solanum tuberosum (St)-As); E: Ps-Fallow (F)-As)] were conducted in major A. sinensis producing areas located in Weiyuan County, Gansu Province. The physicochemical properties and relative abundance of bacterial genomic DNA in rhizosphere soil under different cropping patterns were measured during A. sinensis harvest period to investigate the effects of different cropping patterns on physicochemical properties, bacterial community diversity, and metabolic pathways. The results showed that: 1) the physicochemical properties in A. sinensis rhizosphere soil varied among different cropping patterns. Compared with the control, soil electrical conductivity under C pattern was significantly higher, and lower under B, D and E, CO2 respiration rate for B, C, D and E were significantly increased. 2) Soil bacteria of A. sinensis rhizosphere soil in the five cropping patterns belonged to 26 phyla and 368 genera. The dominant genera were Gemmatimonas from Gemmatimonadetes, Sphingomonas from Proteobacteria, and Subgroup_6 from Acidobacteria. Compared with the control, the relative abundance of Proteobacteria and Actinobacteria under B and C patterns was significantly higher, Acidobacteria in D pattern was significantly lower, while Proteobacteria, Acidobacteria, and Actinobacteria in E pattern was significantly higher. 3) There were significantly negative relations between soil pH, electrical conductivity, contents of organic matter, available nitrogen, phosphorus and potassium with the relative abundance of Proteobacteria in A. sinensis rhizosphere soil across the five cropping patterns. 4) There was significant difference in relative abundance for bacteria of six metabolic pathways under the five cropping patterns. In conclusion, C pattern had a regulating effect on physicochemical properties and bacterial communities in A.sinensis rhizosphere soil, which could be taken as a major practice to overcome the continuous cropping obstacles.
    Effects of continuous cotton monocropping on soil physicochemical properties and nematode community in Xinjiang, China
    CHEN Hong, YANG Lei, ZHANG Feng-hua
    2021, 32(12):  4263-4271.  doi:10.13287/j.1001-9332.202112.026
    Asbtract ( 242 )   PDF (1150KB) ( 51 )  
    The community characteristics of soil nematodes are an important basis for evaluating and indicating soil health. In this study, cotton fields with different continuous monocropping years (5, 10, 15, 20, and 25 years) were selected. High-throughput sequencing technology was used to explore the responses of soil properties and nematode communities to long-term continuous monocropping in cotton fields. The results showed that after 10-15 years, soil pH and soil conductivity increased significantly, and the contents of soil organic carbon, total nitrogen, available phosphorus (AP), available potassium, nitrate and soil microbial biomass carbon (MBC) decreased significantly. A total of 25 genera, 18 families, 7 orders and 3 classes of soil nematodes were identified. Among them, Helicotylenchus was the dominant genus. Parasitic nematodes were the dominant trophic groups, showing a trend of first decreasing and then increasing. Compared with other continuous monocropping years, plant parasitic nematodes increased by 9.1%-208.6% and Helicotylenchus increased by 392.0% under continuous monocropping for 25 years. With the increases of continuous monocropping duration, plant parasitic nematodes such as Tylenchorhynchus, Ditylenchus, Discopersicus, Mesocriconema, and Criconemoides were detected. After continuous monocropping for 15 years, soil nematode richness index and maturity index of free-living nematodes (MI) decreased significantly, PPI/MI increased significantly, and diversity index (Shannon) and Wasilewska index were the lowest. Soil AP and MBC were the main environmental factors affecting the changes of soil nematode communities. Our results suggested that, after 10-15 years continuous monocropping, soil nutrients were unbalanced, the diversity of soil nematodes was reduced, the stability of soil food network was deteriorated, and the parasitic nematodes were increased, which resulted in continuous monocropping obstacles in cotton.
    Scaling laws of altitudinal pattern of soil fauna diversity in Dongling Mountain, Beijing, China
    DING Zhang-qi, XU Guo-rui, ZHANG Shuang, ZHANG Yu-xin, MA Ke-ming
    2021, 32(12):  4272-4278.  doi:10.13287/j.1001-9332.202112.030
    Asbtract ( 196 )   PDF (1539KB) ( 46 )  
    Identifying the complexity of diversity pattern of various taxa within a community is a challenge for ecologist. Scaling law is one of the suitable ways to detecting the complex ecological structure. In this study, we explored the scaling laws of soil fauna diversity pattern along an altitudinal gradient by multifractal analysis, and compared the difference of multifractal spectra between the litter and the soil layers. Consistent with results from plant communities in previous studies, there was power law scaling law for soil fauna diversity, i.e., richness, the exponential of Shannon’s Diversity Index, and the inverse Simpson’s Diversity Index. Moreover, power law scaling law also existed for the richness changes of different relative abundance species in both litter and soil layers. Although multifractal characteristics existed for both litter layer and soil layer of soil fauna diversity, the fractal structure of the diversity in the litter layer was more even than that in the soil layer, and the scaling properties of dominant and rare species showed different patterns in multifractal spectra between litter layer and soil layer. In conclusion, there were power law scaling laws for soil fauna diversity which had high richness and abundance along the altitudinal gradient, which would help us uncovering the spatial distribution mechanism of belowground biodiversity.
    Effects of soil meso- and micro-fauna on litter decomposition under nitrogen deposition and rainfall changes
    LIANG Zhi-wei, HONG Mei, DE Hai-shan, YE He, ZHANG Yuo-chen, YAN Jin, LI Jing
    2021, 32(12):  4279-4288.  doi:10.13287/j.1001-9332.202112.041
    Asbtract ( 300 )   PDF (2586KB) ( 67 )  
    Soil fauna are important regulators of litter decomposition and nutrient transformation. Nitrogen deposition and rainfall changes driven by global changes could affect litter decomposition by changing environment and soil faunal community. Different mesh size (2 mm and 0.01 mm) litter bags were used to explore how soil meso- and micro-fauna contribute to decomposition of Stipa breviflora litter under nitrogen deposition and rainfall changes. The experiment followed split-plot design, with rainfall change (natural rainfall, CK; rainfall addition 30%, W; rainfall reduction 30%, R) as the main trement and nitrogen addition (0, N0; 30, N30; 50, N50; 100 kg·hm-2·a-1, N100) as the sub-treatment. The results showed that: 1) Rainfall change significantly affected litter decomposition rate, which was increased by rainfall addition. Moreover, litter decomposition rate was accelerated with increasing nitrogen addition rates. Litter residual rate decreased gradually with increasing N addition, and got to the highest in N100. Litter decomposition rate decreased first and then increased, and peaked in N50 in rain reduction and natural rainfall treatment. There was no significant interactions between rainfall change and nitrogen addition in affecting litter decomposition. 2) During the whole decomposition process, a total of 1577 soil meso- and micro-fauna were captured, belonged to 1 phyla, 3 classes, 13 orders (including suborders) and 49 families. The dominant groups were Acarina, Coleoptera larvae, and Collembola. Nitrogen addition significantly increased abundance and group numbers of soil meso- and micro-fauna. 3) The litter mass residue rate was significantly negatively correlated with abundance and group numbers of soil meso- and micro-fauna. The contribution rate of soil meso- and micro-fauna to litter decomposition increased with increasing rainfall. In summary, soil meso- and micro-fauna had a positive effect on decomposition of Stipa breviflora litter in desert steppe. Their contribution to litter was promoted by the enhancement of soil mesofauna abundance and group under rainfall and nitrogen addition. Excessive nitrogen would inhibit soil meso- and micro-fauna community and group density when water was insufficient, and would thus weaken the function of soil mesofauna to litter decomposition.
    Fungal community diversity and driving factors in rhizosphere soil of Caragana species across semi-arid regions
    LI Yuan-yuan, XU Ting-ting, AI Zhe, MA Fei
    2021, 32(12):  4289-4297.  doi:10.13287/j.1001-9332.202112.037
    Asbtract ( 210 )   PDF (1236KB) ( 60 )  
    To investigate fungal community diversity in rhizosphere soil of Caragana species, and their relationships with ecological factors in the semi-arid regions of China, we collected soil samples from Caragana liouana and C. microphylla in different habitats. By using Illumina MiSeq high-throughput sequencing technology, we analyzed fungal diversity, community composition, and the driving factors. Results showed that 7 phyla, 20 classes, 43 orders, 66 families, and 78 genera were recorded from the rhizosphere soil. The dominant phyla in all sample sites were Ascomycota (37.7%), Basidiomycota (13.7%), and Zygomycota (4.3%). At the genus level, Penicillium, Geomyces, and Mortierella were dominant, and Rhizophagus and Glomus were also found in rhizosphere soil. The fungal Chao1 index, ACE index, and Simpson index of C. microphylla were significantly higher than those of C. liouana. There was significant difference between the two species in the composition of fungal communities. Results of redundancy analysis showed that soil organic carbon, total nitrogen, electrical conductivity, available potassium, altitude, total phosphorus and aridity index exhibited significant impacts on soil fungal diversity. Overall, those results enhanced our understanding of the relationships between rhizospheric microbial community diversity of Caragana species and ecological factors, and provided important information on the adaptive mechanisms of desert plants in semi-arid regions.
    Accumulation, subcellular distribution, and chemical forms of zinc in three tree species
    TANG Min, ZHANG Xin, TAN Xin-rui, LIU Yan, WANG Mei-xian
    2021, 32(12):  4298-4306.  doi:10.13287/j.1001-9332.202112.033
    Asbtract ( 225 )   PDF (2705KB) ( 38 )  
    In order to explore the mechanism underlying zinc (Zn) accumulation and tolerance in woody garden species, the effects of different Zn concentrations (0, 250, 500, 1000, 2000 mg·kg-1) on leaf, branch, root biomass and leaf ultrastructure of Koelreuteria paniculata, Ailanthus altissima, and Ginkgo biloba were studied in a pot pollution simulation experiment. The concentration of Zn in plant organs, the subcellular distribution and chemical forms of Zn in leaves and roots were further analyzed. The results showed that all the three species could survive under diffe-rent Zn concentrations, but the biomass of leaves, stems and roots decreased compared with the control. Excessive Zn could lead to cell deformation, cell wall rupture and organelle disintegration of leaves in K. paniculata and A. altissima, while the cells in leaves of G. biloba could maintain normal morphology, indicating that G. biloba had a better tolerance to Zn than K. paniculata and A. altissima. With the increases of Zn concentration, Zn concentration in the organs of the three species showed an increasing trend, and the Zn concentration in K. paniculata and A. altissima was significantly higher than that in G. biloba, indicating that the Zn accumulation ability of K. paniculata and A. altissima was stronger than that of G. biloba. Zn was mainly distributed in the cell walls of leaves and roots, accounting for 26.9%-71.8% and 28.1%-82.6%, respectively. Under the treatment with the highest Zn concentration (2000 mg·kg-1), Zn concentration in the soluble components (mainly vacuoles) could be higher than that in the cell walls. In addition, Zn mainly existed in NaCl-, HAc- and HCl-extracted forms in leaves, accounting for 57.4%-82.7%, and Zn mainly existed in NaCl- and HAc-extracted forms in roots, accounting for 42.8%-67.2%, all of which were forms with relatively low activity. Therefore, cell wall retention, vacuoles segregation and accumulating Zn in less active forms might be important mechanisms underlying Zn accumulation and tolerance in the three trees.
    Potential geographical distribution of Rana hanluica in China under climate change
    XIA Xin, LI Yuan, YANG Dao-de, PI Yang-yan
    2021, 32(12):  4307-4314.  doi:10.13287/j.1001-9332.202112.003
    Asbtract ( 378 )   PDF (2638KB) ( 153 )  
    Global warming in the last few decades had strong impacts on biodiversity and geographi-cal distribution of different animal species worldwide, especially amphibians. Rana hanluica, a frog species endemic in China, is still classified as Least Concerned in the Red List of Threatened Species because few studies have been conducted on this species. To understand the survival of Rana hanluica population, we used maximum entropy models (MaxEnt) to analyze its distribution across regions under current climatic conditions based on 47 distribution records and 20 environmental factors. We investigated the changes in distribution of this species under different climate scenarios in China (2050s and 2070s). Finally, current and future suitable habitats for R. hanluica were mode-led, and the impacts of environmental factors in shaping its distribution were evaluated. The results showed that the prediction accuracy of the MaxEnt model was high, and AUC value of the receiver operating curve was 0.993. The total suitable habitat area for R. hanluica was 36.36×104 km2, mainly located in Hunan and Guizhou provinces in China. The major environmental factors influencing the geographic distribution of R. hanluica were precipitation of dryest month and altitude. Under the future climate scenario (2050 and 2070) with two representative concentration pathways (RCPs, SSP1-2.5, SSP5-8.5), the suitable habitat of R. hanluica was reduced in different degrees, resulting in a decreasing trend of the total suitable habitat area. The center of gravity in highly suitable habitat of R. hanluica shifted to high-latitude regions, with the core distribution area in Hunan Province.
    Applicability of multiple remotely sensed vegetation indices for extracting key phenological metrics of Tamarix chinensis shrubs based on CO2 flux observation and Sentinel-2 data
    ZHOU Hao-qiang, BAO Gang, JIN Hugejiletu, DU Ling-tong, ZHANG Si-lian, XU Zi-wei, BAO Yu-hai
    2021, 32(12):  4315-4326.  doi:10.13287/j.1001-9332.202112.007
    Asbtract ( 216 )   PDF (6746KB) ( 35 )  
    We analyzed the relationship between gross primary productivity (GPP) and environmental factors at Sidaoqiao Superstation of the Ejina Oasis in China’s Gobi Desert, by combining eddy flux and meteorological data from 2018 to 2019 and Sentinel-2 remote sensing images from 2017 to 2020. We evaluated the applicability of 12 remote sensing vegetation indices to simulate the growth of Tamarix chinensis and extract key phenological metrics. A seven-parameter double-logistic function (DL-7) + global model function (GMF) was used to fit the growth curves of GPP and vegetation indices. Three key phenological metrics, i.e., the start of the growing season (SOS), the peak of the growing season (POS), and the end of the growing season (EOS), were extracted for each year. Growing season degree days (GDD) and soil water content were the main environmental factors affecting the phenological dynamics of T. chinensis. Compared with 2018, the lower temperatures in 2019 resulted in slower accumulation rate of accumulated temperature before the SOS. T. chinensis required longer heat accumulation to enter growing season, which might cause later SOS in 2019. The hydrothermal conditions between SOS and POS were similar for 2018 and 2019. Howe-ver, the POS in 2019 was 8 days later than that in 2018, because of the late SOS in 2019. Following the POS in 2019, high GDD and low soil water content caused the T. chinensis to suffer from water stress, resulting in a shortened late growing season. The linear regression between the standardized Sentinel-2 vegetation index and the average value of GPP between 10:00 and 14:00 indicated that the enhanced vegetation index of the broadband vegetation index and the chlorophyll red edge index, inverted red edge chlorophyll index, and red-edge normalized difference vegetation index (NDVI705) of the narrowband vegetation index were highly consistent with the GPP of T. chinensis. Remote sensing extraction of SOS and POS of T. chinensis suggested that the Sentinel-2 narrowband vegetation index was more accurate than the broadband vegetation index. The modified chlorophyll absorption in reflectance index provided the most accurate extraction of SOS, while the MERIS terrestrial chlorophyll index provided the most accurate extraction of EOS. Conversely, the Sentinel-2 broadband vegetation index was the most accurate for extracting POS, especially the 2-band enhanced vegetation index and the near-infrared reflectance of vegetation. Overall, NDVI705 was the best index to estimate phenological metrics.
    Modeling spatial distributions of hydrogen and oxygen stable isotopes in surface soils of the upper reaches of Minjiang River based on machine languages
    QIN Wen-yi, CHEN Guo, LI Xiao-zhen, WANG Xiang, WANG Peng
    2021, 32(12):  4327-4338.  doi:10.13287/j.1001-9332.202112.002
    Asbtract ( 221 )   PDF (4221KB) ( 48 )  
    To study the feasibility of simulating the spatial distribution of hydrogen and oxygen stable isotopes composition (δ2H and δ18O) in the surface soil based on the machine learning method and to investigate large-scale distribution of δ2H and δ18O in the upper reaches of Minjiang River, 183 soil samples were collected from the 0-10 cm soil layer. After variable selection, back propagation (BP) neural network, random forests (RF) and support vector machine (SVM) were used to model the δ2H and δ18O of the study area, with the accuracies being evaluated. The structural equation model (SEM) was used to reveal the mechanism between the auxiliary variables and the δ2H and δ18O of soil water. The results showed that the RF model had the highest prediction accuracy, and could explain 75.0% and 64.0% of the variations of δ2H and δ18O in the surface soil, respectively. In this model, soil water content was the most important auxiliary variable, contributing 48.9% and 37.4% to δ2H and δ18O. Vegetation factors had stronger influence on δ2H and δ18O in the surface soil than climate factors, and the influence of climate factors on δ2H and δ18O was media-ted by vegetation factors. Among all the auxiliary variables, hydrogen/oxygen isotope of precipitation had the lowest effect on δ2H and δ18O due to the fractionation. The δ2H and δ18O in the surface soil of the upper reaches of the Minjiang River changed significantly across different months during the growing season. The increases of δ2H and δ18O in the early growing season and the decreases in the late growing season were mainly affected by vegetation, while climate change led to a small fluctuation in the middle growing season.
    Identification of important biodiversity areas by InVEST model considering opographic relief: A case study of Yunnan Province, China
    YANG Wen-xian, LI Shi-hua, PENG Shuang-yun, LI Ying-xin, ZHAO Shou-lou, QIU Li-dan
    2021, 32(12):  4339-4348.  doi:10.13287/j.1001-9332.202112.004
    Asbtract ( 301 )   PDF (2902KB) ( 51 )  
    Accurately identifying important areas of biodiversity is one of the key issues in ecology and biodiversity research, as well as an important basis for the delineation of the red line for ecologi-cal protection and territorial spatial planning. With China’s typical plateau mountainous area (Yunnan Province) as a research case, we used the net primary productivity (NPP) quantitative index method, InVEST model and InVEST model focusing on topographic relief to identify biodiversity important areas. The results showed that NPP quantitative index method was not suitable for the plateau mountainous areas with obvious vertical zonal development. The identified area contained only 26.1% of the protected areas. The InVEST model had higher identification accuracy than the NPP quantitative index method in Yunnan Province. The identified area covered 49.4% of the protected natural areas. Fragmentation was obvious in northwest Yunnan. The InVEST model focusing on topographic relief improved the identification accuracy of important areas of biodiversity, including 71.7% of nature reserves. The deficiency of NPP quantitative index method in water area identification was made up and the fragmentation problem of InVEST model was solved. The area of biodiversity important areas was 119466.94 km2, accounting for 30.3% of the total land area of Yunnan Province. The spatial distribution showed a pattern of “three barriers, two zones and one region for multi-point development”.
    Effect of restoration types on the community structure of microbes harboring nifH and chiA genes in alpine meadow
    ZHANG Lu, WANG Jie, WANG Xiang-tao, LIAO Li-rong, WAN Qian, LIU Guo-bin, ZHANG Chao
    2021, 32(12):  4349-4358.  doi:10.13287/j.1001-9332.202112.040
    Asbtract ( 215 )   PDF (1771KB) ( 137 )  
    Biological nitrogen (N) fixation and organic N degradation are the main sources of soil available N, while microorganisms driving such processes play an important role in soil N supply and the maintenance of soil fertility. In this study, real-time quantitative PCR and amplicon sequencing technology were used to examine the effects of restoration types on the community structure of N2-fixing and chitin-degrading bacteria harboring nifH and chiA genes, respectively, and the gene abundance under four meadows (undisturbed, grazing, fencing, and fencing + reseeding mea-dows) in Qinghai-Tibet Plateau. The results showed that the abundance of nifH and chiA in the four meadows followed the order of undisturbed meadow > grazed meadow > fencing meadow > fencing + reseeding meadow. The abundance of nifH and chiA in the undisturbed meadow was 3.4-6.3 times and 3.3-8.3 times of that in the other three meadows. The α diversity of N2-fixing bacteria in gra-zing, fencing, and fencing + reseeding meadows was significantly higher than that in the undisturbed meadow, while the α diversity of chitin-degrading bacteria was higher in the undisturbed and grazing meadows. Grazing significantly increased the relative abundance of Proteobacteria, but decreased the relative abundance of Cyanobacteria and Actinobacteria. The abundance of nifH and chiA was significantly affected by soil moisture, nutrients, and vegetation characteristics, while the community structure of nifH and chiA was affected by soil moisture, soil organic carbon content, and soil pH. Compared with undisturbed meadow, grazing reduced the potential of N fixation and organic N degradation.The improvement of 10 years grazing prohibition with fencing and reseeding measures on the function of N fixation and organic N degradation was not obvious. The characteristics of functional microbes and their influencing factors should be comprehensively considered during meadow restoration, which might take longer time or take reasonable management measures to restore grazing meadow to undisturbed level.
    Effects of water limiting and nitrogen reduction on nitrogen use and apparent balance of winter wheat in the Guanzhong Plain, Northwest China
    GUO Zeng-hui, LIU Peng-zhao, LUO Wen-he, WANG Rui, LI Jun
    2021, 32(12):  4359-4369.  doi:10.13287/j.1001-9332.202112.022
    Asbtract ( 171 )   PDF (2762KB) ( 28 )  
    Effects of water limiting and nitrogen reduction on yield, nitrogen use efficiency and nitrogen apparent balance of wheat were investigated to explore whether it would be feasible to restrict water and reduce nitrogen in wheat production of the Guanzhong Plain and thus to provide scientific supports for yield-stable, high-efficiency, and environment-friendly developments in the irrigated production of winter wheat. Following a split-plot design with two water regimes as the main plots and four N addition rates as sub-plot factors, a field experiment (2017-2019) was conducted in Yangling, Shaanxi. The two water regimes were conventionally irrigating at the rate of 60 mm during the overwinter period and at the jointing stage, respectively (W2, a conventional practice) and irrigating at a rate of 60 mm during the overwintering period (W1, a restrictive irrigation practice). The four nitrogen addition rates were 300 kg·hm-2(N300, a conventional N rate), 225 kg·hm-2 (N225, a nitrogen rate of 25% less than the convention), 150 kg·hm-2(N150, a nitrogen rate 50% of less than the convention), and 0 kg·hm-2(N0, no nitrogen applied). The decreased irrigation rate and nitrogen rate significantly increased nitrogen content in the plants and grains, yield, N output, nitrogen use efficiency, nitrogen harvest index, nitrogen recovery efficiency, and nitrogen agronomic efficiency, reduced nitrate leaching and N surpluses, and maintained nitrogen balance. With both W1 and N150 adopted, the increased irrigation rate and nitrogen rate did not affect yield and N output of winter wheat in 2017-2019. Plant nitrogen content with both W1 and N150 adopted increased by 0.1%-25.5% and 14.0%-31.6% and the grain nitrogen content increased by 0.1% and 4.6%, compared with those with both W2 and N300 adopted in 2017-2018 and 2018-2019, respectively. Nitrogen use efficiency, nitrogen harvest index, nitrogen recovery efficiency, and nitrogen agronomic efficiency were averagely increased by 95.3%, 4.2%, 81.7% and 33.0% respectively. The N surplus was decreased by 97.2% and 95.1%, which effectively alleviated soil nitrate leaching. Considering all the indicators, irrigating at 600 m3·hm-2 during the overwintering period plus applying nitrogen at 150 kg·hm-2 could achieve high yield, high efficiency, and environment friendly development of winter wheat in the Guanzhong Plain of Shaanxi.
    Effects of elevated CO2 concentration on photosynthetic acclimation of winter wheat under drought condition
    ZONG Yu-zheng, YANG Qi, CHANG Cui-cui, GOU Jun-ying, ZHANG Dong-sheng, HAO Xing-yu, GAO Zhi-qiang
    2021, 32(12):  4370-4380.  doi:10.13287/j.1001-9332.202112.014
    Asbtract ( 210 )   PDF (3673KB) ( 44 )  
    Mechanisms underlying leaf photosynthetic acclimation in winter wheat under elevation of CO2 concentration ([CO2]) remain unclear. The aim of the study was to investigate the effects of source-sink variation on photosynthetic acclimation induced by drought under elevated [CO2]. A winter wheat (Triticum aestivum L. ‘Zhengmai 9023’) pot experiment was conducted in open top climate chambers with [CO2] of 400μmol·mol-1 or 600 μmol·mol-1 and soil water content at 80%±5% or 55%±5% of field capacity. The parameters of chlorophyll fluorescence, electron transport rate, photosynthetic curve, leaf nitrogen content, and grain yield were measured at the elongation and heading stages. Under drought condition, leaf PSⅡ photochemical efficiency was not affected by elevated [CO2], but the maximum electron transport rate and the ratio of electron partitioned to carboxylation reaction in Calvin cycle was increased at the elongation stage, and thus the Rubisco carboxylation rate and maximum photosynthetic rate were increased. Although the maximum electron transportation rate and partitioning ratio of electron to carboxylation reaction in Calvin cycle remained high at the heading stage, the PSⅡ photochemical efficiency, Rubisco carboxylation rate, and triose phosphate utilization rate were decreased by elevated [CO2], which consequently reduced the maximum photosynthetic rate for plant under drought stress. Under drought condition, elevated [CO2] increased wheat tiller biomass, kernel number, and kernel weight per ear, but decreased infertile kernel number, resulting in an overall increase in grain weight. In conclusion, the elevated [CO2]-induced increase in wheat grain yield per tiller under drought condition was mainly caused by enhanced photosynthetic performance at the elongation stage. The photosynthetic acclimation in source leaves during the heading stage under elevated [CO2] was mainly attributed to the reduction in PSⅡ photochemical efficiency and triose phosphate utilization rate, but not to the maximum electron transportation rate, ratio of electron partitioned to carboxylation in Calvin cycle or sink leaf strength.
    Osmotic regulation and chlorophyll fluorescence characteristics in leaves of wheat seedlings under different salt stresses
    ZHANG Ya, SHI Shu-qian, LI Ya-ping, GAO Tian-peng, YANG Ying-li
    2021, 32(12):  4381-4390.  doi:10.13287/j.1001-9332.202112.023
    Asbtract ( 258 )   PDF (2179KB) ( 72 )  
    The damage mechanism of salt stress on plants has attracted much attention. In order to reveal the damage mechanism of different salt stresses, we compared osmotic regulation and photosynthetic characteristics of seedlings of wheat cultivar Xianhan 3 under sodium salt (150 mmol·L-1) and calcium salt (5, 30 mmol·L-1) treatments alone or in combination. The results showed that sodium salt or calcium salt stress alone significantly inhibited the growth of roots and stems, but increased the amount of soluble sugar and proline, regulatory energy-dissipated electron yield, non-photochemical quenching and relative content of zeaxanthin contents in leaves. In contrast, salt treatments alone significantly decreased the levels of chlorophyll a and chlorophyll b, maximum photochemical efficiency, PSⅡ photochemical efficiency, photochemical quenching and photosynthetic electron transport efficiency. Furthermore, the inhibition of wheat seedling growth was more sensitive to calcium salt than to sodium salt stress, whereas the decreases of chlorophyll content and chlorophyll fluorescence parameters were more prominent in response to sodium salt stress. Except for the amount of soluble protein, lutein and the relative level of zeaxanthin, the changes of other parameters in the leaves due to sodium salt stress were effectively blocked by the application of low calcium concentration, but further increased by the presence of high calcium salt concentration. Taken together, sodium or calcium salt stress alone significantly inhibited seedling growth. The toxicity of sodium salt to wheat seedlings was effectively alleviated by low calcium concentration, but was aggravated by high calcium concentration, which were associated with the changes of photosynthetic pigment content, light energy capture, and photosynthetic electron transport process in the leaves of wheat seedlings. Moreover, osmotic regulators played an important role in enhancing the resistance of wheat seedlings to sodium or/and calcium environment.
    Effects of nitrogen reduction combined with organic materials on crop yield, photosynthetic characteristics, and product quality of corn-cabbage rotation system
    CHEN Yun-mei, ZHAO Huan, XIAO Hou-jun, XIE Ting-ting, QIN Song, HU Gang
    2021, 32(12):  4391-4400.  doi:10.13287/j.1001-9332.202112.015
    Asbtract ( 232 )   PDF (2329KB) ( 29 )  
    Nitrogen reduction combined with organic materials is an important measure to achieve or even increase crop yield retention at the background of fertilizer reduction. We conducted a pot experiment to explore the effects of nitrogen reduction combined with organic materials on yield, photosynthetic characteristics, and product quality of agricultural products of maize-cabbage rotation system in yellow soil area of Guizhou. There were five treatments, including no fertilizer (CK), conventional fertilizer (CF), nitrogen reduction (20%, the same below) combined with biochar (RF+B), nitrogen reduction combined with rapeseed cake (RF+O), and nitrogen reduction combined with both biochar and rapeseed cake (RF+BO). Leaf photosynthetic characteristics were measured in maize (seedling stage, jointing stage, heading stage, and mature stage) and cabbage (seedling stage, growing stage and harvest stage). The biological characters, yield and quality indices were investigated in the harvest period. Compared with CF, RF+BO significantly enhanced the yield of corn and cabbage by 9.7% and 39.2%, respectively, while RF+O had no effect, and RF+B did not affet maize yield. RF+BO improved the biological properties of maize and cabbage, including the 100-kernel weight of maize, and plant height, maximum leaf length and total biomass of cabbage. Furthermore, the green holding period and high photosynthetic duration of maize and cabbage were prolonged, among which, maize leaf SPAD was increased respectively by 42.7%, 11.0%, 12.8%, and 30.3% at seedling, jointing, heading, and mature stages, the cabbage leaf SPAD was increased by 13.5%, 9.2%, and 30.3% in seedling, growing and harvest stages, respectively. The net photosynthetic rate (Pn) of maize was increased by 11.1%, 10.9%, and 119.8% in seedling, jointing, and mature stages, while that of cabbage was increased by 12.7% and 14.6% in growing and harvest stages, respectively. The stomatal conductance (gs) of maize was increased by 58.3% and 41.7% in jointing and harvest stages, while that of cabbage was increased by 10%, 64.7%, and 19.2% in seedling, growing, and harvest stages, respectively. The transpiration rate (Tr) of maize was increased by 55.0%, 10.6%, 14.0%, and 143.9% in seedling, jointing, heading, and mature stages, respectively, while that of cabbage was increased by 26.1% in growing stage. The nutritional quality of maize and cabbage was significantly improved. The contents of reducing sugar, starch, and crude protein in maize were increased by 16.2%, 3.5% and 20.3%. The contents of Vc, amino acid, and reducing sugar in cabbage were increased by 26.3%, 21.0% and 27.8%, separately. In conclusion, 20% nitrogen reduction combined with biochar and rapeseed cake had positive effects on crop growth, yield increase, green retention period, high photosynthetic duration, and agricultural product quality improvement in Guizhou yellow soil maize-cabbage rotation system, the overall effect of which was the best. Nitrogen reduction combined with single organic material overally did not affect crop yield, photosynthetic characteristics, and quality.
    Chemical structure characteristics of organic carbon in saline soil aggregates under straw returning condition
    PEI Zhi-fu, HONG Mei, XING An, ZHANG Yue-xian, WEN Xin, ZHAO Hui-xin, SHEN Qin-guo
    2021, 32(12):  4401-4410.  doi:10.13287/j.1001-9332.202112.021
    Asbtract ( 338 )   PDF (2820KB) ( 66 )  
    We examined the regularity of distribution and chemical structure characteristics of organic carbon in soda alkaline fluvo-aquic soil aggregates after straw returning. We set up six different straw returning treatments in 2020, including 0 (CK), 2100 (ST1), 4200 (ST2), 6300 (ST3), 8400 (ST4) and 10500 kg·hm-2(full straw returning, ST5). We measured organic carbon (OC) content and infrared spectroscopy characteristics of aggregates and internal different components through physical fractionation method and infrared spectroscopy technology. The results showed that: 1) the OC content of soil and all aggregates increased with the increasing amount of returned straw; 2) different straw returning treatments significantly increased the content of light organic carbon (LOC) in 53-250 μm aggregates. Compared with CK, ST3 and ST4 treatments significantly increased the content of mineral-bound organic carbon (MOC) in 250-2000 μm aggregates and the content of fine particulate organic carbon (fPOC) in 53-250 μm aggregates. The OC content of different components in aggregates followed the order of LOC>MOC>POC. The fPOC content in 250-2000 μm aggregates was higher than that of coarse particulate organic carbon (cPOC); 3) the results of principal component analysis showed that OC chemical structure of different components in aggregates was seldom affected by the straw returning, but was mainly affected by particle size; 4) the OCs in >250 μm aggregates were mainly derived from aromatic carbon and polysaccharides. The OCs in 53-250 μm aggregates were mainly derived from carbohydrates, such as monosaccharides and polysaccharides, while the OC in <53 μm aggregates was mainly derived from aliphatic carbon, alkyl carbon, aromatic carbon and phenolic alcohols. Within different aggregates, LOC was mainly derived from aliphatic carbon, aromatic carbon and phenolic alcohols. Particulate organic carbon (POC) was mainly derived from carbohydrates. MOC was mainly derived from alkyl carbon. In summary, straw returning increased organic carbon content in soil aggregates in short term, but did not alter organic carbon chemical structure. The organic carbon chemical structures of the same particle size fractions in different aggregates were similar. The organic carbon content increased with the decreases of particle size, and the chemical structure tended to be stable. Therefore, straw returning promoted the fixation of organic carbon by saline soil aggregates in short term, but did not change their chemistry structural characteristics, indicating that the location and protection degree of soil organic carbon in aggregates were the main factors affecting the chemical structure of organic carbon.
    Effects of different amounts of organic acid soil conditioners on soil nutrients and crop growth
    CUI Heng, ZHANG Jiu-dong, BAO Lin, HAN Jie-rong, CHE Zong-xian, BAO Xing-guo, YANG Rui-ju
    2021, 32(12):  4411-4418.  doi:10.13287/j.1001-9332.202112.016
    Asbtract ( 246 )   PDF (1330KB) ( 31 )  
    Northwest China is burdened by declining soil fertility and poor capacity of water and nutrient retention. A pot experiment was conducted to examine the effects of organic acid conditio-ners (OASC) with four application rates (0, 20, 40, and 60 g·kg-1) on soil nutrients and crop growth. Maize and common vetch were the focus crops and loessial soil and irrigated desert soil were the soil types. The results showed that OASC application reduced water evaporation loss and significantly improved soil available nutrient content, with the most obvious effects on soil phosphorus. Available phosphorus content and proportion were increased by 256.5% and 227.4%, respectively, compared with no OASC treatment. The shoot dry weights of maize and common vetch on loessial soil were initially increased with the increasing OASC application rate and were highest at the application rate of 20 g·kg-1. The values progressively decreased with increasing OASC addition rates. Total salt content was significantly increased and the rate of emergence of common vetch decreased at OASC rate exceeding 20 g·kg-1. For irrigated desert soil, the OASC application rate did not affect total salt content when maize was planted. There was significant increase in soil total salt contents when common vetch was planted and at the OASC rate of 60 g·kg-1. The shoot dry weight of common vetch and maize was highest with the OASC application rates at 40 g·kg-1 and 60 g·kg-1, respectively. The optimal OASC rate for planting common vetch and maize on loess soil was 20 g·kg-1. The application rates of 40 g·kg-1 and 60 g·kg-1 were optimal for planting common vetch and maize on irrigated desert, respectively.
    Application effect of high efficiency and stability urea added with biochemical inhibitors and humic acid in loess
    XIAO Fu-rong, LI Dong-po, WU Zhi-jie, XUE Yan, CUI Lei, ZHANG Ke, LI Yong-hua, ZHENG Ye
    2021, 32(12):  4419-4428.  doi:10.13287/j.1001-9332.202112.010
    Asbtract ( 169 )   PDF (1158KB) ( 23 )  
    We carried out pot experiment to investigate nitrogen transformation characteristics, yield increasing effect, and apparent utilization rate of nitrogen fertilizer in loess soils by combining chemi-cal inhibitor and biostimulant humic acid to investigate the application effect and provide a theoretical basis for new type highly efficient and stable urea in loess soil. In this study, blank (CK) and urea (N) were set as controls, and humic acid alone (F), N-butyl thiophosphate-triamine (NBPT), 3,4-dimethyl-pyrazolate phosphate (DMPP), 2-chloro-6-trimethyl-pyridine (CP) and humic acid were respectively combined with three biochemical inhibitors to urea. The results showed that compared with N treatment, F, NBPT+F, DMPP+F, CP+F treatments significantly increased maize yield, chlorophyll content, leaf area index and nitrogen uptake, and had obvious effects on soil available nitrogen contents. The addition of humic acid increased chlorophyll content of maize leaves in all cases compared to the application of biochemical inhibitors alone. Compared with CP treatment, CP+F treatment could significantly increase nitrogen uptake, chlorophyll content, and nitrogen adsorption efficiency of maize. Addition of humic acid with NBPT increased nitrification inhibition rate by 10.7% compared with NBPT alone, but decreased yield, leaf area index, nitrogen uptake, nitrogen use efficiency. Compared with DMPP treatment, DMPP+F treatment significantly reduced maize yield, leaf area index, nitrogen uptake, nitrogen use efficiency and nitrification inhibition rate. Considering maize yield, plant N uptake, N fertilizer uptake rate and soil ammonium N and nitrate N contents, the addition of humic acid and CP is recommended for urea application in loess areas to enhance urea performance, yield, and fertilizer utilization.
    Relationships between phosphate solubilizing ability, pH, and organic acids of Penicillium oxalicum YTY and its mutants
    YANG Tian-you, LI Lin-bo, TIAN Jing, ZHANG Shi-shuang, WU Zhong-qi
    2021, 32(12):  4429-4438.  doi:10.13287/j.1001-9332.202112.038
    Asbtract ( 183 )   PDF (3358KB) ( 18 )  
    The high-efficiency phosphate solubilizing mutants of Penicillium oxalicum YTY were screened by mutagenesis of ion beam combined with UV. We analyzed the changes and correlation of phosphate solubilizing ability, pH, and organic acid for YTY and its mutants, and examined the phosphate solubilizing mechanism of P. oxalicum YTY. The results showed that five high-efficiency mutants, P9-8, P9-9, P15-4, P15-6, and P15-7 were screened, and that the phosphate solubili-zing ability of mutants was increased by more than 60% compared with YTY. In the process of pho-sphorus solubilization, both phosphorus solubilizing ability and rate of mutants were higher than that of YTY, and the mutants pH was significantly lower than YTY. The type and content of organic acids secreted by the mutants showed some variations. All mutants and YTY could secrete lactic acid, acetic acid and oxalic acid, while P9-8 also produced citric acid. The pH and the phosphate solubilizing ability of YTY and its mutants had a significant negative correlation. Phosphate solubilizing ability with organic acid and pH were all significantly correlated for YTY and the mutants, except P15-4. Organic acids and low environmental pH reduced by organic acids were the probable mechanism for P. oxalicum to dissolve phosphorus. Radiation of ion beam combined with UV could change the type and content of organic acids of P. oxalicum YTY, and initiate other H+ releasing pathways to lower pH, and participate phosphorus dissolution. The study provided biological mate-rials and theoretical basis for the research and development of high-efficiency phosphate solubilizing P. oxalicum and understanding the phosphate solubilizing mechanism of P. oxalicum.
    Isolation, identification, and degradation characteristics of three effective PAHs degradation strains
    YAN Shuang-dui, LIU Li-jun, CAO Yan-zhuan, YAN Qiu-yan, DONG Xin-yu, HONG Jian-ping
    2021, 32(12):  4439-4446.  doi:10.13287/j.1001-9332.202112.028
    Asbtract ( 209 )   PDF (974KB) ( 108 )  
    Isolating dominant strains for the degradation of polycyclic aromatic hydrocarbons (PAHs) is of great practical significance for the restoration of ecosystem polluted by PAHs. A total of 11 strains with capacity of degrading PAHs were obtained from soil polluted by PAHs around a coking plant, by enrichment culture, acclimation, and plate isolation. Three of them with effective PAH-degrading capability were identified and screened out by morphological observation, physiobiochemical characterization, and 16S rRNA gene sequencing, and respectively, named as DJ-3, DJ-8 and DJ-10. Based on the results of 16S rRNA gene sequencing, DJ-3, DJ-8, and DJ-10 were identified as Pseudomonas sp. Klebsiella sp., and Bacillus sp. The degradation rate of phenanthrene (200 mg·L-1), pyrene (200 mg·L-1), and naphthol (160 mg·L-1) by three strains (DJ-3, DJ-8 and DJ-10) after seven-day incubation were 48.9%-65.9%, 38.9%-43.1%, and 57.6%-64.9%, respectively. The degradation rates of mixed PAHs sample (1200 mg·L-1) by three strains were 49.1%, 44.5%, and 53.9%, which were significantly higher than other eight strains, indicating that they were highly effective in PAHs degradation. There was no antagonistic relationship among the three strains. This study would lay a foundation for building efficient PAHs degrading strains and improve the in situ bioremediation of PAHs contaminated soil.
    Effects of heavy metal cadmium on Caulerpa lentillifera based on transcriptome analysis
    PANG Mei-xia, HUANG Zhi-li, LIU Xiao-long, TANG Yong-jun, LI Xiao-dong, JIN Gang
    2021, 32(12):  4447-4456.  doi:10.13287/j.1001-9332.202112.035
    Asbtract ( 185 )   PDF (1727KB) ( 29 )  
    With the acceleration of industrialization, the toxic effect of heavy metal cadmium (Cd) pollution has become prominent. In order to explore the molecular mechanism of the physiological regulation of Caulerpa lentillifera under Cd stress, we analyzed the transcriptome of Cd-stressed (Hcd2+) algae tissues using RNA-Seq. A total of 702 differentially expressed genes (DEGs) were screened between the control and Hcd2+ groups, out of which 257 genes were up-regulated and 445 genes were down-regulated in the Hcd2+ group. We conducted functional annotation and enrichment analysis of the obtained DEGs. The results showed that various biological functions of C. lentillifera were affected under Cd2+stress, which eventually showed growth inhibition. Results of GO enrichment analysis showed that the production and removal of reactive oxygen species (ROS) in C. lentillifera were out of balance and caused oxidative damage such as DNA damage. Results of KEGG enrichment analysis showed that many photosynthesis-related pathways were inhibited, indicating that Cd2+ stress led to disorder of photosynthetic reaction of C. lentillifera.
    Evolution mechanism of ecosystem service value at the township-scale in Anxi County of Fujian Province, China
    JIANG Wen-jia, LIN Man-hong, WENG Pei-ying, WEI Dao-zhi, LIN Wen-xiong, SU Kai
    2021, 32(12):  4457-4466.  doi:10.13287/j.1001-9332.202110.011
    Asbtract ( 221 )   PDF (1180KB) ( 36 )  
    To clarify the interaction between land-use change and ecosystem service, the traditional ecosystem service valuation model was modified using the adjustment coefficients of biomass factor and socio-economic factor to evaluate ecosystem service value (ESV) of 24 towns in Anxi County from 1999 to 2019, aiming to understand the ecological response to land-use change. The results showed that ESV of 24 towns in Anxi County decreased gradually during the study period, with a decrease of 0.004 to 0.295 million yuan. Forestland and cultivated land made the largest contribution to ESV, followed by grassland and water bodies, whereas the contribution of other land types were negative. ESV could be increased by converting other land types into forestland, but be decreased by converting other land types into orchard land area (mainly tea plantation). The chemical fertilizer-monoculture management pattern of ratooning tea plantation reduced ESV, while the ecological cultivation of tea plantation could effectively improve soil microbial diversity and soil fertility, and finally enhance ESV.
    Decomposition dynamics and driving factors of manure in reclaimed soils in coal mining area
    LI Ran, DUAN Ying-hua, SUN Nan, WU Lei, ZHANG Qiang, JIN Dong-sheng, GAO Chun-hua, HONG Jian-ping, XU Ming-gang
    2021, 32(12):  4467-4474.  doi:10.13287/j.1001-9332.202112.013
    Asbtract ( 178 )   PDF (1865KB) ( 49 )  
    Understanding the decomposition dynamics and driving factors of manure in the soil subjected to different reclaimed years could provide theoretical basis to rational utilization of manure and soil fertility improvement in coal mining area. Cattle manure and pig manure were mixed with soils subjected to different reclaimed years (one year, R1; 10 years, R10; and 30 years, R30) at the ratio of manure carbon to soil mass of 4 to 100, so as to examine manure decomposition characteristics using the nylon mesh bag (15 cm deep of soil buried) in the Shanxi coal mine reclamation area, with no manure addition as control (CK). Soil samples were collected at day 12, 23, 55, 218, 281, and 365 to measure the contents of soil manure residual, soil microbial biomass carbon (MBC), and dissolved organic carbon (DOC). The contributions of soil properties, manure properties, and hydrothermal condition to manure decomposition were quantified. The results showed that the decomposition rates of pig manure were significantly higher than cattle manure. The humification coefficient of pig manure (average 46.3%) was lower than that of cattle manure (average 71.7%). The humification coefficient of pig manure was significantly lower in the 30-year reclaimed soil (44.5%) compared to the 1-year and 10-year reclaimed soil (average 47.2%). There was no significant difference in the humification coefficient of cattle manure among the three reclaimed soils. The proportion and decomposition rate constant of labile carbon pool of pig manure and cattle manure were significantly different, with values of 52% and 26%, and 0.00085 and 0.00074 ℃-1, respectively. The positive effect of pig manure on MBC and DOC in reclaimed soil was significantly higher than that of cattle manure over 0-218 days, but no difference over 281-365 days. The magnitude of the enhancement of MBC and DOC in those three reclaimed soils after manure amendments showed a similar trend of R1 >R10 ≈ R30. Results of variance partitioning analysis showed that manure decomposition was mainly controlled by manure properties (17.9%) when considering soil properties, manure properties, and hydrothermal condition. In conclusion, the decomposition of pig manure but not cattle manure was regulated by reclamation year. Cattle manure, with higher humification coefficient than pig manure, was recommended for reclaimed mining area to improve soil fertility.
    Nitrogen flow and polluted nitrogen footprint accounting of rural systems at the municipal scale: The case of Shenyang
    ZHAI Jia-ning, LI Hong-qing, REN Wan-xia, TANG Cheng-rui, XUE Bing
    2021, 32(12):  4475-4487.  doi:10.13287/j.1001-9332.202112.019
    Asbtract ( 153 )   PDF (2221KB) ( 17 )  
    Understanding the process of nitrogen flow, emission, and utilization in rural systems is of great importance to reduce pollutant emissions caused by agricultural activities and to promote the sustainable utilization of nutrient resources. Based on the NUFER (nutrient flow in food chain, environment and resources use) model and nitrogen footprint approach, we proposed a nitrogen accounting framework for municipal-scale rural systems, and conducted a quantitative and time series-based comparative study on nitrogen flow, recycling, and footprint from the perspective of three subsystems, namely planting, animal husbandry, and rural human settlement in Shenyang from 1998 to 2018. The results showed that: 1) nitrogen utilization rates of the planting subsystem, animal husbandry subsystem, and rural human settlement subsystem were 36.1%, 59.7%, and 70.1%, respectively in 2018, with a growth rate of 15.9%, 9.1%, and 0.7% respectively compared to 1998; 2) The total polluted nitrogen footprint in Shenyang grew rapidly from 1998 to 2014, but declined from 2014 to 2018. In 2018, the total footprint of nitrogen pollution was 123.5 Gg, increasing by 21.6% compared with 1998; 3) In 2018, the unused nitrogen in the planting subsystem, animal husbandry subsystem, and rural human settlement subsystem were 129.5, 62.2, and 8.7 Gg, respectively, which were equivalent to 420.4, 202.1, and 8.7 Gg of nutrient resources from nitrogen fertilizer, respectively. In general, the nitrogen use efficiency of rural system increased gradually from the production end to the consumption end, but the temporal and spatial variation of nitrogen in rural systems need further studies.
    Evaluation of eutrophication level changes in Baiyangdian Lake based on multiple biological groups
    LI Dai-kui, HE Ping, LIU Cun-qi, WANG Jun-xia, ZHANG Ya-juan, REN Ying, WANG De-wang
    2021, 32(12):  4488-4498.  doi:10.13287/j.1001-9332.202112.027
    Asbtract ( 202 )   PDF (3345KB) ( 38 )  
    Baiyangdian Lake is the largest shallow lake in the North China Plain. Due to the increases in upstream water consumption in recent decades, the amount of natural water entering the lake has decreased, resulting in a significant drop in the water level of Baiyangdian Lake. Severe eutrophication has occurred in Baiyangdian Lake owing to the discharge of domestic sewage around the sub-lakes. With the establishment of the Xiong’an New Area in 2017, the importance of environmental governance and protection of Baiyangdian Lake has been unprecedentedly increased. The implementation of ecological water supplement from upstream reservoirs and other basins has significantly increased water level of Baiyangdian Lake. Moreover, domestic sewage from surrounding rural areas was collected. In order to understand the current state of water quality in Baiyangdian Lake, and to evaluate the effects of the previous water replenishment and pollution control projects, we investigated water physicochemical variables and aquatic organisms of Baiyangdian Lake in August 2019. We evaluated water quality status of Baiyangdian Lake based on water quality, phytoplankton, zooplankton and macrozoobenthos, by comparing with the evaluation based on the survey data in 2010. In addition, submerged macrophytes were used to evaluate the water ecological status. The results showed that the water quality state of Baiyangdian Lake in 2019 was greatly improved compared with 2010. Among them, total phosphorus concentration decreased by 88.6%, total nitrogen concentration decreased by 83.9%, chlorophyll a concentration decreased by 47.8%, and the Seechi depth increased by 43.4%. The diversity of zooplankton and macrozoobenthos were significantly increased, the density of phytoplankton was significantly reduced, the dominant species of phytoplankton changed from a general highly eutrophic type to a eutrophic type, and the distribution of the clean-type submerged macrophytes was narrowed. The evaluation methods based on submerged macrophytes species and based on zooplankton diversity were not suitable for water quality evaluation in Baiyangdian Lake. In the early stage, the pollution source of Baiyangdian Lake was concentrated in the west, which led to lower water quality. Nowadays, due to the water replenishment from the west, north and south, the hydrodynamic conditions changed significantly. In the west, pollution control projects were carried out, and the spatial distribution of water quality in Baiyangdian Lake are showing a characteristics of homogenization.
    Impacts of Spartina alterniflora invasion on the benthic food web in the Yellow River Delta during autumn
    JIANG Shao-yu, CHEN Lin-lin, YAN Lang, LIU Chun-yun, PENG Zi-rui, ZHANG Chuan-xin, LI Bao-quan
    2021, 32(12):  4499-4507.  doi:10.13287/j.1001-9332.202112.009
    Asbtract ( 280 )   PDF (1291KB) ( 52 )  
    Spartina alterniflora was introduced into the Yellow River Delta (YRD) in 1990 with the purpose of shore protection and siltation accretion. However, it spread rapidly and became a severe threat to the local coastal wetland ecosystem. To assess the impacts of S. alterniflora invasion on the benthic food web, we sampled the potential food sources of macrobenthos in November 2020, analyzed the trophic level and the benthic food web structure based on stable isotope technique. Results showed that the average δ13C values of macrobenthic food sources followed an order: sediment organic matter (SOM) > S. alterniflora > benthic microalgae > particulate organic matter (POM) > Suaeda salsa. The average δ15N values significantly differed among food sources, ranging from 1.24‰ to 9.03‰. The trophic levels of different macrobenthos ranged from 1.73 to 4.19, of which the bivalve species was the lowest one. S. alterniflora and the decayed debris were the most important food sources for macrobenthos, but without any impact on the trophic level structure of macro-benthos. In conclusion, Spartina alterniflora invasion distinctly changed the composition of food sources of macrobenthos through a “bottom-up” effect, which would probably impact the local food web structure in the YRD wetland.
    Patterns of maternal transfer of trace elements in the pelagic thresher shark (Alopias pelagicus) of the Eastern Pacific Ocean
    LI Ze-zheng, LI Yun-kai
    2021, 32(12):  4508-4514.  doi:10.13287/j.1001-9332.202112.039
    Asbtract ( 143 )   PDF (1185KB) ( 24 )  
    Pelagic sharks are vulnerable to overfishing due to slow growth rates, late-at-maturity and low fecundity, 90% of which are Near Threatened with an elevated risk of extinction according to IUCN Red List Criteria. Trace elements can be accumulated by marine predators and may have detrimental effects on population dynamics. In this study, we analyzed the concentrations of 11 trace elements (Zn, Cu, Cr, Ni, Mn, Se, Co, Hg, Cd, Pb, and As) in muscle and liver tissues of 10 pregnant pelagic thresher sharks (Alopias pelagicus) and their 18 embryos. The results showed that four essential elements (Cu, Cr, Mn, and Se) were accumulated in both tissue types of embryos. Ni and Zn concentrations were higher in embryonic muscle than that in the liver. For nonessential elements, concentrations of As, Cd and Hg in both embryonic tissues were lower than those of their mothers. Though maternal and embryonic tissues had high levels of Hg, the Se/Hg molar ratios in both tissues of the embryo were above 1 with larger values in the embryos, indicating that Se played a protective role against Hg toxicity in embryonic tissues. Liver is the primary energy resource of embryo development. There was no correlation for element concentrations between embryonic and maternal liver tissues, indicating there is a regulatory mechanism to maintain the stability of element contents during maternal transfer in pelagic thresher shark.
    Preliminary analysis on the composition of fish and cephalopod in diet of Dosidicus gigas in the high sea of Eastern Pacific Ocean
    CAO Liang-liang, LI Qing-ying, LIU Bi-lin
    2021, 32(12):  4515-4522.  doi:10.13287/j.1001-9332.202112.042
    Asbtract ( 203 )   PDF (1445KB) ( 22 )  
    Dosidicus gigas is widely distributed in the Eastern Pacific Ocean and plays an important role in the marine ecosystem. In this study, we identified and analyzed 4131 otoliths and 75 beaks from the stomach content residuals of 62 D. gigas samples, which were collected by Chinese squid jigging vessels in the high sea of Eastern Pacific Ocean from June to December 2019. The results showed that the preys of D. gigas included 10 fish species and 4 cephalopod species. Judging from the frequency of appearance and percentage of amount, Vinciguerria lucetia, Diogenichthys laternatus, and Triphoturus mexicanus were the dominant myctophidae prey of D. gigas. Among those preys, V. lucetia was the most important one, as it was found in all D. gigas with different mantle lengths. D. gigas, Sthenoteuthis oualaniensis, and Onychoteuthis banksii were the dominant cephalopod prey. With the growth of mantle length of D. gigas, the number of prey species, the percentage of cephalopod and myctophidae preys, and the size of preys increased, and thus the trophic level of preys increased as the trophic pattern of preys changed. Our results could provide basic information for evaluating the contribution of different preys in the preying transformation of D. gigas.
    Sampling design optimization based on estimation accuracy and survey cost for fishery-independent surveys
    ZHANG Guo-sheng, ZHANG Chong-liang, XUE Ying, JI Yu-peng, REN Yi-ping, XU Bin-duo
    2021, 32(12):  4523-4531.  doi:10.13287/j.1001-9332.202112.034
    Asbtract ( 194 )   PDF (2175KB) ( 35 )  
    Optimization of sampling design can reduce survey cost, ensure the accuracy of survey data, and get the maximum benefit of survey design. In this study, survey voyage was added into sampling survey design optimization as the evaluation index to quantify the survey cost. Computer simulation and resampling technique were used to simulate the survey plan. Simulation annealing algorithm was used to find the survey design with the shortest voyage. Based on the survey data collected from the bottom trawl survey conducted in the Haizhou Bay and its adjacent waters in spring and autumn from 2013 to 2015, Kriging interpolation was used to simulate the relative abundance distribution of Conger myriaster, Enedrias fangi, and Hexagrammos otakii in the bay in two seasons as their ‘true’ values. Resampling was conducted using stratified random sampling with different sample sizes for simulation study, and the abundance indices of each target species were estimated based on the simulated data. The relative bias (RB) was used to evaluate the accuracy of estimation of abundance index. The average survey cost and probability of budget overshoot (P) were used to measure the cost of survey. Integrated evaluation index (IEI) was developed to measure the survey cost and estimation accuracy comprehensively. The results showed RB values of all target species decreased with sample size. Different target species had different RB values with the same sample size. The optimal numbers of station determined by different evaluation indices were different. The optimal sample size determined by RB was relatively high, while the cost exceeded the budget of survey. The optimal sample size determined by probability of exceeding survey budget was relatively low but the accuracy and precision of estimation was low. IEI balanced the survey cost and estimation accuracy, and the optimal sample size defined by IEI was at intermediate level. The sampling design considering survey cost and route planning could not only determine the optimal sample size, but also record the potential station location, the survey sequence and the corresponding survey cost of simulations. According to the actual cost budget, the desired survey design with low cost and low RB of estimation from simulations could be chosen for fishery-independent surveys.
    Heterogeneity of fork length-weight relationship for juvenile Engraulis japonius based on linear mixed-effects models
    ZHU Wen-bin, ZHU Hai-chen, WANG Ya-li, ZHANG Ya-zhou, LU Zhan-hui, CUI Guo-chen
    2021, 32(12):  4532-4538.  doi:10.13287/j.1001-9332.202112.032
    Asbtract ( 171 )   PDF (1103KB) ( 23 )  
    In order to examine the growth heterogeneity of juvenile Engraulis japonicus, we conducted heterogeneity analysis of fork length-weight relationships of the species based on the samples of juvenile E. japonius from April to June 2019 collected from the coastal waters of Zhejiang Province by using the generalized linear model and 9 linear mixed-effect models. The results showed that the fork length of juvenile E. japonius was between 14-74 mm, with an average of 33 mm. The dominant fork length group was 21-50 mm. The weight range was 0.01-2.96 g, with an average of 0.28 g, and the dominant body weight group was 0.01-0.50 g. According to Akaike information criterion (AIC), the growth-fitting model with both months and regions random effects on the growth parameters a and b of juvenile E. japonius fitted best.The reliable prediction results was testified by the cross-validation. In the best model, the fixed value of growth parameter a was 0.24×10-5 and its estimated value did not fluctuate significantly, while the fixed value of b was 3.246 and the estimated value ranged from 3.206 to 3.272, indicating that juvenile E. japonius was under positive allometric growth. Our results suggested that month and region significantly affected the fork length-weight relationships of juvenile E. japonius.
    Reviews
    Convergent cross mapping method and its application in ecology
    WANG Dan-yu, ZHU Yuan-jun, YANG Xiao-hui
    2021, 32(12):  4539-4548.  doi:10.13287/j.1001-9332.202112.036
    Asbtract ( 867 )   PDF (1899KB) ( 326 )  
    The convergent cross mapping (CCM) is a method to analyze causality of nonlinear time series variables. Different from the traditional linear system analysis method, CCM gets historical information based on their state space reconstruction. The presence of causality can be confirmed when the estimated values perform convergent with time series extension. Here, we introduced the develop-ment history of CCM and its advantages over the traditional Granger causality test, and elaborated the principle, algorithm process, and implementation approach. As a system analysis method aiming at the coupling relationship between variables from weak to moderate, CCM can effectively solve the complex causality among nonlinear multivariable in ecosystems. When it is applied to the causality analysis of multi-point time series variables with spatial information, the spatial autocorrelation among points should be fully considered and combined with the method that can remove the spatial correlation between variables and sequences, so as to ensure more accurate causality analysis using CCM and more convincing results.