起伏地表采集数据的三维直接叠前时间偏移方法

提出一种可对起伏地表采集的三维地震资料直接进行偏移成像的叠前时间偏移方法和流程.它用两个等效速度描述近地表和上覆层对地震波传播的影响,可对炮、检点不在同一水平面的三维地震资料直接进行叠前时间偏移处理.该方法不对近地表地震波传播做垂直出、入射假定,因此可适应高速层出露等不存在明显低、降速带情况.描述近地表和上覆层的两个等效速度参数可依据偏移道集的同相轴是否平直来确定,避免了确定近地表速度的困难;而对已知近地表速度的情况,则可进一步修正近地表速度,获得更好的成像效果.用三维起伏地表的理论数据和中国东部某工区实际数据验证了所发展方法和处理流程的有效性和实用性.     

基于有限差分正演的带地形三维大地电磁反演方法

本研究实现了一套基于有限差分（FD）方法的大地电磁测深数据带地形三维反演算法及代码.其中,在大地电磁场正演数值模拟方面,开发了起伏地形条件下基于交错网格剖分、有限差分方法的大地电磁测深三维正演代码;在满足平面波场假设的前提下,使用长方体网格剖分模拟三维起伏地形,实现了带地形三维正演计算;并设计理论模型进行试算,经试算结果与前人的有限元法计算结果对比,验证了所研发的带地形三维正演计算的正确性与可靠性.在反演方面,本研究基于非线性共轭梯度方法编写了大地电磁测深带地形三维反演代码,试验了不同的共轭梯度搜索因子β,避免了目标函数对海森矩阵（参数二次导数矩阵）的显式计算和存储,初步实现了大地电磁资料的带地形三维反演.最后,对一系列理论模型进行正演计算,利用其生成的合成数据模拟实测数据进行反演,并与现有的不带地形大地电磁测深三维反演结果比较,检验了所研发的带地形三维反演计算的可靠性与稳定性.     

新近月球重力场模型与地形模型的局部导纳和相关性分析

本文利用三个高阶重力场模型LP150Q、GLGM-3和SGM150j以及嫦娥地形模型CLTM-s01,在频率域内使用固定窗口的方法,研究了不同重力场模型的重力/地形局部导纳谱与局部相关谱的全球分布,以及典型质量瘤盆地重力/地形局部导纳谱与局部相关谱分布的特征.结果表明:加入远月面重力场信息有助于增强重力场模型在中高频段的信号强度;三个重力场模型LP150Q、GLGM-3和SGM150j均适合于作近月面重力/地形局部导纳与局部相关性的分析和近月面地球物理参数的估计;重力场模型LP150Q和GLGM-3较适合作远月面大范围的研究,不适合于作远月面重力/地形局部导纳谱与局部相关谱的分析以及作远月面局部区域地球物理参数的估计,SGM150j较适合于作远月面局部区域地球物理参数的估计;近月面大型质量瘤盆地异常质量的尺度比远月面质量瘤盆地异常质量的大,而深度也比远月面的深.     

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.     

The role of topography in controlling throughflow generation

Previous studies of the downslope movement of soil moisture have been largely confined to two-dimensional slope profiles—the influence of contour curvature (three-dimensional) on moisture movement has received little attention. An automatic tensiometer system has been used to monitor soil moisture status continuously in a single hollow and adjacent spurs. Maximum saturated hillslope flow is shown to coincide with the stream discharge ‘throughflow peak’. Relative to the hollow, little saturation is evident at the base of the spurs. This subcatchment model is applied to a major portion of the catchment, showing the general importance of topography in the control of stream discharge.     

Spectral analysis of quasi-stationary sea surface topography from GRACE mission

The mean sea surface height (MSSH) refers to the average of the long-term sea height. The quasi-sea surface topography (QSST) is usually defined as the height difference between the MSSH and the geoid. As to 100 years of time yardstick of geodesy, the time that satellite altimetry data sets spanned is relatively shorter, in this paper, the QSST refers to the residual sea surface height (RSSH) that shows the height dif-ference between MSSH derived from altimetry and the geoid[1]. As w…     

What drives high flow events in the Swiss Alps? Recent developments in wavelet spectral analysis and their application to hydrology

Extreme hydrological events are often triggered by exceptional co-variations of the relevant hydrometeorological processes and in particular by exceptional co-oscillations at various temporal scales. Wavelet and cross wavelet spectral analysis offers promising time-scale resolved analysis methods to detect and analyze such exceptional co-oscillations. This paper presents the state-of-the-art methods of wavelet spectral analysis, discusses related subtleties, potential pitfalls and recently developed solutions to overcome them and shows how wavelet spectral analysis, if combined to a rigorous significance test, can lead to reliable new insights into hydrometeorological processes for real-world applications. The presented methods are applied to detect potentially flood triggering situations in a high Alpine catchment for which a recent re-estimation of design floods encountered significant problems simulating the observed high flows. For this case study, wavelet spectral analysis of precipitation, temperature and discharge offers a powerful tool to help detecting potentially flood producing meteorological situations and to distinguish between different types of floods with respect to the prevailing critical hydrometeorological conditions. This opens very new perspectives for the analysis of model performances focusing on the occurrence and non-occurrence of different types of high flow events. Based on the obtained results, the paper summarizes important recommendations for future applications of wavelet spectral analysis in hydrology.     

基于矢量有限元的带地形大地电磁三维反演研究

研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法（BICG）计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,进而将PARDISO的快速正演方案应用于"拟正演"问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.     

Spatial patterns of water surface topography at a river confluence

Understanding flow structures in river confluences has largely been the product of interpretations made from measured flow velocity data. Here, we turn the attention to the investigation of the patterns of both the average and standard deviations of the micro‐topography of the water surface at an asymmetrical natural discordant confluence for different flow conditions. Water surface topography is measured using a total station to survey the position of a reflector mounted on a custom‐built raft. To limit error problems related to changes in the water level, measurements are taken and analysed by cross‐stream transects where five water surface profiles are taken before moving to the next transect. Three‐dimensional numerical simulations of the flow dynamics at the field site are used to examine predicted water surface topography for a steady‐state situation. The patterns are interpreted with respect to flow structure dynamics, visual observations of boils, and bed topography. Results indicate that coherent patterns emerge at the water surface of a discordant bed confluence for different flow conditions. The zone of stagnation and the mixing layer are characterized by super‐elevation, a lateral tilt is present at the edge of the mixing layer, and a zone of super‐elevation is present on the tributary side at the downstream junction corner. The latter seems associated with periodical upwelling and is not present in the numerical simulations that do not take into account instantaneous velocity fluctuations. Planform curvature, topographic steering related to the tributary mouth bar, and turbulent structures associated with the mixing layer all play a key role in the pattern of both the average and standard deviation of the water surface topography at confluences. Copyright © 2002 John Wiley & Sons, Ltd.     

The impact of land use change on soil water holding capacity and river flow modelling in the Nakambe River, Burkina-Faso

The annual hydrological regime of the Nakambe River shows substantial changes during the period 1955–1998 with a shift occurring around 1970. From 1970 to the mid-1990s, despite a reduction in rainfall and an increase in the number of dams in the basin, average runoff and maximum daily discharges increased. This paper reviews the hydrological behaviour of the Nakambe River from 1955 to 1998 and examines the potential role of land use change on soil water holding capacity (WHC) in producing the counter-intuitive change in runoff observed after 1970. We compare the results of two monthly hydrological models using different rainfall, potential evapotranspiration and WHC data sets. Model simulations with soil WHC values modified over time based upon historical maps of land use, are compared against simulations with a constant value for WHC. The extent of natural vegetation declined from 43 to 13% of the total basin area between 1965 and 1995, whilst the cultivated areas increased from 53 to 76% and the area of bare soil nearly tripled from 4 to 11%. The total reduction in WHC is estimated to range from 33 to 62% depending on the method used, either considering that the WHC values given by the FAO stand for the environmental situation in 1965 or before. There is a marked improvement in river flow simulation using the time-varying values of soil WHC. The paper ends with a discussion of the role of other factors such as surface runoff processes and groundwater trends in explaining the hydrological behaviour of the Nakambe River.     

Quantitative assessment on basin-scale hydrological services of wetlands

Despite recognizing the importance of hydrological function of wetlands, basin-scale wetlands services have rarely been investigated. The PHYSITEL/HYDROTEL modelling platform was used to quantitatively assess the impact of wetlands on quickflow and baseflow with paired simulation scenarios in Duobukuli River Basin, namely with wetlands and without wetlands.Simulation results showed that wetlands exert significant impact on basin hydrological processes by decreasing streamflow and altering streamflow regime(magnitude, frequency, duration and time of flow events). The intensity(significant or not) of wetlands influences on quickflow had daily, monthly and annual variation. Wetlands significantly attenuated quickflow during flood season while slightly support daily, monthly and annual baseflow. The average quickflow attenuation and baseflow support of wetlands were 5.89% and 0.83%, respectively. Although the intensity and effect(mitigation or augment) of wetlands on streamflow temporally varied at daily, monthly, seasonal and annual scales, wetland overall mitigated quickflow and augment baseflow in Duobukuli River Basin. Our results could provide insights for future decision-making for rehabilitation and conservation of wetlands as well as integrated basin water resources management.     

模型扩展和地形平化起伏地形处理方案在地震走时层析成像中的对比研究

中国大陆中西部普遍具有强烈的地形起伏,起伏地形会对地震资料的处理分析产生严重干扰.精细处理起伏地形成为高精度地震成像的必然要求.传统方法通过填充低速介质将不规则模型扩展为规则模型来处理起伏地形.近年来,借助坐标变换将物理空间不规则模型转换为计算空间规则模型的地形平化方法,为解决起伏地形问题提供了新思路.本文基于经典的模型扩展和新发展的地形平化方法分别处理起伏地形,从走时正演、射线追踪和反演成像三个方面,全面细致地评判了两种地形处理方法在起伏地形层析成像中的适用性和有效性.结果表明,模型扩展中阶梯状近似和填充介质速度参与计算,会造成起伏地形走时计算精度损失,出现虚假射线路径和错误出射角,导致反演分辨率降低,成像结果模糊甚至失真;地形平化中采用贴体网格参数化,能够保证离散模型完全匹配起伏地形,并且保持起伏地形在物理空间和计算空间中均为自由表面.在此基础上发展的层析成像技术具有高度的保真性,有效地处理了地形起伏效应,为起伏地形区域精细速度成像提供了有力的技术保障.

     

Effects of surface topography on seismic ground response in the Egion (Greece) 15 June 1995 earthquake

The Greek coastal town of Egion on 15 June 1995 was shaken by a strong, small epicentral distance, earthquake that caused heavy damages to buildings and loss of life. The damages were concentrated in the central elevated part of the town whereas the flat coastal region remained almost intact. This non-uniform distribution of damage is studied in this article in terms of surface topography effects by conducting seismic response analyses of a simplified 2-D profile of the town. A dynamic finite element code implementing the equivalent-linear soil behavior (FLUSHPLUS) was used for the analyses and it was found that the step-like topography amplified greatly the intensity of motion without affecting its frequency content. The analyses showed that the motion recorded by an accelerograph installed at the center of the town is in agreement with the computed values; they also indicated a particularly intense amplification close to the crest of the steep slope, where a multi-story RC residential building partially collapsed. In contrast, the level of motion was found to be low at the flat coastal zone of the town where the earthquake damages were insignificant. It is concluded that the characteristic surface topography of the town played an important role in modifying the intensity of base motion.     

Directional statics in common reflection surface stack for rugged surface topography

Seismic data acquired along rugged topographic surfaces present well‐known problems in seismic imaging. In conventional seismic data processing, datum statics are approximated by the surface consistence assumption, which states that all seismic rays travel vertically in the top layer. Hence, the datum static for each single trace is constant. In case this assumption does not apply, non‐constant statics are required. The common reflection surface (CRS) stack for rugged surface topography provides the capability to deal with this non‐vertical static issue. It handles the surface elevation as a coordinate component and treats the elevation variation in the sense of directional datuming. In this paper I apply the CRS stack method to a synthetic data set that simulates the acquisition along an irregular surface topography. After the CRS stack, by means of the wavefield attributes, a simple algorithm for redatuming the CRS stack section to an arbitrarily chosen planar surface is performed. The redatumed section simulates a stack section whose acquisition surface is the chosen planar surface.     

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免, 其对直流电阻率法勘探影响巨大.近年来, 电阻率三维正演取得诸多进展, 特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟, 但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格, 并考虑到应用GPS/GNSS时, 区域地球物理调查中可非规则布设测网的实际特点, 实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性, 可应用于复杂野外环境下的三维电法勘探.     

Dynamic mechanisms controlling the topography of Longmenshan area

The Longmenshan fault, which defines the eastern edge of the Tibetan Plateau, is one of the steepest margins of the plateau with a sharp elevation drop of about 4 km over a distance less than 100 km across the Longmenshan fault. The mechanism which is responsible for controlling and maintaining the elevation difference is highly debated. Using multiple observations including seismic velocity model, Moho depth, effective elastic thickness of the lithosphere, we conducted a quantitative study for elucidating the contributions from crust and lithospheric mantle by an integrated analysis of lithospheric isostasy and flexure. It is shown that the topography of the Longmenshan fault is supported by both lithospheric isostasy and flexure statically, and lower crustal channel flow and mantle convection dynamically. Different mechanisms have different weights for contribution to the topography of the Songpan-Ganzi block and the Sichuan Basin. The static and dynamic support contribute roughly the same to the topographic difference of ~4 km between the two sides of the Longmenshan fault. The static topographic difference of ~2 km is mainly resulted from the lithospheric isostasy, while the dynamic one of ~2 km is contributed by the uprising of the accumulated material in the lower crust beneath the Songpan-Ganzi block and the downward drag force caused by the upper mantle convection under the Sichuan Basin. It is thus suggested that the lower crustal flow and upper mantle convection are dynamic forces which should be taken into account in the studies on the dynamics in the Longmenshan and surrounding regions.     

Numerical evaluation of slope topography effects on seismic ground motion

This paper presents results of numerical analyses for the seismic response of step-like ground slopes in uniform visco-elastic soil, under vertically propagating SV seismic waves. The aim of the analyses is to explore the effects of slope geometry, predominant excitation frequency and duration, as well as of the dynamic soil properties on seismic ground motion in a parametric manner, and provide qualitative as well as quantitative insight to the phenomenon. Among the main conclusions of this study is that this kind of topography may lead to intense amplification or de-amplification variability at neighboring (within a few tens of meters) points behind the crest of the slope, especially for high frequency excitations. Nevertheless, a general trend of amplification near the crest and de-amplification near the toe of the slope seems to hold for the horizontal motion. As a result of these two findings, it becomes evident that reliable field evidence of slope topography aggravation is extremely difficult to establish. Furthermore, this study highlights the generation of a parasitic vertical component of motion in the vicinity of the slope, due to wave reflections at the slope surface, that under certain preconditions may become as large as the horizontal. Criteria are established for deciding on the importance of topography effects, while approximate relations are provided for the preliminary evaluation of the topographic aggravation of seismic ground motion and the width of the affected zone behind the crest.     

Prediction of landslide occurrence based on slope-instability analysis and hydrological model simulation

Hillslope failure usually occurs as soil resistance deteriorates in the presence of the acting stress developed by a rising groundwater level during rainstorms. The present study adopted a slope-instability analysis and a hydrological model for landslide prediction during heavy rainstorms. Variation of the groundwater table on hillslope was simulated by using the hydrological model and then the temporal groundwater level at each grid was substituted into the slope-instability analysis to determine the instability of the grids in watersheds for prediction of massive landslides.Hydrological records from two landslide-prone areas in northern Taiwan were collected. Digital elevation model was adopted to obtain the geomorphologic factors required for the slope-instability analysis and the hydrological model. The spatial distribution of soil thickness required for performing the infinite slope model was estimated by using a wetness index. Results showed that the temporal variation of the percentage of unstable grids in the study watersheds basically followed the variation of rainfall hyetographs. The percentage of the unstable grids reached a maximum value when the centroid of the hyetograph passed. A comparison between the landslide records and the model analytical results revealed that a massive landslide might occur if more than 50% of the grids in the subwatershed were classified as unstable in the study areas. The predicted time and location of landslide occurrence were consistent with those obtained from field investigations. It is therefore considered promising to apply the developed analytical method for landslide warning to alleviate the loss of lives and property.     

Quantitative analysis of agricultural land use change in China

This article reviews the potential impacts of climate change on land use change in China. Crop sown area is used as index to quantitatively analyze the temporal–spatial changes and the utilization of the agricultural land. A new concept is defined as potential multiple cropping index to reflect the potential sowing ability. The impacting mechanism, land use status and its surplus capacity are investigated as well. The main conclusions are as following;

• 1.During 1949–2010, the agricultural land was the greatest in amount in the middle of China, followed by that in the country’s eastern and western regions. The most rapid increase and decrease of agricultural land were observed in Xinjiang and North China respectively, Northwest China and South China is also changed rapid. The variation trend before 1980 differed significantly from that after 1980.
• 2.Agricultural land was affected by both natural and social factors, such as regional climate and environmental changes, population growth, economic development, and implementation of policies. In this paper, the effects of temperature and urbanization on the coverage of agriculture land are evaluated, and the results show that the urbanization can greatly affects the amount of agriculture land in South China, Northeast China, Xinjiang and Southwest China.
• 3.From 1980 to 2009, the extent of agricultural land use had increased as the surplus capacity had decreased. Still, large remaining potential space is available, but the future utilization of agricultural land should be carried out with scientific planning and management for the sustainable development.

     

Impacts of land use and land cover changes on hydrological processes and sediment yield determined using the SWAT model

Land use and land cover (LULC) changes strongly affect local hydrology and sediment yields.The current study focused on a basin in the Brazilian Amazon and had the following three objectives:(1) to perform an effective diagnosis of flow and sediment yield,(2) to evaluate the impacts of LULC changes over the last 40 years on the hydro-sedimentological variables,and (3) to investigate the impacts of the possible trends or breaking points in the flow,surface runoff,and sediment yield series.The Soi...