Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

开发旅游观光资源与四川西部之发展

川西甘孜州旅游观光资源丰富多彩，遍布的温泉、热泉与雪山、草地、森林相配置成的景点、景区比比皆是．加之悠久的人文、宗教建筑和文物，都具有很高的研究和观赏价值。开发这些资源以引导甘孜州经济的发展，将是大有可为的。     

太湖表层沉积物多环芳烃的空间分布、来源及其风险

随着太湖流域社会与经济的发展,多环芳烃(PAHs)在各种环境介质中逐渐累积,污染日益严重,可能对太湖生态环境及周边人体健康构成威胁。为探究太湖沉积物PAHs的来源及生态风险,于2021年12月在太湖采集30个表层沉积物样品,利用气相色谱-质谱联用仪(GC-MS)检测样品中16种PAHs含量;利用受体模型和苯并[a]芘(BaP)毒性当量法进行来源解析及生态风险评估,并将各来源贡献与毒性当量浓度相结合,量化源风险。结果表明,太湖表层沉积物中16种PAHs总含量介于124~592 ng/g之间,平均值为294 ng/g,中值为279 ng/g;高环多环芳烃(HMW PAHs)为主要组分,占∑PAHs的67%。高含量区域位于竺山湾、梅梁湾、贡湖湾和西太湖,与国内外其他湖泊沉积物相比,太湖沉积物PAHs含量处于较低水平。源解析的结果表明,太湖表层沉积物中PAHs交通排放源贡献率为29.1%、煤炭燃烧源贡献率为26.7%、生物质燃烧源贡献率为28.7%、石油源贡献率为15.6%。生态风险评价结果表明,交通排放源、生物质燃烧源、煤炭燃烧源和石油源的BaP毒性当量含量(TEQ_BaP)均值分别为19.34、17.81、16.33和9.1 ng/g,均小于70 ng/g,几乎处于无风险水平。西太湖、贡湖湾和梅梁湾的部分区域ΣTEQ_BaP大于70 ng/g属于潜在风险区,具有一定潜在毒性。在后续的污染治理中应重点关注太湖西北部地区污染物的排放。本研究可为沉积物中PAHs污染的研究提供数据支撑,为地方政府精准、高效地管控PAHs污染提供理论依据。     

Correspondence Between Astronomical Periods and Sedimentary Cycles

The term of geological cycle was used firstly inHutton's" The Earth Theory" in 1788 (Williams,1980). The geological cycles including tectonic cycles and sedimentary cycles have been describedwidely since than. Some researchers have perceivedthat the geolo…     

Risk Management for Petroleum Reservoir Production: A Simulation-Based Study of Prediction

We consider numerical solutions of the Darcy and Buckley–Leverett equations for flow in porous media. These solutions depend on a realization of a random field that describes the reservoir permeability. The main content of this paper is to formulate and analyze a probability model for the numerical coarse grid solution error. We explore the extent to which the coarse grid oil production rate is sufficient to predict future oil production rates. We find that very early oil production data is sufficient to reduce the prediction error in oil production by about 30%, relative to the prior probability prediction.     

Massive sedimentation of fine sediment with organic matter and enhanced benthic-pelagic coupling by an artificial dyke in semi-enclosed Chonsu Bay, Korea

To assess the impact of an artificial dyke in Chonsu Bay (CBD) on the organic carbon (C(org)) cycle, we measured excess (210)Pb activities, C(org) and nitrogen content in sediment cores. The C(org) oxidation rates (C(ox)) on the surface sediment and benthic nutrient fluxes were also quantified with an in situ benthic chamber. The higher excess (210)Pb inventory, C(org) and nitrogen in cores near the CBD indicated lateral transport and local, massive deposition of particulate matter due to tidal circulation altered by artificial dyke construction. The C(ox) in sediment near the CBD was about twice as high as that out of the bay, suggesting the importance of benthic remineralization of organic matter. The benthic fluxes of dissolved inorganic nitrogen and phosphate were four to six times higher than those outside the bay, corresponding to 141% and 131% respectively, of the requirements for primary production.     

鲜水河断裂带炉霍县卡扎村大探槽古地震研究

本文测定分析了特殊地区的两种不同饱和状态的干压实土体在不同等级荷栽作用下产生的沉陷量，探讨了施荷时间对沉陷量的影响从而提示了两种干压实土体在荷载作用下沉陷量变化规律，并分析研究了这两种干压实土体的最大压缩模量与饱和度的关系。这一研究成果对干旱特殊地区的公路设计施工具有一定的参考价值。     

P-T Paths Derived from Garnet Growth Zoning in Danba Domal Metamorphic Terrain, Sichuan Province, West China

Danba (丹巴) domal metamorphic terrain belongs to Songpan (松潘)-Ganze (甘孜) orogenic belt, where typical Barrovian and Buchan metamorphic zones are preserved. The former included chlorite, biotite, garnet, staurolite, kyanite and sillimanite zones, while the latter only developed silimanite+muscovite and sillimanite+K-feldspar zones. Integrated study has been carried on metamorphic reactions of garnet production and consumption, P-T paths and P-T-X-M phase relation and thermal tectonic model for Danba metamorphic zones. Petrological textures in thin sections show that garnet production and consumption in kyanite-sillimanite zone is mainly attributed to ChI+Ms+PI+Q=Grt+Bt+H2O and kyanite=sillimanite respectively. Based on mineral compositions, the geothermobarometry gives an average P, T condition of (4.9±0.3)×108 Pa, 543±30℃ for the first growth stage of the garnet and (5.8±0.3)±108 Pa, 534±29 ℃ for the second stage of garnet growth respectively. Anti-counter clockwise P-T paths were drawn using Gibbs method by NCMnKFMASH system for sample G98686 in the kyanite zone. The P-T-X-M modeling for the first mineral assemblages shows that the prediction is similar to the measured values in gossular, almandine and spessartine but mole fraction of pyrope and Fe/(Fe+Mg) deviated far from the contours; while that for the second mineral assemblages exhibits that the prediction is consistent with the measured value of pyrope, grossular content and Fe/(Fe+Mg) of garnet. A thermal tectonic model that there are at least three structure levels across the thrnst-decollement zones is presented according to the P-T paths, metamorphic grades and deformation styles for the staurolite-kyanite zone of the Barroviau type metamorphism, which will provide some constraints for the evolution of the nappe complex.