在重庆市北碚区西南大学柑桔研究所沃柑园(29°45′ N，106°22′ E)于2022-2024年连续开展两年重复试验。果园土壤为碱性(pH值8.2)，有机质17.9 g/kg，碱解氮、有效磷、速效钾分别为60、25、168 mg/kg。试验以枳橙为基砧、福本为中间砧高接8年生的沃柑为试材。供试肥料为含N 47%的尿素、含P₂O₅ 12%的过磷酸钙、含K₂O 51%的硫酸钾、复合肥(N-P₂O₅-K₂O在春、夏、秋三季施用比例分别为21∶8∶11、15∶5∶25、15∶15∶15)，土壤改良剂为纯度70.0%~85.0%脱硫石膏和99.0%硫磺粉，生产厂家分别为河南省巩义市元亨净水材料厂和山西省牧颜宸科技有限公司。根据李玉波等^[28]的研究结果，确定脱硫石膏和硫磺粉适宜施用量为15 000 kg/hm²。烟气脱硫石膏样品重金属含量满足《土壤环境质量农用地土壤污染风险管控标准(试行)》(GB 15618—2018)中pH值≥7.53的二级标准要求。

试验分别设置复合肥(CK)、配方肥(T1)、配方肥配施硫磺粉(T2)、配方肥配施脱硫石膏(T3)、复合肥配施硫磺粉(T4)、复合肥配施脱硫石膏(T5)等6个处理，具体改良施肥材料及其施用量如表 1所示，各处理N、P₂O₅、K₂O养分施用量一致，年施肥量分别为0.42 kg/株、0.22 kg/株、0.44 kg/株，每处理6株。每年于3月、7月、10月进行施肥，T1、T2、T3处理3次施肥比例均为2∶5∶3，CK、T4、T5处理3次施肥比例均为8∶10∶7。于树体行间两侧滴水线外侧位置各挖一条长度为100 cm、宽度为20 cm、深度为30 cm的施肥沟，将准备好的肥料、土壤改良材料与从沟内挖出的土壤充分搅拌混合均匀后填回施肥沟内，且回填土比周边土面高10 cm左右。果园植保、树体修剪等其他田间管理按常规技术统一进行。

分别于2022年8月至2023年1月以及2023年8月至2024年1月，即从果实膨大期开始，每月中旬在每株树上沿树冠中部东、南、西、北4个方向各采4个果作为1个混合样，单株为1次重复，每处理共6个重复样，用于测定果实内外品质指标。分别于2023年1月和2024年1月果实成熟期按单株全部采收称重测产量，2024年1月对每株树采收的所有果实进行分级。分级标准为：果实横径55.0~60.0 mm为小果型(S)，60.0~65.0 mm为中果型(M)，65.0~70.0 mm为大果型(L)，70.0~80.0 mm为超大果型(XL)^[29]。

果实经去离子水清洗、擦干后，以游标卡尺测量果皮厚度，采用GY-4果实硬度计测定果皮硬度。榨汁后，使用PAL-1数显糖度仪测定可溶性固形物(TSS)含量，采用NaOH中和滴定法、2，6-二氯靛酚钠滴定法测定可滴定酸(TA)含量和维生素C(V_C)含量，比值法计算固酸比(TSS/TA)。

采用Microsoft Excel 2019统计试验数据，利用R语言进行方差分析，使用Origin Pro 2025软件绘图，并运用配套插件进行主成分分析和综合评价。

由图 1分析可知，2024年、2023年各处理沃柑单株平均产量分别为51.42 kg、34.51 kg，前者比后者提高了49.0%；相较于CK处理，各碱性土改良施肥处理的沃柑产量均有所提升。2023年各碱性土改良施肥处理的沃柑单株平均产量较CK增产了7.8%~60.2%，其中T5处理的产量最高，与CK处理差异显著；2024年各碱性土改良施肥处理的沃柑单株平均产量较CK增产了3.5%~60.2%，以T3最高，与CK处理差异显著。

从8月至翌年1月沃柑果皮厚度检测结果表明(图 2)，两年试验的动态变化规律一致，即随着时间的推移，呈先增加后降低再增加最后缓慢降低的趋势。2023年1月采收时，T3处理的果皮厚度最薄，为2.51 mm；T1处理最厚，为2.89 mm；但各处理差异均不显著。2024年1月采收时，T3处理果皮厚度最薄，为2.15 mm；T1处理最厚，为2.88 mm；T3处理的果皮厚度比T1降低了25.3%，且两个处理间差异显著。

从9月至翌年1月沃柑果皮硬度监测结果表明(图 3)，2023年9月至10月果皮硬度迅速降低，10月至12月缓慢降低至平稳，至2023年1月采收时又略有降低；2023年1月果实成熟采收时各处理果皮硬度为16.63~19.44 N，其中T4处理最小，T1处理最大，T4比T1处理显著降低了14.5%。2023年9月至11月果皮硬度快速降低，11月至12月缓慢降低，至2024年1月采收时又略有升高；2024年1月果实成熟采收时各处理果皮硬度为19.21~22.69 N，以T2处理相对最小，CK处理最大，T2比CK处理显著降低了15.4%。

由表 2可知，2024年各处理大果型(L级和XL级)级别分布频率以T3、T5、T4处理相对较高，均显著高于CK处理；各处理M级中型果分布频率以T2、T3、T4、T5处理较高，显著高于T1处理；各处理S级小型果分布频率以CK、T1处理显著高于其他处理。

从9月至翌年1月沃柑果实TSS含量检测结果表明(图 4)，2022年9月至11月果实TSS含量缓慢升高，11月至12月快速升高，12月至翌年1月缓慢升高至最大值。2023年1月采收时，T3处理的果实TSS含量最大，为12.98%；T1处理最小，为11.70%；T3比T1处理显著提高了11.0%。2024年9月至11月果实TSS含量缓慢升高，11月至翌年1月快速升高至最大值。2023年12月各处理的果实TSS含量，以T3处理相对最高，为12.27%；CK和T1处理最小，为11.23%；T3处理比CK、T1处理均显著增加了9.2%。2024年1月采收时，T2处理的TSS含量最大，比CK处理提高了2.7%，但差异未达显著水平。

从9月至翌年1月沃柑果实TA含量检测结果表明(图 5)，两年试验的动态变化规律一致，即9月至11月果实TA含量快速下降，11月至翌年1月缓慢下降，至采收时达相对最低。2023年1月采收时，T3和T4处理的果实TA含量最小，为0.51%；T1处理最大，为0.54%；T3和T4分别比T1处理降低了5.0%、5.0%，但各处理间无显著性差异。2023年12月，T5处理的果实TA含量相对最小，为0.56%；CK处理最大，为0.62%；T5比CK处理显著降低了10.1%。2024年1月采收时，T3处理的果实TA含量最小，为0.41%，T3比CK处理降低了11.8%，但各处理间无显著差异。

从9月至翌年1月沃柑果实TSS/TA检测结果表明(图 6)，2022年9月至10月沃柑果实TSS/TA缓慢升高，10月至11月快速升高，11月至翌年1月缓慢升高至最大值。2023年1月采收时，T3处理的TSS/TA最大，为25.42；T1处理最小，为21.87；T3比TI显著提高了16.2%。2023年9月至10月沃柑果实TSS/TA缓慢升高，10月至翌年1月快速升高至最大值。2023年12月，T3处理的果实TSS/TA最大，为21.02；CK处理最小，为18.20；T3处理比CK显著增加了15.5%。2024年1月采收时，T3处理的TSS/TA最大，为30.94，但各处理间差异不显著。

从9月至翌年1月沃柑果实V_C含量检测结果表明(图 7)，2022年9月至10月沃柑果实V_C含量缓慢升高，10月至12月快速升高至最高值，至2023年1月采收时有所降低。2023年1月采收时，T3处理的果实V_C含量最大，为227.6 mg/L；CK处理最小，为205.3 mg/L；T3比CK提高10.9%，但不存在显著性差异。2023年9月至10月沃柑果实V_C含量迅速升高，10月至12月缓慢升高至最高值，至2024年1月采收时有所下降。2024年1月采收时，T4处理的果实V_C含量最大，为236.6 mg/L；T1处理最小，为221.8 mg/L；T4与T1处理差异显著。

将2023年和2024年1月中旬果实成熟采收期各处理指标的数据取平均值，通过多项式拟合，探究果皮硬度、果皮厚度分别与果实TSS、TA、V_C含量之间的相关性。数据分析结果表明，果皮硬度与果实TA含量、果实V_C含量之间存在显著的相关性，呈现出较高的曲线拟合度(R²＞0.8)。从图 8a可知，果皮硬度与TA含量呈正相关，其变化趋势在T3处理时出现转折。从图 8b可见，果皮硬度与V_C含量呈负相关，于T2处理时发生明显转折。由上可见，当果实变软时，TA下降，V_C升高。各碱性土改良施肥处理均能有效提高沃柑果实外观内质。

各碱性土改良施肥处理产量及品质指标两年平均值的主成分分析载荷图表明(图 9)，主成分1和主成分2贡献率分别为58.9%、16.5%，累计贡献率达75.4%，主成分1中起决定作用的指标有产量、果皮厚度、果皮硬度、TSS/TA，主成分2主要由TSS、TA、V_C含量决定。

由表 4所示，选取特征值大于1的两个主成分，对不同碱性土改良施肥处理进行综合评价。不同碱性土改良施肥处理的综合得分，由大到小的顺序为T3、T2、T5、T4、CK、T1。

碱性土改良施肥在保证氮、磷、钾等养分的供应下，配施脱硫石膏降低了土壤含盐量，配施硫磺粉降低了土壤的pH值，为作物创造了适宜其生长的土壤环境，有利于作物的生长发育。2023年和2024年各碱性土改良施肥处理的沃柑产量分别较CK增加了7.8%~60.2%和3.5%~60.2%，且分别以T5和T3处理最高，这与前人^{[14-15, 30-31]}研究结果一致，可能与配施脱硫石膏置换土壤中的Na⁺，降低土壤含盐量，提高根系活力，有效延长叶片功能期，增加干物质积累量，促进果实生长发育有关。

果实中的钙以果胶酸钙的形式参与细胞壁和胞间层的构成，从而使果实器官或组织具有一定的机械强度，进而提高了果皮硬度^[32]。本试验结果表明，果实成熟采收时各碱性土改良施肥处理的果皮厚度及果皮硬度比单施化肥均有一定程度的改善，其中2024年度果皮厚度以T3处理最薄，以T1处理最厚，T3比T1处理显著降低了25.3%，这可能与脱硫石膏与土壤中的碳酸盐反应生成中性盐和碳酸钙，降低了土壤pH值，提高了土壤有效磷含量，吸收利用的磷促进果肉生长、降低果皮厚度有关。果皮硬度以T2处理相对最小，比CK处理显著降低了15.4%，这可能与硫磺粉氧化生成硫酸氧化物，置换出一定量的氢离子，与土壤中的碳酸钙发生反应生成可溶性钙盐，易随雨水或灌溉水淋失，造成了土壤中钙的流失，间接影响果实中的钙含量，从而降低果皮硬度有关^[33]。果实大果型(L级和XL级)比例变化范围为7.00%~42.00%，其中以T3处理的比例最高，与CK相比显著提升，可能与脱硫石膏降低了土壤pH值，提高了土壤有效氮、磷、钾等含量和促进了养分的吸收利用，有利于果实的膨大有关。盐碱地改良试验结果表明，施用脱硫石膏可提高葡萄的叶面积和叶绿素含量，有利于叶片的光合作用，从而促进了果实中光合产物的积累，最终提高了果实品质^[34]。连续两年的试验结果表明，果实成熟采收时各碱性土改良施肥处理的沃柑果实内在品质与单施化肥相比均有一定程度的提升，其中2023年和2024年果实TSS含量均以T3处理最高，比T1处理显著增加了11.0%和9.2%，可能与脱硫石膏通过置换Na⁺，降低了土壤pH值，缓解了碱性土对沃柑的胁迫，提高了土壤中磷、铁、锌等元素的有效性，有利于沃柑根系对土壤矿质元素的吸收利用，促进了叶片光合产物的合成与转运，进而提高了沃柑果实TSS含量有关。T3处理的沃柑果实TA含量分别比T1、CK降低了5.0%、11.8%，可能与脱硫石膏提供了一定量的钙，钙通过稳定细胞膜的结构，减少柠檬酸、苹果酸等有机酸在液泡中的积累，从而降低了沃柑果实TA含量有关。V_C含量分别以T3和T4处理相对最大，T4比T1处理显著增加了6.7%，可能与硫磺粉提供了硫，可提高叶片光合效率，增加糖分供应，而糖是V_C合成的碳骨架来源，进而提高V_C含量有关，这与前人^{[16-17, 34]}的研究结果一致。

可见，施用酸性肥料及脱硫石膏、硫磺粉等碱性土改良剂，不仅能改善果实外观品质，还显著增加了沃柑产量，提高了果实TSS、V_C含量，以及降低了果实TA含量。对沃柑成熟期果实外观内质指标进行曲线拟合，分析发现果皮硬度与果实TA、V_C含量之间有较好的相关性，当硬度下降时，TA下降，V_C升高，分别以T3和T2处理为最佳。借助果皮硬度与果实内在品质指标的良好相关性量化模型，结合果皮硬度快速监测技术，甚至可以通过触摸果皮“手感”硬度的经验性判断，今后可通过果皮硬度的“触感”来快速、实时、无损估测果实内在品质与成熟度情况。利用主成分分析对成熟期沃柑产量及外观内质指标进行综合评价，分析不同碱性土改良施肥处理对果实品质形成的影响，以T3处理即配方肥(尿素+硫酸钾+过磷酸钙)配施脱硫石膏的综合得分最高。

本试验条件下，与常规单施复合肥相比，酸性肥料(过磷酸钙和硫酸钾)和土壤改良剂(硫磺粉和脱硫石膏)配施，有利于果园碱性土壤的改良培肥，进而提高了沃柑产量，改善了果实外观内质。基于主成分分析的沃柑产量、品质评价得出，T3处理即配方肥(尿素+硫酸钾+过磷酸钙)配施脱硫石膏的综合得分最高，亦即西南地区石灰性紫色土柑橘园通过脱硫石膏结合施用酸性磷、钾肥的综合改良效果为最佳。