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.     

Silicon limitation on primary production and its destiny in Jiaozhou Bay, China——Ⅳ：Study on cross-bay transect from estuary to ocean

The authors analyzed the data collected in the Ecological Station Jiaozhou Bay from May 1991 to November 1994, including 12 seasonal investigations, to determine the characteristics, dynamic cycles and variation trends of the silicate in the bay. The results indicated that the rivers around Jiaozhou Bay provided abundant supply of silicate to the bay. The silicate concentration there depended on river flow variation. The horizontal variation of silicate concentration on the transect showed that the silicate concentration decreased with distance from shorelines. The vertical variation of it showed that silicate sank and deposited on the sea bottom by phytoplankton uptake and death, and zooplankton excretion. In this way, silicon would endlessly be transferred from terrestrial sources to the sea bottom. The silicon took up by phytoplankton and by other biogeochemical processes led to insufficient silicon supply for phytoplankton growth. In this paper, a 2D dynamic model of river flow versus silicate concentration was established by which silicate concentrations of 0.028–0.062 μmol/L in seawater was yielded by inputting certain seasonal unit river flows (m³/s), or in other words, the silicate supply rate; and when the unit river flow was set to zero, meaning no river input, the silicate concentrations were between 0.05–0.69 μmol/L in the bay. In terms of the silicate supply rate, Jiaozhou Bay was divided into three parts. The division shows a given river flow could generate several different silicon levels in corresponding regions, so as to the silicon-limitation levels to the phytoplankton in these regions. Another dynamic model of river flow versus primary production was set up by which the phytoplankton primary production of 5.21–15.55 (mgC/m²·d)/(m³/s) were obtained in our case at unit river flow values via silicate concentration or primary production conversion rate. Similarly, the values of primary production of 121.98–195.33 (mgC/m²·d) were achieved at zero unit river flow condition. A primary production conversion rate reflects the sensitivity to silicon depletion so as to different phytoplankton primary production and silicon requirements by different phytoplankton assemblages in different marine areas. In addition, the authors differentiated two equations (Eqs. 1 and 2) in the models to obtain the river flow variation that determines the silicate concentration variation, and in turn, the variation of primary production. These results proved further that nutrient silicon is a limiting factor for phytoplankton growth. This study was funded by NSFC (No. 40036010), and the Director's Fund of the Beihai Sea Monitoring Center, the State Oceanic Administration.     

复杂模型重力位场快速正反演方法及应用

提要 本文详细讨论了一种三维重力位场快速正反演方法。作者在前人工作的基础上,对算法作了行之有效的改进,通过对反演中的不稳定因素进行各种理论模型试算,得出保证迭代反演稳定收敛的准则,编制出可在微型机IBM—PC上运行的人机对话式自动正反演程序。本文还对各种不均质模型进行了模似计算并将该方法应用于某含油气沉积盆地的双层界面构造研究,揭示出了储油有利地段。     

Design of stable drainage networks in upland forestry plantations

Erosion of man-made, forestry drainage channels occurs when the plough cuts through the surface peat layer into the underlying erodible mineral soil. A procedure is developed, based on hydraulic considerations, which will allow the drainage engineer to design stable drainage networks in upland forestry plantations. An example design chart is given for an erodible sandy loam type soil.     

区域GPS网实时计算可降水量的若干问题

目前地基GPS气象学测得的可降水量 (PWV )精度好于 2mm ,但在利用区域GPS网实时计算每个测站上空的PWV时 ,要涉及到很多常规GPS资料处理时所忽略的问题 ,如需考虑数据处理软件和计算方式的选择、站坐标的确定和约束、轨道的使用方法、网外辅助站最佳数量的确定、海潮对实时计算PWV的影响以及实时应用于气象服务时的端部效应等问题。利用上海GPS综合应用网获取的 2 0 0 2年 6、7月份长江三角洲地区入梅前后的数据 ,分析了利用区域性的GPS网实时计算高精度的PWV时要解决的各种问题 ,探讨了其数据处理方案     

Dissipative-acoustic instability in accretion disks at a nonlinear stage

Nonstationary hydrodynamic models of a viscous accretion disk around a central compact object were constructed. Two different numerical methods (TVD and SPH) are used to study the dynamics of dissipatively unstable acoustic perturbations at the nonlinear stage in terms of the standard α-disk model. The standard disk accretion in the Shakura-Sunyaev model is unstable against acoustic waves for various parameters of the system. If the α parameter, which specifies the level of turbulent viscosity, exceeds α?0.03, then a complex nonstationary system of small-scale weak shock waves is formed. The growth rate of the perturbations is higher in the central disk region. For α?0.2, the relative shock amplitude can exceed 50% of the equilibrium disk parameters. The reflection of waves from the disk boundaries and their nonlinear interaction are important factors that can produce unsteady accretion. The luminosity of such a disk undergoes quasi-periodic oscillations at a level of several percent (?5%) of the equilibrium level.     

Deep Processes of Basin Evolution:A Case Study of the Southern Songliao Basin and Adjacent Area,NE China in Mesozoic-Cenozoic

1IntroductionOne of the key factors related to basin geody-namics is deep process controlling formation and evolution of sedimentary basin. Depth and tempera-ture of asthenosphere,existence of mantle plume,occurrence of mantle melting,and amounts of melts under depressurization during thinning of lithosphere are controlling factors influenced formation and characteristics of extensional basin directly (Li,1994). Rifting is probably governed by frictional forces exerted on the base of litho…     

用层次分析法决策地下水库回灌方案的排序问题

本文引用了一种新的系统分析方法(层次分析法)对地下水库人工回灌的复杂系统进行了择优排序研究,取得了满意的结果。实践证明,层次分析法对于地下水库人工回灌方案的择优是一种简捷实用的决策方法。     

REVIEWS OF GEOGRAPHIC SOFTWARE

Reviews of geographic software in this article: DEMO-GRAPHICS: WORLD POPULATIONS AND PROJECTIONS. ESP GAUSS. CEMODEL S. Damus LIMDEP. William H. Greene MICROSTAT 4.1 OTIS PCIPS. (Personal Computer Image Processing System) . H.J. Meyers and R. Bernstein. REGRESSION ANALYSIS OF TIME SERIES (RATS) SPSS/PC+ URBAN DATA MANAGEMENT SOFTWARE (UDMS)     

Application of stress release models to historical earthquakes from North China

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.