CN113111531A - 面向分布式水文模型的季节性冻土区冻土层厚度模拟方法 - Google Patents

面向分布式水文模型的季节性冻土区冻土层厚度模拟方法 Download PDF

Info

Publication number
CN113111531A
CN113111531A CN202110439329.8A CN202110439329A CN113111531A CN 113111531 A CN113111531 A CN 113111531A CN 202110439329 A CN202110439329 A CN 202110439329A CN 113111531 A CN113111531 A CN 113111531A
Authority
CN
China
Prior art keywords
soil
temperature
frozen soil
layer
frozen
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110439329.8A
Other languages
English (en)
Other versions
CN113111531B (zh
Inventor
孙青言
陆垂裕
郭辉
严聆嘉
何鑫
秦韬
吴初
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Heilongjiang Prov Investigation Design And Research Institute Of Water Conservancy And Hydropower
China Institute of Water Resources and Hydropower Research
Original Assignee
Heilongjiang Prov Investigation Design And Research Institute Of Water Conservancy And Hydropower
China Institute of Water Resources and Hydropower Research
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Heilongjiang Prov Investigation Design And Research Institute Of Water Conservancy And Hydropower, China Institute of Water Resources and Hydropower Research filed Critical Heilongjiang Prov Investigation Design And Research Institute Of Water Conservancy And Hydropower
Priority to CN202110439329.8A priority Critical patent/CN113111531B/zh
Publication of CN113111531A publication Critical patent/CN113111531A/zh
Application granted granted Critical
Publication of CN113111531B publication Critical patent/CN113111531B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,包括:S1、建立季节性冻土区冻土层中土壤温度与其影响因素之间的关系式;S2、基于土壤温度与其影响因素之间的关系式,建立冻土层厚度与土壤温度之间的关系,进而实现季节性冻土区冻土层厚度模拟。本发明方法通过先建立土壤温度与相关影响因素之间的关系,然后再建立冻土厚度与土壤温度之间的关系,提高了分布式水文模型在气候寒冷的季节性冻土区的适用性,同时,利用冻土区厚度监测数据率定模型中的参数变量,进一步提高模拟精度。

Description

面向分布式水文模型的季节性冻土区冻土层厚度模拟方法
技术领域
本发明涉及季节性冻土技术领域,具体涉及一种面向分布式水文模型的季节性冻土区冻土层厚度模拟方法。
背景技术
分布式水文模型用于季节性冻土区水文循环模拟时,往往面临着土壤冻融的影响,与暖热地区相比,寒区降雪量的比例更高,土壤冻结的时间更长,气温回升后融雪水因冻土的存在不能及时下渗补给地下水,在冻土区以上形成高含水率土壤层,甚至形成饱和土壤层,造成规模较大的春汛。水文模型如果不能较好地刻画冻土层的冻融过程,将影响地下水补给、土壤水含量、春汛径流的模拟精度。
季节性冻土的形成、发展及消融受多种因素的影响,包括气温、土温、土壤类型、土壤含水率、地表覆盖度等,在时间跨度上和空间分布上存在很大的差异。目前的分布式水文模型对冻土的动态一般都是简化处理甚至不予考虑,这在低纬度地区的水文模拟中尚能接收,但是在气候寒冷的高海拔、高纬度地区,则难以适用。
发明内容
针对现有技术中的上述不足,本发明提供的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法解决了现有的分布式水文模型在气候寒冷的季节性冻土区的适用性和模拟精度都不高的问题。
为了达到上述发明目的,本发明采用的技术方案为:面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,包括以下步骤:
S1、建立季节性冻土区冻土层中土壤温度与其影响因素之间的关系式;
S2、基于土壤温度与其影响因素之间的关系式,建立冻土层厚度与土壤温度之间的关系,进而实现季节性冻土区冻土层厚度模拟。
进一步地,所述步骤S1中,以日为尺度,建立的土壤温度与其影响因素之间的关系式为:
Figure BDA0003034407590000021
式中,Tsoil(z,n)和Tsoil(z,n-1)分别为土壤深度z处、年内第n日和第n-1日的土壤温度,单位为℃;fl为第n-1日土壤温度对n日土壤温度的延迟影响因子,其取值范围为0~1;fd为深度对土壤温度的影响因子;
Figure BDA0003034407590000022
为年平均气温,由气温站观测数据统计获得,单位为℃;Tssurf为土表温度,单位为℃。
进一步地,所述影响因子fd的计算公式为:
Figure BDA0003034407590000023
式中,zs为土壤深度z与土壤温度稳定深度ds的比值;
其中,土壤温度稳定深度ds为:
Figure BDA0003034407590000024
式中,
Figure BDA0003034407590000025
为土壤含水量对土壤温度稳定埋深的影响因子,与土壤容量、土壤剖面厚度及其含水量有关,其计算公式为:
Figure BDA0003034407590000026
dsmax为最大土壤温度稳定深度,与土壤容量有关,其计算公式为:
Figure BDA0003034407590000027
式中,ρb为土壤容重,单位为Mg/m3;SW为水文模型模拟的土壤剖面总厚度内的总含水量,单位为mm。
进一步地,所述土表温度Tssurf为:
Tssurf=fcv·Tsoil(10,n-1)+(1-fcv)·Tbare
式中,fcv为地表覆盖物影响土表温度的权重因子;Tsoil(10,n-1)为深度为10mm处第n-1日的土壤温度,用于表示影响土表温度的前一日表层土壤温度,单位为℃;Tbare为裸土土表温度℃;
其中,地表覆盖物影响土表温度的权重因子fcv为:
fcv=max(fveg,fsno)
式中,fveg为植被生物量和残余量覆盖影响土表温度的权重;fsno为积雪覆盖地面影响土表温度的权重;
其中,fveg的计算公式为:
Figure BDA0003034407590000031
fsno的计算公式为:
Figure BDA0003034407590000032
式中,veg为植被生物量和残余量,单位为kg/ha;sno为积雪当量,单位为mm;
裸土土表温度Tbare的计算公式为:
Figure BDA0003034407590000033
式中,
Figure BDA0003034407590000034
Tmax和Tmin分别为日平均气温、日最高气温和日最低气温,单位为℃;εsr为辐射项,其计算公式为:
Figure BDA0003034407590000035
式中,Hday为当天到达地表的太阳辐射量,单位为MJ m-2d-1;α为当天的短波反射率或反照率。
进一步地,所述步骤S2具体为:
S21、将土壤剖面分为M层,基于土壤温度与其影响因素之间的关系式,通过线性插值法计算冻土层的上边界和下边界;
S22、基于冻土层的上边界和下边界,计算冻土层的厚度,进而实现季节性冻土区冻土层厚度模拟。
进一步地,所述步骤S21中,冻土层的上边界zup的计算公式为:
Figure BDA0003034407590000041
式中,zup为冻土层的上边界的深度,单位为mm;下标m为土壤剖面的分层序号,从地表z=0处到土壤剖面底部z=ztot处共分为M层,序号m分别为1,2,3,...,M,m取1~M-1;zm和zm+1分别为土壤第m层和第m+1层中心的深度,单位为mm;T(z1,n)、T(zm,n)和T(zm+1,n)分别为土壤第1层、第m层和第m+1层中心的深度,单位为℃;
冻土层的下边界zdown的计算公式为:
Figure BDA0003034407590000042
T(zm,n)≤0and T(zm+1,n)>0
所述步骤S22中,冻土层的厚度Δz为:
Δz=zdown-zup
本发明的有益效果为:
本发明方法通过先建立土壤温度与相关影响因素之间的关系,然后再建立冻土厚度与土壤温度之间的关系,提高了分布式水文模型在气候寒冷的季节性冻土区的适用性,同时,利用冻土区厚度监测数据率定模型中的参数变量,进一步提高模拟精度。
附图说明
图1为本发明提供的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法流程图。
图2为本发明提供的季节性冻土区冻土发展过程示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
本发明实施例中,冻土的冻结与融化直接受到土壤温度的影响,而土壤温度与气温、地表覆盖物、土壤含水率等相关,因此,如图1所示,本发明中的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,包括以下步骤:
S1、建立季节性冻土区冻土层中土壤温度与其影响因素之间的关系式;
S2、基于土壤温度与其影响因素之间的关系式,建立冻土层厚度与土壤温度之间的关系,进而实现季节性冻土区冻土层厚度模拟。
上述步骤S1中,根据水文模型开展日尺度水循环模拟的需求,以日尺度计算土壤日平均温度和冻土日平均厚度(以下无特殊说明,均以土壤温度和冻土厚度简称),某一年第n日、深度z处的土壤温度与年平均气温、土表温度、前一天的土壤温度有关,基于此,以日为尺度,建立的土壤温度与其影响因素之间的关系式为:
Figure BDA0003034407590000051
式中,Tsoil(z,n)和Tsoil(z,n-1)分别为土壤深度z处、年内第n日和第n-1日的土壤温度,单位为℃;fl为第n-1日土壤温度对第n日土壤温度的延迟影响因子,其取值范围为0~1;fd为深度对土壤温度的影响因子;
Figure BDA0003034407590000061
为年平均气温,由气温站观测数据统计获得,单位为℃;Tssurf为土表温度,单位为℃。
其中,初始土壤温度Tsoil(z,1)和延迟影响因子fl可以由用户赋初值,土壤深度影响因子fd和土表温度Tssurf需要进一步计算;
具体地,影响因子fd的计算公式为:
Figure BDA0003034407590000062
式中,zs为土壤深度z与土壤温度稳定深度ds的比值,即:
Figure BDA0003034407590000063
从上式(1)~(3)中可以看出,土壤深度越接近地表,zs和fd趋近于0,土壤温度与土表温度关系越密切;随着土壤深度增大,越接近土壤温度稳定深度ds,zs趋近于1,fd趋近于0.95,土壤温度与年平均气温越密切。由于土壤温度稳定深度ds受到土壤含水量的影响,每天都随着土壤含水量的变化而变化,因此,ds需要逐日计算,土壤温度稳定深度ds的计算公式为:
Figure BDA0003034407590000064
式中,
Figure BDA0003034407590000065
为土壤含水量对土壤温度稳定埋深的影响因子,与土壤容量、土壤剖面厚度及其含水量有关,其计算公式为:
Figure BDA0003034407590000066
dsmax为最大土壤温度稳定深度,与土壤容量有关,其计算公式为:
Figure BDA0003034407590000071
式中,ρb为土壤容重,可分类检测获得,也可以在“土壤科学数据库”获得,单位为Mg/m3;SW为水文模型模拟的土壤剖面总厚度内的总含水量,单位为mm。
具体地,上述土表温度Tssurf的计算公式为:
Tssurf=fcv·Tsoil(10,n-1)+(1-fcv)·Tbare (7)
式中,fcv为地表覆盖物影响土表温度的权重因子;Tsoil(10,n-1)为深度为10mm处第n-1日的土壤温度,用于表示影响土表温度的前一日表层土壤温度,初始值可以由用户给定,单位为℃;Tbare为裸土土表温度℃;
地面覆盖物一般分为两类,一是植被及其残余,二是积雪;前者不断影响土壤温度,而且影响蒸散发,因此分布式水文模型一般都具有模拟地表植被生长、成熟、休眠、收割等过程的功能,能够计算地表植被生物量及其脱落地表的残余量。模拟降雪-积雪-融雪过程是分布式水文模型的基本功能之一。基于此,地表覆盖物影响土表温度的权重因子fcv为:
fcv=max(fveg,fsno) (8)
式中,fveg为植被生物量和残余量覆盖影响土表温度的权重;fsno为积雪覆盖地面影响土表温度的权重;
其中,fveg的计算公式为:
Figure BDA0003034407590000072
fsno的计算公式为:
Figure BDA0003034407590000073
式中,veg为植被生物量和残余量,可在水文模型中逐日计算,单位为kg/ha;sno为积雪当量,可在水文模型中逐日模拟,单位为mm;
裸土的土表温度Tbare不但与气温有关,还受到太阳辐射的影响,计算公式为:
Figure BDA0003034407590000081
式中,
Figure BDA0003034407590000082
Tmax和Tmin分别为日平均气温、日最高气温和日最低气温,单位为℃;εsr为辐射项,其计算公式为:
Figure BDA0003034407590000083
式中,Hday为当天到达地表的太阳辐射量,由水文模型根据辐射/日照实测数据计算得到,单位为MJ m-2d-1;α为当天的短波反射率或反照率,为土壤的重要参数,可由用户赋值。
在季节性冻土区,随着气温降低,表层土壤首先冻结,然后冻土层逐日增厚,气温回升到0度以上,表层土壤开始融化,直至冻土层融通(如图2所示);
一般而言,土壤温度沿深度垂向连续变化,这样模型可以将土壤剖面分为若干层,分别计算每层土壤中心处的土壤温度,其他深度处的温度可以插值计算,以减少模型运算的负担。如果相邻两个土壤层中心处的温度,一个大于0℃,一个小于0℃,则通过插值可以确定温度为0℃的深度,由此确定冻土层边界的深度。在土壤剖面分层方面,可等厚度平分,也可在水文模型的土壤分层基础上进一步细分,温度变化比较剧烈的表层土壤可以细分,总之,用户可根据需求灵活分层。因此,上述步骤S2具体为:
S21、将土壤剖面分为M层,基于土壤温度与其影响因素之间的关系式,通过线性插值法计算冻土层的上边界和下边界;
S22、基于冻土层的上边界和下边界,计算冻土层的厚度,进而实现季节性冻土区冻土层厚度模拟。
上述步骤S21中,冻土层的上边界zup的计算公式为:
Figure BDA0003034407590000091
式中,zup为冻土层的上边界的深度,单位为mm;下标m为土壤剖面的分层序号,从地表z=0处到土壤剖面底部z=ztot处共分为M层,序号m分别为1,2,3,...,M,m取1~M-1;zm和zm+1分别为土壤第m层和第m+1层中心的深度,单位为mm;T(z1,n)、T(zm,n)和T(zm+1,n)分别为土壤第1层、第m层和第m+1层中心的深度,单位为℃;土表温度和土壤第一层中心温度都低于0℃时,说明冻土上边界已经达到土表;土表温度大于0℃,土壤第一层中心温度低于0℃,冻土上边界位于两者之间,也采用线性插值方法计算。
冻土层的下边界zdown的计算公式为:
Figure BDA0003034407590000092
步骤S22中,冻土层的厚度Δz为:
Δz=zdown-zup (15)
季节性冻土区冻土层在垂向上一般分为单层,局部气候异常可能导致局部区域、个别时段出现多层,可以采用上述方法分别计算各冻土层的厚度;本实施例中的季节性冻土区冻土层厚度模拟需要嵌入到分布式水文模型中,需要用户对部分参数变量赋初值,还可以利用冻土厚度监测值率定模型中的这些参数变量,提高模拟精度。

Claims (6)

1.面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,其特征在于,包括以下步骤:
S1、建立季节性冻土区冻土层中土壤温度与其影响因素之间的关系式;
S2、基于土壤温度与其影响因素之间的关系式,建立冻土层厚度与土壤温度之间的关系,进而实现季节性冻土区冻土层厚度模拟。
2.根据权利要求1所述的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,其特征在于,所述步骤S1中,以日为尺度,建立的土壤温度与其影响因素之间的关系式为:
Figure FDA0003034407580000011
式中,Tsoil(z,n)和Tsoil(z,n-1)分别为土壤深度z处、年内第n日和第n-1日的土壤温度,单位为℃;fl为第n-1日土壤温度对第n日土壤温度的延迟影响因子,其取值范围为0~1;fd为深度对土壤温度的影响因子;
Figure FDA0003034407580000012
为年平均气温,由气温站观测数据统计获得,单位为℃;Tssurf为土表温度,单位为℃。
3.根据权利要求2所述的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,其特征在于,所述影响因子fd的计算公式为:
Figure FDA0003034407580000013
式中,zs为土壤深度z与土壤温度稳定深度ds的比值;
其中,土壤温度稳定深度ds为:
Figure FDA0003034407580000014
式中,
Figure FDA0003034407580000015
为土壤含水量对土壤温度稳定埋深的影响因子,与土壤容量、土壤剖面厚度及其含水量有关,其计算公式为:
Figure FDA0003034407580000016
dsmax为最大土壤温度稳定深度,与土壤容量有关,其计算公式为:
Figure FDA0003034407580000021
式中,ρb为土壤容重,单位为Mg/m3;SW为水文模型模拟的土壤剖面总厚度内的总含水量,单位为mm。
4.根据权利要求2所述的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,其特征在于,所述土表温度Tssurf的计算公式为:
Tssurf=fcv·Tsoil(10,n-1)+(1-fcv)·Tbare
式中,fcv为地表覆盖物影响土表温度的权重因子;Tsoil(10,n-1)为深度为10mm处第n-1日的土壤温度,用于表示影响土表温度的前一日表层土壤温度,单位为℃;Tbare为裸土土表温度℃;
其中,地表覆盖物影响土表温度的权重因子fcv为:
fcv=max(fveg,fsno)
式中,fveg为植被生物量和残余量覆盖影响土表温度的权重;fsno为积雪覆盖地面影响土表温度的权重;
其中,fveg的计算公式为:
Figure FDA0003034407580000022
fsno的计算公式为:
Figure FDA0003034407580000023
式中,veg为植被生物量和残余量,单位为kg/ha;sno为积雪当量,单位为mm;
裸土土表温度Tbare的计算公式为:
Figure FDA0003034407580000031
式中,
Figure FDA0003034407580000032
Tmax和Tmin分别为日平均气温、日最高气温和日最低气温,单位为℃;εsr为辐射项,其计算公式为:
Figure FDA0003034407580000033
式中,Hday为当天到达地表的太阳辐射量,单位为MJ m-2d-1;α为当天的短波反射率或反照率。
5.根据权利要求2所述的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,其特征在于,所述步骤S2具体为:
S21、将土壤剖面分为M层,基于土壤温度与其影响因素之间的关系式,通过线性插值法计算冻土层的上边界和下边界;
S22、基于冻土层的上边界和下边界,计算冻土层的厚度,进而实现季节性冻土区冻土层厚度模拟。
6.根据权利要求5所述的面向分布式水文模型的季节性冻土区冻土层厚度模拟方法,其特征在于,所述步骤S21中,冻土层的上边界zup的计算公式为:
Figure FDA0003034407580000034
式中,zup为冻土层的上边界的深度,单位为mm;下标m为土壤剖面的分层序号,从地表z=0处到土壤剖面底部z=ztot处共分为M层,序号m分别为1,2,3,...,M,m取1~M-1;zm和zm+1分别为土壤第m层和第m+1层中心的深度,单位为mm;T(z1,n)、T(zm,n)和T(zm+1,n)分别为土壤第1层、第m层和第m+1层中心的深度,单位为℃;
冻土层的下边界zdown的计算公式为:
Figure FDA0003034407580000041
所述步骤S22中,冻土层的厚度Δz为:
Δz=zdown-zup
CN202110439329.8A 2021-04-23 2021-04-23 面向分布式水文模型的季节性冻土区冻土层厚度模拟方法 Active CN113111531B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110439329.8A CN113111531B (zh) 2021-04-23 2021-04-23 面向分布式水文模型的季节性冻土区冻土层厚度模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110439329.8A CN113111531B (zh) 2021-04-23 2021-04-23 面向分布式水文模型的季节性冻土区冻土层厚度模拟方法

Publications (2)

Publication Number Publication Date
CN113111531A true CN113111531A (zh) 2021-07-13
CN113111531B CN113111531B (zh) 2021-10-22

Family

ID=76719548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110439329.8A Active CN113111531B (zh) 2021-04-23 2021-04-23 面向分布式水文模型的季节性冻土区冻土层厚度模拟方法

Country Status (1)

Country Link
CN (1) CN113111531B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114707344A (zh) * 2022-04-15 2022-07-05 西南交通大学 一种基于系统动力学的多年冻土活动层厚度计算方法
CN114818367A (zh) * 2022-05-18 2022-07-29 珠江水利委员会珠江水利科学研究院 适用于高寒地区的径流模拟分析方法、系统及介质
CN114970143A (zh) * 2022-05-24 2022-08-30 武汉大学 水库消落带温室气体净通量的评估方法
CN115540788A (zh) * 2022-11-08 2022-12-30 中南大学 多年冻土活动层厚度估算方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4943918A (en) * 1985-01-09 1990-07-24 Phillips Petroleum Company Seismic data processing method
CN104036134A (zh) * 2014-06-10 2014-09-10 北京林业大学 冻土含水量和温度的计算方法
CN206002050U (zh) * 2016-09-23 2017-03-08 国网四川省电力公司经济技术研究院 一种季节性冻土厚度监测装置
CN106951612A (zh) * 2017-03-06 2017-07-14 河海大学 土壤冻融过程中动态蓄水容量产流计算方法
CN109871618A (zh) * 2019-02-22 2019-06-11 中国水利水电科学研究院 一种在土壤水热耦合模拟中确定底层温度边界的方法
CN110442985A (zh) * 2019-08-09 2019-11-12 中国地质科学院探矿工艺研究所 一种考虑冻融循环效应的泥石流防治参数设计优化方法
CN112016192A (zh) * 2020-08-07 2020-12-01 武汉大学 盐渍化季节性冻融土壤蒸发的模拟方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4943918A (en) * 1985-01-09 1990-07-24 Phillips Petroleum Company Seismic data processing method
CN104036134A (zh) * 2014-06-10 2014-09-10 北京林业大学 冻土含水量和温度的计算方法
CN206002050U (zh) * 2016-09-23 2017-03-08 国网四川省电力公司经济技术研究院 一种季节性冻土厚度监测装置
CN106951612A (zh) * 2017-03-06 2017-07-14 河海大学 土壤冻融过程中动态蓄水容量产流计算方法
CN109871618A (zh) * 2019-02-22 2019-06-11 中国水利水电科学研究院 一种在土壤水热耦合模拟中确定底层温度边界的方法
CN110442985A (zh) * 2019-08-09 2019-11-12 中国地质科学院探矿工艺研究所 一种考虑冻融循环效应的泥石流防治参数设计优化方法
CN112016192A (zh) * 2020-08-07 2020-12-01 武汉大学 盐渍化季节性冻融土壤蒸发的模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李保琦 等: ""寒区水文循环模型研究进展"", 《西南民族大学学报( 自然科学版)》 *
蔡汉成 等: ""青藏铁路沿线多年冻土区气温和多年冻土"", 《岩石力学与工程学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114707344A (zh) * 2022-04-15 2022-07-05 西南交通大学 一种基于系统动力学的多年冻土活动层厚度计算方法
CN114707344B (zh) * 2022-04-15 2023-09-19 西南交通大学 一种基于系统动力学的多年冻土活动层厚度计算方法
CN114818367A (zh) * 2022-05-18 2022-07-29 珠江水利委员会珠江水利科学研究院 适用于高寒地区的径流模拟分析方法、系统及介质
CN114970143A (zh) * 2022-05-24 2022-08-30 武汉大学 水库消落带温室气体净通量的评估方法
CN115540788A (zh) * 2022-11-08 2022-12-30 中南大学 多年冻土活动层厚度估算方法
CN115540788B (zh) * 2022-11-08 2023-08-29 中南大学 联合多轨道InSAR形变观测和未冻水含量的活动层厚度估算方法

Also Published As

Publication number Publication date
CN113111531B (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN113111531B (zh) 面向分布式水文模型的季节性冻土区冻土层厚度模拟方法
López‐Moreno et al. Response of snow processes to climate change: spatial variability in a small basin in the Spanish Pyrenees
Zhang et al. Simulating potential response of hydrology, soil erosion, and crop productivity to climate change in Changwu tableland region on the Loess Plateau of China
Zhou et al. Measurement and modelling of evapotranspiration over a reed (Phragmites australis) marsh in Northeast China
Fang et al. Snowmelt runoff sensitivity analysis to drought on the Canadian prairies
Jarrett et al. Regional interdisciplinary paleoflood approach to assess extreme flood potential
Jin et al. The long–term field–scale hydrology of subsurface drainage systems in a cold climate
Timmis et al. Frost injury prediction model for Douglas-fir seedlings in the Pacific Northwest
KR20160057911A (ko) 농업가뭄지수 도출 시스템 및 이를 이용한 농업가뭄지수 도출 방법
Pnevmatikos et al. The changing rainfall regime in Greece and its impact on climatological means
Zhou et al. Modeling of subsurface tile drainage using MIKE SHE
Sands et al. Comparing the subsurface drainage flow prediction of the DRAINMOD and ADAPT models for a cold climate
Flerchinger et al. Modeling soil freezing and thawing on a rangeland watershed
Meng et al. Spring flood forecasting based on the WRF-TSRM mode
Shankar et al. Irrigation scheduling for maize and Indian-mustard based on daily crop water requirement in a semi-arid region
Fulton et al. Evaluation of crop coefficients and evapotranspiration in English walnut
Sabbagh et al. Comparison of EPIC-WT and DRAINMOD simulated performance of land drainage systems
Bruland et al. Glacial mass balance of Austre Brøggerbreen (Spitsbergen), 1971–1999, modelled with a precipitation‐run‐off model
Anderson The response of Ka Roimata o Hine Hukatere Franz Josef Glacier to climate change
Aubry‐Wake et al. Predicting hydrological change in an alpine glacierized basin and its sensitivity to landscape evolution and meteorological forcings
Saxton et al. Climatic impacts on dryland winter wheat yields by daily soilwater and crop stress simulations
Bengtson et al. A fluctuating watertable model for the management of a controlled-drainage/subirrigation system
Moeletsi Agroclymatic characterization of Lesotho for dryland maize production
Gustafson Simulation of water discharge rates from a clay-till soil over a ten year period
Temesgen et al. Assessment of irrigation water potential and water requirements of selected crops in the Wabe-Shebelle River Basin, Ethiopia

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant