植物生态学报 ›› 2010, Vol. 34 ›› Issue (7): 792-799.DOI: 10.3773/j.issn.1005-264x.2010.07.004 cstr: 32100.14.j.issn.1005-264x.2010.07.004
收稿日期:2009-07-18
接受日期:2010-01-17
出版日期:2010-07-18
发布日期:2010-07-01
作者简介:* E-mail: gszhou@ibcas.ac.cn
JI Yu-He1,2, ZHOU Guang-Sheng1,*(
)
Received:2009-07-18
Accepted:2010-01-17
Online:2010-07-18
Published:2010-07-01
摘要:
基于辽河三角洲14年的水稻(Oryza sativa)与玉米(Zea mays)产量、42年的芦苇(Phragmites australis)产量及相应的气候、农业人口、灌溉和施肥等资料, 主要利用主成分分析、简单相关分析和偏相关分析, 阐述了辽河三角洲3种主要植被生产功能的年际变化及其驱动因子。结果表明: 玉米和水稻的单产呈弱上升趋势, 但年际间波动较大; 芦苇单产则呈显著上升趋势, 但年际间波动不大。气候因子是玉米和水稻产量年际间变化的主要驱动因子, 而人为因子(农业人口、灌溉和施肥)的作用不明显; 气候因子和人为因子(灌溉)均明显影响芦苇产量的年际变化。但驱动辽河三角洲3种主要植被(玉米、水稻和芦苇)生产功能年际波动的主导气候因子并不相同, 分别是年日照数、年降水量和年蒸发量。
汲玉河, 周广胜. 辽河三角洲植被生产功能的年际变化及其驱动 因子分析. 植物生态学报, 2010, 34(7): 792-799. DOI: 10.3773/j.issn.1005-264x.2010.07.004
JI Yu-He, ZHOU Guang-Sheng. Inter-annual change of vegetation productivity and its driving factors in the Liaohe Delta, northeastern China. Chinese Journal of Plant Ecology, 2010, 34(7): 792-799. DOI: 10.3773/j.issn.1005-264x.2010.07.004
| 地区 Region | 盘锦 Panjin | 锦州 Jinzhou | 营口 Yingkou | 鞍山 Anshan | 沈阳 Shenyang | 辽阳 Liaoyang |
|---|---|---|---|---|---|---|
| 百分比 Percentage | 100% | 35.4% | 27.7% | 15.4% | 0.5% | 0.2% |
表1 各地区属于辽河三角洲的面积与该地区的面积比
Table 1 The area percentage of each region falling into the Liaohe Delta
| 地区 Region | 盘锦 Panjin | 锦州 Jinzhou | 营口 Yingkou | 鞍山 Anshan | 沈阳 Shenyang | 辽阳 Liaoyang |
|---|---|---|---|---|---|---|
| 百分比 Percentage | 100% | 35.4% | 27.7% | 15.4% | 0.5% | 0.2% |
| 主成分 Principal component | 特征值 Eigenvalue | 贡献率 Contribution (%) | 累计贡献率 Accumulative contribution (%) |
|---|---|---|---|
| 1 | 1.796 | 59.844 | 59.844 |
| 2 | 1.047 | 34.894 | 94.739 |
| 3 | 0.158 | 5.261 | 100 |
表2 从影响玉米产量的人为因子中提取的主成分特征值和贡献率
Table 2 The eigenvalue and contribution rate of principal component extracting from anthropic factors which effect the yield of maize
| 主成分 Principal component | 特征值 Eigenvalue | 贡献率 Contribution (%) | 累计贡献率 Accumulative contribution (%) |
|---|---|---|---|
| 1 | 1.796 | 59.844 | 59.844 |
| 2 | 1.047 | 34.894 | 94.739 |
| 3 | 0.158 | 5.261 | 100 |
| 主成分载荷 Principal component loading | 特征向量 Eigenvector | |||
|---|---|---|---|---|
| L1 | L2 | T1 | T2 | |
| 农业人口 Agricultural population | -0.961 | -0.005 | -0.717 | -0.005 |
| 玉米灌溉面积 Irrigation area of maize | 0.630 | 0.752 | 0.470 | 0.735 |
| 玉米施肥量 Fertilizer for maize | 0.689 | -0.694 | 0.514 | -0.678 |
表3 从影响玉米产量的人为因子中提取的主成分载荷和主成分特征向量
Table 3 The loading and eigenvector of principal component extracting from anthropic factors which effect the yield of maize
| 主成分载荷 Principal component loading | 特征向量 Eigenvector | |||
|---|---|---|---|---|
| L1 | L2 | T1 | T2 | |
| 农业人口 Agricultural population | -0.961 | -0.005 | -0.717 | -0.005 |
| 玉米灌溉面积 Irrigation area of maize | 0.630 | 0.752 | 0.470 | 0.735 |
| 玉米施肥量 Fertilizer for maize | 0.689 | -0.694 | 0.514 | -0.678 |
图1 1992-2005年间辽河三角洲玉米粮食单产(kg·hm-2)的年际变化(辽宁省统计局, 1993-2006)。
Fig. 1 The inter-annual change of maize yield per unit (kg·hm-2) in the Liaohe Delta during 1992-2005 (Liaoning Provincial Bureau of Statistics, 1993-2006).
图2 1992-2005年间辽河三角洲水稻粮食单产(kg·hm-2)的年际变化(辽宁省统计局, 1993-2006)。
Fig. 2 The inter-annual change of rice yield per unit (kg·hm-2) in the Liaohe Delta during 1993-2005 (Liaoning Provincial Bureau of Statistics, 1993-2006).
图3 1949-1990年辽河三角洲5个农场的芦苇单产的年际变化(kg·hm-2) (盘锦市人民政府地方志办公室, 1998)。
Fig. 3 The inter-annual change of reed yield per unit in the five reed farms in the Liaohe Delta during 1949-1990 (kg·hm-2) (Panjin Municipal Government Office of local Chorography, 1998).
| 皮尔逊相关 Pearson correlation | p | N | |
|---|---|---|---|
| 人为综合因子得分 Synthesized score of human factors | -0.019 | 0.947 | 14 |
| 气候因子综合得分 Synthesized score of climatic factors | -0.279 | 0.334 | 14 |
表4 1992-2005年玉米单产与影响玉米产量的人为因素综合得分、气候因子综合得分的皮尔逊相关关系
Table 4 The Pearson correlation between maize yield per unit and the synthesized score of climatic factors, and the synthesized score of human factors effected maize during 1992-2005
| 皮尔逊相关 Pearson correlation | p | N | |
|---|---|---|---|
| 人为综合因子得分 Synthesized score of human factors | -0.019 | 0.947 | 14 |
| 气候因子综合得分 Synthesized score of climatic factors | -0.279 | 0.334 | 14 |
| 皮尔逊相关 Pearson correlation | p | N | |
|---|---|---|---|
| 人为因子综合得分 Synthesized score of human factors | -0.046 | 0.876 | 14 |
| 气候因子综合得分 Synthesized score of climatic factors | -0.537* | 0.048 | 14 |
表5 1992-2005年水稻单产与影响水稻产量的人为因素综合得分、气候因子综合得分的皮尔逊相关关系
Table 5 The Pearson correlation between rice yield per unit and the synthesized score of climatic factors, and the synthesized score of human factors effected rice during 1992-2005
| 皮尔逊相关 Pearson correlation | p | N | |
|---|---|---|---|
| 人为因子综合得分 Synthesized score of human factors | -0.046 | 0.876 | 14 |
| 气候因子综合得分 Synthesized score of climatic factors | -0.537* | 0.048 | 14 |
| 相关系数 r | 显著度 p | 自由度 df | 控制变量 Control variable | |
|---|---|---|---|---|
| 年日照数 Annual hours of sunshine (Ann_shine) | -0.555 | 0.096 | 8 | Ann_preci, Frost free, Ann_evap, Ann_temp |
| 年降水量 Annual precipitation (Ann_preci) | -0.431 | 0.213 | 8 | Ann_shine, Frost free, Ann_evap, Ann_temp |
| 无霜期 Frost free | -0.075 | 0.837 | 8 | Ann_shine, Ann_preci, Ann_evap, Ann_temp |
| 年蒸发量 Annual evaporation (Ann_evap) | 0.024 | 0.948 | 8 | Ann_shine, Ann_preci, Frost free, Ann_temp |
| 年平均温度 Mean annual temperature (Ann_temp) | 0.109 | 0.764 | 8 | Ann_shine, Ann_preci, Frost free, Ann_evap |
表6 1992-2005年玉米单产与气候因子的偏相关关系
Table 6 The partial correlation between climatic factors and maize yield per unit during 1992-2005
| 相关系数 r | 显著度 p | 自由度 df | 控制变量 Control variable | |
|---|---|---|---|---|
| 年日照数 Annual hours of sunshine (Ann_shine) | -0.555 | 0.096 | 8 | Ann_preci, Frost free, Ann_evap, Ann_temp |
| 年降水量 Annual precipitation (Ann_preci) | -0.431 | 0.213 | 8 | Ann_shine, Frost free, Ann_evap, Ann_temp |
| 无霜期 Frost free | -0.075 | 0.837 | 8 | Ann_shine, Ann_preci, Ann_evap, Ann_temp |
| 年蒸发量 Annual evaporation (Ann_evap) | 0.024 | 0.948 | 8 | Ann_shine, Ann_preci, Frost free, Ann_temp |
| 年平均温度 Mean annual temperature (Ann_temp) | 0.109 | 0.764 | 8 | Ann_shine, Ann_preci, Frost free, Ann_evap |
| 相关系数 r | 显著度 p | 自由度 df | 控制变量 Control variable | |
|---|---|---|---|---|
| 年降水量 Annual precipitation (Ann_preci) | -0.666 | 0.035 | 8 | Ann_shine, Frost free, Ann_evap, Ann_temp |
| 年日照数 Annual hours of sunshine (Ann_shine) | -0.430 | 0.215 | 8 | Ann_preci, Frost free, Ann_evap, Ann_temp |
| 年平均温度 Mean annual temperature (Ann_temp) | 0.058 | 0.874 | 8 | Ann_shine, Ann_preci, Frost free, Ann_evap |
| 年蒸发量 Annual evaporation (Ann_evap) | 0.077 | 0.832 | 8 | Ann_shine, Ann_preci, Frost free, Ann_temp |
| 无霜期 Frost free | 0.211 | 0.559 | 8 | Ann_shine, Ann_preci, Ann_evap, Ann_temp |
表7 1992-2005年水稻单产与气候因子的偏相关关系
Table 7 The partial correlation between climatic factors and rice yield per unit during 1992-2005
| 相关系数 r | 显著度 p | 自由度 df | 控制变量 Control variable | |
|---|---|---|---|---|
| 年降水量 Annual precipitation (Ann_preci) | -0.666 | 0.035 | 8 | Ann_shine, Frost free, Ann_evap, Ann_temp |
| 年日照数 Annual hours of sunshine (Ann_shine) | -0.430 | 0.215 | 8 | Ann_preci, Frost free, Ann_evap, Ann_temp |
| 年平均温度 Mean annual temperature (Ann_temp) | 0.058 | 0.874 | 8 | Ann_shine, Ann_preci, Frost free, Ann_evap |
| 年蒸发量 Annual evaporation (Ann_evap) | 0.077 | 0.832 | 8 | Ann_shine, Ann_preci, Frost free, Ann_temp |
| 无霜期 Frost free | 0.211 | 0.559 | 8 | Ann_shine, Ann_preci, Ann_evap, Ann_temp |
图4 1949-1990年辽河三角洲5个芦苇农场的芦苇总产量和灌溉面积、苇田面积(盘锦市人民政府地方志办公室, 1998)。
Fig. 4 The irrigation area, reed area and total yield of the five reed farms in the Liaohe Delta during 1949-1990 (Panjin Municipal Government Office of local Chorography, 1998).
| 年平均温度 Ann_temp (℃) | 年降水量 Ann_preci (mm) | 年蒸发量 Ann_evap (mm) | 年日照数 Ann_shine (h) | 无霜期 Frost free (d) | |
|---|---|---|---|---|---|
| 皮尔逊相关(r) Pearson correlation | 0.348 | 0.79 | -0.979 | -0.775 | 0.364 |
| p | 0.5 | 0.062 | 0.001 | 0.07 | 0.478 |
| N | 6 | 6 | 6 | 6 | 6 |
表8 1958-1990年5个芦苇农场的芦苇单产与大洼县气候因子的相关关系
Table 8 The correlation between climatic factors in the Dawa county and reed yield per unit at the five reed farms during 1958-1990
| 年平均温度 Ann_temp (℃) | 年降水量 Ann_preci (mm) | 年蒸发量 Ann_evap (mm) | 年日照数 Ann_shine (h) | 无霜期 Frost free (d) | |
|---|---|---|---|---|---|
| 皮尔逊相关(r) Pearson correlation | 0.348 | 0.79 | -0.979 | -0.775 | 0.364 |
| p | 0.5 | 0.062 | 0.001 | 0.07 | 0.478 |
| N | 6 | 6 | 6 | 6 | 6 |
| [1] | Alexandrov VA, Hoogenboom G (2001). Climate variation and crop production in Georgia, USA, during the twentieth century. Climate Research, 17, 33-34. |
| [2] | Chen ZM, Babiker IS, Chen ZX, Komaki KM, Mohamed AA, Kato K (2004). Estimation of inter annual variation in productivity of global vegetation using NDVI data. International Journal of Remote Sensing, 25, 3139-3159. |
| [3] | Dilley M (1997). Climatic factors affecting annual maize yields in the valley of Oaxaca, Mexico. International Journal of Climatology, 17, 1549-1557. |
| [4] | Donohue RJ, Mcvicar T, Roderick ML (2009). Climate related trends in Australian vegetation cover as inferred from satellite observations, 1981-2006. Global Change Biology, 15, 1025-1039. |
| [5] | Ewe SML, Gaiser EE, Childers DL, Iwaniec D, Rivera- Monroy VH, Twilley RR (2006). Spatial and temporal patterns of aboveground net primary productivity (ANPP) along two freshwater estuarine transects in the Florida coastal everglades. Hydrobiologia, 569, 459-474. |
| [6] | Hao YP (郝永萍), Chen YF (陈育峰), Zhang XY (张兴有) (1998). Progress in estimation of net primary productivity and its responses to climate change. Advance in Earth Sciences (地球科学进展), 13, 564-571. (in Chinese with English abstract) |
| [7] | Huenneke LF, Anderson JP, Remmenga M, Schlesinger WH (2002). Desertification alters patterns of aboveground net primary production in Chihuahuan ecosystems. Global Change Biology, 8, 247-264. |
| [8] | Kahabka JE, Vanes HM, Clenahan EJ, Cox WJ (2004). Spatial analysis of maize response to nitrogen fertilizer in central New York. Precision Agriculture, 5, 463-476. |
| [9] | Kar G, Verma HN, Singh R (2006). Effects of winter crop and supplemental irrigation on crop yield, water use efficiency and profitability in rain-feed rice based cropping system of eastern India. Agricultural Water Management, 79, 280-292. |
| [10] | Kumar KK, Kumar KR, Ashrit RG, Deshpande NR, Hansen JW (2004). Climate impacts on Indian agriculture. International Journal of Climatology, 24, 1375-1393. |
| [11] | Li JL (李加林), Zhao HB (赵寒冰), Cao YG (曹云刚), Liu C (刘闯) (2006). Analysis on change of landscape pattern of wetland in Liaohe Delta. Urban Environment Urban Ecology (城市环境与城市生态), 19, 5-7. (in Chinese with English abstract) |
| [12] | Li XW (李晓文), Xiao DN (肖笃宁), Hu YM (胡远满) (2001). The landscape planning scenarios designing and the measures identification in the Liaohe River Delta wetland. Acta Ecologica Sinica (生态学报), 21, 353-364. (in Chinese with English abstract) |
| [13] | Liaoning Provincial Bureau of Statistics (辽宁省统计局) (1993-2006). Liaoning Statistical Yearbook (辽宁统计年鉴). China Statistics Press, Beijing. (in Chinese) |
| [14] | Lobell DB, Field CB (2007). Global scale climate crop yield relationships and the impacts of recent warming. Environmental Research Letters, 2, 1-7. |
| [15] | Meng JH (蒙继华) (2006). Research to Crop Growth Monitoring Indicators with Remote Sensing (农作物长势遥感监测指标研究). PhD dissertation, Institute of Remote Sensing Applications, Chinese Academy of Sciences, Beijing. (in Chinese with English abstract) |
| [16] | Nagy J, Huzsvai L (2005). Effect of weather on maize yields and the efficiency of fertilization. Acta Agronomica Hungarica, 53, 31-39. |
| [17] | Panjin Municipal Government Office of Local Chorography (盘锦市人民政府地方志办公室) (1998). Panjin Municipal Chorography (盘锦市志). Chorography Press, Panjin, China. (in Chinese) |
| [18] |
Peng SB, Huang JL, Sheehy JE, Laza RC, Visperas RM, Zhong XH, Centeno GS, Khush GS, Cassman KG (2004). Rice yield decline with higher night temperature from global warming. Proceedings of the National Academy of Sciences of the United States of America, 101, 9971-9975.
DOI URL PMID |
| [19] | Pennington DD, Collins SL (2007). Response of an arid-land ecosystem to inter-annual climate variability and prolonged drought. Landscape Ecology, 22, 897-910. |
| [20] | Piovesan G, Biondi F, Filippo A, Alessandrini A, Maugeri M (2008). Drought driven growth reduction in old beech (Fagus sylvatica L.) forests of the central Apennines, Italy. Global Change Biology, 14, 1-17. |
| [21] | Potgieter AB, Hammer GL, Butler D (2002). Spatial and temporal patterns in Australian wheat yield and their relationship with ENSO. Australia Journal of Agricultural Research, 53, 77-89. |
| [22] |
Rathore AL, Pal AR, Sahu RK, Chaudhary JL (1996). On farm rainwater and crop management for improving productivity of rain fed areas. Agricultural Water Management, 31, 253-267.
DOI URL |
| [23] | Sang WG, Su HX (2009). Interannual NPP variation and trend of Picea schrenkiana in the Tianshan Mountain, Xinjiang, China. Ecology Research, 24, 441-452. |
| [24] |
Stige LC, Stave J, Chan KS, Ciannelli L, Pettorelli N, Glantz M, Herren HR, Stenseth NC (2006). The effect of climate variation on agro-pastoral production in Africa. Proceedings of the National Academy of Sciences of the United States of America, 103, 3049-3053.
DOI URL PMID |
| [25] | Veron SR, Paruelo JM, Sala OE, Lauenroth WK (2002). Environmental controls of primary in agricultural systems of the Argentine Pampas. Ecosystems, 5, 625-635. |
| [26] | Walker NJ, Schulze RE (2008). Climate change impacts on agro-ecosystem sustainability across three climate regions in the maize belt of South Africa. Agriculture, Ecosystem and Environment, 124, 114-124. |
| [27] | Wang EL, Yu Q, Wu DR, Xia J (2008). Climate, agricultural production and hydrological balance in the North China Plain. International Journal of Climatology, 28, 1959-1970. |
| [28] | Wang L (王凌), Li XZ (李秀珍), Guo DF (郭笃发), Wang XG (王绪高), Shu Y (舒莹), Hu YM (胡远满) (2004). Scenarios study on returning farmland to wetlands for protecting red-crowned cranes in the Liaohe Delta. Acta Ecologica Sinica (生态学报), 24, 1770-1719. (in Chinese with English abstract) |
| [29] | Wang XL (王宪礼), Hu YM (胡远满), Bu RC (布仁仓) (1996). Analysis of wetland landscape changes in Liaohe Delta. Scientia Geographica Sinica (地理科学), 16, 260-265. (in Chinese with English abstract) |
| [30] | Yan M (晏明), Liu ZM (刘志明), Yan XY (晏晓英) (2005). Estimating of main crop yields in Jinlin Province using meteorological satellite data. Meteorological Science and Technology (气象科技), 23, 350-354. (in Chinese with English abstract) |
| [31] | Yan Y (阎雨), Chen SB (陈圣波), Tian J (田静), Tan XP (谭笑平) (2004). Development and prospects of production estimation by satellite remote sensing. Journal of Jinlin Agricultural University (吉林农业大学学报), 26, 187-191. (in Chinese with English abstract) |
| [32] | Zhou GS (周广胜), Zhou L (周莉), Guan EK (关恩凯) (2006). Liaohe Delta wetland and global change. Journal of Meteorology and Environment (气象与环境学报), 22, 7-12. (in Chinese with English abstract) |
| [1] | 吴光进, 郭垚鑫, 任成杰, 王俊, 岳明, 赵发珠. 秦岭北麓不同植被类型土壤有机碳含量分布及其影响因素[J]. 植物生态学报, 2026, 50(预发表): 1-. |
| [2] | 张静, 陈洁, 李艳朋, 盘李军, 许涵, 李意德, 何海生. 南亚热带针阔混交人工林植物生物量比较及其影响因子分析[J]. 植物生态学报, 2026, 50(2): 400-416. |
| [3] | 徐恩相, 周蕾, 章晓炜, 张国萍, 仲杜伟, 黄智, 刘派, 迟永刚. 基于不同生育阶段冠层光谱和碳通量的水稻产量估算[J]. 植物生态学报, 2026, 50(1): 82-93. |
| [4] | 张健, 孙慕梵, 宋坤, 赖丽茗, 宋永昌. 基于生态外貌-植物区系分类系统的中国常绿阔叶林分类修订[J]. 植物生态学报, 2026, 50(1): 107-122. |
| [5] | 张琳, 袁伟影, 宋创业, 吴冬秀. 1998-2010年中国典型生态系统长期监测样地环境要素、物种丰富度和生物量动态数据集[J]. 植物生态学报, 2025, 49(8): 1182-1190. |
| [6] | 严文秀, 赵诗晗, 郑春燕, 张萍, 沈海花, 常锦峰, 徐亢. 基于多物候指标的人工饲草长势监测及产量估测[J]. 植物生态学报, 2025, 49(7): 1096-1109. |
| [7] | 范亚冉, 夏少攀, 于冰冰, 朱紫琪, 杨威, 范豫川, 刘晓雨, 张旭辉, 郑聚锋. 大气CO2浓度升高和增温对土壤有机碳库积累、分子组成和结构稳定性的影响[J]. 植物生态学报, 2025, 49(7): 1053-1069. |
| [8] | 杜燕, 刘鑫, 张瀚曰, 马少伟, 包维楷. 中国高山松群系的群落特征[J]. 植物生态学报, 2025, 49(6): 816-824. |
| [9] | 熊高明, 申国珍, 赵常明, 徐文婷, 王杨, 谢宗强, 李家湘, 徐耀粘, 李跃林, 陈芳清. 中国檵木灌丛植被类型及群落特征[J]. 植物生态学报, 2025, 49(6): 875-887. |
| [10] | 惠城阳, 章巧依, 刘腾腾, 刘维勇, 周丽娜, 金鑫杰, 张永华, 刘金亮. 温州大罗山主要植被类型及物种组成特征[J]. 植物生态学报, 2025, 49(6): 990-998. |
| [11] | 曹毅, 张松林, 王旭峰, 杨安昌, 任敏慧, 杨浩, 韩超. 兰州市南北两山植物群落数据集[J]. 植物生态学报, 2025, 49(6): 975-989. |
| [12] | 熊高明, 申国珍, 徐文婷, 谢宗强, 李跃林, 徐耀粘, 陈芳清, 李家湘. 中国低山丘陵热性常绿阔叶灌丛主要类型及群落特征[J]. 植物生态学报, 2025, 49(6): 865-874. |
| [13] | 李家湘, 刘文倩, 蒋国平, 赵丽娟, 徐永福, 吴磊, 喻勋林. 植被科学画符号系统设计和标准化绘图范式——以森林群落剖面图为例[J]. 植物生态学报, 2025, 49(6): 897-910. |
| [14] | 郝杰, 刁华杰, 苏原, 武帅楷, 高阳阳, 梁雯君, 牛慧敏, 杨倩雯, 常婕, 王袼, 许雯丽, 马腾飞, 董宽虎, $\boxed{\hbox{王常慧}}$. 降水调控农牧交错带盐渍化草地净初级生产力对氮添加及刈割的响应[J]. 植物生态学报, 2025, 49(5): 710-719. |
| [15] | 王贝贝, 吴苏, 王苗苗, 胡锦涛. 日光诱导叶绿素荧光不同组分在作物总初级生产力估算中的贡献比例: 多时间尺度分析[J]. 植物生态学报, 2025, 49(4): 562-572. |
| 阅读次数 | ||||||
|
全文 |
|
|||||
|
摘要 |
|
|||||
Copyright © 2026 版权所有 《植物生态学报》编辑部
地址: 北京香山南辛村20号, 邮编: 100093
Tel.: 010-62836134, 62836138; Fax: 010-82599431; E-mail: apes@ibcas.ac.cn, cjpe@ibcas.ac.cn
备案号: 京ICP备16067583号-19