跳到主要內容

發表文章

目前顯示的是有「Soybean」標籤的文章

農業溫室氣體排放估算 - 水田甲烷(參數篇)

基本概念可參考上一篇文章 https://learncropmodel.blogspot.com/2025/05/blog-post.html 上一篇文章提及碳排係數在Tier1當中的拆解方式如下,我們來看看文獻裡的參數吧 $EF_i = EF_c \cdot SF_w \cdot SF_p \cdot SF_o$ EF i : i條件的每日排放係數 EF c : 無施用有機肥料下持續灌溉水田的排放係數 SF w : 不同水分管理區域的調整因子(表 5.12) SF p : 不同種植前水分管理調整因子(表 5.13) SF o : 有機質肥料型態或施用量之調整因子(Eq. 5.3, 表5.14) Baseline emission factor 也就是EF c 的預設係數來自於表5.1 (page 5.53) 表5.1 預設水田甲烷排放係數,假設水稻栽培時間小於180天、連續淹水、未施用有機質 World Regional Emission factor (kg CH 4 ha -1 d -1 ) Error range (kg CH 4 ha -1 d -1 ) Region Emission factor (kg CH 4 ha -1 d -1 ) Error range (kg CH 4 ha -1 d -1 ) 1.19 0.80 - 1.76 Africa 1.19 0.80 - 1.76 East Asia 1.32 0.89 - 1.96 Southeast Asia 1.22 0.83 - 1.81 South Asia 0.85 0.58 - 1.26 Europe 1.56 1.06 - 2.31 ...

GLYCIM模式介紹-植物參數

GLYCIM模式最早在1982由Dr. Acock等人建立,經過Dr. V.R. Reddy, Dr. Yakov Pachevsky, Dr. Dennis Timlin 等人測試與修訂,根據Dr. Reddy口述,Ken Boote最早也有參與GLYCIM模式開發,後來依據GLYCIM模式的基礎建立CROPGRO模式,因此CROPCRO與GLYCIM模式有一定的相似度。原本的GLYCIM模式使用Penman monteith進行水平衡模擬、Hyperbolic equation模擬光合作用,該模式相當令人驚豔之處在於能夠反映淹水情境下,根部活性降低、生長速率下降與根瘤死亡,可以從上述的現象探討淹水情境下,大豆的生理變化以及元素平衡,對於最終產量的影響。 Dr. Timlin在20年以前就開始協助GLYCIM模式土壤模型的建立,這幾年重新將GLYCIM模式與2DSOIL結合成為2DGLYCIM,並且使用FvCB取代Hyperbolic equation進行光合作用模擬,可想像是另一個昇級版的大豆模式。目前模式的原始碼仍然是公開的,只要有興趣的人都可以在Github上找到GLYCIM,Dr. Reddy的團隊一直歡迎大家取用模式,這個氣度永遠是我無法企及的。 GLYCIM模式所需的作物參數檔下表,可以看到參數檔當中有相當大的部分用來模擬作物的發育,模式的時步(time step)是小時,在這個模式既可以觀察到Dr. Acock盡量使用生長/發育速率的取代累積溫度,很大的程度影響到MAIZSIM模式。GLYCIM模式假設的基礎溫度是0℃,並無optimum temperature,值得觀察高溫環境下是否能使用beta function來取代現行的linear equation。 代碼 Definition 定義 備註 MG maturity group number (group 00 = 1) 成熟群 SEEDLB number of seeds per pound weight typical for cultivar 每磅籽粒數 Sink FILL Seed f...

GLYCIM模式待討論事項

  1.        株高 :印象中這邊是原始碼有點問題 2.       模擬的莢數 >> 實際莢數: 3.       模擬的分支數 >> 實際分支數: rosetta 程式需要有visual studio才能開啟(還要有C++套件),有可能是當初building的時候有問題 1. 株高 GLYCIM程式中,與株高相關的參數為PARM(21)與PARM(22),模式假設株高是Vstage的指數關係,也就是 株高 = Parm(21)*VSTAGE^PARM(22) 為了表示增加速率,我們把上式進行微分可獲得 d株高/dvstage = Parm(21)*Parm(22)*Vstage^(Parm(22)-1) 程式碼在2750行 !YA PDMH=(PARM(21)+(PARM(22)*(VSTAGE+PDV/2.0)**1.37))*PDV*SLOW*HTFATR IF (VSTAGE.GT.0.0) THEN F1=PARM(21)*PARM(22)*(VSTAGE)**(PARM(22)-1.) ELSE F1=0.0 END IF IF (VSTAGE+PDV.GT.0.0) THEN F2=PARM(21)*PARM(22)*(VSTAGE+PDV)**(PARM(22)-1.) ELSE F2=0.0 END IF PDMH=0.5*(F1+F2)*PDV 這裡的計算方法很值得學習,它分別計算 Vstage (這一個時間點)和 Vstage + dv (下一個時間點所對應的株高增長速度,再計算平均植以取得較準確的生長速率值。 值得注意的是,在第4270行的時候開始累計株高,這邊使用Maturity group和 R stage作為條件 ! CALCULATE CHANGES IN MAINSTEM AND B...