留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

大亚湾夏季浮游植物群落结构的长期变化及其与环境因子的关系

杨熙 李开枝 谭烨辉 吕意华

杨熙,李开枝,谭烨辉,等. 大亚湾夏季浮游植物群落结构的长期变化及其与环境因子的关系[J]. 海洋学报,2022,44(8):110–122 doi: 10.12284/hyxb2022130
引用本文: 杨熙,李开枝,谭烨辉,等. 大亚湾夏季浮游植物群落结构的长期变化及其与环境因子的关系[J]. 海洋学报,2022,44(8):110–122 doi: 10.12284/hyxb2022130
Yang Xi,Li Kaizhi,Tan Yehui, et al. Long-term changes of phytoplankton community structure with relation to environmental factors in the Daya Bay in summer[J]. Haiyang Xuebao,2022, 44(8):110–122 doi: 10.12284/hyxb2022130
Citation: Yang Xi,Li Kaizhi,Tan Yehui, et al. Long-term changes of phytoplankton community structure with relation to environmental factors in the Daya Bay in summer[J]. Haiyang Xuebao,2022, 44(8):110–122 doi: 10.12284/hyxb2022130

大亚湾夏季浮游植物群落结构的长期变化及其与环境因子的关系

doi: 10.12284/hyxb2022130
基金项目: 自然资源部海洋环境探测技术与应用重点实验室2020年度自主设立课题(MESTA-2020-C006);广东省平台基地及科技基础条件建设项目(2021B1212050025);中国海洋发展基金会项目(CODF-002-ZX-2021)。
详细信息
    作者简介:

    杨熙(1989-),男,湖北省监利市人,博士,主要从事海草生态学、海洋浮游生物生态学研究。E-mail:yang1209xi@163.com

  • 中图分类号: Q178.53

Long-term changes of phytoplankton community structure with relation to environmental factors in the Daya Bay in summer

  • 摘要: 本文分析了1999−2017年大亚湾夏季浮游植物群落结构的长期变化及其与环境因子的关系,结果显示,大亚湾海域海水温度呈显著下降趋势,盐度呈显著上升趋势;溶解无机氮浓度出现较大幅度提升,2008−2017年间大亚湾溶解无机氮浓度平均值比1999−2007年提升了72.73%;大亚湾浮游植物种类数变化趋势不明显,主要优势种没有发生明显变化,柔弱伪菱形藻(Pseudonitzschia delicatissima)为区域第一优势种,其次为中肋骨条藻(Skeletonema costatum);浮游植物总丰度、硅藻丰度、甲藻丰度以及主要种类中的柔弱伪菱形藻、中肋骨条藻和叉角藻(Ceratium furca)丰度均呈现显著上升趋势;浮游植物生物多样性指数(H′)和均匀度(J)均呈下降趋势。人类活动所引起的溶解无机氮浓度大幅升高以及外海水入侵加强所引起的海水温度降低和盐度上升导致了浮游植物丰度的上升、优势种的单一化和生物多样性指数的下降。
  • 图  1  调查区域和站位分布

    “+”为中国科学院大亚湾海洋生物综合实验站调查站位;“▲”为国家海洋局南海环境监测中心调查站位“+” represents the stations of Daya Bay Marine Biology Research Station, Chinese Academy of Sciences; “▲” represents the stations of South China Sea Environmental Monitoring Center, State Oceanic Administration

    Fig.  1  Distribution of study area and monitoring stations

    图  2  1999−2017年大亚湾夏季环境因子的长期变化

    Fig.  2  Long-term changes of environmental variables in the Daya Bay in summer between 1999 and 2017

    图  3  1999−2017年大亚湾夏季浮游植物种类数长期变化

    Fig.  3  Long-term changes of phytoplankton species number in the Daya Bay in summer between 1999 and 2017

    图  4  1999−2017年大亚湾夏季浮游植物生物量的长期变化

    Fig.  4  Long-term changes of phytoplankton groups abundance in the Daya Bay in summer between 1999 and 2017

    图  5  1999−2017年大亚湾夏季几种主要浮游植物种类丰度和占比的长期变化

    Fig.  5  Long-term changes of phytoplankton dominant species abundance and their abundance percentages in the Daya Bay in summer between 1999 and 2017

    图  6  1999−2017年大亚湾夏季浮游植物生物多样性指数(H′)和均匀度指数(J)的长期变化

    Fig.  6  Long-term changes of species diversity indes (H) and evenness indes (J) in the Daya Bay in summer between 1999 and 2017

    图  7  冗余分析图

    a. 各年份与环境因子的相关性;b. 不同浮游植物种类丰度(小向量)与环境因子(大向量)的相关性;图中代码所代表的环境因子和浮游植物种类见表2;其中环境因子对浮游植物丰度变化的解释率为29.30%,前两轴的特征值分别用λ1和λ2表示

    Fig.  7  Correlation plots of the redundancy analysis

    a. The relationship between the environmental factors and sample years; b. the relationship between the abundance of phytoplankton species (small vectors) and environmental factors (large vectors); codes of the phytoplankton species are as in Table 2; the plots display 29.30% of the variance in the phytoplankton abundance, and eigenvalues of the first two axes are indicated by λ1 and λ2

    表  1  调查数据基本信息

    Tab.  1  Basic information of survey data

    来源调查时间完成站位数完成调查的具体站位
    中国科学院大亚湾海洋生物综合实验站1999年8月4日10S3、S4、S5、S6、S7、S8、S9、S10、S11、S12
    2000年7月22日10S3、S4、S5、S6、S7、S8、S9、S10、S11、S12
    2001年8月9日12S1、S2、S3、S4、S5、S6、S7、S8、S9、S10、S11、S12
    2002年8月9日12S1、S2、S3、S4、S5、S6、S7、S8、S9、S10、S11、S12
    2003年8月27日12S1、S2、S3、S4、S5、S6、S7、S8、S9、S10、S11、S12
    国家海洋局南海
    环境监测中心
    2004年7月18 −19日8D1、D2、D4、D7、D9、D12、D14、D16
    2005年7月29日至8月3日12D1、D2、D4、D5、D7、D9、D11、D12、D14、D15、D16、D17
    2006年8月29日至9月1日12D1、D2、D4、D5、D7、D9、D11、D12、D14、D15、D16、D17
    2007年8月12−13日12D1、D2、D4、D5、D7、D9、D11、D12、D14、D15、D16、D17
    2008年8月28 −31日12D1、D2、D4、D5、D7、D9、D11、D12、D14、D15、D16、D17
    2009年8月2−6日12D1、D2、D4、D5、D7、D9、D11、D12、D14、D15、D16、D17
    2010年8月29日至9月1日12D1、D2、D4、D5、D7、D9、D11、D12、D14、D15、D16、D17
    2011年8月24 −27日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    2012年8月14 −19日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    2013年8月27−31日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    2014年8月19−21日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    2015年8月21 −24日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    2016年8月24 −28日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    2017年8月9 −13日18D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D11、D12、D13、D14、D15、D16、D17、D18
    注:数据包括表底层温度、表底层盐度、表底层营养盐(NO3-N、NO2-N、NH4-N和PO4-P)浓度、网采浮游植物丰度。
    下载: 导出CSV

    表  2  环境因子和浮游植物群落结构长期变化趋势的Mann-Kendall检验结果

    Tab.  2  Results of the Mann-Kendall test for detection of long-term changes in the environmental factors and the phytoplankton community structure

    参数(缩写)单位最小值至最大值平均值年份数统计值显著性(p
    温度(T26.25~30.5028.3319−2.10<0.05
    盐度(S27.82~33.8731.53192.17<0.05
    无机氮浓度(DIN)µg/L26.54~203.2173.16190.98>0.05
    无机磷浓度(DIP)µg/L1.86~17.798.00190.42>0.05
    氮磷比(DIN∶DIP)5.85~85.4630.75190.21>0.05
    总丰度(Total)105 cells/m37.22~3605.27471.85192.24<0.05
    硅藻丰度(DIAT)105 cells/m35.46~3581.99468.24192.17<0.05
    角毛藻属丰度(Cha-sp.)105 cells/m30.01~259.5430.29190.56>0.05
    圆筛藻属丰度(Cos-sp.)105 cells/m30.03~2.470.34192.31<0.05
    北方劳德藻丰度(Lau-bo)105 cells/m30~129.247.23161.76<0.05
    柔弱伪菱形藻丰度(Pse-de)105 cells/m30.01~1579.22201.91192.03<0.05
    尖刺伪菱形藻丰度(Pse-pu)105 cells/m30~14.612.4417−0.87>0.05
    根管藻属丰度(Rhi-sp.)105 cells/m30.08~11.792.27190.56>0.05
    中肋骨条藻丰度(Ske-co)105 cells/m30.17~2692.95186.33192.03<0.05
    掌状冠盖藻丰度(Ste-pa)105 cells/m30~0.350.0410−1.61>0.05
    扭鞘藻丰度(Str-th)105 cells/m30~1.700.17113.58<0.01
    菱形海线藻丰度(Tha-ni)105 cells/m30~178.4111.00191.87<0.05
    伏氏海毛藻丰度(Tha-fr)105 cells/m30.23~210.0517.11180.30>0.05
    甲藻丰度(DINO)104 cells/m33.30~218.0534.71192.73<0.05
    短角角藻丰度(Cer-br)104 cells/m30~1.550.199 −2.04 <0.05
    叉角藻丰度(Cer-fur)104 cells/m30.14~115.1315.38193.15<0.01
    纺锤角藻丰度(Cer-fus)104 cells/m30.10~6.011.6019−0.63>0.05
    大角角藻丰度(Cer-ma)104 cells/m30~8.291.74182.12<0.05
    三叉角藻丰度(Cer-tr)104 cells/m30~7.411.4711−0.93>0.05
    具尾鳍藻丰度(Din-ca)104 cells/m30~6.590.80170.78>0.05
    夜光藻丰度(Noc-sc)104 cells/m30~2.670.7114−0.66>0.05
    原多甲藻属丰度(Pro-sp.)104 cells/m30.17~13.651.88191.26>0.05
    H′0.79~3.882.3619−2.10<0.05
    J0.17~0.750.5019−1.89<0.05
    下载: 导出CSV

    表  3  1999−2017年大亚湾夏季浮游植物优势种优势度的长期变化

    Tab.  3  Dominance changes of phytoplankton dominant species in the Daya Bay in summer between 1999 and 2017

    中文名拉丁名类群1999年2000年2001年2002年2003年2004年2005年2006年2007年2008年2009年2010年2011年2012年2013年2014年2015年2016年2017年
    日本星杆藻Asterionella japonica硅藻0.035
    变异辐杆藻Bacteriastrum varians硅藻0.044
    海洋角管藻Cerataulina pelagica硅藻0.021
    窄隙角毛藻Chaetoceros affinis硅藻0.0520.0580.0540.085
    扁面角毛藻Chaetoceros compressus硅藻0.036
    旋链角毛藻Chaetoceros curvisetus硅藻0.0780.0320.0320.0230.026
    柔弱角毛藻Chaetoceros debilis硅藻0.0610.0210.034
    并基角毛藻Chaetoceros decipiens硅藻0.041
    远距角毛藻Chaetoceros distans硅藻0.022
    异角角毛藻Chaetoceros diversus硅藻0.1200.0310.140
    罗氏角毛藻Chaetoceros lauderi硅藻0.022
    洛氏角毛藻Chaetoceros lorenzianus硅藻0.0540.0390.027
    暹罗角毛藻Chaetoceros siamense硅藻0.167
    中心圆筛藻Coscinodiscus centralis硅藻0.255
    长角弯角藻Eucampia cornuta硅藻0.047
    薄壁半管藻Hemiaulus membranaceus硅藻0.020
    北方劳德藻Lauderia borealis硅藻0.030
    柔弱伪菱形藻Pseudo-nitzschia delicatissima硅藻0.1290.2920.2590.0840.9580.7410.8130.6360.5260.1800.5240.8580.3830.154
    尖刺伪菱形藻Pseudonitzschia pungens硅藻0.1620.0230.1040.0360.0430.071
    翼根管藻纤细变型Rhizosolenia alata f. gracillima硅藻0.634
    笔尖根管藻Rhizosolenia styliformis硅藻0.022
    中肋骨条藻Skeletonema costatum硅藻0.2660.0260.2420.0420.5010.0850.0420.5320.0250.1470.0260.0750.110
    热带骨条藻Skeletonema tropicum硅藻0.0230.637
    菱形海线藻Thalassionema nitzschioides硅藻0.0290.3330.0480.2330.3220.0480.0430.0640.0220.1190.0630.0380.037
    伏氏海毛藻Thalassiothrix fraenfeldii硅藻0.0610.2830.0240.212
    叉角藻Ceratium furca甲藻0.0210.0230.0240.266
    纺锤角藻Ceratium fusus甲藻0.029
    大角角藻Ceratium macroceros甲藻0.024
    三叉角藻Ceratium trichoceors甲藻0.0260.021
    锥状斯克里普藻Scrippsiella trochoidea甲藻0.027
    注:空白表示不占优势。
    下载: 导出CSV
  • [1] Bricker S B, Longstaff B, Dennison W, et al. Effects of nutrient enrichment in the nation’s estuaries: a decade of change[J]. Harmful Algae, 2008, 8(1): 21−32. doi: 10.1016/j.hal.2008.08.028
    [2] Boyce D G, Lewis M R, Worm B. Global phytoplankton decline over the past century[J]. Nature, 2010, 466(7306): 591−596. doi: 10.1038/nature09268
    [3] Field C B, Behrenfeld M J, Randerson J T, et al. Primary production of the biosphere: integrating terrestrial and oceanic components[J]. Science, 1998, 281(5374): 237−240. doi: 10.1126/science.281.5374.237
    [4] Murtugudde R, Beauchamp J, McClain C R, et al. Effects of penetrative radiation on the upper tropical ocean circulation[J]. Journal of Climate, 2002, 15(5): 470−486. doi: 10.1175/1520-0442(2002)015<0470:EOPROT>2.0.CO;2
    [5] Sabine C L, Feely R A, Gruber N, et al. The oceanic sink for anthropogenic CO2[J]. Science, 2004, 305(5682): 367−371. doi: 10.1126/science.1097403
    [6] Harding Jr L W, Gallegos C L, Perry E S, et al. Long-term trends of nutrients and phytoplankton in Chesapeake Bay[J]. Estuaries and Coasts, 2016, 39(3): 664−681. doi: 10.1007/s12237-015-0023-7
    [7] Yoshida K, Chiba S, Ishimaru T. Long-term variation in the wintertime diatom community structure in Tokyo Bay, Japan (1981–2000)[J]. Plankton and Benthos Research, 2011, 6(4): 195−205. doi: 10.3800/pbr.6.195
    [8] Li Yaqin, Smayda T J. Temporal variability of chlorophyll in Narragansett Bay, 1973−1990[J]. ICES Journal of Marine Science, 1998, 55(4): 661−667. doi: 10.1006/jmsc.1998.0383
    [9] 丘耀文, 王肇鼎, 朱良生. 大亚湾海域营养盐与叶绿素含量的变化趋势及其对生态环境的影响[J]. 台湾海峡, 2005, 24(2): 131−139.

    Qiu Yaowen, Wang Zhaoding, Zhu Liangsheng. Variation trend of nutrient and chlorophyll contents and their effects on ecological environment in Daya Bay[J]. Journal of Oceanography in Taiwan Strait, 2005, 24(2): 131−139.
    [10] Ke Zhixin, Tan Yehui, Huang Liangmin, et al. Significantly depleted 15N in suspended particulate organic matter indicating a strong influence of sewage loading in Daya Bay, China[J]. Science of the Total Environment, 2019, 650(Pt 1): 759−768.
    [11] Wang Youshao, Lou Zhiping, Sun Cuici, et al. Ecological environment changes in Daya Bay, China, from 1982 to 2004[J]. Marine Pollution Bulletin, 2008, 56(11): 1871−1879. doi: 10.1016/j.marpolbul.2008.07.017
    [12] 王友绍, 王肇鼎, 黄良民. 近20年来大亚湾生态环境的变化及其发展趋势[J]. 热带海洋学报, 2004, 23(5): 85−95. doi: 10.3969/j.issn.1009-5470.2004.05.012

    Wang Youshao, Wang Zhaoding, Huang Liangmin. Environment changes and trends in Daya Bay in recent 20 years[J]. Journal of Tropical Oceanography, 2004, 23(5): 85−95. doi: 10.3969/j.issn.1009-5470.2004.05.012
    [13] Wu Meilin, Wang Youshao, Wang Yutu, et al. Scenarios of nutrient alterations and responses of phytoplankton in a changing Daya Bay, South China Sea[J]. Journal of Marine Systems, 2017, 165: 1−12. doi: 10.1016/j.jmarsys.2016.09.004
    [14] Hao Yanju, Tang Danling, Boicenco L, et al. Environmental ecological response to increasing water temperature in the Daya Bay, Southern China in 1982−2012[J]. Natural Resources, 2016, 7(4): 184−192. doi: 10.4236/nr.2016.74017
    [15] Song Xingyu, Huang Liangmin, Zhang Jianlin, et al. Harmful algal blooms (HABs) in Daya Bay, China: an in situ study of primary production and environmental impacts[J]. Marine Pollution Bulletin, 2009, 58(9): 1310−1318. doi: 10.1016/j.marpolbul.2009.04.030
    [16] 徐恭昭. 大亚湾环境与资源[M]. 合肥: 安徽科学技术出版社, 1989.

    Xu Gongzhao. Environments and Resources of Daya Bay[M]. Hefei: Anhui Science and Technology Press, 1989.
    [17] Shannon C E. A mathematical theory of communication[J]. Bell System Technical Journal, 1948, 27(3): 379−423. doi: 10.1002/j.1538-7305.1948.tb01338.x
    [18] Whittaker R H. Evolution and measurement of species diversity[J]. Taxon, 1972, 21(2/3): 213−251.
    [19] 许金电, 蔡尚湛, 宣莉莉, 等. 粤东至闽南沿岸海域夏季上升流的调查研究[J]. 热带海洋学报, 2014, 33(2): 1−9. doi: 10.3969/j.issn.1009-5470.2014.02.001

    Xu Jindian, Cai Shangzhan, Xuan Lili, et al. Observational study on summertime upwelling in coastal seas between eastern Guangdong and southern Fujian[J]. Journal of Tropical Oceanography, 2014, 33(2): 1−9. doi: 10.3969/j.issn.1009-5470.2014.02.001
    [20] 韩舞鹰, 马克美. 大亚湾海水混合交换特征[J]. 海洋科学, 1991, 15(2): 64−67.

    Han Wuying, Ma Kemei. Study on the process of sea water exchange in Daya Bay[J]. Marine Sciences, 1991, 15(2): 64−67.
    [21] 杨熙, 谭烨辉. 夏季外海水入侵对大亚湾浮游植物群落结构的影响[J]. 海洋科学, 2019, 43(7): 96−105. doi: 10.11759/hykx20190321003

    Yang Xi, Tan Yehui. Effects of shelf seawater intrusion on phytoplankton community structure in Daya Bay in the summer[J]. Marine Sciences, 2019, 43(7): 96−105. doi: 10.11759/hykx20190321003
    [22] 柯伟雄. 大亚湾统计年鉴(2018)[M]. 惠州: 惠州大亚湾经济技术开发区统计局, 2018.

    Ke Weixiong. Daya Bay Statistical Yearbook (2018)[M]. Huizhou: Statistics Bureau of Huizhou Daya Bay Economic and Technological Development Zone, 2018.
    [23] 黄小平, 黄良民, 宋金明, 等. 营养物质对海湾生态环境影响的过程与机理[M]. 北京: 科学出版社, 2019.

    Huang Xiaoping, Huang Liangmin, Song Jinming, et al. Process and Mechanism of the Eco-Environmental Effects of Nutrients on the Semi-Enclosed Bay[M]. Beijing: Science Press, 2019.
    [24] 彭云辉, 王肇鼎, 高红莲, 等. 大亚湾大鹏澳养殖网箱水体无机氮的生物地球化学[J]. 海洋通报, 2001, 20(2): 16−24. doi: 10.3969/j.issn.1001-6392.2001.02.003

    Peng Yunhui, Wang Zhaoding, Gao Honglian, et al. Biogeochemistry of inorganic nitrogen in the water of aquicultural boxes at Dapeng’ao, Daya Bay[J]. Marine Science Bulletin, 2001, 20(2): 16−24. doi: 10.3969/j.issn.1001-6392.2001.02.003
    [25] Justić D, Rabalais N N, Turner R E. Stoichiometric nutrient balance and origin of coastal eutrophication[J]. Marine Pollution Bulletin, 1995, 30(1): 41−46. doi: 10.1016/0025-326X(94)00105-I
    [26] 刘胜, 黄晖, 黄良民, 等. 大亚湾核电站对海湾浮游植物群落的生态效应[J]. 海洋环境科学, 2006, 25(2): 9−12, 25. doi: 10.3969/j.issn.1007-6336.2006.02.003

    Liu Sheng, Huang Hui, Huang Liangmin, et al. Ecological response of phytoplankton to the operation of Daya Bay nuclear power station[J]. Marine Environmental Science, 2006, 25(2): 9−12, 25. doi: 10.3969/j.issn.1007-6336.2006.02.003
    [27] 韦桂峰, 王肇鼎, 练健生. 大亚湾大鹏澳水域春季浮游植物优势种的演替[J]. 生态学报, 2003, 23(11): 2285−2292. doi: 10.3321/j.issn:1000-0933.2003.11.011

    Wei Guifeng, Wang Zhaoding, Lian Jiansheng. Succession of dominant phytoplankton species in spring 2002 at Dapeng’ao Cove, Daya Bay[J]. Acta Ecologica Sinica, 2003, 23(11): 2285−2292. doi: 10.3321/j.issn:1000-0933.2003.11.011
    [28] Baek S H, Shimode S, Han M S, et al. Growth of dinoflagellates, Ceratium furca and Ceratium fusus in Sagami Bay, Japan: the role of nutrients[J]. Harmful Algae, 2008, 7(6): 729−739. doi: 10.1016/j.hal.2008.02.007
    [29] Baek S H, Shimode S, Kikuchi T. Growth of dinoflagellates, Ceratium furca and Ceratium fusus in Sagami Bay, Japan: the role of temperature, light intensity and photoperiod[J]. Harmful Algae, 2008, 7(2): 163−173. doi: 10.1016/j.hal.2007.06.006
    [30] 孙晓霞, 孙松, 赵增霞, 等. 胶州湾营养盐浓度与结构的长期变化[J]. 海洋与湖沼, 2011, 42(5): 662−669. doi: 10.11693/hyhz201105006006

    Sun Xiaoxia, Sun Song, Zhao Zengxia, et al. Long-term changes in nutrient concentration and structure in the Jiaozhou Bay[J]. Oceanologia et Limnologia Sinica, 2011, 42(5): 662−669. doi: 10.11693/hyhz201105006006
    [31] 吴玉霖, 孙松, 张永山. 环境长期变化对胶州湾浮游植物群落结构的影响[J]. 海洋与湖沼, 2005, 36(6): 487−498. doi: 10.3321/j.issn:0029-814X.2005.06.002

    Wu Yulin, Sun Song, Zhang Yongshan. Long-term change of environment and it’s influence on phytoplankton community structure in Jiaozhou Bay[J]. Oceanologia et Limnologia Sinica, 2005, 36(6): 487−498. doi: 10.3321/j.issn:0029-814X.2005.06.002
  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  450
  • HTML全文浏览量:  186
  • PDF下载量:  83
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-01-06
  • 录用日期:  2022-03-09
  • 修回日期:  2022-02-28
  • 网络出版日期:  2022-08-15
  • 刊出日期:  2022-08-15

目录

    /

    返回文章
    返回