微分算子展开与耦合势三维有限体积法快速计算随钻超深前视电阻率测井响应非线性Born逼近

本文将非均质各向异性地层中电磁场耦合势Helmhotz方程分解为各向异性背景场方程与散射场方程,并应用算子展开技术推导出耦合势Born级数解,建立一套三维非均质各向异性地层中随钻超深前视电阻率测井响应各阶散射电磁场与非线性Born逼近的有效计算方法.首先,利用传输线法与二维插值技术确定背景电磁场空间分布,得到零阶Born逼近解,然后根据电磁场耦合势的Born级数解推导出各阶散射电磁场方程与递推关系,通过三维有限体积法实现各阶散射场耦合势方程的离散,在此基础上结合MKL PARDISO并行技术,通过递推方式逐步计算各阶散射电磁场,并根据Born级数部分和确定不同阶次的非线性Born逼近.最后利用数值模拟结果对正演算法加以检验,并深入研究分析各阶散射场的变化特征与Born级数的收敛性质.

     

A 3D multi-threshold,region-growing algorithm for identifying dust storm features from model simulations

Dust storms cause significant damage to health, property, and the environment worldwide every year. To help mitigate the damage, dust forecast models simulate and predict upcoming dust events, providing valuable information to scientists, decision makers, and the public. These simulation outputs are in four-dimensions (i.e., latitude, longitude, elevation, and time) and represent spatially heterogeneous dust storm features and their evolution over space and time. This research investigates and proposes an automatic multi-threshold, region-growing-based algorithm to identify critical dust storm features from 3D dust storm simulations. A multi-threshold scheme is defined for the identification of dust storm features with different dust concentrations. Based on the multi-thresholds, dust storm features are iteratively identified by developing a region-growing algorithm that splits a clustered dust storm feature into multiple sub-features. The proposed approach is compared with three commonly used methods in image processing and thunderstorm identification. The proposed approach outperforms the other three methods in sensitivity and quantitative/qualitative accuracy. This research approach may also be slightly adjusted to identify critical 3D features from simulation outputs for other severe weather and geographical phenomena.     

Comparison of ecophysiological characteristics between introduced and indigenous mangrove species in China

Due to its rapid growth, the introduced mangrove species Sonneratia apetala from Bangladesh has been widely used in mangrove restoration in southeastern China since 1985. As an indigenous mangrove species in Hainan, China, Sonneratia caseolaris was also planted in Guangdong Province for afforestation purposes. Both species have developed well in their new habitats, but their ecophysiological differences with the native mangrove species have not been studied. In this study, leaf gas exchange, water and nitrogen use efficiencies of two Sonneratia species were compared with those of selected native mangrove species (Avicennia marina, Aegiceras corniculatum, Kandelia candel, and Excoecaria agallocha) in Hainan and Shenzhen. The introduced S. apetala maintained lower carbon assimilation rate (A) and photosynthetic nitrogen use efficiency (PNUE) than the indigenous S. caseolaris. In Shenzhen, the two introduced Sonneratia had comparable photosynthetic rates and water use efficiency (WUE) with the native mangrove species, except that PNUE in S. caseolaris was significantly higher than in the native mangrove species. The two Sonneratia species showed significant overlap in PNUE and long-term WUE. Photosynthetic parameters derived from leaf photosynthetic light–response curves and A–C_i curves also suggested lower carbon assimilation capacities for the introduced Sonneratia than for the native mangrove species in both study sites. The lower light compensation point (LCP) of two introduced Sonneratia in both study sites also indicated a better adaptation to a low light regime than the native mangrove species. The results of photosynthetic capacities indicated that the introduced mangrove species have little competitive advantage over local native mangrove species in their respective new habitats.     

地震孕育、发生、发展动态过程的三维有限元数值模拟

为探讨用数值方法模拟地震孕育、发生、发展动态过程的途径,从线性流变体介质内制约质点运动的运动方程出发,导出了模拟地震孕育、发生、发展动态过程的三维有限元方程及程序。还给出了模拟进程中地震孕育、发生、发展过程的约束条件,使得可以用同一个程序完整地模拟地震孕育、发生、发展的全过程,为地震过程本质的认识以及物理预报的实现提供了一个潜在的新手段。     

地磁场能量在地球内部的分布及其长期变化

用国际参考地磁场模型(IGRF)分析了地磁场能量在地球内部的分布及其长期变化.结果表明,从1900年到2005年,地核以外地磁场总能量由6.818×1018J减少到6.594×1018J,减小了3.3%,地表以外地磁场总能量由8.658×101J减小到.63×101J,减小了11.4%.分析地球内部不同圈层地磁场能量的变化表明,地壳(A层)、上地幔(B层)、转换带(C层)、下地幔D′层的地磁场总能量在减小,但是下地幔"层的地磁场总能量却在快速增加.磁能密度随时间的变化更清楚地显示出磁能增加和减小的分界面在r=3840km处.上述结果表明,地核和地表以外地磁场总能量在趋势性减小的同时,也在进行重新分配.进一步分析表明,下地幔D"层磁能快速增长,主要是由高阶磁多极子的增强引起的.在地磁场倒转前,偶极矩减小而多极性相对增强在能量分布上的表现就是磁能向下地幔底部(特别是D"层)集中.     

BIBLIOGRAPHIE

Abstract

Time series modelling approaches are useful tools for simulating and forecasting hydrological variables and their change through time. Although linear time series models are common in hydrology, the nonlinear time series model, the Generalized Autoregressive Conditional Heteroscedasticity (GARCH) model, has rarely been used in hydrology and water resources engineering. The GARCH model considers the conditional variance remaining in the residuals of the linear time series models, such as an ARMA or an ARIMA model. In the present study, the advantages of a GARCH model against a linear ARIMA model are investigated using three classes of the GARCH approach, namely Power GARCH, Threshold GARCH and Exponential GARCH models. A daily streamflow time series of the Matapedia River, Quebec, Canada, is selected for this study. It is shown that the ARIMA (13,1,4) model is adequate for modelling streamflow time series of Matapedia River, but the Engle test shows the existence of heteroscedasticity in the residuals of the ARIMA model. Therefore, an ARIMA (13,1,4)-GARCH (3,1) error model is fitted to the data. The residuals of this model are examined for the existence of heteroscedasticity. The Engle test indicates that the GARCH model has considerably reduced the heteroscedasticity of the residuals. However, the Exponential GARCH model seems to completely remove the heteroscedasticity from the residuals. The multi-criteria evaluation for model performance also proves that the Exponential GARCH model is the best model among ARIMA and GARCH models. Therefore, the application of a GARCH model is strongly suggested for hydrological time series modelling as the conditional variance of the residuals of the linear models can be removed and the efficiency of the model will be improved.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Modarres, R. and Ouarda, T.B.M.J., 2013. Modelling heteroscedasticty of streamflow times series. Hydrological Sciences Journal, 58 (1), 1–11.     

地磁低点位移法的改进

中低纬度地磁场平静太阳日变化Z分量的极小值出现的时刻常被用于地震预报,即"低点位移"法,但现有算法的时间格式和人工绘制突变界线的模式会产生较大的误差,本文在第一部分针对这些误差改进了原有算法.第一项是改进了地磁低点时间的格式,第二项是引入了Voronoi剖分技术替代人工画线.前者过滤了假的突变界线,消除了系统误差;后者建立了突变界线与地磁台网低点数据之间的解析关系,使它们一一对应,消除了随机误差.本文第二部分内容是提出了新的时间格式与数字化突变界线的计算公式,并用程序实现计算任务.本文第三部分内容是在对地磁台网的低点时间跟踪的过程中,捕捉到了甘肃省周边几次地震前的地磁低点位移现象,并分析这些案例,在结论中证明该方法的可靠性.     

基于地质记录用时域组合模型预测气候变化趋势的初步研究

用时域分析组合模型建立了１００万年来６０°Ｎ６月份太阳辐射量时间序列、宝鸡黄土粒度时间序列、渭南夏季风指数时间序列的动态模型．研究结果表明,时域分析组合模型较好地提取了时间序列的信息,得到的几个显著周期Ｔ＝１３３,１００,８９,４１,２３,１９ｋａ,与地球轨道三要素的变化周期接近．用组合模型拟合实测数据,精度是高的;用其预测未来气候替代性指标时间序列的变化情况,发现未来气候有向干冷方向发展的趋势．时域分析组合模型为研究气候变化趋势提供了一种定量分析、预测的方法．     

地球重力场、重力、重力梯度在三维直角坐标系中的表达式

现代卫星重力测量主要利用星载GPS接收机、加速度计、星载测距仪等来确定重力卫星的轨道,削弱非保守力的干扰,由此根据卫星的位置、速度及其变率来确定地球重力场。而上述GPS等星载仪器所提供的数据,包括卫星轨道坐标及其速率、扰动加速度、星间距离及其变率,都是以三维直角坐标(x,y,z)的形式表示的,因此,地球重力场、重力和重力梯度在三维直角坐标系中的表达式在卫星重力解算中具有实际意义。     

基于时间序列统计特性的森林变化监测

森林动态变化分析对揭示生态系统环境变化及植被恢复和布局重建等具有重要意义,时间序列的遥感数据为森林监测提供了基础数据。本文根据森林植被的统计学特性,在暗目标法的基础上,利用归一化植被指数NDVI实现森林样本自动选择;并融合NDVI构建了新的综合森林特征指数(Integrated Forest Z-Score,IFZ);以时间序列的IFZ分析森林动态信息,实现森林变化动态监测。以三峡大坝及周边区域森林为研究区,利用2001年至2012年每年生长季节(5月—10月)的Landsat TM影像检验本文算法。基于2002年、2006年和2010年三期7月—9月的Quick Bird影像的精度分析结果发现:研究区森林变化检测的总体精度可达96.53%,Kappa系数为0.9512。在添加NDVI指数后构建的IFZ提高了总体监测精度。其中,毁林类别的检测精度提高显著,漏检率和误检率分别为2.74%和3.64%;干扰后重建的森林类别的检测精度有一定提高,其漏检率和误检率分别为10.79%和10.51%。研究结果表明,改进暗目标法能提高森林样本的选样效率,添加NDVI的IFZ能提高森林动态变化的识别度。此外,本算法不仅能定性识别森林变化,而且能定量提供森林干扰发生时间和干扰强度。