基于貝葉斯方法的土壤溶質(zhì)遷移轉(zhuǎn)化反演研究為了更好地認(rèn)識(shí)養(yǎng)分和農(nóng)藥等土壤溶質(zhì)的環(huán)境行為,數(shù)值模型常用來(lái)模擬土壤溶質(zhì)的遷移轉(zhuǎn)化過(guò)程。溶質(zhì)在土壤中的遷移轉(zhuǎn)化本質(zhì)上是一個(gè)確定性過(guò)程,但由于觀測(cè)數(shù)據(jù)的稀缺性以及土壤物理、化學(xué)和生物性質(zhì)的空間非均質(zhì)性,遷移轉(zhuǎn)化過(guò)程必然存在不確定性。面對(duì)復(fù)雜的自然環(huán)境狀況,遷移轉(zhuǎn)化模擬不僅需要識(shí)別正確的模型形式,即處理模型結(jié)構(gòu)的不確定性,也需要確定模型參數(shù),即處理參數(shù)的不確定性。如何從包含有測(cè)量誤差的觀測(cè)數(shù)據(jù)中提取有效信息,反演識(shí)別出正確的模型結(jié)構(gòu)及相應(yīng)的最優(yōu)參數(shù),減小模擬過(guò)程中的不確定性,對(duì)于準(zhǔn)確模擬土壤溶質(zhì)遷移轉(zhuǎn)化過(guò)程以及理解其內(nèi)在機(jī)理具有重要的理論和實(shí)踐意義。貝葉斯推斷是一種概率方法,已經(jīng)廣泛應(yīng)用在數(shù)值模擬的反演研究中。基于貝葉斯所進(jìn)行的反演估計(jì)是利用先驗(yàn)認(rèn)識(shí)和觀測(cè)數(shù)據(jù)推斷出目標(biāo)量的后驗(yàn)概率。本文將貝葉斯方法與土壤溶質(zhì)遷移數(shù)值模型耦合,分別對(duì)氮和農(nóng)藥五氯酚(Pentachlorophenol,PCP)的遷移轉(zhuǎn)化過(guò)程進(jìn)行反演識(shí)別,定量解析各個(gè)過(guò)程對(duì)溶質(zhì)消減的貢獻(xiàn)率,具體工作如下:(1)基于廣義速率方程的動(dòng)力學(xué)模型貝葉斯推斷方法在土壤溶質(zhì)遷移轉(zhuǎn)化模擬中,反應(yīng)動(dòng)力學(xué)類型的選擇是很重要的。但在很多情況下,我們卻無(wú)法確定反應(yīng)動(dòng)力學(xué)的形式。為了解決這個(gè)難題,本章發(fā)展了一種新的確定反應(yīng)動(dòng)力學(xué)形式的貝葉斯估計(jì)方法。該方法通過(guò)定義一個(gè)含有離散和連續(xù)參數(shù)的廣義速率方程,利用可以實(shí)現(xiàn)對(duì)離散分布進(jìn)行采樣的馬爾科夫鏈蒙特卡羅算法(MarkovChainMonteCarlo,MCMC)同時(shí)進(jìn)行動(dòng)力學(xué)模型的選擇和動(dòng)力學(xué)參數(shù)的估計(jì),并通過(guò)數(shù)值研究和氮轉(zhuǎn)化實(shí)驗(yàn)驗(yàn)證方法的有效性。數(shù)值算例表明,此方法可以準(zhǔn)確識(shí)別真實(shí)的動(dòng)力學(xué)模型及其動(dòng)力學(xué)參數(shù),在不同的觀測(cè)誤差水平和不同的觀測(cè)數(shù)據(jù)缺失情況下也表現(xiàn)良好。在氮轉(zhuǎn)化的實(shí)驗(yàn)研究中,本方法能夠從多達(dá)19683(39)個(gè)候選氮轉(zhuǎn)化模型中選擇出最佳模型,最佳模型的模擬結(jié)果能夠準(zhǔn)確再現(xiàn)觀測(cè)數(shù)據(jù)。(2)不同溫度和土壤孔隙含水率(Water-filledporespace,WFPS)條件下水稻土中產(chǎn)生N2O主導(dǎo)反應(yīng)的識(shí)別在上述研究的基礎(chǔ)上,我們利用氮轉(zhuǎn)化模型對(duì)不同溫度和WFPS條件下水稻土中的氮轉(zhuǎn)化濃度數(shù)據(jù)進(jìn)行模擬研究,利用MCMC對(duì)不同處理下氮轉(zhuǎn)化進(jìn)行動(dòng)力學(xué)模型選擇和動(dòng)力學(xué)參數(shù)估計(jì)。各處理組的模擬結(jié)果能夠準(zhǔn)確匹配NH4+、N03-和NO2-濃度數(shù)據(jù),也能夠準(zhǔn)確預(yù)測(cè)20T90W(20℃,90%WFPS,下同)、20T120W、30T90W和30T120W的N2O釋放速率數(shù)據(jù)。通過(guò)分析不同溫度和WFPS條件下各反應(yīng)路徑的速率水平及其不確定性,我們推斷在30T120W處理組下RNO2-(an)-N2O(s)(即反硝化菌主導(dǎo)的反硝化反應(yīng)或硝酸鹽異化還原成銨(Dissimilatorynitratereductiontoammonium,DNRA)與ONH4-(ae)-N2O(s)(即自養(yǎng)硝化菌主導(dǎo)的硝化反應(yīng))是N20主要的產(chǎn)生路徑。在30T90W、20T90W和20T120W處理組下,ONH4-(a+,、N2O(s)是N2O主要的產(chǎn)生路徑。在高溫高濕條件下,NH4(ae)-N2O(s)N2O能夠從土壤中逸出,其他處理中土壤中的N2O基本上被還原成N2。(3)基于貝葉斯模型平均的水稻土氮轉(zhuǎn)化路徑反演在上述研究中,我們利用了較復(fù)雜的氮模型,但模擬結(jié)果仍無(wú)法匹配低溫低濕條件下5T和60W處理組的N20釋放速率。為了全面考慮模型結(jié)構(gòu)中的不確定性,避免使用單一模型帶來(lái)的偏差,我們構(gòu)建了12個(gè)不同路徑組合的氮轉(zhuǎn)化模型,應(yīng)用貝葉斯模型平均(Bayesianmodelaveraging,BMA)方法,通過(guò)加權(quán)平均整合12個(gè)模型的模擬結(jié)果進(jìn)行預(yù)測(cè),找到了符合低溫低濕條件下最為合理的模型結(jié)構(gòu)并完成了較為全面的不確定性分析。模擬結(jié)果表明BMA模擬能夠整合各候選模型對(duì)不同氮庫(kù)的模擬能力,全面考慮模型結(jié)構(gòu)的不確定性,從而獲得更好的預(yù)測(cè)效果。(4)PCP在水稻土好氧厭氧界面中消減過(guò)程的反演模擬為了驗(yàn)證反演模型在更復(fù)雜場(chǎng)景中的適用性,我們利用遷移轉(zhuǎn)化模型模擬了水稻土土水好氧-厭氧界面中五氯酚的消減行為。我們主要考慮PCP及其降解產(chǎn)物在土水界面中的擴(kuò)散、吸附和降解過(guò)程,利用MCMC反演了未知的擴(kuò)散反應(yīng)模型參數(shù),采用最大后驗(yàn)概率估計(jì)的參數(shù)反演結(jié)果(maximum-a-posteriori,MAP)模擬再現(xiàn)了PCP在水稻土好氧厭氧界面中的消減,定量化解析了各過(guò)程對(duì)于消減的貢獻(xiàn)率。根據(jù)PCP和Cl-濃度數(shù)據(jù)得出的反演結(jié)果表明,0-6mm土層是整個(gè)土壤剖面中PCP消減的活躍區(qū)域,其中PCP的消減主要是由固相中的降解作用導(dǎo)致的,同時(shí)擴(kuò)散作用將較深層(2.4-4.8mm)的PCP傳輸?shù)綔\層以促進(jìn)PCP的完全消減;積水層的脫氯水平最高,且不同深度土層的PCP脫氯水平相似,2,3,4,5-TeCP和3,4,5-TCP是PCP的主要脫氯產(chǎn)物。綜上所述,本文發(fā)展的方