Utilizing LiDAR technology to estimate forest aboveground biomass in Qianjiangyuan National Park, Jiangxi Province of eastern China
-
摘要:目的 采用激光雷达结合地面调查的方法，拟合并预测钱江源国家公园森林的地上生物量，为国家公园森林生态系统保护和动物生境质量评估提供基础数据。方法 采用冠层结构指数量化亚热带常绿阔叶林的冠层高度和垂直复杂度，并结合地面调查数据和地形数据，用模型拟合和预测钱江源国家公园的地上生物量。结果 冠层结构指数中，与地上生物量的相关程度从大到小依次为冠层平均高度、冠层垂直复杂度、冠层最大高度。钱江源国家公园生物量分布主要在27.24 ~ 210.31 Mg/hm2（0.05 ~ 0.95分位数），平均为111.21 Mg/hm2，总的生物量约为2.57 × 106 Mg。结论 结合冠层指数和地形变量能较好地模拟和预测钱江源国家公园地上生物量。钱江源国家公园有大面积低海拔集中分布的常绿阔叶林，冠层平均高度和冠层垂直复杂度都比较大，地上生物量高；同时，也有大面积人工林和次生林，这些人工林和次生林主要分布在国家公园成立后才开始保护的两个片区，冠层平均高度和冠层垂直复杂度小，地上生物量较低，有待进一步恢复。Abstract:Objective Utilizing technologies of airborne lidar together with ground survey, we modelled and predicted aboveground biomass of subtropical forests in Qianjiangyuan National Park, Jiangxi Province of eastern China to provide basis for assessing conservation effectiveness and quality of animal habitats in national parks.Method We used canopy structural indices to quantify canopy structural variation and maximal height of subtropical forests, then modelled forest aboveground biomass with canopy structural indices and topographic variables.Result Mean canopy height, vertical complexity and maximal height were closely correlated with aboveground biomass. Aboveground biomass was mainly between 27.24−210.31 Mg/ha (quantile 0.05−0.95), mean value was 111.21 Mg/ha, and total aboveground biomass was2.57 × 106 Mg.Conclusion The forest aboveground biomass is well modelled and predicted using canopy structural indices and topographic variables. Large area of zonal evergreen broadleaved forests is distributed in the national park with complex canopy structure and high aboveground biomass, while large area plantation and secondary forests are also distributed in two new parts of the national park with simple canopy structure and low aboveground biomass.
-
板栗（Castanea mollissima）属于壳斗科（Fagaceae）栗属（Castanea）植物，是我国重要的木本粮食树种之一，其果实营养价值高，除富含淀粉外，还有蛋白质、脂肪、矿物质和维生素等营养物质，深受消费者喜爱，具有较高的经济价值[1]。目前我国板栗种植规模和总产量均居世界第一，但由于板栗喜光，枝条顶端优势强，结果部位外移快且自交不亲和，其结实率低，效益低下[2]。因此，培育具有优质、高产、短枝型等优良性状的板栗新品种具有重要的现实意义。
现阶段分子辅助育种已被广泛应用到植物育种工作中，但板栗新品种的培育仍以杂交育种为主[3]，如：新品种‘明丰2号’（C. mollissima ‘Mingfeng 2’）是以‘燕明’（C. mollissima ‘Yanming’）为母本、‘燕山早丰’（C. mollissima ‘Yanshanzaofeng’）为父本杂交培育而来[4]；‘南垂5号’（C. mollissima ‘Nanchui5’）则是以母本‘垂枝’（C. mollissima ‘Chuizhi’）与父本‘南沟1号’（C. mollissima ‘Nangou1’）杂交培育得来[5]。在杂交授粉过程中，如何提高成功率，以获取丰富的杂种资源是板栗育种亟待解决的问题。较高活性的花粉和柱头是授粉成功的关键[6]，准确掌握花粉与柱头活性检测方法，探明其日变化规律，明确最适杂交授粉时间，可有效提高板栗杂交育种成功率。
目前，植物花粉活力检测方法多样，如2，3，5-三苯基氯化四氮唑（TTC）染色法、碘–碘化钾（I2-KI）法、MTT染色法、离体萌发法、墨水染色法等[7−9]。在柱头可授性检测方面，多数学者采用联苯胺–过氧化氢法[10−11]检测柱头可授性，部分研究采用MTT染色法[12]、直接授粉法等[13]。前期板栗花粉活力研究发现，不同板栗品种采用不同花粉活力测定方法，结果存在较大差异。陈可可[14]采用TTC染色法检测5个板栗品种花粉活力发现，‘泰安薄壳栗’（C. mollissima ‘Taianbaoke’）花粉活力最高可达72.18%；覃腾[15]以广西4个板栗品种为研究对象，采用离体萌发法检测其花粉活力，发现‘封开2号’（C. mollissima ‘Fengkai2’）萌发率最高，但仅为15.76%。在板栗柱头可授性研究方面，王倩[16]采用联苯胺–过氧化氢法检测了不同时期燕山板栗柱头可授性，发现盛花期柱头可授性最强；周晶[17]采用联苯胺–过氧化氢法检测板栗‘燕山红’（C. mollissima ‘Yanshanhong’）柱头可授性发现，柱头可授性高峰值出现于盛花期单日08:00—12:00。以上研究对花粉活力和柱头可授性的测定方法和日变化规律进行了初步探索，但并未将花粉活力、柱头可授性日变化规律与其坐果率、结实率、空苞率等进行综合分析。
为准确把握板栗最适授粉时间，提高板栗杂交育种成功率，为后期的良种培育提供充足的杂种资源，本研究首次以河北省迁西县主栽国审板栗品种‘燕山早丰’（果肉质地细腻、味香甜、熟食品质上等、丰产、质感效应较强）为母本[18−19]，父本选择具有多粒、矮化等特异性状的‘王钱大栗’（C. mollissima ‘Wangqiandali’）、‘短枝芽变’（C. mollissima ‘Duanzhiyabian’）、‘三屯红栗’（C. mollissima ‘Santunhongli’），采用联苯胺–过氧化氢法、MTT法检测母本柱头可授性，采用TTC染色法、I2-KI法、醋酸洋红染色法、离体萌发法等检测花粉活力，通过研究盛花期花粉活力、柱头可授性、不同时间授粉对结实情况的影响，筛选‘燕山早丰’柱头可授性和杂交花粉活力的适宜检测方法，探明两者的日变化规律，筛选出最适杂交授粉时间。研究结果将为板栗人工杂交育种的成功提供理论指导，同时对解析板栗生殖生物学特性及繁殖机制具有重要意义。
1. 材料与方法
1.1 试验地概况
试验地位于河北省唐山市迁西县北京林业大学经济林（板栗）育种与栽培试验实践基地（40°18′ N，118°54′ E），属暖温带大陆性季风气候，年平均气温10.1 ℃，最低月平均气温−6.5 ℃，最高月平均气温25.4 ℃。年平均降水量817 mm，其中5—10月降雨量657.6 mm，占全年降水量的88%。年平均相对湿度59%，年平均无霜期176 d，全年日照2 581.5 h，日照充足。该地属于低山丘陵区，土壤质地为沙壤土，pH 6.44，有机质含量2.89 mg/g，土壤肥力为中等水平[20]。
1.2 试验材料
材料为11 年生板栗品种‘燕山早丰’‘王钱大栗’‘短枝芽变’‘三屯红栗’，长势优良，无病虫害，株行距4 m × 3 m，平均树高4.3 m，冠幅4.0 m × 4.5 m，树体为三主枝，开张角度 > 60°，树形为自然开心形。
1.3 研究方法
1.3.1 试验设计
参照周晶[17]研究结果，选择盛花期进行人工授粉。选取长势一致且健壮的板栗试验树60株，包括‘燕山早丰’45株，‘王钱大栗’‘短枝芽变’‘三屯红栗’各5株。40株‘燕山早丰’作为母本进行柱头可授性检测与授粉试验，5株‘燕山早丰’用于花粉活力检测；‘王钱大栗’‘短枝芽变’‘三屯红栗为父本，用于花粉活力日变化规律研究与授粉试验。各杂交授粉组合为 ‘燕山早丰’♀ × ‘王钱大栗’♂、‘燕山早丰’♀ × ‘短枝芽变’♂、‘燕山早丰’♀ × ‘三屯红栗’♂和1个自交组合‘燕山早丰’♀ × ‘燕山早丰’♂。
1.3.2 人工授粉
每授粉组合选择10株母树，每棵母树去雄套袋后，于盛花期单日08:00—18:00每隔2 h进行杂交授粉。采用随机区组设计，以单株为试验小区，共10个重复。各授粉组合每株母树每授粉时段在冠层中部东、南、西、北4个方向外围各均匀选取3个长势相近且健壮的雌花簇，1个授粉时段每母株共选取12个雌花簇，每个授粉组合每时段授粉120个雌花簇，4个组合共选取2 400个雌花簇。授粉时去掉套袋，用毛笔蘸取花粉轻轻点涂雌花柱头上，授粉完成后立即套袋。
1.3.3 花粉活力测定方法筛选
在上午08:00—10:00，采集‘燕山早丰’雄花盛花期的雄花序，每株树均从东、南、西、北4个方向取样，混合后分别采用TTC染色法[21]、醋酸洋红染色法[22]、I2-KI染色法[23]、离体萌发培养法[24]测定花粉活力，并以清水作为对照，筛选板栗花粉活力适宜检测方法。
1.3.4 花粉活力日变化规律研究
在08:00—18:00每隔2 h采集‘燕山早丰’‘王钱大栗’‘短枝芽变’‘三屯红栗’的雄花盛花期的花粉，选择1.3.3筛选出的花粉活力检测方法进行花粉活力日变化规律研究。
1.3.5 柱头可授性检测方法筛选
从40株‘燕山早丰’授粉母株中随机选取15株，在雌花盛花期每株树冠层中部东、南、西、北4个方向各均匀选取3个雌花簇，在上午08:00—10:00采集柱头，分别利用联苯胺–过氧化氢法[25]和MTT染色法[26]测定柱头可授性，并以清水作为对照，筛选板栗柱头可授性适宜检测方法。联苯胺–过氧化氢法判定标准：若柱头周围产生气泡数量多且体积较大，则柱头可授性很强（+++）；若柱头周围产生大量小气泡，则柱头可授性较强（++）；若柱头周围产生少量气泡，则柱头可授性较弱（+）；若柱头无气泡产生，则柱头无可授性（−）。MTT染色法判定标准：若柱头具有可授性，则柱头被染成深蓝色或黑色，颜色由浅变深，柱头可授性由弱变强；若没有可授性，则柱头没有颜色反应。
1.3.6 柱头可授性日变化规律研究
在08:00—18:00每隔2 h采集‘燕山早丰’盛花期的柱头，选择1.3.5筛选出的柱头可授性检测方法进行柱头可授性日变化规律研究。
1.3.7 结实情况统计
授粉20 d后去除纸袋，使球苞自然生长，统计坐果率、落果率，9月上旬果实成熟后，计算结实率、空苞率。坐果率为未凋落雌花簇数与授粉雌花簇数的比值，落果率为凋落雌花簇数与授粉雌花簇数的比值，空苞率为空苞数与授粉雌花簇数的比值，结实率为果实成熟时的结实苞数与授粉雌花簇数的比值[15]。
1.4 数据分析
试验数据采用Excel软件进行整理，SPSS 26.0进行方差分析。
2. 结果与分析
2.1 花粉活力适宜检测方法筛选
花粉在清水处理的情况下，颜色均为淡黄色，接近半透明，分为形状规则的卵圆形与形状不规则两种形态（图1a）。0.5%和2%的TTC染色法的染色效果均不明显，延长染色时间后花粉依旧呈淡黄色（图1b，1c），故此方法不适合板栗花粉活力检测。醋酸洋红染色法效果明显，有活力的花粉呈鲜红色，形状呈规则的卵圆形，而无活力花粉未着色，体积小且形状不规则（图1d）；该方法操作简单，染色时间短，效果好，有活力花粉与无活力花粉较易辨别，可作为板栗花粉活力检测方法。I2-KI染色法染色效果明显，有活力花粉颜色呈深棕色，无活力花粉不着色，但两者之间存在色浅、黄褐色的花粉，花粉活力较难判别，故此方法也不适合板栗花粉活力检测（图1e）。离体萌发培养法可以直接观察到花粉的萌发情况（图1f），可判断花粉是否具有活力，但花粉离体培养耗时较长，且培养条件较为严苛，易染菌，需要在无菌且温度恒定的环境下培养，花粉整体萌发率较低。
![]()
图 1 不同方法测定‘燕山早丰’花粉活力效果a. 清水对照；b. 0.5%TTC染色法；c. 2%TTC染色法；d. 醋酸洋红染色法；e. I2-KI染色法；f. 离体萌发法Figure 1. Effects of different methods on measuring pollen viability of ‘Yanshanzaofeng’综合反应条件、反应时间、判断标准等方面，发现醋酸洋红染色法更适合进行板栗盛花期单日花粉活力规律研究。
2.2 不同板栗品种盛花期单日花粉活力变化规律
采用醋酸洋红染色法测定不同板栗品种单日内花粉活力变化规律，结果如表1所示。各品种盛花期单日内花粉活力均保持在75%以上，但不同品种板栗花粉在单日内不同时间活力存在差异。‘短枝芽变’和‘三屯红栗’单日内花粉活力波动较小，最高值均出现在10:00—12:00。‘燕山早丰’‘王钱大栗’单日内花粉活力波动相对较大，但花粉活力最高值也出现在10:00—12:00，最低值出现时段存在差异。其中，‘燕山早丰’10:00—12:00花粉活力显著高于16:00—18:00，‘王钱大栗’10:00—12:00花粉活力显著高于其余时间花粉活力。综上，以这4个品种为父本时，应该选择上午10:00—12:00采集花粉，此时花粉活力最高。
表 1 不同板栗品种盛花期花粉活力变化Table 1. Changes in pollen viability of different chestnut cultivars at full flowering stage% 时段 ‘燕山早丰’ ‘短枝芽变’ ‘王钱大栗’ ‘三屯红栗’ 08:00—10:00 84.78ab 90.47a 83.71b 88.58a 10:00—12:00 90.71a 91.32a 87.26a 88.97a 12:00—14:00 80.65ab 88.93a 75.34c 82.23a 14:00—16:00 81.14ab 90.63a 79.49c 83.76a 16:00—18:00 78.63b 86.89a 78.26c 80.53a 注：不同小写字母表示同一品种单日内不同时间花粉活力存在显著差异（P < 0.05）。下同。 2.3 柱头可授性适宜检测方法筛选
采用联苯胺–过氧化氢法检测‘燕山早丰’柱头可授性，由于板栗柱头形状结构限制，染色效果不佳，仅有柱头尖端可观察到轻微蓝色反应；但其表面气泡反应明显，区分度较高。柱头表面只产生少量气泡，且气泡体积较小，柱头可授性较弱，标记为+，记1分（图2a-1）；柱头表面产生较多小气泡，柱头可授性较强，标记为++，记2分（图2a-2）；柱头表面迅速产生大量气泡，且均匀附着整个柱头，柱头可授性最强，标记为+++，记3分（图2a-3）。
![]()
图 2 联苯胺–过氧化氢法（a）和MTT染色法（b）检测‘燕山早丰’柱头可授性Figure 2. Detecting stigma receptivity of ‘Yanshanzaofeng’ according to benzidine hydrogen peroxide method (a) and MTT method (b)采用MTT染色法检测柱头可授性，短时间内柱头表面无明显颜色变化，避光染色1 h后，染色情况较为明显，有活力的柱头顶部变为黑色，且随着柱头可授性增强，染色长度逐渐增大。具体而言，仅有柱头尖端部分染色，染色面积微小，呈点状，此柱头可授性较弱（图2b-1）；柱头尖端面积染色面积增大，其柱头底部也部分染色，即柱头可授性较强（图2b-2）；柱头尖端染色部分延伸至柱头中部，且花柱上端也有部分染色，此柱头可授性极强（图2b-3）。
综上，联苯胺–过氧化氢法和MTT染色法都可以判断板栗柱头可授性强弱。综合考虑染色时长、可操作性等因素，在实际生产中，采用联苯胺–过氧化氢法检测板栗柱头可授性更适宜。因而，选择联苯胺–过氧化氢法进行板栗盛花期单日柱头可授性变化规律研究。
2.4 ‘燕山早丰’盛花期单日柱头可授性变化规律
表2显示：‘燕山早丰’盛花期单日内不同时段柱头均具有可授性，整体呈现出“上升—下降—上升—下降”的变化趋势。上午08:00—10:00柱头就已经具有较强可授性，10:00—12:00可授性达到最高，柱头表面具有大量气泡；随后急剧下降，12:00—14:00可授性要明显低于其他时段，之后可授性又重新升高，于14:00—16:00达到日变化第2个高峰，16:00—18:00可授性降低（图3）。
表 2 联苯胺–过氧化氢法检测‘燕山早丰’盛花期单日不同时间柱头可授性Table 2. Detecting single day stigma receptivity during full flowering period of ‘Yanshanzaofeng’ by benzidine-hydrogen peroxide method时段 柱头反应 柱头可授性 08:00—10:00 浸泡在溶液中3 ~ 5 min后，柱头表面出现大量气泡，且体积较大 +++ 10:00—12:00 浸泡溶液3 min左右，柱头表面出现大量气泡，气泡密集且体积较大 +++ 12:00—14:00 浸泡溶液5 min后，柱头表面出现少量气泡，体积较小 + 14:00—16:00 浸泡溶液5 min后，柱头表面出现气泡，气泡密集 +++ 16:00—18:00 浸泡溶液5 min后，柱头表面出现大体积气泡，但数量较少 ++ 注：+++. 柱头可授性极强；++. 柱头可授性较强；+. 柱头具可授性。 ![]()
图 3 联苯胺–过氧化氢法检测‘燕山早丰’盛花期单日柱头可授性a. 08:00—10:00；b. 10:00—12:00；c. 12:00—14:00；d. 14:00—16:00；e. 16:00—18:00Figure 3. Detecting single day stigma receptivity during blooming period of ‘Yanshanzaofeng’ by benzidine-hydrogen peroxide method2.5 盛花期单日不同授粉时间对板栗结实情况的影响
各杂交组合中，不同时间授粉落果率与坐果率均无显著差异，但结实率与空苞率在不同授粉时间存在显著差异（图4）。 ‘燕山早丰’ × ‘王钱大栗’杂交组合在盛花期08:00—10:00与10:00—12:00时段授粉，其结实率分别高达79.30%与80.36%，显著高于12:00—14:00与16:00—18:00两个授粉时段（图4）；且在10:00—12:00时段授粉，3粒坚果球苞结实率相较于其他授粉时段最高，达55.10%（表3）；空苞率在10:00—12:00授粉时段最低，仅为14.74%，而16:00—18:00授粉时段的空苞率高达37.88%，二者差异达到显著水平。
![]()
图 4 各授粉组合盛花期单日不同时段授粉结果a.‘燕山早丰’♀ × ‘王钱大栗’♂；b.‘燕山早丰’♀ × ‘短枝芽变’♂；c.‘燕山早丰’♀ × ‘三屯红栗’♂；d.‘燕山早丰’♀ × ‘燕山早丰’♂。不同小写字母表示同一授粉组合单日内不同时间授粉后结实情况存在显著差异（P < 0.05）。Figure 4. Pollination results of different combinations at varied times in a single day during blooming period表 3 各授粉组合盛花期单日不同时间授粉结实情况Table 3. Pollination and seed setting at different times of a single day during blooming period of each pollination combination授粉组合 授粉时段 结实率（1粒）/% 结实率（2粒）/% 结实率（3粒）/% 总结实率/% ‘燕山早丰’♀ × ‘王钱大栗’♂ 08:00—10:00 10.76a 14.71a 53.83a 79.30a 10:00—12:00 9.64a 15.62a 55.10a 80.36a 12:00—14:00 12.24a 17.78a 30.58b 60.60bc 14:00—16:00 5.21a 20.96a 44.24ab 70.42ab 16:00—18:00 12.35a 12.92a 29.86b 55.12c ‘燕山早丰’♀ × ‘短枝芽变’♂ 08:00—10:00 10.59a 8.07ab 40.69a 59.35a 10:00—12:00 4.98a 7.40ab 43.72a 56.10a 12:00—14:00 7.13a 3.25b 26.20ab 36.59b 14:00—16:00 5.67a 12.61a 34.56a 52.85a 16:00—18:00 7.28a 8.01ab 23.62ab 38.90ab ‘燕山早丰’♀ × ‘三屯红栗’♂ 08:00—10:00 3.25b 19.76a 53.41a 76.42a 10:00—12:00 10.60a 17.28ab 50.98a 78.86a 12:00—14:00 14.65a 17.31ab 29.02b 60.98ab 14:00—16:00 13.82a 17.32ab 39.06ab 70.20ab 16:00—18:00 11.38a 16.50b 29.02b 56.91b ‘燕山早丰’♀ × ‘燕山早丰’♂ 08:00—10:00 23.80a 4.70a 2.82a 31.32a 10:00—12:00 21.12a 4.17ab 2.50ab 27.79ab 12:00—14:00 17.63a 3.25ab 1.95ab 22.83ab 14:00—16:00 20.07a 3.96ab 2.38ab 26.40ab 16:00—18:00 16.84a 2.73b 1.64b 21.21b 注：不同小写字母表示同一授粉组合单日内不同时段授粉后结实情况存在显著差异（P < 0.05）。 在‘燕山早丰’ × ‘短枝芽变’杂交组合中，结实率在08:00—10:00授粉时段最高（59.35%），10:00—12:00时段的结实率与08:00—10:00时段相比略有下降（56.10%），并显著高于12:00—14:00授粉时段（图4）；3粒坚果球苞结实率以10:00—12:00时段最高（43.72%）（表3）；空苞率则以10:00—12:00授粉时段最低（27.64%），12:00—14:00与16:00—18:00的空苞率则高达46.34%与44.59%，均显著高于10:00—12:00（图4）。
在‘燕山早丰’ × ‘三屯红栗’杂交组合中，08:00—10:00、10:00—12:00授粉时段的结实率较高，达76.42%、78.86%，16:00—18:00结实率最低（56.91%），与上午差异达到显著水平（图4）；3粒坚果球苞结实率以08:00—10:00时段最高（53.41%），10:00—12:00时段略有下降（50.98%）（表3）；空苞率则在10:00—12:00授粉时段最低（8.94%），显著低于12:00—14:00和16:00—18:00授粉时段。
在‘燕山早丰’ × ‘燕山早丰’自交组合中，结实率均在30%左右（图4）。在08:00—10:00授粉时段，结实率最高（31.32%），空苞率最低（30.38%）；在16:00—18:00授粉时段结实率最低（21.21%）（图4）。结实的球苞大多只含有1粒坚果，2粒坚果球苞与3粒坚果球苞结实率均以08:00—10:00最高，16:00—18:00最低；空苞率则以12:00—14:00最高，达53.54%（表3）。相比其他杂交组合，自交组合中的结实率较低，落果率与空苞率较高；但是结实率最高值出现在上午授粉时段，这与其他杂交组合相同（图4）。
综合不同板栗品种授粉组合在不同时时段授粉后的结实指标可以得出，结实率和3粒坚果球苞结实率最高值均出现在上午。因此，为提高板栗杂交授粉结实率，应选择上午进行板栗杂交授粉。
3. 讨 论
3.1 花粉活力检测
前人系统研究花粉活力发现，不同植物有不同的适宜检测方法。如：花生（Arachis hypogaea）花粉活力的最适检测方法为离体萌发法与红墨水染色法[7]，薄壳山核桃（Carya illinoinensis）的花粉采用MTT染色法着色效果最佳[8]，多花黄精（Polygonatum cyrtonema）的花粉活力最适检测方法为TTC染色法[27]。TTC试剂是依靠呼吸作用产生的NADH2或NADPH2对花粉进行染色。本试验中TTC试剂无法使板栗花粉着色，这可能是由于板栗自身呼吸作用产生的脱氢酶较少，无法达到理想的染色状态。此外，利用I2-KI对花粉进行染色时，染色效果会因直链淀粉与支链淀粉的比例不同而产生差异[28]。当花粉直链淀粉含量高时，会被试剂染为蓝色，如芝麻（Sesamum indicum）[23]、油茶（Camellia oleifera）[29]；反之则为紫红色，如花生[7]、百合（Lilium lancifolium）[30]等。醋酸洋红染液则是根据花粉内物质积累量对花粉进行染色，正常花粉物质积累较多，对醋酸洋红的吸附作用与亲和力强，易被染为深红色，反正，呈淡红色或无色，如油菜（Brassica campestris）[31]。本研究中，醋酸洋红染色法、I2-KI染色法均可将‘燕山早丰’花粉成功染色，但是醋酸洋红染色法操作简单，染色快速，有活力花粉与无活力花粉易辨别，结果准确；而I2-KI染色法中存在浅色、黄褐色的介于有活力和无活力之间的花粉，活力界定具有主观性，花粉活力较难准确判别。因而，建议生产上使用醋酸洋红染色法检测‘燕山早丰’花粉活力。
3.2 柱头可授性检测
利用联苯胺–过氧化氢法检测柱头可授性，主要是依靠柱头表面有活力的过氧化氢酶，试剂中的过氧化氢在酶的作用下产生强氧化性的原子氧，将联苯胺氧化变色，并产生气泡。故当柱头表面气泡体积较大，数量较多时，柱头具有较强可授性，反之则柱头可授性较弱。此方法已广泛用于植物柱头可授性检验，如川贝母（Fritillaria cirrhosa）[32]、兜兰（Paphiopedilum）[33]等。MTT染色法检测柱头可授性主要依据柱头表面活性酶与试剂反应变色原理，可授性较强的柱头被染为深色，且染色面积较大，反之柱头颜色较浅[34−35]。本试验采用联苯胺–过氧化氢法和MTT染色法检验‘燕山早丰’的柱头可授性，其反应结果与其他植物基本相同，但MTT染色法检验板栗柱头可授性耗时较长，可能导致错过最适授粉时间。因此，在生产上可选取联苯胺–过氧化氢法快速检测柱头可授性，为杂交授粉提供指导。
3.3 花粉活力与柱头可授性日变化规律
花粉活力是决定授粉成功的重要因素，柱头是植物承接花粉的部位，二者对于准确把握植物授粉时间有重要意义，共同决定植物育种成功率[36]。不同植物的花粉活力与柱头可授性变化规律不同[23,37]，因此，植物育种授粉前，应准确掌握花粉活力和柱头的变化规律，由此得到最适授粉时间，提高育种成功率。
本研究发现，上午08:00—12:00柱头可授性与花粉活力均较强，而在正午12:00—14:00柱头可授性较弱，且此时花粉活力显著降低，这可能与空气温度、湿度变化有关。板栗柱头为针状形，在雌花成熟过程中可观察到柱头尖端有蛋白质、脂类、多糖等分泌物溢出[38]，中午温度升高，空气湿度降低，导致柱头分泌物内酶活性降低，柱头可授性降低。这与温度对甜樱桃（Prunus avium）柱头可授性的影响结果[39]相一致。此外，正午高温导致板栗花粉活力降低，这与辣椒（Capsicum annuum）的研究结果[40]相似。正午过后，气温降低，日照强度减弱，柱头内酶活性增强，柱头可授性与花粉活力回升；傍晚气温继续下降导致柱头可授性与花粉活力又重新降低。本研究得出的柱头可授性变化规律与周晶[17]对板栗柱头可授性的研究结果基本一致，且本研究中板栗花粉活力检测方法和变化规律方面的内容也弥补了板栗相关研究的缺失。
3.4 杂交授粉时间与结实情况
授粉时间显著影响授粉受精完成情况，如周建辉等[41]对厚皮甜瓜（Cucumis melo）进行不同时间杂交授粉，发现06:00—12:00授粉种子发芽率较高。对西瓜（Citrullus lanatus）[42]、西葫芦（Cucurbita pepo）[43]的研究也证实上午授粉比下午效果好，主要表现为坐果率较高，结籽数量多。此外，对玉米（Zea mays）进行不同时间授粉处理，发现不同自交系均在09:00—11:00和16:00—18:00授粉果穗的结实性状较好[44]。这与本研究得到的板栗08:00—12:00与14:00—16:00杂交授粉结实情况较好的结果相似，即最适授粉时间均避开了正午时段。除可以根据授粉结实情况筛选植物授粉时间外，也可根据柱头可授性和花粉活力筛选植物的适宜授粉时间。对青杨派树种进行柱头可授性检测，发现最佳授粉期为花序顶部露出芽鳞84 ~ 96 h[45]。本研究中，上午10:00—12:00进行授粉结实效果优良，这与前期柱头可授性与花粉活力变化规律的结果相符合，进一步验证了此时段为板栗最适杂交授粉时间。
通过将3个杂交授粉组合与自交授粉结实情况相对比，可以看出‘燕山早丰’存在明显的自交不亲和现象，其自交落果率与空苞率均较高，而结实率较低。不同杂交授粉组合间结实情况也存在差异，以‘王钱大栗’‘三屯红栗’为父本的授粉组合的结实率和3粒坚果球苞结实率与以‘短枝芽变’为父本的授粉组合存在较大差距，即存在明显的花粉直感效应，这与前人的研究结果[19]一致。后续在进行板栗杂交育种和种质优化时，宜选取性状优良的品种为父本，以丰富子代板栗优良性状。
4. 结 论
本研究首次以国审板栗品种‘燕山早丰’为试验对象，筛选其柱头可授性与花粉活力的高效便捷检测方法，并在此基础上以其为母本进行了多组合的分时段授粉，从而筛选出板栗适宜的杂交授粉时段。研究得到以下主要结论。
（1）板栗‘燕山早丰’柱头可授性最适检测方法为联苯胺–过氧化氢法，花粉活力最适检测方法为醋酸洋红染色法。
（2）盛花期单日内柱头可授性呈“上升—下降—上升—下降”趋势，且以10:00—12:00时间段柱头可授性最强，此时段不同父本板栗品种的花粉活力也最高。
（3）不同杂交组合中，总结实率与3粒坚果球苞结实率最高值均出现在08:00—12:00。
综合上述结果，实际生产中建议在上午10:00—12:00对‘燕山早丰’杂交授粉。研究结果不仅为优化板栗人工杂交授粉技术，提高板栗杂交授粉成功率，以获得充足的杂种资源提供了理论指导，同时有助于丰富板栗生殖生物学内容，揭示其繁殖机制。本研究设置的授粉试验时间梯度为2 h，在后续研究中可将时间梯度细化，以期筛选出更加精确的最适授粉时间，从而获得更理想的结实效果。
-
![]()
图 1 钱江源国家公园地形和森林实测样地分布图
Figure 1. Digital elevational model and forest survey sample plots in Qianjiangyuan National Park
![]()
图 2 冠层结构参数、海拔与地上生物量观测值的关系以及回归方程（8）对地上生物量（AGB）的拟合度
Figure 2. Relationship between canopy structural parameters, elevation and observed aboveground biomass, and the goodness offit for aboveground biomass (AGB) model using equation (8)
![]()
图 3 钱江源国家公园的冠层结构指数和地上生物量的分布
Figure 3. Canopy structural indices and aboveground biomass distribution in Qianjiangyuan National Park
表 1 航飞设计参数
Table 1 Aeronautical flight design parameters
航飞参数 Flight parameter 值 Value 飞行速度 Flight speed/（km·h−1） 250 飞行高度 Flight height/m 1 960 相对航高 Relative flight height/m 1 660 脉冲频率 Pulse repetition rate/kHz 700 视场角 Field angle of view/（°） 58.52 扫描速率 Scanning rate/lps 246 扫描幅宽 Scanning width/m 1 860 点密度 Point density/（number·m−2） > 4 旁向重叠度 Side overlap/% 20 
下载: 导出CSV
表 2 钱江源国家公园森林地上生物量模型解释变量的系数值
Table 2 Coefficients of explanatory variables in forest aboveground biomass model of Qianjiangyuan National Park
模型解释变量 Model explanatory variable 系数 Coefficient t值 t value p值 p value 1 m × 1 m网格冠层平均高度
Mean canopy height of 1 m × 1 m grid0.424 5.11 < 0.001 2 m × 2 m网格冠层平均高度
Mean canopy height of 2 m × 2 m grid−0.360 −4.09 < 0.001 叶面积垂直分布的标准差
Standard deviation of vertical foliage profile−0.308 −4.55 < 0.001 叶面积垂直分布的香农指数
Shannon index of vertical foliage profile1.274 3.85 < 0.001 冠层最大高度 Maximal canopy height 0.117 4.12 < 0.001 凹凸度 Convexity 0.034 7.53 < 0.001 海拔 Elevation 0.001 7.27 < 0.001
下载: 导出CSV
-
[1] Lohbeck M, Poorter L, Martínez-Ramos M, et al. Biomass is the main driver of changes in ecosystem process rates during tropical forest succession[J]. Ecology, 2015, 96(5): 1242−1252. doi: 10.1890/14-0472.1
[2] Poorter L, Bongers F, Aide T M, et al. Biomass resilience of Neotropical secondary forests[J]. Nature, 2016, 530: 211−214. doi: 10.1038/nature16512
[3] 马克平. 国家公园首先是自然保护基地[J]. 生物多样性, 2014, 22(4): 415−417. Ma K P. Nature conservation is the first priority for a national park[J]. Biodiversity Science, 2014, 22(4): 415−417.
[4] Coops N C, Tompalski P, Goodbody T R H, et al. Modelling lidar-derived estimates of forest attributes over space and time: a review of approaches and future trends[J]. Remote Sensing of Environment, 2021, 260: 112477. doi: 10.1016/j.rse.2021.112477
[5] Pandit S, Tsuyuki S, Dube T. Estimating above-ground biomass in sub-tropical buffer zone community forests, nepal, using sentinel 2 data[J]. Remote Sensing, 2018, 10(4): 601. doi: 10.3390/rs10040601
[6] 周俊宏, 王子芝, 廖声熙, 等. 基于GF-1影像的普达措国家公园森林地上生物量遥感估算[J]. 农业工程学报, 2021, 37(4): 216−223. doi: 10.11975/j.issn.1002-6819.2021.4.026 Zhou J H, Wang Z Z, Liao S X, et al. Remote sensing estimation of forest aboveground biomass in Potatso National Park using GF-1 images[J]. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(4): 216−223. doi: 10.11975/j.issn.1002-6819.2021.4.026
[7] Rodríguez-Veiga P, Wheeler J, Louis V, et al. Quantifying forest biomass carbon stocks from space[J]. Current Forestry Reports, 2017, 3(1): 1−18. doi: 10.1007/s40725-017-0052-5
[8] Jucker T, Caspersen J, Chave J, et al. Allometric equations for integrating remote sensing imagery into forest monitoring programmes[J]. Global Change Biology, 2017, 23(1): 177−190. doi: 10.1111/gcb.13388
[9] Fischer R, Knapp N, Bohn F, et al. The relevance of forest structure for biomass and productivity in temperate forests: new perspectives for remote sensing[J]. Surveys in Geophysics, 2019, 40(4): 709−734. doi: 10.1007/s10712-019-09519-x
[10] Coomes D A, Dalponte M, Jucker T, et al. Area-based vs tree-centric approaches to mapping forest carbon in Southeast Asian forests from airborne laser scanning data[J]. Remote Sensing of Environment, 2017, 194: 77−88. doi: 10.1016/j.rse.2017.03.017
[11] 于明坚, 胡正华, 余建平, 等. 浙江古田山自然保护区森林植被类型[J]. 浙江大学学报(农业与生命科学版), 2001, 27(4): 23−28. Yu M J, Hu Z H, Yu J P, et al. Forest vegetation types in Gutianshan Natural Reserve in Zhejiang[J]. Journal of Zhejiang University (Agriculture & Life Science), 2001, 27(4): 23−28.
[12] 张田田, 王璇, 任海保, 等. 浙江古田山次生与老龄常绿阔叶林群落特征的比较[J]. 生物多样性, 2019, 27(10): 1069−1080. doi: 10.17520/biods.2019059 Zhang T T, Wang X, Ren H B, et al. A comparative study on the community characteristics of secondary and old-growth evergreen broad-leaved forests in Gutianshan, Zhejiang Province[J]. Biodiversity Science, 2019, 27(10): 1069−1080. doi: 10.17520/biods.2019059
[13] Guo Q, Su Y, Hu T, et al. An integrated UAV-borne lidar system for 3D habitat mapping in three forest ecosystems across China[J]. International Journal of Remote Sensing, 2017, 38(8−10): 2954−2972.
[14] Harms K E, Condit R, Hubbell S P, et al. Habitat associations of trees and shrubs in a 50-ha neotropical forest plot[J]. Journal of Ecology, 2001, 89: 947−959. doi: 10.1111/j.1365-2745.2001.00615.x
[15] Yamakura T, Kanzaki A, Itoh T, et al. Topography of a large-scale research plot established within a tropical rain forest at Lambir, Sarawak[J]. Tropics, 1995, 5: 41−56. doi: 10.3759/tropics.5.41
[16] Getzin S, Wiegand K, Schöning I. Assessing biodiversity in forests using very high-resolution images and unmanned aerial vehicles[J]. Methods in Ecology and Evolution, 2012, 3(2): 397−404. doi: 10.1111/j.2041-210X.2011.00158.x
[17] Knapp N, Fischer R, Cazcarra-Bes V, et al. Structure metrics to generalize biomass estimation from lidar across forest types from different continents[J]. Remote Sensing of Environment, 2020, 237: 111597. doi: 10.1016/j.rse.2019.111597
[18] Roussel J R, Auty D, Coops N C, et al. LidR: an R package for analysis of Airborne Laser Scanning (ALS) data[J]. Remote Sensing of Environment, 2020, 251: 112061. doi: 10.1016/j.rse.2020.112061
[19] Hijmans R J, Bivand R, Forner K, et al. Terra: spatial data analysis. R package version 1.6-7[Z/OL]. [2022−07−21]. https://CRAN.R-project.org/package=terra.
[20] R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, 2022, Vienna, Austria[Z/OL]. [2022−07−21]. https://www.R-project.org/.
[21] Gough C M, Atkins J W, Fahey R T, et al. High rates of primary production in structurally complex forests[J]. Ecology, 2019, 100(10): e02864.
[22] Ali A. Forest stand structure and functioning: current knowledge and future challenges[J]. Ecological Indicators, 2019, 98: 665−677. doi: 10.1016/j.ecolind.2018.11.017
[23] 王晓莉, 常禹, 陈宏伟, 等. 黑龙江省大兴安岭森林生物量空间格局及其影响因素[J]. 应用生态学报, 2014, 25(4): 974−982. doi: 10.13287/j.1001-9332.2014.0082 Wang X L, Chang Y, Chen H W, et al. Spatial pattern of forest biomass and its influencing factors in the Great Xing’an Mountains, Heilongjiang Province, China[J]. Chinese Journal of Applied Ecology, 2014, 25(4): 974−982. doi: 10.13287/j.1001-9332.2014.0082
[24] Lin D, Lai J, Yang B, et al. Forest biomass recovery after different anthropogenic disturbances: relative importance of changes in stand structure and wood density[J]. European Journal of Forest Research, 2015, 134: 769−780. doi: 10.1007/s10342-015-0888-9
[25] 余建平, 伊晓霞, 余顺海, 等. 基于景观格局指数的钱江源国家公园生态系统完整性评价分析[J]. 浙江林业科技, 2020, 40(4): 30−36. doi: 10.3969/j.issn.1001-3776.2020.04.005 Yu J P, Yi X X, Yu S H, et al. Analysis on ecological integrity of Qianjiangyuan National Park based on landscape pattern index[J]. Zhejiang Forest Science and Technology, 2020, 40(4): 30−36. doi: 10.3969/j.issn.1001-3776.2020.04.005
计量
- 文章访问数: 756
- HTML全文浏览量: 384
- PDF下载量: 161
