高效优化非分裂PML边界二阶标量波方程数值模拟方法

卷积完全匹配层（convolution perfectly matched layer,CPML）吸收边界是一种高效处理波动方程数值模拟中人工边界反射波的方法。本文基于传统的一阶系统CPML吸收边界条件推广并推导了新的二阶系统CPML边界条件（NCPML）。与常规二阶系统CPML边界条件不同,新边界条件推导的核心思想是在复数-频率域中忽略部分衰减因子空间变化特性,避免其在时间域产生复杂卷积算子,然后反变换至时间域得到基于CPML吸收条件的二阶标量波方程,并应用于二阶标量波方程数值模拟。均匀介质模型测试验证了NCPML吸收条件在内存使用上相对于常规二阶系统CPML与SPML（split PML）吸收条件更少。在对人工边界反射的吸收效果上,NCPML稍逊色于常规二阶系统CPML,但二者均相对于SPML优势明显。最后通过层状模型和Marmousi模型测试验证了NCPML的稳定性及其在效率上的优势。     

基于两种插值算法的三维地质建模对比

三维地质建模过程中,插值算法对模型准确性具有显著影响。为评价不同插值算法对三维建模准确性的影响,本文选取反距离权重插值法和自然邻域插值法开展对比研究。通过理论分析和案例研究,从统计学原理、插值误差和可视化效果等三方面进行了对比分析。结果表明：反距离权重插值法在建模中精度较高,适应面更广。与自然邻域插值法相比,反距离权重插值法更加适用于地层缺失严重的层位,能够更好地保留地层缺失的特征;同时,反距离权重插值法能够更好地处理断层构造,对于地层的错断起伏情况表现效果更好;反距离权重插值法在沉积地层中误差更小,与实际情况更接近。     

城市规划对大气环境变化及空气质量的影响

采用北京大学三维复杂地形中尺度数值模式,利用Landsat-TM卫星影像图和"北京海淀区北部地区总体规划--用地规划图",分别模拟了北京海淀区北部规划完成后和目前情况下的气象场和大气环境,并分析了规划区的建成对本地和整个海淀区的影响.通过将观测的直径小于10 μm的粒子即PM10总量加权分配得到源强分布,计算海淀区的PM10浓度分布,在目前源强分布调查很困难的情况下不失为一种较为合理可行的方法.试验表明规划区的建成会增加当地和其下风方向的热岛强度和PM10浓度.对规划区进行较为合理的布局,如相对增加树木减少建筑物的面积将会缓和上述影响.     

Integration schemes with substepping algorithms for creep analysis in geomaterials

This paper describes a robust and efficient methodology for predicting displacements, deformations, and stresses in geomaterials that are susceptible to creep. The methodology is based on two integration schemes, which consider substepping algorithms. The first scheme is used for integrating space-time relations in a global sense, whereas the second scheme is used for integrating stress-stain relations in a local sense. Different from previous studies, both integration schemes are easy to implement and general in the sense that they can be applied to any type of creep law. Through an in-house finite element simulator, several numerical tests are performed. They include triaxial and wellbore closure analyses considering soft soils and salt rocks. The results show that the combination of both schemes leads to stable and accurate solutions with reduced computational time.     

Response to the paper by H. Zia,B. Lecampion “Explicit versus implicit front advancing schemes for the simulation of hydraulic fracture growth” (Int. J. Numer. Anal. Methods Geomech., 2019, 43 (6), 1300–1315)

The response presents three comments, which (i) explain theoretical background for the success of the explicit front advancing scheme by Zia and Lecampion (comment 1), (ii) suggest possible extensions (comment 2), and (iii) pay attention on expected limitations of the scheme to be checked by the authors in numerical experiments (comment 3). The response does not intend to present novel results: it is entirely based on known results published in the references cited.     

Towards Regional Synergy: Reconciling Rangeland Ecological Functioning with Forage Production of Cultivated Pasture

Animal husbandry and crop farming are specialized for development in separate areas on the Tibetan Plateau. Such a pattern of isolation has led to current concerns of rangeland and farming system degradation due to intensive land use. The crop-livestock integration, however, has been proven to increase food and feed productivity thorough niche complementarity, and is thereby especially effective for promoting ecosystem resilience. Regional synergy has emerged as an integrated approach to reconcile rangeland livestock with forage crop production. It moves beyond the specialized sectors of animal husbandry and intensive agriculture to coordinate them through regional coupling. Therefore, crop-livestock integration (CLI) has been suggested as one of the effective solutions to forage deficit and livestock production in grazing systems. But it is imperative that CLI moves forward from the farm level to the regional scale, in order to secure regional synergism during agro-pastoral development. The national key R & D program, Technology and Demonstration of Recovery and Restoration of Degraded Alpine Ecosystems on the Tibetan Plateau, aims to solve the problems of alpine grassland degradation by building up a grass-based animal husbandry technology system that includes synergizing forage production and ecological functioning, reconciling the relationship between ecology, forage production and animal husbandry, and achieving the win-win goals of curbing grassland degradation and changing the development mode of animal husbandry. It is imperative to call for regional synergy through integrating ecological functioning with ecosystem services, given the alarming threat of rangeland degradation on the Tibetan Plateau. The series of papers in this issue, together with those published previously, provide a collection of rangeland ecology and management studies in an effort to ensure the sustainable use and management of the alpine ecosystems.     

Estimating groundwater recharge for a freshwater lens in an arid region: Formative and stability assessment

The formation of subsurface freshwater lenses on top of brackish groundwater is a fascinating hydrologic phenomenon that creates groundwater supplies of great potential value in arid regions. Information on the recharge quantity and mechanism of these lenses is both scarce and uncertain. This study examines the formation and macroscale stability of the Rawdatain freshwater lens in Kuwait, for which significant pre-development data are available. The Rawdatain is a large (150 million m³) subsurface freshwater lens overlying brackish groundwater compared to the other freshwater lenses in the Arabian Peninsula. In this study, a three-dimensional (3D) density-dependent groundwater flow model is tested against the following data targets to estimate long-term diffuse and focused groundwater recharge: (i) groundwater head, (ii) total dissolved solids (TDS) groundwater concentration, (iii) volume and vertical thickness of stored groundwater of three different water quality TDS ranges (0–700, 700–1000 and 1000–2000 mg/L) and (iv) geometrical shape features of the lens along cross-sections. To better represent the spatial variation in TDS, six different recharge zones were assigned to allocate diffuse and focused recharge conditions. Twelve recharge rate scenarios, encompassing a wide range of feasible long-term average annual recharge values (200,000–5,000,000 m³/year), were tested against the multiple targets and compared with the groundwater age of the Rawdatain lens. Based on comparison with data targets, the long-term average annual recharge is estimated to be 500,000 m³/year. Scenarios of reduced recharge, which may occur due to changes in land-use or climate, demonstrate the extremely slow response of the lens, which is in agreement with the slow development and formation of the lens (>2,000 years). Within a 100-year time frame, a 50% reduction in annual recharge reduces the lens volumes by 21, 17 and 9% for the three water quality categories, respectively. This study demonstrates the stability of freshwater lenses in arid regions and also provides methodology for similar focused rainfall recharge freshwater lenses.     

Regional extreme‐dry‐spell frequency analysis using the L‐moments method in the middle reaches of the Yellow River Basin,China

In this research, the regional extreme‐dry‐spell frequency in the middle reaches of the Yellow River Basin (YRB) is studied by the L‐moments method. The research area has been divided into three subregions (regions 1, 2 and 3), which have been identified as homogenous regions. The results of a goodness‐of‐fit test indicate that a generalized normal distribution is the optimal regional model for regions 1 and 2 whereas a generalized Pareto distribution is the optimal regional model for region 3. The return period analysis figures out that the maximum length‐of‐dry‐spell (MxDS) values increase from south to north in the southern part and increase from northeast to southwest in the northern part of the middle reaches of the YRB under different return periods. The increments of quantiles of dry spell under different return levels indicate that drought risk in region 1 is higher than that in regions 2 and 3. The analysis of the occurrence day of MxDS shows that MxDS mostly occurred during winter of 1998 and spring of 1999 in most stations during the considered period. By comparing summer MxDS events, it can be found that mean MxDS values have slightly increased in regions 1 and 2 during the last five decades. The maximum mean MxDS values appeared in the 2000s for regions 1 and 2 and in the 1990s for region 3. The atmospheric circulation shows that the positive anomaly centre in the west of North China, negative anomaly centre in the east of North China and the strong western Pacific subtropical high led to the decrease of precipitation in North China during the summer of 1997. Copyright © 2013 John Wiley & Sons, Ltd.     

Analysis of uncertainties in the hydrological response of a model‐based climate change impact assessment in a subcatchment of the Spree River,Germany

Climate change impact assessments form the basis for the development of suitable climate change adaptation strategies. For this purpose, ensembles consisting of stepwise coupled models are generally used [emission scenario → global circulation model → downscaling approach (DA) → bias correction → impact model (hydrological model)], in which every item is affected by considerable uncertainty. The aim of the current study is (1) to analyse the uncertainty related to the choice of the DA as well as the hydrological model and its parameterization and (2) to evaluate the vulnerability of the studied catchment, a subcatchment of the highly anthropogenically impacted Spree River catchment, to hydrological change. Four different DAs are used to drive four different model configurations of two conceptually different hydrological models (Water Balance Simulation Model developed at ETH Zürich and HBV‐light). In total, 452 simulations are carried out. The results show that all simulations compute an increase in air temperature and potential evapotranspiration. For precipitation, runoff and actual evapotranspiration, opposing trends are computed depending on the DA used to drive the hydrological models. Overall, the largest source of uncertainty can be attributed to the choice of the DA, especially regarding whether it is statistical or dynamical. The choice of the hydrological model and its parameterization is of less importance when long‐term mean annual changes are compared. The large bandwidth at the end of the modelling chain may exacerbate the formulation of suitable climate change adaption strategies on the regional scale. Copyright © 2013 John Wiley & Sons, Ltd.     

城市创意产业空间动态集聚演化的计算与可视优化方法

基于人文地理视角下的城市创意产业图像可视化分析对城市深层次空间综合和区域创新发展具有重大意义。但Swarm群智能动态时空建模难以满足创意产业空间集聚的可视化发展。本文研究目标是,从城市区域创意产业空间聚类影响因素指标出发,创新性地提出区域空间动态集聚轨迹算法（Density-Based Interest Spatial Clustering of Path,DBICP）,并与计算机浏览器共建聚类可视化图像,为城市管理提供决策依据。首先,根据影响因素指标体系,利用2014—2018年空间卡口流量数据和产业指标数据进行预处理,构建空间标准聚类算法DBSCAN。然后,对其进行聚类密度分级优化形成全新DBICP算法并得出初步轨迹图像。最后,通过源码转译实现了浏览器界面下空间动态集聚轨迹图像的输出。结果表明：以上海市为例,普陀区、浦东新区、徐汇地区的创意产业空间分布形成了3种不同的聚类模式,并相应提出了分摊、均布、虹吸的管控策略。此方法克服了传统图像的聚类分级和轨迹测量的缺失,可以有效地从指标数据中发现图像轨迹聚类信息,体现了地理信息科学和人文社会学科的交叉融合。也为大数据动态图像的集聚方法提供了全新视角和借鉴价值。