Numerical modelling of permafrost in bedrock in northern Fennoscandia during the Holocene

The occurrence of permafrost in bedrock in northern Fennoscandia and its dependence on past and presently ongoing climatic variations was investigated with one- (1D) and two-dimensional (2D) numerical models by solving the transient heat conduction equation with latent heat effects included. The study area is characterized by discontinuous permafrost occurrences such as palsa mires and local mountain permafrost. The ground temperature changes during the Holocene were constructed using climatic proxy data. This variation was used as a forcing function at the ground surface in the calculations. Several versions of the present ground temperature were applied, resulting in different subsurface freezing–thawing conditions in the past depending on the assumed porosity and geothermal conditions.Our results suggest that in high altitude areas with a cold climate (present mean annual ground temperature between 0°C and −3°C), there may have been considerable variations in permafrost thickness (ranging from 0 to 150 m), as well as periods of no permafrost at all. The higher is the porosity of bedrock filled with ice, the stronger is the retarding effect of permafrost against climatic variations.Two-dimensional models including topographic effects with altitude-dependent ground temperatures and slope orientation and inclination dependent solar radiation were applied to a case of mountain permafrost in Ylläs, western Finnish Lapland, where bedrock permafrost is known to occur in boreholes to a depth of about 60 m. Modelling suggests complicated changes in permafrost thickness with time as well as contrasting situations on southern and northern slopes of the mountain.Extrapolating the climatic warming of the last 200 years to the end of the next century when the anticipated increase in the annual average air temperature is expected to be about 2 K indicates that the permafrost occurrences in bedrock in northern Fennoscandia would be thawing rapidly in low-porosity formations. However, already a porosity of 5% filled with ice would retard the thawing considerably.     

The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

鲫鱼(Carassius auratus)在微冻保鲜过程中的质量变化

该文研究鲫鱼 ( Carassius auratus)在微冻 ( - 3℃ )保鲜过程中的质量变化规律。以感官评价、细菌总数、T- VBN值、p H值、ATPase活性、K值等作为质量指标。实验结果表明 ,微冻可以明显抑制细菌总数的增长 ,维持较低的 T- VBN和 K值。但是 ,微冻条件下鲫鱼的 ATPase活性下降速度较快     

青藏高原土壤水热分布特征及冻融过程在季节转换中的作用

利用GAME-Tibet期间所取得的高分辩率土壤温度和含水量资料，对青藏高原（主要是藏北高原）土壤水热分布特征及冻融过程在季节转换中的作用进行了分析。指出藏北高原4cm学深处土壤在10月份开始冻结，次年4-5月份开始消融，冻结持续时间长达5-7个月。冻结过程有利于土壤维持其水分，因此，在刚刚开始消融时土壤含水量仍然很高。从而为夏季风爆发前土壤通过蒸发向大气提供水分打下了基础。指出土壤冻融过程可能在高原季节转换中起着重要作用。     

地下非均匀非饱和带中地下洞室的渗流问题数值模拟--介质参数的灵敏度分析

影响任何一个地下废物处置场所长期行为的一个重要因素是进入地下废物贮置洞室内的渗流量的大小，预测地下洞室中的渗流量是困难的，特别是当地下废物贮置洞室位于非均质非饱和带中，三维数值模型用于研究地下非均质中非饱和流动及流入地下洞室中的渗流量，讨论了非均质模型与均质模型的比较以及数值剖分尺度对计算结果的影响，进入洞室中的渗漏率随着引入系统中入渗量的增加而增大，选用4个参数来衡量渗流场的非均质程度：(1)介质的平均渗透率K0；(2)Van Genuchten参数α；(3)渗透率空间分布相关尺度；(4)渗透率空间分布变化的标准差σ，根据一个随机实现的渗透率分布，通过改变平均渗透率来研究其对流入洞室中渗透量的影响，对一个固定的入渗率而言，流入洞室中的渗漏率将随着Van Genuchten参数α的减小而减小，模拟结果表明流入洞室的渗漏率与介质的平均渗透率相关，即随平均渗透率的增大，流入洞室的渗漏率亦增大，流入洞室的渗漏率还高度依赖于非均质渗透率场的空间分布相关长度与标准差，一个大的相关长度或高的标准差均能导致流入洞室的渗漏率增大。     

冻土路基表面的融化指数与冻结指数

在冻土层之上筑路，由于会改变地-气界面的热物理特性，进而影响冻土层的热力→动力稳定性，故而修筑一定高度的路基成为保护冻土层所采取的一种常规措施．在修筑路基之后，与路基边坡的朝向有关的热效应是冻土路基工程保护措施必须考虑的问题．在数理分析与数值模拟分析的基础上，给出了可根据气温的年最大和最小月平均值计算路基表面的融化指数与冻结指数以及有关热状况参数的方法，并以青藏铁路北麓河段2002年为例进行了计算分析．实例分析表明，即便是没有修筑道路，北麓河地区的冻土也已经处于临界状态；路基相对的两个坡面，由于朝向不同会造成温度分布的强非均匀性，其中南和偏南方向与北和偏北方向的路基坡面热状况差异最大，有必要对路基相对的两个坡面采用不同的防护措施，一方面改善就地取土修筑路基对其下伏冻土层的直接不良影响，同时也尽可能减小路基表面温度分布的非均匀性，以避免纵向裂缝的发生。     

南京地铁联络通道冻结法施工措施分析

对南京地铁一期工程TA4标联络通道冻结法施工的成功经验进行了总结。指出地铁联络通道冻结法施工中,必须采取必要的施工技术措施,才能保证施工安全、顺利的进行。本工程的成功经验可供其他工程参考。     

Generalized Coarse Graining Procedures for Flow in Porous Media

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.