2016-2020年黄河三角洲尾闾河道悬浮泥沙时空分布

三角洲河道悬浮泥沙的时空分布与河口流路稳定性及水下三角洲发育密切相关。基于Sentinel-2高分辨率卫星影像反演悬沙浓度,结合河道形态和水动力数值模拟,分析了近年来黄河尾闾河道悬沙的时空分布特征和影响因素。结果表明：2016-2020年,黄河尾闾河道平均悬沙浓度具有显著的洪枯季差异和不同的沿程分布;洪季悬沙浓度受上游来沙和底沙起动共同控制,王家庄和清加2断面是底沙起动再悬浮的主要断面;径流流速的增大对不同河段悬沙浓度的影响不尽相同,河口段悬沙浓度更易受径流流速的影响;以王家庄断面为界,上下游悬沙浓度分别呈现“深泓低两岸高”和“深泓高两岸低”两种截然不同的分布模式;尾闾河道悬沙浓度对调水调沙响应敏感,2018年调水调沙期间悬沙浓度达到10.75 g/L,但随着河床底沙粗化,尾闾河道悬沙浓度对调水调沙的响应敏感性下降。     

Modified state–space procedures for pseudodynamic testing

The existing on‐line numerical integration algorithms are derived from the Newmark method, which is based on an approximation of derivatives in the differential equation. The state–space procedure (SSP), based on an interpolation of the discrete excitation signals for piecewise convolution integral, has been confirmed as more reliable than the Newmark method in terms of numerical accuracy and stability. In an attempt to enhance the pseudodynamic test, this study presents an on‐line integration algorithm (referred to as the OS–SSP method) via an integration of the state–space procedure with Nakashima's operator‐splitting concept. Numerical stability and accuracy assessment of the proposed algorithm in addition to the explicit Newmark method and the OS method were investigated via an eigenvalue, frequency‐domain and time‐domain analysis. Of the on‐line integration algorithms investigated, the OS–SSP method is demonstrated as the most accurate method with an acceptable stability (although not unconditionally stable) characteristic. Therefore, the OS–SSP method is the most desirable method for pseudodynamic testing if the numerical stability criterion (Δt/T⩽0.5) is ensured for every vibration mode involved. Copyright © 2001 John Wiley & Sons, Ltd.     

Threshold criterion for debris flow initiation in seasonal gullies

A series of flume experiments were done to investigate the effect of grain composition on the critical gradient and discharge of debris flows initiated in seasonal gullies. The results indicated that the critical gradient and discharge for debris flow initiation decrease initially, and then increase as the mass content of fine particles (<2 mm) increases. As the mass content of fine particles increases, the angle of repose, permeability of widely graded gravel soils, and the incipient motion conditions of the coarse grains in non-uniform sediments decrease at first, and then increase. The mass content of fine particles of all inflection points is the same. The theoretical model based on the combination of hydrodynamic force and shear stress is more applicable to the prediction of the critical gradient for debris flow initiation. The critical discharge model considering the effect of non-homogeneity of the soil and the size of coarse and fine grains provides a more accurate prediction of debris flow initiation than other models based on the mean diameter.     

Complexity analysis of precipitation using the Lempel–Ziv algorithm and a multi-scaling approach: a case study in Jilin province,China

Precipitation is an important part of the hydrologic cycle, and its complexity is closely related to surface runoff and changing groundwater dynamics, which in turn influences the accuracy of precipitation forecasts. In this study, we used the Lempel–Ziv algorithm (LZA) and a multi-scaling approach to assess precipitation complexity for 1958–2011 by analyzing time series data from 28 gauging stations located throughout Jilin province, China. The spatial distribution of normalized precipitation complexity was measured by LZA, a symbolic dynamics algorithm, and by a multi-scaling approach, which is described by fractals. In addition, the advantages and limitations of these two methods were investigated. The results indicate that both methods are applicable and consistent for calculating precipitation complexity, and that the degree of relief is a primary factor controlling precipitation complexity in the mountainous area; in the plain terrain, however, the prominent influencing factor is climate.     

Comparison of parameter uncertainty analysis techniques for a TOPMODEL application

Parameter uncertainty in hydrologic modeling is crucial to the flood simulation and forecasting. The Bayesian approach allows one to estimate parameters according to prior expert knowledge as well as observational data about model parameter values. This study assesses the performance of two popular uncertainty analysis (UA) techniques, i.e., generalized likelihood uncertainty estimation (GLUE) and Bayesian method implemented with the Markov chain Monte Carlo sampling algorithm, in evaluating model parameter uncertainty in flood simulations. These two methods were applied to the semi-distributed Topographic hydrologic model (TOPMODEL) that includes five parameters. A case study was carried out for a small humid catchment in the southeastern China. The performance assessment of the GLUE and Bayesian methods were conducted with advanced tools suited for probabilistic simulations of continuous variables such as streamflow. Graphical tools and scalar metrics were used to test several attributes of the simulation quality of selected flood events: deterministic accuracy and the accuracy of 95 % prediction probability uncertainty band (95PPU). Sensitivity analysis was conducted to identify sensitive parameters that largely affect the model output results. Subsequently, the GLUE and Bayesian methods were used to analyze the uncertainty of sensitive parameters and further to produce their posterior distributions. Based on their posterior parameter samples, TOPMODEL’s simulations and the corresponding UA results were conducted. Results show that the form of exponential decline in conductivity and the overland flow routing velocity were sensitive parameters in TOPMODEL in our case. Small changes in these two parameters would lead to large differences in flood simulation results. Results also suggest that, for both UA techniques, most of streamflow observations were bracketed by 95PPU with the containing ratio value larger than 80 %. In comparison, GLUE gave narrower prediction uncertainty bands than the Bayesian method. It was found that the mode estimates of parameter posterior distributions are suitable to result in better performance of deterministic outputs than the 50 % percentiles for both the GLUE and Bayesian analyses. In addition, the simulation results calibrated with Rosenbrock optimization algorithm show a better agreement with the observations than the UA’s 50 % percentiles but slightly worse than the hydrographs from the mode estimates. The results clearly emphasize the importance of using model uncertainty diagnostic approaches in flood simulations.     

Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space

The scattering of plane SH-waves by topographic features in a layered transversely isotropic (TI) half-space is investigated by using an indirect boundary element method (IBEM). Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green’s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green’s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.     

汶川地震震前与同震断层滑移的比较分析

本文利用GPS观测的1999-2007年汶川震前3期地表变形数据和2008年汶川同震地表变形数据,结合地震位错理论,通过高斯变换和坐标旋转建立断层模型,运用遗传算法,反演了龙门山断裂带断层震前3期和同震滑动参数。结果表明龙门山断层震前3期平均走滑位移为-5.39mm,倾向位移为2.66mm,与同震断层滑移相比较,发现震前断层的滑移趋势与同震断层滑移一致,均为逆冲兼右旋的挤压运动。比较震前3期逆冲方向的滑移量,发现逆冲滑移有加速的现象。并根据震前和同震的断层滑动量估算了汶川地震复发周期。     

Retrieval of middle and upper atmospheric wind based on non-full circular fringe recorded by Fabry-Perot Inteferometer

Fabry-Perot Interferometer(FPI) has been used widely for wind measurements of the middle and upper atmosphere.To date, most of FPIs have been based on full-closed circular fringe, which needs 15–25 min to obtain a group of wind velocity(zonal and meridional). However, it is hard to improve the temporal resolution because full-closed circular fringe in several directions cannot be easily imaged onto the same Charge-Coupled Device(CCD) with enough airglow intensity. In this paper, a data processing method is proposed for non-full circular fringe of FPI, which can support CCD with enough area of observations in several directions simultaneously. The method is focused on the center determination of non-full fringe. It includes radial cross-section, peak coordinate determination, and center calculation. Based on the calculated center, the fringe is annular summed. Then its radius is determined subsequently using Gaussian fitting. Finally, the wind is retrieved from the fringe radius. For validation, fringes from two ground-based FPIs were used, which are deployed in Kelan(38.71°N, 111.58°E) and Xinglong(40.40°N, 117.59°E) in China. The results retrieved from non-full fringes of FPIs were compared with that from full-closed circular fringe. The averaged wind deviation between them demonstrates reasonable difference with 5.38 ms~-(1) for 892.0 nm airglow emission, 5.81 ms~-(1) for 630.0 nm emission, and 3.03 ms~-(1) for 557.7 nm emission. Besides, wind results of Xinglong FPI are compared roughly with measurements of meteor radar which is deployed in Ming Tombs of Beijing(40.3°N,116.2°E). Good agreement demonstrates that this method is robust enough for FPI wind retrieval of mesosphere and thermosphere.     

利用背景噪声成像研究合肥-金华剖面地壳速度结构及径向各向异性的东西差异

背景噪声成像能够为地壳内部结构提供地震波速的观测证据,能够增进对地壳结构与物性的认识.本文利用背景噪声成像方法对合肥-金华流动地震台阵55个台站数据进行了分析反演,获得了华北克拉通与华南块体东部构造边界附近区域的地壳S波速度与径向各向异性结构.成像结果表明,以东经118°-118.5°附近为界,西北部表现为地壳低速异常和中下地壳正的径向各向异性,东南部表现为地壳高速异常,下部地壳存在负的径向各向异性.由此,本文推测地震波速结构差异反映了地壳内部物质与温度的差异,这种差异与两个地区自中生代晚期以来经历了截然不同的岩浆活动相关.     

宽频带超导重力仪的研制进展

本文报道一种构建宽频带高分辨率超导重力仪（SG）的方法.与已商业化的GWR超导重力仪类似,新型超导重力仪同样使用磁悬浮超导检验质量来构建弹簧振子结构,不同的是弹簧振子的固有频率更高,可以有效覆盖地脉动频带.同时,利用基于超导量子干涉器件（SQUID）的位移传感技术,大幅提高了位移传感灵敏度,以保证仪器在测量带宽扩展之后仍具有与GWR仪器相当的测量分辨率,从而实现对时变重力信号与地球背景噪声信号的同时高分辨观测.宽频带超导重力仪核心部件垂向超导加速度计（VSA）的设计、组装和初步测试已经完成.结果表明,加速度计的噪声本底为8 μGal/Hz^1/2@0.1 Hz,实现了对地脉动信号的高分辨率测量.而为了将仪器噪声本底降低到GWR仪器水平,加速度计的固有频率需要进一步减小,该部分参数优化工作正在进行.本文将对垂向超导加速度计的工作原理、结构和测试结果进行详细介绍.