Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

长江三峡大宁河流域3000年来的环境演变与人类活动

长江三峡大宁河流域张家湾东周-汉代古遗址的孢粉学和沉积学等分析结果表明3000年来该区环境演变与人类活动的互动影响，连续的孢粉和粒度记录证实，该区古人类活动初期处于暖湿的气候，东周-汉代时期转变为温凉略干的温带气候，东周-汉代后，该区经历了一次规模较大的大宁河洪水泛滥期，导致汉代文化层中断。此后不久，由于人类活动剧烈，自然植被严重破坏，水土流失比较严重，山洪灾害频繁，导致汉代之后该区没有连续的文化层堆积，而由于人类活动剧烈，自然植被严重破坏，水土流失比较严重，山洪灾害频繁，导致汉代之后该区没有连续的文化层堆积，而由于短暂洪水暴发引起大面积坡面片流，导致坡积物的大量堆积，汉代之后该区频繁的兴涝灾害和山洪灾害体现了人地关系的复杂性。     

福州市植被信息提取及其景观格局空间变化

ASTER(高级空间热辐射热反射探测仪)是高光谱遥感影像时代到来的标志,与一般遥感影像相比能为对地观测提供更高质量的信息源。根据福福州市1988年的TM影像和2001年ASTER影像,在NDVI和分形几何理论方法的基础上,求取两个时期植被斑块的分形维数,结果表明:植被景观整体分形维数上升,说明植被斑块空间形态趋向复杂化,稳定性降低,其空间格局变化表现出2种基本规律:一是破碎化,孔隙化,主要发生在城区内部,二是空间萎缩,主要表现在城市边缘地区。     

科尔沁沙地固定沙丘土壤水体积分数变化

利用2006-2009年5-9月科尔沁沙地固定沙丘0~150 cm土壤水体积分数监测数据及降水资料,采用偏度、峰度检验法检验了固定沙丘不同月及不同土层土壤水体积分数的正态性,并得到了其正态分布相关参数的估计值及其置信区间;同时,对不同土层、不同月及不同年间的土壤水体积分数变化和差异进行了分析。结果表明:(1)固定沙丘同一月的土壤水体积分数以显著性水平α=0.1通过正态分布的假设检验;固定沙丘土壤水体积分数的平均值范围1.82%~4.95%,主要集中在2%~4%;(2)受降水影响,固定沙丘7月土壤水体积分数较高,且与其他月的差异显著;7月土壤水体积分数均值和方差的点估计及置信水平为0.95的置信区间、统计特征均与其他月的有较大差异;(3)固定沙丘0~10 cm土壤水体积分数较低,且与其他土层有显著差异;而10~30 cm土壤水体积分数较高;(4)2008年土壤水体积分数较低,而2009年土壤水体积分数较高;除2007年与2008年土壤水体积分数没有显著差异外,其余年份间的土壤水体积分数均有显著差异。     

2000-2012年陕甘宁黄土高原区地表蒸散时空分布及影响因素

分析陕甘宁黄土高原区地表蒸散变化特征及其影响因素,可以为区域水资源规划、生态环境改善提供依据。本文利用MOD16蒸散数据,统计分析了陕甘宁黄土高原区2000-2012年地表实际蒸散量的时空变化特征,并结合国家气象站点观测数据和基于像元的相关分析法探讨了其影响因素。结果表明:(1) 2000-2002年蒸散量迅速上升,在2003年达到最高值378.6 mm, 2003-2006年蒸散量呈下降趋势,2006年之后蒸散量呈现缓慢上升趋势。(2) 近13年来,陕甘宁黄土高原区多年平均蒸散量具有明显的空间差异,蒸散量自西北至东南递增,最南部的六盘山、子午岭、黄龙山地是3个主要的高值区;年蒸散量以夏季最多,其次是春季,秋季和冬季最少,且季节蒸散的分布与年蒸散的空间分布格局基本一致。(3) 陕甘宁黄土高原区蒸散量草地和耕地的贡献率最高,密灌丛、疏灌丛次之,常绿针叶林、森林草原贡献率则较小。(4) 陕甘宁黄土高原区动力因素对地表蒸散量影响以正相关为主,风速对该区影响较大;热力因素对地表蒸散量影响以负相关为主,其中气温与蒸散在空间上呈负相关的区域较大,日照时数与蒸散在空间上的负相关区域的面积次之;水分条件(降水量、相对湿度)对蒸散的影响也以正相关为主。     

中国大气本底观测发展概述

中国早期的大气本底观测网建设起步于上世纪80年代初，并于10a后基本建成。目前，中国气象局、中国科学院在原有观测站的基础上各自建设和发展着自己下属的大气本底观测网络。全球大气观测网给出了一整套开展大气本底观测的技术标准，各国均予以严格遵循。20多年来，中国大气本底观测工作在基础数据积累以及观测标准制定等方面取得了较大进展，与此相关的科学研究也取得了长足进步。最后对中国大气本底观测发展过程中存在的一些问题进行了探讨。     

中国城市价值链功能分工及其影响因素

价值链的空间重组正在深刻的改变着城市体系的经济景观,建立在价值链分工基础上的城市功能结构的研究已经成为经济地理学的重要课题。将中国上市公司500强企业网络划分为公司总部、商务服务、研究开发、传统制造、现代制造、物流仓储和批发零售七种类型功能区块,研究了中国城市价值链功能分工及其影响因素。结果发现：沿着价值链的功能分工已经成为中国城市体系经济景观的显著特征,功能多样化城市和功能专业化城市并存于中国城市体系,东部地区和城市密集地区的城市在价值链分工中占据了更好的地位;中国城市按照价值链中的优势功能可以划分为九种类型,少数城市转变为承载公司总部、研究开发、商务服务等多样化功能的高等级中心城市,而大量中小城市则转变为传统制造专业化基地;市场潜力、关键资源、区位条件、营商环境等城市属性特征是城市功能分工的重要影响因素,城市资源、区位可达性等属性特征的增强将提高城市成为总部基地、商务中心和研究基地的概率,而降低城市成为传统制造基地的概率。     

Variation characteristics of evaporation in the Gulang River Basin during 1959-2013

Based on monthly evaporation of two meteorological stations in the Gulang River Basin of China, the inter-annual variation of evaporation during 1959-2013 were analyzed using Mann-Kendall and wavelet analysis. The results demonstrated that the annual evaporation show a fluctuating increase over the past 50 years approximately, with an average increase rate of 4.26 mm per decade. The overall trend was decrease-increase-decrease. According to the cumulative anomaly curve,the turning point of the annual evaporation occurred in 1979, in which the evaporation increased in the early stage and decreased in the later stage. Meanwhile, the seasonal variation of the evaporation shows that it decreased in Spring and Autumn, and increased in Summer and Winter, especially obvious for the later. The evaporation abruptly changed in Spring and Summer in 2008 and in Winter in 1994. In addition, all evaporation increased after the changes. However,the evaporation in Autumn abruptly changed in 1986 and 1999, which show a trend of increase-decrease-increase.Wavelet analysis shows that evaporation in Summer and wet season would decrease in the next few years, and in the other seasons would increase. Based on the aforementioned analysis, it can be concluded that increased evaporation is mainly induced by increase of evaporation in dry season, especially in Winter, and this trend to be continued in the future for the Gulang River Basin.     

Population,urbanization and economic scenarios over the Belt and Road region under the Shared Socioeconomic Pathways

The countries throughout the Belt and Road region account for more than 60% of the world's population and half of the global economy. Future changes in this area will have significant influences on the global economic growth, industrial structure and resource allocation. In this study, the proportion of the urban population to the total population and the gross domestic product were used to represent the levels of urbanization and economic development, respectively. The population, urbanization and economic levels of the Belt and Road countries for 2020–2050 were projected under the framework of the IPCC's shared socioeconomic pathways(SSPs), and the following conclusions are drawn.(1) The population, urbanization and economic levels in the Belt and Road region will likely increase under all five pathways. The population will increase by 2%–8%/10 a during 2020–2050 and reach 5.0–6.0 billion in 2050. Meanwhile, the urbanization rate will increase by 1.4%–7.5%/10 a and reach 49%–75%. The GDP will increase by 17%–34%/10 a and reach 134–243 trillion USD.(2) Large differences will appear under different scenarios. The SSP1 and SSP5 pathways demonstrate relatively high urbanization and economic levels, but the population size is comparatively smaller; SSP3 shows the opposite trend. Meanwhile, the economy develops slowly under SSP4, but it has a relatively high urbanization level, while SSP2 exhibits an intermediate trend.(3) In 2050, the population will increase relative to 2016 in most countries, and population size in the fastest growing country in Central Asia and the Middle East countries will be more than double. Urbanization will develop rapidly in South Asia, West Asia and Central Asia, and will increase by more than 150% in the fastest growing countries. The economy will grow fastest in South Asia, Southeast Asia and West Asia, and increase by more than 10 times in some counties with rapid economic development.     

Effects of land use and cultivation time on soil organic and inorganic carbon storage in deep soils

The vertical distribution and exchange mechanisms of soil organic and inorganic carbon(SOC, SIC) play an important role in assessing carbon(C) cycling and budgets. However, the impact of land use through time for deep soil C(below 100 cm) is not well known. To investigate deep C storage under different land uses and evaluate how it changes with time, we collected soil samples to a depth of 500 cm in a soil profile in the Gutun watershed on the Chinese Loess Plateau(CLP); and determined SOC, SIC, and bulk density. The magnitude of SOC stocks in the 0–500 cm depth range fell into the following ranking: shrubland(17.2 kg m~(-2)) grassland(16.3 kg m~(-2)) forestland(15.2 kg m~(-2)) cropland(14.1 kg m~(-2)) gully land(6.4 kg m~(-2)). The ranking for SIC stocks were: grassland(104.1 kg m~(-2)) forestland(96.2 kg m~(-2)) shrubland(90.6 kg m~(-2)) cropland(82.4 kg m~(-2)) gully land(50.3 kg m~(-2)). Respective SOC and SIC stocks were at least 1.6-and 2.1-fold higher within the 100–500 cm depth range, as compared to the 0–100 cm depth range. Overall SOC and SIC stocks decreased significantly from the 5 th to the 15 th year of cultivation in croplands, and generally increased up to the 70 th year. Both SOC and SIC stocks showed a turning point at 15 years cultivation, which should be considered when evaluating soil C sequestration. Estimates of C stocks greatly depends on soil sampling depth, and understanding the influences of land use and time will improve soil productivity and conservation in regions with deep soils.