原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
BaseflowCalculationinMATSIRO
Base flow Calculation in MATSIRO Kei YOSHIMURA 基底流量（Rb）の算定 平林(2004)修正。Ks０を深度2mの値とする。fは0.33。 ?Sand/GSWP2だと9.9E-3[m/s]=9.9[mm/s] ?Clayloam/GSWP2だと5.8E-4[m/s]=0.58[mm/s] 地下水面の位置で大きく異なる。 ?z=0.5mのとき0.6、z=1mのとき0.13、z=2mのとき0.007 tanβの感度が高い。（Ls＝2√3/tanβ） ?grtans-islscp1 grzsd-gswp2での最大値は0.001程度。 ?GTOPO30/日本域/0.1度メッシュで平均スロープ、高度標準偏差を作った場合は、最大0.0003程度 Clayloam、z=0.5mだと、最大9mm/day。 Clayloam、z=2mだと、最大0.09mm/day。 問題点：地下水面（ｚ）形成 土壌層内における地下水面形成過程： 重力と毛管力の釣り合いを考慮。 （地下水面z＝土壌層上端Zk_upから水分ポテンシャルψk分下がった位置。） ? ψk（飽和時0.1m～0.8m）が土壌層厚を上回るときに問題。 旧： ψkが土壌層厚以内のときのみ、地下水面が形成。?土壌が厚いところにしかできない。地下水面が上がらない。 新：飽和している場合、少なくともその層の上端は地下水面とする。 修正コード （たったの３行） 全陸面シミュレーション結果 ERA15を用いた1979年ラン（除く氷面） 全球平均基底流出／トータル流出(P:962.6[mm/y]) 0.15% ? 31.5% (0.5/343.6 [mm/y] ? 112.1/355.1 [mm/y]) 参考：蒸散/総蒸発散 182.5/618.6 [mm/y] ? 181.6/607.3 [mm/y] 河川流量シミュレーション結果 利根川（140.5N，36.5N） 河川流量シミュレーション結果 ナコンサワン（100.5N，16.5N） 日本域シミュレーション結果 GPV-MSMを用いた2004年ラン（Thanks to 小岩君） 平均基底流出／トータル流出(P:1792.0[mm/y]) 0.3% ? 11.2% (3.0/1004.1 [mm/y] ? 112.9/1009.7 [mm/y]) 参考：蒸散/総蒸発散 222.0/700.1 [mm/y] ? 222.1/695.4 [mm/y] 河川流量シミュレーション結果 利根川栗橋（139.65N，36.15N）, 8939km2 MATSIRO-MSM 06/02/06-07 * ① ② ③ In matroi.F (matrof.F in the original MATSIRO) * mean water table depth * DO 1300 IJ = ISTA, IEND IF ( ILSFC( IJ ).GT.1 .AND. ISAT ( IJ ).GE.0 .AND. K.LT.KFTAB( IJ ) ) THEN WTABDX = WLEVM( K ) - GPSI( IJ,K ) IF ( WTABDX .LE. WLEVM( K+1 ) ) THEN WTABD( IJ ) = WTABDX KWTAB( IJ ) = K *kei ELSEIF ( K.NE.KWMAX ) THEN WTABD( IJ ) = WLEVM( K+1 ) KWTAB( IJ ) = K+1 *kei ENDIF IF ( GLW(IJ,K) .LT. GWS(IJ,K)-EPSGW ) ISAT( IJ ) = -1 ENDIF 1300 CONTINUE 修正後 修正前 *
文档评论（0）