Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Exposure dating and reinterpretation of coarse debris accumulations (‘rock glaciers’) in the Cairngorm Mountains,Scotland

Relict rock glaciers have considerable potential for contributing to palaeoclimatic reconstruction, but this potential is often undermined by lack of dating control and problems of interpretation. Here we reinvestigate and date four proposed ‘rock glaciers’ in the Cairngorm Mountains and show that the morphology of only one of these appears consistent with that of a true rock glacier produced by creep of underlying ice or ice‐rich sediment. All four features comprise rockslide or rock avalanche runout debris, and the possibility that all four represent unmodified runout accumulations cannot be discounted. Surface exposure dating of the four debris accumulations using cosmogenic ¹⁰Be produced uncertainty‐weighted mean ages of 15.4 ± 0.8 ka, 16.2 ± 1.0 ka, 12.1 ± 0.6 ka and 12.7 ± 0.8 ka. All four ages imply emplacement under cold stadial conditions, two prior to the Windermere Interstade of ca. 14.5–12.9 cal. ka BP and two during the Loch Lomond Stade of ca. 12.9–11.5 cal. ka BP. The above ages indicate that paraglacial rock‐slope failure on granite rockwalls occurred within a few millennia after deglaciation. The mean exposure ages obtained for runout debris at two sites – Strath Nethy (16.2 ± 1.0 ka) and Lairig Ghru (15.4 ± 0.8 ka) – are consistent with basal radiocarbon ages from Loch Etteridge, 22 km to the southwest (mean = 15.6 ± 0.3 cal. ka BP) and imply widespread deglaciation of the Cairngorms and adjacent valleys before 15 ka and possibly 16 ka. Copyright © 2008 John Wiley & Sons, Ltd.     

Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

中国兴凯湖北岸平原晚全新世花粉记录及泥炭沼泽形成

通过兴凯湖北岸平原泥炭剖面高分辩率花粉分析研究，对晚全新世花粉划分4个组合带。XKH-4组合时期(1857～1746aB．P．)为沼泽发育前期，这一时期花粉浓度小，陆生草本植物占优势，气候干冷。XKH-3组合时期(1746～1287aB．P．)为沼泽发育早期，这一时期花粉浓度较大，且水生植物花粉含量为剖面最高，喜温落叶阔叶植物大发展时期，气候温和湿润。XKH-2组合时期(1287～602aB．P．)为沼泽发育中期，这一时期花粉浓度最小，但以陆生草本植物为主，木本植物为辅，水生植物急剧减少，针叶植物出现两次高峰，气候向冷干方向发展。XKH-1组合时期(602aB．P．至今)为沼泽发育盛期，这一时期花粉浓度最大，陆生草本植物大发展时期，气候波动较大。     

The Geological Characteristics of Chengkou Type Shale Gas of Lower Cambrian in Dabashan

<正>1 Introduction After the successful implementation of large-scale shale gas exploration、development and utilization in United States,the domestic shale gas exploration and development is in a rapid development.However,the progress of exploration and understanding of the highEvolution shale gas in the complex structure zone of the outer basin are relatively slow,the evaluation method and     

基于卷积神经网络的ICESat-2光子点云去噪分类

ICESat-2（Ice, Cloud, and land Elevation Satellite-2）是美国NASA（National Aeronautics and Space Administration）在2018年发射的激光测高卫星,其上搭载的激光测高系统ATLAS（Advanced Topographic Laser Altimeter System）采用微脉冲多波束光子计数激光雷达系统,因其低能耗、高探测灵敏度、高重复频率的特性极大改善了沿轨采样密度,但也使获取的数据中包含大量的噪声,如何有效实现光子点云去噪分类成为后续应用的关键,也是当前研究的热点和难点,为此本文提出一种基于卷积神经网络的光子点云去噪和分类算法。首先将光子点云按照沿轨和高程方向划分格网,去除明显的噪声光子,并将每个粗信号光子点栅格化为影像;然后基于少量样本构建的卷积神经网络分类模型实现光子点云精去噪和分类;最后利用机载激光雷达数据进行验证,并与ATL08产品的去噪分类结果进行对比。结果表明,对于裸地和森林区域,卷积神经网络算法均能有效去除噪声光子,特别对于森林区域,可同时实现去噪和分类;其中,裸地区域地表计算的R²和RMSE分别为1.0和0.72 m,森林区域地表和树冠计算的R²分别为1.0和0.70, RMSE分别为1.11 m和4.99 m。本文利用深度学习算法实现光子点云去噪分类,在裸地和森林区域均取得了较好的结果,为后续光子点云数据处理提供了参考。     

Long-term statistics of extreme tsunami height at Crescent City

Historically, Crescent City is one of the most vulnerable communities impacted by tsunamis along the west coast of the United States, largely attributed to its offshore geography. Trans-ocean tsunamis usually produce large wave runup at Crescent Harbor resulting in catastrophic damages, property loss and human death. How to determine the return values of tsunami height using relatively short-term observation data is of great significance to assess the tsunami hazards and improve engineering design along the coast of Crescent City. In the present study, the extreme tsunami heights observed along the coast of Crescent City from 1938 to 2015 are fitted using six different probabilistic distributions, namely, the Gumbel distribution, the Weibull distribution, the maximum entropy distribution, the lognormal distribution, the generalized extreme value distribution and the generalized Pareto distribution. The maximum likelihood method is applied to estimate the parameters of all above distributions. Both Kolmogorov-Smirnov test and root mean square error method are utilized for goodness-of-fit test and the better fitting distribution is selected. Assuming that the occurrence frequency of tsunami in each year follows the Poisson distribution, the Poisson compound extreme value distribution can be used to fit the annual maximum tsunami amplitude, and then the point and interval estimations of return tsunami heights are calculated for structural design. The results show that the Poisson compound extreme value distribution fits tsunami heights very well and is suitable to determine the return tsunami heights for coastal disaster prevention.     

麻城1932年6．0级地震区地球物理场与深部环境研究

研究了麻城地区重、磁异常场特征，经反演计算获得了深部结构特征。对重力布格异常利用Ｐａｒｋｅｒ－Ｏｌｄｅｎｂｕｒｇ快速位场反演方法计算出莫霍面的展布，对航磁异常利用对数功率谱法和矩谱法计算出居里面深度。在此基础上，探讨该震区深部蕴震条件。主要结论是：①麻城－团凤断裂带为该区切割较深的主要控震构造；②１９３２年麻城６．０级地震发生在区域重力异常等值线梯度最大处，磁力异常区两局部正异常间的低值带，莫霍面及居里面隆起区边缘；③发震部位有其特定的深部蕴震构造及应力蕴震条件。     

基于Web服务的省级林业电子政务信息共享平台的研究

随着信息化建设的深入发展,计算机及网络应用范围不断拓展,应用系统的规模不断扩大,所面临的一个共同问题是如何在原有的建设基础上,进行信息资源的整合,实现信息资源共享。本文以福建省数字林业政务信息共享平台构建为例,在尽可能利用和保护原有投资存量的前提下,针对分布式异构环境中的林业数据和业务应用系统的信息共享要求,提出了基于.NET平台,运用先进的XML&Web Services/MS BizTalk Server2004核心技术的信息交换和共享的整体解决方案及其实现方法。