scholarly journals Development and Testing of a Frozen Soil Parameterization for Cold Region Studies

2007 ◽  
Vol 8 (4) ◽  
pp. 690-701 ◽  
Author(s):  
Xia Zhang ◽  
Shu Fen Sun ◽  
Yongkang Xue
Keyword(s):  
Phase Change ◽  
Water Flow ◽  
Liquid Water ◽  
Frozen Soil ◽  
Picard Iteration ◽  
Richards Equation ◽  
Cold Region ◽  
Field Station ◽  
Mixed Form ◽  
Ice Content

Abstract Proper simulation of soil freezing and thawing processes is an important issue in cold region climate studies. This paper reports on a frozen soil parameterization scheme for cold region studies that includes water flow and heat transfer in soil with water phase change. The mixed-form Richards’ equation is adopted to describe soil water flow affected by thermal processes in frozen soil. In addition, both liquid water and ice content have been taken into account in the frozen soil hydrologic and thermal property parameterization. To solve the complex nonlinear equation set and to ensure water conservation during simulation of complex phase change processes, efficient computational procedures have been designed and a new modified Picard iteration scheme is extended to solve the mixed-form Richards’ equation with phase change. The frozen soil model was evaluated using observational data from the field station at Rosemount, Minnesota, and the Tibet D66 site. The results show that the model is capable of providing good simulations of the evolution of temperature and liquid water content in frozen soil. Comparisons of simulation results from sensitivity studies indicate that there is a maximum difference of about 50 W m−2 in sensible and ground heat fluxes with and without the inclusion of the effect of ice content on matric potential and that using the exponential relationship between hydraulic conductivity and ice content produces realistic results.

scholarly journals Development of physically based liquid water schemes for Greenland firn-densification models

2019 ◽  
Vol 13 (7) ◽  
pp. 1819-1842 ◽  
Author(s):  
Vincent Verjans ◽  
Amber A. Leeson ◽  
C. Max Stevens ◽  
Michael MacFerrin ◽  
Brice Noël ◽  
...  
Keyword(s):  
Water Flow ◽  
Liquid Water ◽  
Preferential Flow ◽  
Greenland Ice Sheet ◽  
Single Domain ◽  
Richards Equation ◽  
Retention Capacity ◽  
Water Percolation ◽  
Physically Based ◽  
Dual Domain

Abstract. As surface melt is increasing on the Greenland Ice Sheet (GrIS), quantifying the retention capacity of the firn layer is critical to linking meltwater production to meltwater runoff. Firn-densification models have so far relied on empirical approaches to account for the percolation–refreezing process, and more physically based representations of liquid water flow might bring improvements to model performance. Here we implement three types of water percolation schemes into the Community Firn Model: the bucket approach, the Richards equation in a single domain and the Richards equation in a dual domain, which accounts for partitioning between matrix and fast preferential flow. We investigate their impact on firn densification at four locations on the GrIS and compare model results with observations. We find that for all of the flow schemes, significant discrepancies remain with respect to observed firn density, particularly the density variability in depth, and that inter-model differences are large (porosity of the upper 15 m firn varies by up to 47 %). The simple bucket scheme is as efficient in replicating observed density profiles as the single-domain Richards equation, and the most physically detailed dual-domain scheme does not necessarily reach best agreement with observed data. However, we find that the implementation of preferential flow simulates ice-layer formation more reliably and allows for deeper percolation. We also find that the firn model is more sensitive to the choice of densification scheme than to the choice of water percolation scheme. The disagreements with observations and the spread in model results demonstrate that progress towards an accurate description of water flow in firn is necessary. The numerous uncertainties about firn structure (e.g. grain size and shape, presence of ice layers) and about its hydraulic properties, as well as the one-dimensionality of firn models, render the implementation of physically based percolation schemes difficult. Additionally, the performance of firn models is still affected by the various effects affecting the densification process such as microstructural effects, wet snow metamorphism and temperature sensitivity when meltwater is present.

scholarly journals Effects of Frozen Soil on Snowmelt Runoff and Soil Water Storage at a Continental Scale

2006 ◽  
Vol 7 (5) ◽  
pp. 937-952 ◽  
Author(s):  
Guo-Yue Niu ◽  
Zong-Liang Yang
Keyword(s):  
Thermal Properties ◽  
Soil Water ◽  
Liquid Water ◽  
Water Storage ◽  
Frozen Soil ◽  
Soil Water Storage ◽  
Continental Scale ◽  
Community Land Model ◽  
Wide Range ◽  
Ice Content

Abstract The presence of ice in soil dramatically alters soil hydrologic and thermal properties. Despite this important role, many recent studies show that explicitly including the hydrologic effects of soil ice in land surface models degrades the simulation of runoff in cold regions. This paper addresses this dilemma by employing the Community Land Model version 2.0 (CLM2.0) developed at the National Center for Atmospheric Research (NCAR) and a simple TOPMODEL-based runoff scheme (SIMTOP). CLM2.0/SIMTOP explicitly computes soil ice content and its modifications to soil hydrologic and thermal properties. However, the frozen soil scheme has a tendency to produce a completely frozen soil (100% ice content) whenever the soil temperature is below 0°C. The frozen ground prevents infiltration of snowmelt or rainfall, thereby resulting in earlier- and higher-than-observed springtime runoff. This paper presents modifications to the above-mentioned frozen soil scheme that produce more accurate magnitude and seasonality of runoff and soil water storage. These modifications include 1) allowing liquid water to coexist with ice in the soil over a wide range of temperatures below 0°C by using the freezing-point depression equation, 2) computing the vertical water fluxes by introducing the concept of a fractional permeable area, which partitions the model grid into an impermeable part (no vertical water flow) and a permeable part, and 3) using the total soil moisture (liquid water and ice) to calculate the soil matric potential and hydraulic conductivity. The performance of CLM2.0/SIMTOP with these changes has been tested using observed data in cold-region river basins of various spatial scales. Compared to the CLM2.0/SIMTOP frozen soil scheme, the modified scheme produces monthly runoff that compares more favorably with that estimated by the University of New Hampshire–Global Runoff Data Center and a terrestrial water storage change that is in closer agreement with that measured by the Gravity Recovery and Climate Experiment (GRACE) satellites.

scholarly journals Development of physically based liquid water schemes for Greenland firn-densification models

2019 ◽  
Author(s):  
Vincent Verjans ◽  
Amber Leeson ◽  
C. Max Stevens ◽  
Michael MacFerrin ◽  
Brice Noël ◽  
...  
Keyword(s):  
Water Flow ◽  
Liquid Water ◽  
Preferential Flow ◽  
Greenland Ice Sheet ◽  
Single Domain ◽  
Richards Equation ◽  
Retention Capacity ◽  
Water Percolation ◽  
Physically Based ◽  
Dual Domain

Abstract. As surface melt is increasing on the Greenland ice sheet (GrIS), quantifying the retention capacity of the firn layer is critical to link meltwater production to meltwater runoff. Firn-densification models have so far relied on empirical approaches to account for the percolation-refreezing process, and more physically based representations of liquid water flow might therefore bring improvements to model performance. Here we implement three types of water percolation schemes into the Community Firn Model: the tipping bucket approach, the Richards Equation in a single-domain and the Richards Equation in a dual-domain, which accounts for partitioning between matrix and fast preferential flow. We investigate their impact on firn densification at four locations on the GrIS and compare model results with observations. We find that for all of the flow schemes, significant discrepancies remain with respect to observed firn density, particularly the density variability in depth, and that inter-model differences are large. The simple bucket scheme is as efficient in replicating observed density profiles as the single-domain Richards Equation. The most physically detailed dual-domain scheme does not necessarily reach best agreement with observed data. However, we find that the implementation of preferential flow does allow for more frequent ice layer formation and for deeper percolation. We also find that the firn model is more sensitive to the choice of densification scheme than to the choice of water percolation scheme. The disagreements with observations and the spread in model results demonstrate that progress towards an accurate description of water flow in firn is necessary. The numerous uncertainties surrounding firn micro- and macro-structure, its hydraulic properties, and the one dimensionality of firn models render the implementation of physically based percolation schemes difficult. An improved understanding of the parameters affecting evolution of polar firn, of the effects of the climatic forcing on the densification process and more accurate treatment of liquid water would benefit further developments.

scholarly journals Finite analytic method based on mixed-form Richards’ equation for simulating water flow in vadose zone

2016 ◽  
Vol 537 ◽  
pp. 146-156 ◽  
Author(s):  
Zaiyong Zhang ◽  
Wenke Wang ◽  
Tian-chyi Jim Yeh ◽  
Li Chen ◽  
Zhoufeng Wang ◽  
...  
Keyword(s):  
Water Flow ◽  
Vadose Zone ◽  
Richards Equation ◽  
Analytic Method ◽  
Mixed Form ◽  
Finite Analytic Method

Experimental Study on Thermal Characteristics of Finned Coil LHSU Using Paraffin as Phase Change Material

2017 ◽  
Vol 139 (4) ◽  
Author(s):  
Guansheng Chen ◽  
Nanshuo Li ◽  
Huanhuan Xiang ◽  
Fan Li
Keyword(s):  
Heat Transfer ◽  
Experimental Study ◽  
Phase Change ◽  
Phase Change Material ◽  
Water Flow ◽  
Water Flow Rate ◽  
Thermal Characteristics ◽  
Heat Transfer Fluid ◽  
Latent Heat Storage ◽  
Change Material

It is well known that attaching fins on the tubes surfaces can enhance the heat transfer into and out from the phase change materials (PCMs). This paper presents the results of an experimental study on the thermal characteristics of finned coil latent heat storage unit (LHSU) using paraffin as the phase change material (PCM). The paraffin LHSU is a rectangular cube consists of continuous horizontal multibended tubes attached vertical fins at the pitches of 2.5, 5.0, and 7.5 mm that creates the heat transfer surface. The shell side along with the space around the tubes and fins is filled with the material RT54 allocated to store energy of water, which flows inside the tubes as heat transfer fluid (HTF). The measurement is carried out under four different water flow rates: 1.01, 1.30, 1.50, and 1.70 L/min in the charging and discharging process, respectively. The temperature of paraffin and water, charging and discharging wattage, and heat transfer coefficient are plotted in relation to the working time and water flow rate.

scholarly journals Modelling Unsaturated Flow in Porous Media Using an Improved Picard Iteration Scheme

2021 ◽  
Author(s):  
S.R. Zhu ◽  
L.Z. Wu ◽  
T. Ma ◽  
S.H. Li
Keyword(s):  
Porous Media ◽  
Convergence Rate ◽  
Unsaturated Flow ◽  
Control Volume ◽  
Linear Equations ◽  
Solution Process ◽  
Coefficient Matrix ◽  
Picard Iteration ◽  
Flow In Porous Media ◽  
Richards Equation

Abstract The numerical solution of various systems of linear equations describing fluid infiltration uses the Picard iteration (PI). However, because many such systems are ill-conditioned, the solution process often has a poor convergence rate, making it very time-consuming. In this study, a control volume method based on non-uniform nodes is used to discretize the Richards equation, and adaptive relaxation is combined with a multistep preconditioner to improve the convergence rate of PI. The resulting adaptive relaxed PI with multistep preconditioner (MP(m)-ARPI) is used to simulate unsaturated flow in porous media. Three examples are used to verify the proposed schemes. The results show that MP(m)-ARPI can effectively reduce the condition number of the coefficient matrix for the system of linear equations. Compared with conventional PI, MP(m)-ARPI achieves faster convergence, higher computational efficiency, and enhanced robustness. These results demonstrate that improved scheme is an excellent prospect for simulating unsaturated flow in porous media.

Ice or rock matrix? Improved quantitative imaging of Alpine permafrost evolution through time-lapse petrophysical joint inversion

2021 ◽  
Author(s):  
Johanna Klahold ◽  
Christian Hauck ◽  
Florian Wagner
Keyword(s):  
Liquid Water ◽  
Joint Inversion ◽  
Time Lapse ◽  
Swiss Alps ◽  
Permafrost Degradation ◽  
Rock Matrix ◽  
Validation Data ◽  
Water Ice ◽  
Borehole Temperature ◽  
Ice Content

<p>Quantitative estimation of pore fractions filled with liquid water, ice and air is one of the prerequisites in many permafrost studies and forms the basis for a process-based understanding of permafrost and the hazard potential of its degradation in the context of global warming. The volumetric ice content is however difficult to retrieve, since standard borehole temperature monitoring is unable to provide any ice content estimation. Geophysical methods offer opportunities to image distributions of permafrost constituents in a non-invasive manner. A petrophysical joint inversion was recently developed to determine volumetric water, ice, air and rock contents from seismic refraction and electrical resistivity data. This approach benefits from the complementary sensitivities of seismic and electrical data to the phase change between ice and liquid water. A remaining weak point was the unresolved petrophysical ambiguity between ice and rock matrix. Within this study, the petrophysical joint inversion approach is extended along the time axis and respective temporal constraints are introduced. If the porosity (and other time-invariant properties like pore water resistivity or Archie exponents) can be assumed invariant over the considered time period, water, ice and air contents can be estimated together with a temporally constant (but spatially variable) porosity distribution. It is hypothesized that including multiple time steps in the inverse problem increases the ratio of data and parameters and leads to a more accurate distinction between ice and rock content. Based on a synthetic example and a field data set from an Alpine permafrost site (Schilthorn, Swiss Alps) it is demonstrated that the developed time-lapse petrophysical joint inversion provides physically plausible solutions, in particular improved estimates for the volumetric fractions of ice and rock. The field application is evaluated with independent validation data including thaw depths derived from borehole temperature measurements and shows generally good agreement. As opposed to the conventional petrophysical joint inversion, its time-lapse extension succeeds in providing reasonable estimates of permafrost degradation at the Schilthorn monitoring site without <em>a priori </em>constraints on the porosity model.</p>

Vadose Zone, Part I : Characterization and Flow Processes

2003 ◽  
Author(s):  
Yoram Rubin
Keyword(s):  
Porous Media ◽  
Water Flow ◽  
Governing Equation ◽  
Vadose Zone ◽  
Stochastic Analysis ◽  
Transport Processes ◽  
Richards Equation ◽  
Saturated Porous Media ◽  
Flow And Transport ◽  
Flow Processes

Many of the principles guiding stochastic analysis of flow and transport processes in the vadose zone are those which we also employ in the saturated zone, and which we have explored in earlier chapters. However, there are important considerations and simplifications to be made, given the nature of the flow and of the governing equations, which we explore here and in chapter 12. The governing equation for water flow in variably saturated porous media at the smallest scale where Darcy’s law is applicable (i.e., no need for upscaling of parameters) is Richards’ equation (cf. Yeh, 1998)

Sign in / Sign up

Export Citation Format

Share Document