基于碳、氮同位素和礦物元素的青蘿卜產(chǎn)地判別

姜浩琛，王世成，李波，郝淋淋，孫辭，崔杰華，李國琛*，徐志祥

1(山東農(nóng)業(yè)大學(xué) 食品科學(xué)與工程學(xué)院，山東 泰安，271018)2(中國科學(xué)院沈陽應(yīng)用生態(tài)研究所，遼寧 沈陽，110016)

蘿卜一詞最早在中國古代專著《爾雅》被提及，稱其為“蘆菔”。青蘿卜是十字花科蘿卜屬一種常見的根莖類蔬菜，種植歷史悠久，可生食、做湯、干腌和鹽漬等，不僅富含人體所需的營養(yǎng)物質(zhì)，而且具有一定藥用價(jià)值[1-3]。但是，由于不同地理來源的青蘿卜維生素C含量、外形、口感等存在差異，例如天津的沙窩蘿卜通常以水果蘿卜出名[4]，消化功能顯著，而山東濰坊的浮橋蘿卜鈣含量豐富。這類特色農(nóng)產(chǎn)品因其特殊的文化和自然因素、較高的口碑和經(jīng)濟(jì)價(jià)值，易遭遇摻假和假冒。因此，研究建立精準(zhǔn)的青蘿卜產(chǎn)地溯源技術(shù)體系，對(duì)判別青蘿卜真實(shí)來源、維護(hù)消費(fèi)者權(quán)益、規(guī)范市場(chǎng)秩序尤為重要。

第三，在教學(xué)組織形式上應(yīng)充分考慮到學(xué)生的特點(diǎn)來進(jìn)行改進(jìn)。英語視聽說是一門不同于傳統(tǒng)聽力課的課程，因此在教學(xué)組織形式上也應(yīng)該進(jìn)行一定的創(chuàng)新，教師應(yīng)創(chuàng)造性地運(yùn)用個(gè)人展示、小組討論、辯論等形式。

目前在農(nóng)產(chǎn)品地理溯源方面常用到的產(chǎn)地判別技術(shù)包括礦物元素分析技術(shù)、穩(wěn)定同位素質(zhì)譜技術(shù)、近紅外光譜技術(shù)及電子鼻技術(shù)等[5-7]。生物在生長(zhǎng)的過程中受氣候、海拔和土壤性質(zhì)等環(huán)境因素的影響，不同產(chǎn)地的農(nóng)產(chǎn)品體內(nèi)穩(wěn)定同位素比值和礦物元素含量呈現(xiàn)出獨(dú)特性和差異性，具有明顯的地域特征[8-9]，因此，礦物元素分析技術(shù)、穩(wěn)定同位素質(zhì)譜技術(shù)在農(nóng)產(chǎn)品產(chǎn)地溯源中具備較高的應(yīng)用價(jià)值，并且這兩項(xiàng)技術(shù)具有前處理操作簡(jiǎn)單、判別模型精準(zhǔn)度高等優(yōu)點(diǎn)，在眾多特色地標(biāo)農(nóng)產(chǎn)品中被廣泛應(yīng)用[10]。連思雨等[11]和BERTOLDI等[12]通過測(cè)定枸杞中的穩(wěn)定同位素、礦物元素等成分，結(jié)合主成分分析(principal components analysis，PCA)、偏最小二程判別分析(partial least squares discrimination analysis，PLS-DA)等化學(xué)計(jì)量方法建立判別模型，準(zhǔn)確率高達(dá)95%以上；MA等[13]和王潔[14]基于礦物元素、礦物元素與穩(wěn)定同位素比率結(jié)合等方式對(duì)扁形茶、龍井茶的產(chǎn)地溯源技術(shù)進(jìn)行探究，建立高精準(zhǔn)度的鑒別模型，進(jìn)一步肯定這2項(xiàng)技術(shù)在產(chǎn)地判別方面的作用；此外，礦物元素分析技術(shù)、穩(wěn)定同位素質(zhì)譜技術(shù)等在谷物[15-17]、果蔬[18-19]以及藥材[20-21]、水產(chǎn)品[22-24]等眾多領(lǐng)域有著廣泛且高效的應(yīng)用，但在蘿卜地理來源的鑒別方面鮮有應(yīng)用研究。

本研究結(jié)合當(dāng)前產(chǎn)地研究的現(xiàn)狀，選取山東、天津、河南3個(gè)省份具有地理標(biāo)志認(rèn)證的10個(gè)品種青蘿卜，通過測(cè)定穩(wěn)定同位素比值和礦物元素含量，開展青蘿卜地理信息的相關(guān)性研究與應(yīng)用，并結(jié)合PLS-DA、線性判別分析(linear discriminant analysis，LDA)等化學(xué)計(jì)量學(xué)方法建立判別模型。旨在服務(wù)“三品一標(biāo)”，為保護(hù)青蘿卜產(chǎn)品品種、品質(zhì)、品牌以及提高經(jīng)濟(jì)效益提供有效方法。

1 材料與方法

1.1 材料與試劑

采集地理來源真實(shí)可靠的青蘿卜樣品177份，其中天津57份，河南59份、山東61份，具體信息如表1所示。穩(wěn)定同位素標(biāo)準(zhǔn)物質(zhì)USGS43(δ13C=-21.34‰，δ15N=+8.44‰)，美國地質(zhì)勘探局Reston穩(wěn)定同位素實(shí)驗(yàn)室；元素國家標(biāo)準(zhǔn)物質(zhì)菠菜GSBW10015，地球物理地球化學(xué)勘查研究所；元素標(biāo)液為電感耦合等離子體分析用多元素標(biāo)準(zhǔn)溶液GNM-M311129-2013，國家有色金屬及電子材料分析測(cè)試中心逐級(jí)稀釋配得；濃硝酸(65%)，美國默克制藥。

1.2 儀器與設(shè)備

Flash 2000元素分析儀、Delta V Advantage穩(wěn)定同位素比值質(zhì)譜儀，美國賽默飛世爾公司；NEXION 300X電感耦合等離子體質(zhì)譜儀，美國珀金埃爾默股份有限公司；破壁機(jī)E310，美國Vitamix公司。

1.3 實(shí)驗(yàn)方法

1.3.1 樣品前處理

青蘿卜樣品采集后，在實(shí)驗(yàn)室用自來水將泥土清洗、陰干后，測(cè)定單個(gè)重量、長(zhǎng)度等基礎(chǔ)數(shù)據(jù)。四分法分取蘿卜樣品并切塊，用破壁機(jī)將樣品勻漿粉碎，分取100 g勻漿后的樣品放入冷凍干燥機(jī)進(jìn)行冷凍干燥(-50 ℃)12 h以上，干燥后樣品放入-80 ℃超低溫冰箱保存?zhèn)溆谩?/p>

1.3.3 元素測(cè)定

1.3.2 穩(wěn)定同位素測(cè)定

應(yīng)用Origin 8.0軟件對(duì)碳、氮穩(wěn)定同位素進(jìn)行箱線圖分析，并應(yīng)用SPSS 25.0進(jìn)行單因素ANOVA檢驗(yàn)，最后分別以所測(cè)青蘿卜中18種礦物元素為變量，應(yīng)用SIMCA-P 13.0分析軟件進(jìn)行PLS-DA建模，應(yīng)用MPP軟件結(jié)合LDA算法建立青蘿卜產(chǎn)地判別模型。

儀器工作參數(shù)：氧化管溫度為1 150 ℃，還原管溫度為850 ℃，通氧時(shí)間為70 s，元素分析儀氦氣吹掃流量為230 mL/min，進(jìn)入質(zhì)譜儀氦載氣流量為100 mL/min，高純N2、CO2(純度>99.999%)為參考?xì)狻?/p>

將數(shù)據(jù)從純文本數(shù)據(jù)轉(zhuǎn)換成數(shù)據(jù)庫格式數(shù)據(jù)。通過對(duì)文本數(shù)據(jù)進(jìn)行識(shí)別、讀取和存儲(chǔ)，實(shí)現(xiàn)對(duì)純文本文檔的解析和轉(zhuǎn)換。

用電子天平準(zhǔn)確稱量0.3 g各產(chǎn)地青蘿卜凍干粉末樣品于消解罐中，每個(gè)樣品3個(gè)平行。向消解罐中加入10 mL濃硝酸，過夜，第2天不蓋蓋子在智能石墨消解爐中進(jìn)行120 ℃石墨消解至澄清，約4 h，冷卻至室溫后從微波消解儀中取出，用少量去離子水分3次洗滌消解管，合并洗滌液及消解液，并定容至25 mL，進(jìn)行上機(jī)檢測(cè)。

儀器工作參數(shù)：電感耦合等離子體質(zhì)譜儀(inductively coupled plasma mass spectrometry，ICP-MS)每次測(cè)定前均進(jìn)行自動(dòng)調(diào)諧，根據(jù)儀器的狀態(tài)對(duì)參數(shù)進(jìn)行調(diào)整，以滿足靈敏度、氧化物、穩(wěn)定性等各項(xiàng)指標(biāo)。儀器調(diào)諧完畢后，將內(nèi)標(biāo)進(jìn)樣管插入1 μg/mL 內(nèi)標(biāo)溶液中進(jìn)行矯正，樣品管依次插入各濃度的標(biāo)準(zhǔn)溶液進(jìn)行制作標(biāo)曲，最后測(cè)定樣品。射頻功率1 230 W；載氣流速1.21 L/min；Babington霧化器，霧化室溫度2 ℃；蠕動(dòng)泵轉(zhuǎn)速0.1 r/s。用內(nèi)標(biāo)法進(jìn)行校正。儀器調(diào)諧完畢后，將內(nèi)標(biāo)進(jìn)樣管插入1 μg/mL內(nèi)標(biāo)溶液中(分析過程中始終進(jìn)內(nèi)標(biāo)溶液)，樣品管依次插入各濃度的標(biāo)準(zhǔn)溶液(濃度依次遞增)和待測(cè)樣中進(jìn)行測(cè)定。

1.4 多元統(tǒng)計(jì)分析

碳、氮同位素比值的測(cè)定：取0.2 mg青蘿卜粉末樣品于錫杯(5 mm×8 mm)中壓實(shí)包好后直接放入元素分析儀自動(dòng)進(jìn)樣器中，在高溫下燃燒、還原，最終生成氣體，按N2、CO2的順序先后進(jìn)入同位素比率質(zhì)譜儀進(jìn)行檢測(cè)。進(jìn)樣序列中每隔10個(gè)樣品插入一個(gè)標(biāo)物進(jìn)行質(zhì)量控制，樣品穩(wěn)定同位素比值應(yīng)用標(biāo)物進(jìn)行單點(diǎn)校正。儀器標(biāo)準(zhǔn)偏差小于0.2‰。

2 結(jié)果與分析

2.1 不同產(chǎn)區(qū)青蘿卜多種礦質(zhì)元素含量和穩(wěn)定同位素比值分析

化學(xué)計(jì)量學(xué)方法結(jié)合的方式建立青蘿卜樣品的判別模型。目前，在產(chǎn)地溯源的研究中常用到的化學(xué)計(jì)量學(xué)方法主要包括PCA、PLS-DA、支持向量機(jī)(support Vector Machine，SVM)、人工神經(jīng)網(wǎng)絡(luò)(back propagation-artificial neural network，BP-ANN)以及LDA等。本文對(duì)青蘿卜樣品中的礦物元素含量進(jìn)行分析，結(jié)合PCA-X、PLS-DA、LDA 3種算法建立產(chǎn)地判別模型。

通過三年多在醫(yī)院宣傳管理工作中應(yīng)用PDCA循環(huán)，如對(duì)比圖表結(jié)果顯示：新管理制度體系和方法開展后，全院通訊員稿件數(shù)量和質(zhì)量持續(xù)提升，醫(yī)院自媒體、網(wǎng)站稿件數(shù)量和質(zhì)量快速增加，見圖1。全院職工對(duì)宣傳工作參與的積極性和主動(dòng)性明顯增強(qiáng)，醫(yī)院宣傳隊(duì)伍體系完善，后備力量充足，形成良好的宣傳氛圍。與此同時(shí)，院外媒體期刊發(fā)稿數(shù)量和省級(jí)以上媒體刊發(fā)稿件數(shù)量也有所增加，說明醫(yī)院宣傳影響力和影響面逐漸擴(kuò)大。醫(yī)院標(biāo)識(shí)廣告等真正達(dá)到統(tǒng)一、溫馨、便捷且有醫(yī)院文化內(nèi)涵。由此可見，在醫(yī)院宣傳管理中應(yīng)用PDCA循環(huán)可將醫(yī)院整體形象顯著提升，宣傳引導(dǎo)能力和宣傳影響力也會(huì)明顯提高。

A-不同省份δ13C值范圍；B-不同省份δ15N值范圍；C-不同亞地區(qū)δ13C值范圍；D-不同亞地區(qū)δ15N值范圍圖1 不同地區(qū)青蘿卜樣品C、N同位素比值分析Fig.1 Analysis of C and N isotope ratios of green turnip samples from different regions

其次，應(yīng)用SPSS對(duì)所測(cè)定的18種礦物元素進(jìn)行單因素ANOVA檢驗(yàn)，由表2可知，天津、河南、山東3個(gè)省份的整體元素含量分為3個(gè)梯度，Na、Ca、P、Mg 4種元素含量較高，平均含量均超過1 000 mg/kg，且Na、Ca、Mg 3種元素在不同省份間均存在顯著差異，天津地區(qū)青蘿卜樣品中Na含量最高，達(dá)10 786 mg/kg，河南地區(qū)Ca含量最高為6 182.05 mg/kg，山東地區(qū)Mg含量最高為2 237.28 mg/kg；Fe、Sr、Zn、B、Ba、Mn、Cu相對(duì)較低，但整體元素含量均高于1 mg/kg，F(xiàn)e、Sr、Zn、B、Mn 5種元素在天津、河南、山東3個(gè)省份間均存在明顯差異；Mo、Cr、Ni、Se、Cd、As、Co 7種元素含量最低，平均含量小于1 mg/kg，其中只有As在不同省份間差異顯著(P<0.05)。

首先通過箱線圖對(duì)不同地區(qū)的穩(wěn)定同位素值進(jìn)行分析，如圖1所示。從整體分析，天津、河南、山東3個(gè)省份的青蘿卜樣品中的δ15N值為0.37‰～2.37‰，且均有重合，差異較小，天津和山東兩省樣本δ15N值范圍較廣，河南樣本的δ15N值范圍較為集中；山東、天津、河南3個(gè)省份的青蘿卜樣品δ13C值分布也顯示整體差異不顯著，范圍為-29.68‰～-27.88‰；從不同省份不同地區(qū)樣品細(xì)分發(fā)現(xiàn)，天津省的武清、西青兩個(gè)地區(qū)蘿卜的δ15N值存在差異，河南省固始縣與上蔡縣、固始縣與郟縣間的δ15N值存在差異，山東壽光和濰城、滕州兩個(gè)地區(qū)δ15N值存在差異；但不同地區(qū)樣品間δ13C值并未存在顯著差異。這表示碳、氮同位素在判別不同省份的青蘿卜時(shí)判別貢獻(xiàn)率較低，但δ15N值對(duì)進(jìn)一步精確鑒別小產(chǎn)區(qū)青蘿卜的地理來源時(shí)可能會(huì)具有相對(duì)顯著貢獻(xiàn)，這與當(dāng)?shù)氐耐寥罈l件、施肥狀況以及農(nóng)藥使用情況有關(guān)，可為今后進(jìn)一步研究青蘿卜的產(chǎn)地提供依據(jù)。

2.2 判別分析

植物中穩(wěn)定同位素比值和礦物元素含量的差異反映了其生長(zhǎng)的環(huán)境及氣候條件，為進(jìn)行產(chǎn)地溯源、建立模型提供了依據(jù)。由圖2分析結(jié)果可知，不同地區(qū)青蘿卜中的碳、氮同位素整體差異范圍小，因此，本研究在進(jìn)行農(nóng)產(chǎn)品產(chǎn)地溯源時(shí)，采用礦物元素含量和

2.2.1 PLS-DA

基于梯形機(jī)構(gòu)自身結(jié)構(gòu)特點(diǎn)存在的不足，需從本質(zhì)上消除轉(zhuǎn)向時(shí)車輪側(cè)滑現(xiàn)象的產(chǎn)生，才可提高架車的可操作性。故針對(duì)設(shè)備自身結(jié)構(gòu)特點(diǎn)，對(duì)轉(zhuǎn)向機(jī)構(gòu)進(jìn)行改進(jìn)，其結(jié)構(gòu)示意如圖6所示。圖中的回轉(zhuǎn)支承用來替換原梯形機(jī)構(gòu)以實(shí)現(xiàn)整車轉(zhuǎn)向，與原機(jī)構(gòu)相比其本質(zhì)區(qū)別在于兩前輪并非與架體固連，在架車轉(zhuǎn)向過程中，前后車輪的軸線均交于一點(diǎn)即轉(zhuǎn)向中心處，滿足無側(cè)滑理論轉(zhuǎn)向條件。

表2 不同地區(qū)青蘿卜元素含量分析 單位：mg/kg

本研究共測(cè)定天津、河南、山東3個(gè)省份177份青蘿卜樣品中的碳、氮同位素比值和18種礦物元素的含量。

圖2 三個(gè)產(chǎn)地青蘿卜VIP值Fig.2 The results of VIP in green turnip from different producing areas

音樂是成功創(chuàng)編的基礎(chǔ)，好的音樂可烘托操的氣氛，表現(xiàn)操的完美性，能產(chǎn)生強(qiáng)烈的藝術(shù)感，增強(qiáng)操的藝術(shù)感染力。音樂的選擇要體現(xiàn)獨(dú)特的個(gè)性特點(diǎn)與技術(shù)風(fēng)格，音樂風(fēng)格與編排風(fēng)格要統(tǒng)一，并使技巧、過渡的視覺效果與音樂出現(xiàn)的特效相配合。

從2015年到2018年，短短3年時(shí)間，銀隆估值漲了近10倍。如今，董明珠等外部股東與大股東魏銀倉的矛盾徹底公開，銀隆陷入風(fēng)暴漩渦之中。

PLS-DA是一種有監(jiān)督的識(shí)別方法，它既可以實(shí)現(xiàn)對(duì)數(shù)據(jù)進(jìn)行降維，同時(shí)可以結(jié)合回歸模型利用一定的判別閾值對(duì)回歸結(jié)果進(jìn)行判別分析，同時(shí)因其所建的模型比較穩(wěn)定，且具有很強(qiáng)的抗干擾能力，目前建立產(chǎn)地判別模型時(shí)被廣泛應(yīng)用。對(duì)所有青蘿卜樣品按照省份進(jìn)行分組，分為天津、河南和山東3個(gè)組別，提取青蘿卜樣品中5個(gè)主成分應(yīng)用PLS-DA算法結(jié)合礦物元素含量進(jìn)行數(shù)據(jù)統(tǒng)計(jì)，如圖2元素VIP貢獻(xiàn)值圖看出，Cr、Ca、Na、Mg、Ni、Cd共6種元素在青蘿卜判別時(shí)貢獻(xiàn)率較大。其中Ca、Na、Mg 3種元素與2.1中單因素ANOVA檢驗(yàn)的判別結(jié)果一致。如圖3所示，由3個(gè)省份二維模型圖可以看出，天津樣品、河南樣品以及山東樣品3個(gè)省份的青蘿卜之間均有不同程度的重合，在判別時(shí)，天津僅有2個(gè)被誤判為山東樣品，山東樣品僅有1個(gè)被誤判為天津，不同省份樣品未發(fā)生明顯聚集。因此，基于青蘿卜礦物元素含量結(jié)合PLS-DA算法建立天津、山東、河南3個(gè)省份的青蘿卜樣品判別模型具有一定的實(shí)際應(yīng)用意義。

圖3 蘿卜樣品PLS-DA歸類判別Fig.3 Classification of green turnip from different producing areas by PLS-DA model

2.2.2 LDA

LDA是一種有監(jiān)督判別方法，以投影后類內(nèi)方差最小，類間方差最大的原則對(duì)樣品數(shù)據(jù)進(jìn)行降維，實(shí)現(xiàn)樣品的分類。應(yīng)用LDA算法對(duì)3個(gè)省份的青蘿卜樣品中18種礦物元素進(jìn)行統(tǒng)計(jì)分析，建立判別模型。結(jié)果如圖4所示，在進(jìn)行模型驗(yàn)證時(shí)，河南有4個(gè)樣品被誤判為山東樣品，在山東樣品判別時(shí)，分別有3個(gè)青蘿卜樣品被誤判為河南，2個(gè)青蘿卜樣品被誤判為天津，天津57個(gè)青蘿卜樣品52個(gè)均能被準(zhǔn)確判別，僅有5個(gè)被誤判為山東樣品。河南樣品和天津樣品都容易誤判為山東，這與箱線圖的分析結(jié)果相符。由表3可知，應(yīng)用PLS-DA對(duì)青蘿卜進(jìn)行產(chǎn)地判別模型，天津蘿卜的預(yù)判效果較好，模型驗(yàn)證和交叉訓(xùn)練時(shí)模型的判別率為91.23%，但河南和山東樣品在模型中識(shí)別準(zhǔn)確率較低；采用LDA算法建模，3個(gè)省份模型的驗(yàn)證集和訓(xùn)練集判別率均在90%以上，表明LDA算法在進(jìn)行蘿卜產(chǎn)地溯源時(shí)的判別時(shí)的效果明顯高于PLS-DA，在模型訓(xùn)練時(shí)整體判別率為93.79%，且Ca、Na、Mg元素在建模時(shí)依舊顯示較大貢獻(xiàn)率。

2.創(chuàng)新管理方面。一號(hào)煤礦力爭(zhēng)達(dá)到安全管理的最高境界，強(qiáng)勢(shì)推行自主管理，實(shí)行安全自主管理，充分調(diào)動(dòng)職工的安全積極性、主動(dòng)性，使各單位、各部門和每個(gè)職工各司其職、各負(fù)其責(zé)，解決了安全管理中靠少數(shù)人抓安全、靠監(jiān)督檢查推動(dòng)安全的依賴思想，實(shí)現(xiàn)了安全管理由被動(dòng)向主動(dòng)的轉(zhuǎn)變，很大程度上提高了安全管理質(zhì)量和效率。推行崗位價(jià)值精細(xì)管理，以“崗位”為最小核算單元，核算到崗，使職工由傳統(tǒng)的“勞動(dòng)者”向“經(jīng)營者”轉(zhuǎn)變，以公開透明的查詢平臺(tái)形成激勵(lì)機(jī)制，在生產(chǎn)任務(wù)保證的前提下，礦區(qū)形成了“比、學(xué)、趕、幫、超”的良好氛圍。

圖4 三個(gè)產(chǎn)地青蘿卜樣品LDA歸類判別Fig.4 Classification of green turnip from different producing areas by LDA model

表3 LDA和PLS-DA模型預(yù)測(cè)準(zhǔn)確率 單位：%

目前，王游游等[25]對(duì)絕明子產(chǎn)地溯源時(shí)發(fā)現(xiàn)，δ15N值、δ13C值等穩(wěn)定同位素和礦物元素結(jié)合對(duì)鑒別決明子產(chǎn)區(qū)有更加顯著的效果，篩選出具有顯著差異的19種礦物元素以及其中重要權(quán)重因子，肯定δ15N 值對(duì)農(nóng)產(chǎn)品產(chǎn)地判別時(shí)的貢獻(xiàn)，并提出礦物元素對(duì)產(chǎn)地判別的效果明顯優(yōu)于同位素；楊健等[26]發(fā)現(xiàn)，穩(wěn)定同位素可單獨(dú)達(dá)到鑒別何首烏產(chǎn)地的效果，但穩(wěn)定同位素結(jié)合元素的模型判別率得到進(jìn)一步提高；本文在分析青蘿卜樣品中的礦物元素含量以及同位素比值時(shí)，同樣發(fā)現(xiàn)礦物元素對(duì)產(chǎn)地判別的效果明顯優(yōu)于穩(wěn)定同位素，且單獨(dú)應(yīng)用礦物元素對(duì)天津、山東、河南3個(gè)地區(qū)的青蘿卜樣品的判別率均在90%以上。因此，本文建立青蘿卜樣品的產(chǎn)地判別模型進(jìn)一步肯定了礦物元素在青蘿卜產(chǎn)地鑒別中的實(shí)際應(yīng)用意義，精簡(jiǎn)進(jìn)行青蘿卜溯源需要測(cè)定的變量，這為未來青蘿卜產(chǎn)地溯源提供理論依據(jù)，具有實(shí)際應(yīng)用意義。

3 結(jié)論

本研究通過測(cè)定天津、河南、山東3個(gè)省份不同地區(qū)青蘿卜樣品中的穩(wěn)定同位素比值和礦物元素含量，首次為青蘿卜建立了產(chǎn)地鑒別模型。研究發(fā)現(xiàn)，這3個(gè)產(chǎn)地的樣品有明顯的地域分布特征，箱線圖分析、單因素ANOVA檢驗(yàn)表明，元素在進(jìn)行青蘿卜產(chǎn)地判別的貢獻(xiàn)明顯高于碳、氮穩(wěn)定同位素，18種礦物元素中，Na、Ca、Mg共3種元素在結(jié)合不同的化學(xué)計(jì)量學(xué)方法分析時(shí)均表現(xiàn)出顯著差異性。PLS-DA算法建立判別模型時(shí)，對(duì)山東和天津的樣品識(shí)別率較低，對(duì)天津樣品識(shí)別率較高，模型訓(xùn)練整體判別率僅為86.44%；應(yīng)用LDA算法進(jìn)行產(chǎn)地鑒別時(shí)，模型對(duì)天津、河南、山東青蘿卜樣品的識(shí)別率均較高，模型在驗(yàn)證和訓(xùn)練的總體判別率分別為92.09%和93.79%。2種化學(xué)計(jì)量學(xué)方法進(jìn)行判別分析的結(jié)果顯示，LDA對(duì)不同省份樣品的區(qū)分效果明顯優(yōu)于PLS-DA算法。本文首次對(duì)青蘿卜樣品的礦物元素和穩(wěn)定同位素進(jìn)行分析，并結(jié)合化學(xué)計(jì)量學(xué)方法建立的產(chǎn)地判別模型，LDA算法在建模過程中表現(xiàn)良好。根據(jù)青蘿卜不同產(chǎn)地之間的差異性建模，礦物元素結(jié)合LDA算法的效果理想，可為下一步青蘿卜產(chǎn)地溯源技術(shù)的研究提供新的研究方向和理論依據(jù)。