站點觀測AGPP源于ChinaFLUX觀測數(shù)據(jù)及公開數(shù)據(jù)。
自2002年以來,ChinaFLUX已在62個生態(tài)系統(tǒng)陸續(xù)開展渦度相關觀測,積累了332個站點年的原始觀測結果。基于ChinaFLUX原始觀測數(shù)據(jù),利用ChinaFLUX通用數(shù)據(jù)處理流程對數(shù)據(jù)質量進行控制和插補
[11],獲得各站點年ChinaFLUX觀測數(shù)據(jù)的AGPP值。
公開數(shù)據(jù)包括兩部分:文獻數(shù)據(jù)和通量網(wǎng)絡公開共享數(shù)據(jù)。以“渦度相關”或“eddy covariance”為關鍵詞,本研究自中國知網(wǎng)(CNKI)和Web of Science檢索中國區(qū)域渦度相關觀測的中英文文獻。對于搜索到的條目,本數(shù)據(jù)集逐篇篩選包含AGPP觀測數(shù)值的文獻進行提取,具體篩選標準包括:(1)文獻報道的結果為渦度相關直接觀測的AGPP數(shù)值,排除部分文獻使用模型模擬報道的結果;(2)具有一年以上(含一年)連續(xù)觀測數(shù)據(jù);(3)有總初級生產(chǎn)力年累積量或所有月份月總量的數(shù)值報道
[12]。經(jīng)過篩選,本數(shù)據(jù)集逐篇提取各文獻報道的站點信息(包括經(jīng)緯度、海拔、氣候要素、觀測年份),并以年為單位記錄各AGPP的觀測數(shù)值。除了文獻數(shù)據(jù),本數(shù)據(jù)集同步查詢了已有通量觀測研究網(wǎng)絡如FLUXNET、ChinaFLUX共享的公開數(shù)據(jù),進一步補充觀測數(shù)據(jù)集。公開數(shù)據(jù)共包括128個生態(tài)系統(tǒng),540個站點年。
同時,為了避免同一站點同一觀測年份出現(xiàn)多個數(shù)據(jù)源的重復數(shù)據(jù)(如ChinaFLUX數(shù)據(jù)共享與文獻數(shù)據(jù)均有某一站點某一年的觀測),如果某一站點擁有多個數(shù)據(jù)源的數(shù)值,本研究選用ChinaFLUX使用統(tǒng)一數(shù)據(jù)處理流程所得結果。基于文獻搜索和數(shù)據(jù)整編,共獲得166個站點,641個站點年,跨越2002–2020年的AGPP數(shù)據(jù)。
輔助數(shù)據(jù)是基于觀測AGPP進行尺度擴展的重要工具,但文獻數(shù)據(jù)及通量觀測網(wǎng)絡公開數(shù)據(jù)提供輔助數(shù)據(jù)的一致性也存在差異,如文獻收集的MAT和MAP與站點觀測值的相關系數(shù)分別為0.95和0.88,而海拔可以達到0.99。因此,本研究基于觀測數(shù)據(jù)的地理信息提取相應氣候土壤生物輔助數(shù)據(jù)。氣候要素包括6個變量,分別為:年均氣溫(MAT)、年總降水量(MAP)、年總光合有效輻射(PAR)、年總潛在蒸散(PET)、年均飽和水汽壓差(VPD)、年均CO2密度(ρc,yr )。土壤要素包括3個變量:土壤年均濕度(SM)、土壤有機碳含量(SOC)、土壤總氮含量(STN)。生物要素包括2個變量:年均葉面積指數(shù)(LAI)和年最大葉面積指數(shù)(MLAI)。
各站點年的MAT、MAP及飽和水汽壓差(VPD)數(shù)據(jù)源自英國氣候研究組織(CRU)的長時間序列結果
[13]。基于CRU的0.5°空間分辨率的逐月氣溫、降水量、水氣壓數(shù)據(jù),結合世界氣象組織(www.worldclim.org)的各月長期平均氣溫、降水量和水氣壓結果,本研究利用Delta空間降尺度方法生成中國區(qū)域1 km空間分辨率的數(shù)據(jù)產(chǎn)品
[14]。
基于氣候研究組織(Climate Research Unit, CRU)的長時間序列逐月潛在蒸散(PET)數(shù)據(jù),利用通過樣條插值法將0.5°空間分辨率的逐月PET至30 arcsecond(約1 km),進而求算得到各年PET。
年總光合有效輻射(PAR)數(shù)據(jù)采用全球陸地表面衛(wèi)星數(shù)據(jù)產(chǎn)品(Global Land Surface Satellite, GLASS)的數(shù)據(jù)(0.05°)。基于GLASS的逐日PAR數(shù)據(jù)
[15-16]累加得到年PAR值,并使用樣條插值方法插值到30 arcsecond(約1 km)。
本研究利用全球陸地表面衛(wèi)星(global land surface satellite, GLASS)產(chǎn)品提取8天尺度的吸收光合有效輻射所占比例(FAPAR)和逐日的光合有效輻射,并基于8天尺度的FAPAR及逐日PAR計算得到該年的APAR數(shù)據(jù),并使用樣條插值方法插值到30arcsecond(約1 km)。
年均CO
2密度(
ρc,yr )則基于大氣CO
2濃度(
bc )、CO
2摩爾質量(44 g mol
-1)和當前狀態(tài)的氣體摩爾體積(
V)計算而得到,其中大氣CO
2濃度采用Mauna Roa的觀測值來替代,當前狀態(tài)的氣體摩爾體積采用理想氣體狀態(tài)方程結合當?shù)貧鈮骸AT計算而得到。當?shù)貧鈮簞t利用年均氣溫和當?shù)睾0尾捎脡焊吖接嬎愣玫?sup>[
17]。
本研究利用空間分辨率為0.05°、時間分辨率為年的微波遙感反演的2002–2018年中國土壤水分降尺度數(shù)據(jù)集
[18],表征觀測年份的土壤濕度,并使用樣條插值方法插值到30 arcsecond(約1 km)。
基于全球土壤數(shù)據(jù)庫(Global Soil Dataset for use in Earth System Models, GSDE)的表層5 cm深的土壤有機碳和土壤總氮密度及土壤容重數(shù)據(jù),計算得到各站點的表層土壤有機碳(gC m
-2)和全氮含量(gN m
-2)
[19]。
年均葉面積指數(shù)(LAI)和最大葉面積指數(shù)(MLAI)基于空間分辨率為1 km的MODIS數(shù)據(jù)計算而獲取。本研究利用經(jīng)過驗證的8天尺度MODIS數(shù)據(jù)
[20] (http://globalchange.bnu.edu.cn/research/lai)計算各像元的年均LAI及年內LAI最大值(MLAI)。