Method and system for calculating the liquefaction water content of frozen soil in highland areas.

Abstract

This invention provides a method and system for calculating the liquid water content of permafrost in highland soil. [Solution] The method includes: step S1 acquiring a soil temperature and humidity observation network dataset from the central part of the plateau; step S2 performing texture classification on the soil to obtain soil texture classification results; step S3 combining the soil texture classification results and the soil temperature and humidity observation network dataset to fit the relationship between solid water content and soil temperature for soils with different textures; step S4 determining a formula for calculating the liquid water content of the soil based on the relationship between solid water content and soil temperature; step S5 improving the VIC model based on the formula for calculating the liquid water content; and step S6 calculating the liquid water content of the plateau permafrost using the improved VIC model.

Description

[Technical Field]

[0001] The present invention relates to the technical field of hydrology in the cryosphere, and more particularly to a method and system for calculating the liquid water content of plateau permafrost. [Background technology]

[0002] The calculation method for the liquefied water content of plateau permafrost refers to a specific method for measuring the liquefied water content in plateau permafrost layers. Generally, this involves considering the multilayer structure of the soil and the dynamics of moisture in different soil types in combination with elements such as meteorological data and surface characteristics, and simulating precipitation, evaporation, thawing, and moisture movement within the permafrost layer to calculate the liquefied water content in the permafrost.

[0003] The river headwaters of southwestern China are high-altitude, cold regions with extensive distribution of seasonal permafrost and permafrost. Against the backdrop of climate warming, in addition to the supplementation of river flow by glacial meltwater, the freeze-thaw process of permafrost also affects river runoff. Calculating solid and liquid water in the soil is an important part of calculating runoff and combined flow rates in hydrological models. Variable infiltration capacity (VIC) hydrological models are preferred for analyzing the process of supplementing river flow by permafrost meltwater in high-altitude, cold regions.

[0004] However, when calculating the liquid water content of soil using conventional VIC models, it relies on several parameters with high uncertainty (e.g., saturated water content, latent heat flux, etc.). The physical meaning of these parameters is not clear, and the state transitions between solid and liquid water in soil are complex and frequent. Furthermore, there are phenomena in the model where using different parameters yields the same effect, increasing the uncertainty of the calculation. This makes the calculation process for the permafrost temperature field complicated and time-consuming, reducing the efficiency of the model and making it difficult to achieve efficient and accurate calculation of liquid water content in the complex permafrost environment of plateaus. [Overview of the project] [Problems that the invention aims to solve]

[0005] When calculating the liquid water content of soil using the conventional VIC model, it relies on several parameters with high uncertainty (e.g., saturated water content, latent heat flux, etc.). However, the physical meaning of these parameters is not clear, the state transitions between solid water and liquid water in soil are complex and frequent, and there are phenomena in the model where using different parameters yields the same effect, thus increasing the uncertainty of the calculation. This makes the calculation process for the frozen soil temperature field complicated and time-consuming, reducing the efficiency of the model. To solve the technical problems of achieving efficient and accurate calculation of liquid water content in the complex frozen soil environment of plateaus, the present invention provides a method and system for calculating the liquid water content of plateau frozen soil. [Means for solving the problem]

[0006] The technical solutions provided in the embodiments of the present invention are as follows:

[0007] In the first embodiment, The method for calculating the liquid water content of plateau frozen soil provided in the embodiments of the present invention is: Step S1 involves obtaining a dataset from a soil temperature and humidity observation network in the central part of the plateau. Step S2 involves performing texture classification on the soil and obtaining the soil texture classification result. Step S3 involves combining the soil texture classification results with the soil temperature and humidity observation network dataset to fit the relationship between solid water content and soil temperature for soils with different textures. Step S4 involves determining the formula for calculating the liquid water content of the soil based on the relationship between the solid water content and soil temperature, Step S5 involves improving the VIC model based on the formula for calculating the liquid water content, The improved VIC model includes step S6, which calculates the liquefaction water content of the plateau permafrost.

[0008] In the second embodiment, The calculation system for the liquid water content of plateau frozen soil provided in the embodiment of the present invention is: Processor and The system includes a memory in which computer-readable instructions are stored, and when the computer-readable instructions are executed by the processor, a method for calculating the liquid water content of plateau frozen soil in a first embodiment is realized.

[0009] In the third embodiment, Embodiments of the present invention provide a computer-readable storage medium in which a computer program is stored, and when the program is executed by a processor, the method for calculating the liquid water content of plateau frozen soil described in the first embodiment is realized. [Effects of the Invention]

[0010] The beneficial effects of the technical solutions provided in the embodiments of the present invention include at least the following:

[0011] In this invention, by utilizing actual soil temperature and humidity observation data from the central plateau in combination with soil texture classification results, and by closely integrating regional characteristics with model calculations, the accuracy of liquid water content calculations is significantly improved. Furthermore, by making the calculation of liquid water content more scientific and rigorous based on a fitting equation between solid water content and temperature, and by introducing a formula for calculating liquid water content instead of the formula for unfrozen water content in the conventional VIC model, the number of uncertain parameters on which the formula for calculating unfrozen water content in the VIC model depends is reduced, the algorithm for unfrozen water content is optimized, the accuracy of the calculation results for unfrozen water content is improved, and the hydrological simulation results of the VIC model are effectively improved. [Brief explanation of the drawing]

[0012] To more clearly explain the technical solutions of the embodiments of the present invention, the drawings necessary for describing the embodiments are briefly described below. Naturally, the drawings described below represent only a few embodiments of the present invention, and those skilled in the art can obtain other drawings based on these without requiring any creative effort. [Figure 1] It is a schematic flowchart of a method for calculating the liquid water content of permafrost on the Qinghai-Tibet Plateau provided in an embodiment of the present invention. [Figure 2a] It is a relationship diagram between the solid water content rate of soil and soil temperature provided in an embodiment of the present invention. [Figure 2b] It is a relationship diagram between the solid water content rate of soil and soil temperature provided in an embodiment of the present invention. [Figure 3] It is a structural schematic diagram of a calculation system for the liquid water content of permafrost on the Qinghai-Tibet Plateau provided in an embodiment of the present invention.

Embodiments for Carrying Out the Invention

[0013] Hereinafter, the technical solution of the present invention will be described while referring to the drawings.

[0014] In the embodiments of the present disclosure, terms such as "exemplary", "for example", etc. are used to indicate examples, illustrations or explanations. In the present invention, no embodiment or design described as "exemplary" should be construed as having more preferable or superior advantages than other embodiments or designs. Strictly speaking, the use of the term "exemplary" is intended to present concepts specifically. Also, in the embodiments of the present invention, the meaning represented by "and / or" may be that both exist together or either one may be selected.

[0015] In the embodiments of the present invention, "image" and "picture" may be used interchangeably. It should be understood that when the distinction is not emphasized, the meanings they intend to express are the same. "Of", "corresponding", and "corresponding" may be used interchangeably. It should be understood that when the distinction is not emphasized, the meanings they intend to express are the same.

[0016] In embodiments of the present invention, subscripts such as W1 may be written in a form other than a subscript such as W1, and if the distinction is not emphasized, the meaning they intend to express is the same.

[0017] The technical problems, technical solutions, and advantages of this invention will be described in detail below with reference to the drawings and specific embodiments, in order to clarify them further.

[0018] Referring to Figure 1 of the specification, a schematic flowchart of the method for calculating the liquid water content of the Qinghai-Tibet Plateau permafrost provided in an embodiment of the present invention is shown.

[0019] The method for calculating the liquid water content of the frozen soil in the Qinghai-Tibet Plateau, as provided in the embodiments of the present invention, includes S1 to S6.

[0020] In S1, we obtain the soil temperature and humidity observation network dataset for the central Qinghai-Tibet Plateau.

[0021] Furthermore, by acquiring the soil temperature and humidity observation network dataset for the central Qinghai-Tibet Plateau, we will provide high-resolution spatiotemporal data, accurately reflect the regional soil temperature and humidity dynamics, support comparative studies under different geological and climatic conditions, provide reliable experimental evidence for model improvement and validation of calculation results, and improve calculation accuracy.

[0022] In S2, texture classification is performed on the soil to obtain the soil texture classification results.

[0023] In one possible embodiment, S2 is specifically as follows:

[0024] We adopted the International Soil Texture Classification Standard and considered a combination of clay, silt, sand particles, and porosity to perform texture classification on soils in the central Qinghai-Tibet Plateau and obtain soil texture classification results.

[0025] Furthermore, using international standards ensures the comparability of soil classification results from different regions, making it versatile. The proportions of clay, silt, and sand particles accurately describe the physical properties of the soil, such as hydrology and nutrient retention capacity, providing a scientific basis for fields such as land use and agricultural production. Considering porosity helps evaluate soil permeability and water retention, which is particularly important for hydrological studies in permafrost regions. A detailed soil texture classification allows for better simulation and prediction of changes in soil moisture under different climatic conditions.

[0026] Referring to Figure 2 of the specification, a diagram is shown illustrating the relationship between the solid water content of the soil and the soil temperature in the embodiment of the present invention.

[0027] Figure 2 shows the relationship between solid water content and soil temperature for two types of soil (ALi's loam soil and Maqu's loam soil). Here, the horizontal axis represents soil temperature (°C), and the vertical axis represents the solid water content of the soil. As can be seen from the figure, the solid water content of both types of soil tends to decrease with increasing soil temperature. This indicates that ice in the soil gradually melts into liquid water. The range of solid water content in Ali's loam soil is wide, and the magnitude of change with temperature changes is large, while the change in solid water content in Maqu's loam soil is relatively stable, and the effect of temperature on it is relatively small. This reflects the difference in moisture states of different soil textures during the freeze-thaw process.

[0028] In actual operations, the coefficient of influence of temperature change on solid water content was -0.01536 for aliloam soil and -0.01547 for machuloam soil, showing only a slight difference, indicating that the influence of temperature change on the two types of soil is similar. For the regression constant, it was 0.3573 for aliloam soil and 0.3371 for machuloam soil, showing only a slight difference. For the degree of fitting of the regression model, it was 0.8607 for both, indicating that the model fitting effect is good. For the root mean square error used to evaluate the prediction accuracy of the model, it was 0.02799 for aliloam soil and 0.01289 for machuloam soil, indicating that the prediction accuracy of machuloam soil is higher.

[0029] In S3, the soil texture classification results and the soil temperature and humidity observation network dataset are combined to fit the relationship between solid water content and soil temperature for soils with different textures.

[0030] Furthermore, by combining soil texture data with actual temperature and humidity data, it is possible to accurately fit the relationship between solid water content and soil temperature to different soil types, thereby improving the accuracy of hydrological models for permafrost regions. By combining this with data from temperature and humidity observation networks, it is possible to analyze changes in soil temperature and humidity at a finer spatial scale, providing highly directional environmental change assessments. This effectively supports the analysis of the laws governing changes in moisture content in permafrost layers and provides important support for assessing the impact of climate change on permafrost regions.

[0031] In S4, the formula for calculating the liquid water content of the soil is determined based on the relationship between the solid water content and soil temperature.

[0032] Here, solid water content refers to the amount of solid water (ice) in the soil, and liquid water content refers to the amount of liquid water in the soil.

[0033] In one possible embodiment, S4 specifically includes S401 to S404.

[0034] In S401, the relationship between the total amount of solid and liquid water in the soil and the soil temperature is fitted based on a general-purpose ideal gas law.

[0035] Here, the general-purpose ideal gas law is an equation that describes the relationship between the pressure, volume, temperature, and amount of a gas.

[0036] Furthermore, fitting the relationship between moisture and temperature using physical formulas has relatively strong scientific basis and versatility, allowing for the quantification of the laws governing the change in the total amount of solid and liquid water in soil. This provides a reliable theoretical foundation for calculating liquid water content, enhances the dynamic adaptability of the model, and facilitates the analysis of moisture changes in permafrost under different climatic conditions.

[0037] In S402, the sequence of solid water content in different soil layers is calculated by combining the relationship between solid water content and soil temperature, based on the total amount of solid and liquid water in the soil.

[0038] Here, the sum of solid and liquid water in the soil refers to the total amount of solid and liquid water in the soil, the solid water content is the ratio of the mass of solid water in the soil to the soil mass, and the solid water content sequence is a data sequence that describes the change in water content in a substance and reflects the tendency to change over time or with other variables (e.g., temperature, pressure, or environmental conditions).

[0039] Furthermore, by comprehensively considering the relationship between solid water, liquid water, and temperature, it is possible to accurately calculate the solid water content of different soil layers. This provides quantifiable data for in-depth analysis of the moisture distribution and dynamic changes in permafrost layers, thereby improving the accuracy of soil moisture models.

[0040] In S403, the relationship between soil ice content and soil temperature is fitted based on the solid water content sequence.

[0041] In S404, a calculation formula for the liquid water content rate of the soil is determined by combining the relationship between the total of solid water and liquid water in the soil and the soil temperature, and the relationship between the ice content of the soil and the soil temperature.

[0042] In addition, in combination with the relationship between the moisture and temperature of the soil, the dynamic changes of liquid water in frozen soil can be more accurately described, thereby providing an accurate calculation method for liquid water and being suitable for hydrological analysis under different soil layers and temperature conditions.

[0043] In one possible embodiment, the equation of state for a general gas is specifically PV =nRT, where P represents the pressure of the gas in the soil, V a represents the volume of the gas in the soil, n represents the amount of substance of the gas in the soil, R represents the Avogadro constant, and T represents the soil temperature at the corresponding depth.

[0044] In one possible embodiment, the relationship between the total of solid water and liquid water in the soil and the soil temperature is specifically V a =V - V s - V<00000 or 5>- V i =[(SR) / (gM)]T, V i + V<000000 -> 8>=V - V s -[(SR) / (gM)]T, where V a represents the volume of the gas in the soil, V represents the total volume of the soil, V s represents the volume of the soil particles, V w represents the volume of the liquid water in the soil, V i represents the volume of the solid water in the soil, S represents the cross-sectional area of the soil, T represents the soil temperature at the corresponding depth, and M represents the molar mass of air.

[0045] In one possible embodiment, the relationship between the ice content of the soil and the soil temperature is specifically V i_sim =kT + b, where Vi_sim k represents the ice content of the soil, k represents the slope in the relationship between soil ice content and soil temperature, and b represents the intercept in the relationship between soil ice content and soil temperature.

[0046] In one possible embodiment, the formula for calculating the liquefied water content of the soil is, specifically, V w =VV s -[(S×R) / (g×M)]×TV i_sim、 V w =VV s -[(S×R) / (g×M)]×T-(k*T+b), Here, V represents the total volume of the soil, s V represents the volume of soil particles. w This represents the volume of liquid water in the soil, V i_sim θ represents the ice content of the soil, S represents the cross-sectional area of ​​the soil, R represents Avogadro's constant, T represents the soil temperature at the corresponding depth, g represents the acceleration due to gravity, and M represents the molar mass of air.

[0047] Furthermore, the formula for deriving the liquid water content based on the relationship between solid water content and soil temperature can accurately simulate the moisture state of the soil at different temperatures, thereby improving the adaptability and accuracy of hydrological models for permafrost regions.

[0048] In S5, the VIC model is improved based on the formula for calculating the liquid water content.

[0049] In one possible embodiment, S5 is specifically as follows:

[0050] The VIC model is improved by replacing the formula for calculating the unfrozen water content with the formula for calculating the liquefied water content of the soil.

[0051] Furthermore, by introducing an accurate formula for calculating liquid water content into the VIC model, the accuracy and efficiency of the model in simulating hydrological dynamics in permafrost regions can be significantly improved, allowing the VIC model to more accurately capture changes in soil moisture conditions due to temperature changes.

[0052] In S6, the improved VIC model is used to calculate the liquefaction water content of the Qinghai-Tibet Plateau permafrost.

[0053] Furthermore, by calculating the liquid water content of the Qinghai-Tibet Plateau permafrost using the improved VIC model, it becomes possible to more accurately simulate the water dynamics of the permafrost region and effectively predict changes in liquid water during the soil freeze-thaw process, thereby improving the accuracy and reliability of the regional hydrological model.

[0054] The beneficial effects of the technical solutions provided in the embodiments of the present invention include at least the following:

[0055] In this invention, by utilizing actual soil temperature and humidity observation data from the central Qinghai-Tibet Plateau in combination with soil texture classification results, and by closely integrating regional characteristics with model calculations, the accuracy of liquid water content calculations is significantly improved. Furthermore, by making the calculation of liquid water content more scientific and rigorous based on a fitting equation between solid water content and temperature, and by introducing a formula for calculating liquid water content instead of the formula for unfrozen water content in the conventional VIC model, the number of uncertain parameters on which the formula for calculating unfrozen water content in the VIC model depends is reduced, the algorithm for unfrozen water content is optimized, the accuracy of the calculation results for unfrozen water content is improved, and the hydrological simulation results of the VIC model are effectively improved.

[0056] Referring to Figure 3 of the specification, a schematic diagram of the structure of the calculation system for the liquid water content of the frozen soil in the Qinghai-Tibet Plateau provided by the present invention is shown.

[0057] The present invention further provides a calculation system 20 for the liquid water content of permafrost on the Qinghai-Tibet Plateau, which is applicable to the above-mentioned method for calculating the liquid water content of permafrost on the Qinghai-Tibet Plateau. Processor 201 and, The method for calculating the liquid water content of the Qinghai-Tibet Plateau permafrost described in the embodiment of the method is realized when the computer-readable instructions are executed by the processor 201, which includes a memory 202 in which computer-readable instructions are stored.

[0058] The Qinghai-Tibet Plateau permafrost liquid water content calculation system 20 provided in the present invention can perform the above-described method for calculating the Qinghai-Tibet Plateau permafrost liquid water content and achieves the same or similar technical effects. Therefore, to avoid duplication, the present invention will not describe it again.

[0059] The beneficial effects of the technical solutions provided in the embodiments of the present invention include at least the following:

[0060] In this invention, by utilizing actual soil temperature and humidity observation data from the central Qinghai-Tibet Plateau in combination with soil texture classification results, and by closely integrating regional characteristics with model calculations, the accuracy of liquid water content calculations is significantly improved. Furthermore, by making the calculation of liquid water content more scientific and rigorous based on a fitting equation between solid water content and temperature, and by introducing a formula for calculating liquid water content instead of the formula for unfrozen water content in the conventional VIC model, the number of uncertain parameters on which the formula for calculating unfrozen water content in the VIC model depends is reduced, the algorithm for unfrozen water content is optimized, the accuracy of the calculation results for unfrozen water content is improved, and the hydrological simulation results of the VIC model are effectively improved.

[0061] In embodiments of the present invention, the processor may be a Central Processing Unit (CPU), and it should be understood that the processor may further be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field programmable gate arrays (FPGAs), or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any general-purpose processor, etc.

[0062] It should be understood that the memory in the embodiments of the present invention may be volatile memory or non-volatile memory, or may include both volatile and non-volatile memory. Here, the non-volatile memory may be read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), or flash memory. The volatile memory may be random access memory (RAM), which is used as external cache memory. Exemplary examples include, but are not limited to, many forms of random access memory (RAM), such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchlink dynamic random access memory (SLDRAM), and direct rambus random access memory (DR RAM).

[0063] All or part of the above embodiments can be implemented using software, hardware (e.g., circuitry), firmware, or any other combination. When implemented using software, all or part of the above embodiments can be implemented in the form of a computer program product. A computer program product includes one or more computer instructions or computer programs. When a computer instruction or computer program is loaded or executed on a computer, all or part of a flow or function according to an embodiment of the present invention is generated. The computer may be a general-purpose computer, a dedicated computer, a computer network, or other programmable device. Computer instructions may be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired method (e.g., infrared, wireless, microwave, etc.). The computer-readable storage medium may be any available medium accessible to the computer, or it may be a data storage device such as a server or data center that includes one or more sets of available media. The usable media may be magnetic media (e.g., floppy disks, hard disks, magnetic tapes), optical media (e.g., DVDs), or semiconductor media. Semiconductor media may be solid-state disks.

[0064] The terms "and / or" in this specification merely describe the relationship between related objects, indicating that there are three types of relationships. For example, A and / or B can represent three cases: when A exists alone, when A and B exist together, and when B exists alone, where A and B may be singular or plural. Furthermore, the symbol " / " in this specification generally indicates that the preceding and following related objects are in an "or" relationship, but it can also indicate an "and / or" relationship, and its specific meaning must be understood by referring to the context.

[0065] In this invention, "at least one" means one or more, and "multiple" means two or more. "At least one of the following" or similar expressions means any combination of these items, including any combination of single items (one) or multiple items (multiple). For example, the expression "at least one of a, b, and c" can represent a, b, c, ab, ac, bc, or abc, where a, b, and c may be one or more.

[0066] In each embodiment of the present invention, the magnitude of the number of each process described above does not indicate the order of execution. The execution order of each process should be determined based on its function and inherent logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0067] Those skilled in the art will recognize that the units and algorithmic steps of each example described with reference to the embodiments disclosed herein can be implemented in electronic hardware or in combination of computer software and electronic hardware. Whether these functions are performed in hardware or software will depend on the specific application and design constraints of the technical solution. Those skilled in the art will recognize that the functions described can be implemented in different ways for each specific application, but such implementations should not be considered to be beyond the scope of the present invention.

[0068] As will be obvious to those skilled in the art, for the sake of convenience and simplicity of explanation, the specific operating processes of the devices, apparatus, and units described above can be found by referring to the corresponding processes in the previously described method embodiments, and therefore a detailed explanation is omitted here.

[0069] It should be understood that the devices, apparatus, and methods disclosed in some embodiments of the present invention can also be implemented in other forms. For example, the embodiments of the apparatus described above are merely illustrative, and the division of units is merely a logical functional division, and in practice, they can be divided in other forms, for example, multiple units or components may be combined, integrated into another device, or may not reflect or perform some features. Furthermore, the coupling, direct coupling, or communication connection shown or discussed may be an indirect coupling or communication connection via several interfaces, devices, or units, and may be in an electrical, mechanical, or other form.

[0070] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units; that is, they may be located in the same place or distributed across multiple network units. Depending on the requirements, only some of these units may be used, or all of them may be used to achieve the objectives of this embodiment.

[0071] Furthermore, each functional unit in each embodiment of the present invention may be integrated into a single processing unit, each unit may exist as a separate physical component, or two or more units may be integrated into a single unit.

[0072] The functions may be implemented in the form of software function units and, if sold or used as independent products, may be stored in computer-readable storage media. Based on this understanding, the spirit of the technical solutions of the present invention, in other words, the part that contributes to the prior art or a part of the technical solutions, may be embodied in the form of a computer software product, which is stored on a storage medium and includes several commands for causing a computer device (which may be a personal computer, server, or network device, etc.) to perform all or part of the steps of each embodiment of the present invention. The aforementioned storage mediums include various media capable of storing program code, such as USB flash memory, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or compact disks.

[0073] An embodiment of the present invention provides a computer-readable storage medium in which a computer program is stored, and when the program is executed by a processor, the method for calculating the liquid water content of the Qinghai-Tibet Plateau frozen soil described in the embodiment of the method is realized.

[0074] The computer-readable storage medium provided in this invention can achieve the steps and effects of the method for calculating the liquid water content of the Qinghai-Tibet Plateau frozen soil described in the above-described embodiment of the method, and will not be described again in this invention to avoid duplication.

[0075] The beneficial effects of the technical solutions provided in the embodiments of the present invention include at least the following:

[0076] In this invention, by utilizing actual soil temperature and humidity observation data from the central Qinghai-Tibet Plateau in combination with soil texture classification results, and by closely integrating regional characteristics with model calculations, the accuracy of liquid water content calculations is significantly improved. Furthermore, by making the calculation of liquid water content more scientific and rigorous based on a fitting equation between solid water content and temperature, and by introducing a formula for calculating liquid water content instead of the formula for unfrozen water content in the conventional VIC model, the number of uncertain parameters on which the formula for calculating unfrozen water content in the VIC model depends is reduced, the algorithm for unfrozen water content is optimized, the accuracy of the calculation results for unfrozen water content is improved, and the hydrological simulation results of the VIC model are effectively improved.

[0077] The above are merely specific embodiments of the present invention, and the scope of protection of the present invention is not limited thereto. Any modifications or substitutions that a person skilled in the art can easily conceive within the scope of the art disclosed herein should be included within the scope of protection of the present invention. Accordingly, the scope of protection of the present invention is subject to the content of the claims.

[0078] The following points need to be explained.

[0079] (1) The drawings of embodiments of the present invention show only the structures relating to embodiments of the present invention, and other structures can be referenced to general designs.

[0080] (2) For clarity, in the drawings illustrating embodiments of the present invention, the thickness of layers or regions is enlarged or reduced, meaning these drawings are not drawn to actual scale. When a component such as a layer, film, region, or substrate is said to be "on top of" or "below" another component, it should be understood that the component may be located "directly" on top of or below the other component, or an intermediate component may be present.

[0081] (3) New embodiments can be obtained by combining the embodiments and features of the embodiments, as long as they do not contradict each other.

[0082] The above are merely specific embodiments of the present invention, and the scope of protection of the present invention is not limited thereto, but is subject to the content of the claims.

Claims

1. Step S1 involves obtaining a dataset from a soil temperature and humidity observation network in the central part of the plateau. Step S2 involves performing texture classification on the soil and obtaining the soil texture classification result. Step S3 involves combining the soil texture classification results and the soil temperature and humidity observation network dataset to fit the relationship between the solid water content and soil temperature of soils with different textures. Step S4 involves determining a formula for calculating the liquid water content of the soil based on the relationship between the solid water content and the soil temperature. Step S5 involves improving the VIC model based on the formula for calculating the liquid water content, The improved VIC model includes step S6, which calculates the liquid water content of the plateau permafrost, Step S4 is, Step S401 involves fitting the relationship between the total amount of solid and liquid water in the soil and the soil temperature based on a general-purpose ideal gas law, Step S402 calculates a sequence of solid water content in different soil layers, based on the total amount of solid water and liquid water in the soil, in combination with the relationship between the solid water content and the soil temperature. Step S403 involves fitting the relationship between the ice content of the soil and the soil temperature based on the solid water content sequence, Step S404 includes determining a formula for calculating the liquid water content of the soil by combining the relationship between the total amount of solid water and liquid water in the soil and the soil temperature, and the relationship between the ice content of the soil and the soil temperature. Step S5 is, The VIC model is improved by replacing the formula for calculating the unfrozen water content in the VIC model with the formula for calculating the liquid water content of the soil. A method for calculating the liquid water content of plateau permafrost, characterized by the following features.

2. The method for calculating the liquid water content of frozen soil in a plateau according to claim 1, characterized in that step S2 involves adopting an international soil texture grading standard, considering a combination of clay, silt, sand particles, and porosity, performing a texture classification on the soil in the central part of the plateau, and obtaining a soil texture classification result.

3. The aforementioned general-purpose ideal gas law is: PV ａ = nRT, Here, P represents the pressure of the gas in the soil, and V ａ The method for calculating the liquid water content of plateau frozen soil according to claim 1, characterized in that n represents the volume of gas in the soil, n represents the amount of substance of gas in the soil, R represents the Avogadro constant, and T represents the soil temperature at the corresponding depth.

4. The relationship between the total amount of solid and liquid water in the soil and the soil temperature is as follows: V ａ =V-V ｓ -V ｗ -V ｉ =[(SR) / (gM)]T、 V ｉ +V ｗ =V-V ｓ - Here, V ａ V represents the volume of gas in the soil, and V represents the total volume of the soil. ｓ V represents the volume of soil particles. ｗ V represents the volume of liquid water in the soil. ｉ The method for calculating the liquid water content of plateau frozen soil according to claim 1, characterized in that: represents the volume of solid water in the soil, S represents the cross-sectional area of ​​the soil, T represents the soil temperature at the corresponding depth, R represents the Avogadro constant, g represents the acceleration due to gravity, and M represents the molar mass of air.

5. The relationship between the ice content of the soil and the soil temperature is, V ｉ＿ｓｉｍ = kT + b, Here, V ｉ＿ｓｉｍ The method for calculating the liquid water content of frozen soil in a plateau according to claim 1, characterized in that k represents the ice content of the soil, k represents the slope in the relationship between the ice content of the soil and the soil temperature, b represents the intercept in the relationship between the ice content of the soil and the soil temperature, and T represents the soil temperature at the corresponding depth.

6. The formula for calculating the liquid water content of the aforementioned soil is: V ｗ =V-V ｓ -[(S×R) / (g×M)]×T-V ｉ＿ｓｉｍ 、 V ｗ = V - V ｓ -[(S × R) / (g × M)] × T - (k * T + b), Here, V represents the total volume of the soil, ｓ V represents the volume of soil particles. ｗ V represents the volume of liquid water in the soil. ｉ＿ｓｉｍ The method for calculating the liquid water content of frozen soil in a plateau according to claim 1, characterized in that: is the ice content of the soil, S is the cross-sectional area of ​​the soil, R is the Avogadro constant, T is the soil temperature at the corresponding depth, g is the acceleration due to gravity, M is the molar mass of air, k is the slope in the relationship between the ice content of the soil and the soil temperature, and b is the intercept in the relationship between the ice content of the soil and the soil temperature.

7. Processor and A memory containing computer-readable instructions is stored, and when the computer-readable instructions are executed by the processor, the method for calculating the liquid water content of plateau frozen soil according to any one of claims 1 to 6 is realized. A system for calculating the liquid water content of permafrost in highland areas, characterized by the following features.

8. When a computer program is stored and executed by a processor, the method for calculating the liquid water content of plateau frozen soil according to any one of claims 1 to 6 is realized. A computer-readable storage medium characterized by the following features.

Images