Influence of southern oscillation index on the variability and predictability of Indian monsoon. A reappraisal

The relationship between the monsoon rainfall throughout all India, northwest India and peninsular India as well as the onset dates of the monsoon and two indices of southern oscillation (SOI), namely Isla de Pascua minus Darwin (I-D) and Tahiti minus Darwin (T-D) pressure anomaly have been studied for different periods. The study indicates that the monsoon rainfall shows a strong and significant direct relationship with SOI for the concurrent, succeeding autumn and succeeding winter seasons. The magnitude of the direct correlation coefficient for the SOI using (I-D) is enhanced over all India and peninsular India if the above seasons happen to be associated with an easterly phase of the QBO (Quasi-Biennial Oscillation) at 50 mb. The result indicates that the strength of the monsoon plays an important role in the following southern oscillation events in the Pacific Ocean. The premonsoon tendency of the SOI anomaly spring minus winter SOI shows a significant positive correlation with monsoon rainfall over all India, northwest India and peninsular India. The absolute value of the positive correlation coefficient becomes highly enhanced over all India, northwest India as well as peninsular India if the 6-month period from December to March is associated with the westerly phase of the QBO. Hence, the premonsoon SOI tendency parameter can be a useful predictor of Indian monsoon rainfall especially if it happens to be associated with the westerly QBO. Significant negative association is also found between the anomaly of monsoon onset dates and SOI of the previous spring season, the absolute value being higher for SOI (T-D) than for SOI (I-D). The negative correlation coefficient becomes enhanced if the previous springs are associated with a westerly phase of the QBO. It shows that the previous spring SOI has some predictive value for the onset date of Indian monsoon, a positive SOI followed by an early onset of monsoon, andvice versa, especially if it is associated with a westerly phase of the QBO.     

Evaluating the effect of land development on sediment transport using a probability density function

An important problem in sedimentation analysis is the development of a channel section that preserves, as best as possible, the current sedimentation regime even though the flood frequency tendencies have been altered due to land development within the catchment. In order to accomplish this task, a methodology is needed that estimates sediment transport capacity for various channel configurations. Such a procedure is described which allows the computation of the total sediment transport capacity for each of several T-year return frequency runoff hydrographs. This information is used to obtain an approximate probability distribution for the total sediment transport capacity, and the mean and standard deviation of this distribution are computed.Comparing the results for the catchment in its present state with a future developed state, using a selection of new channel parameters, indicates how to improve the channel to control changes in sedimentation due to development. The analysis procedure provides a basis for estimating a new channel configuration such that the new flow conditions retain, as best as possible, the existing condition sedimentation effects, and hence retain the natural sediment supply and transport trends even though runoff flow rates have changed due to land development within the catchment.The results of Wilson Creek are typical of the several sites examined, see Table 3 below. The T=2, T=5, T=25, and T=100 year values for total sediment transport capacity, in kilotons, are 6.9, 39.4, 61.3, and 96.7 with a mean of 17.1 and standard deviation of 19.3. After development with no change in the channel the respective values increase to: 17.9, 84.6, 128.1, and 258.0 with a mean of 39.1 and standard deviation of 44.3. A new channel can be constructed which will reduce these sediment transport capacity values, after development, to 5.2, 41.0, 62.0, and 124.8 with a mean of 17.4 and standard deviation of 22.0.     

重力资料的三维密度反演

分别论述了在空间域与频率域中由重力资料反演三维密度分布的理论及其实现的方法,并提出了一种简便、实用的频率域三维密度反演算法。     

Sq的经度效应

本文提供了一个事实:同纬度不同经度台站的S_q曲线可以有很大的差异。对拉萨台的S_q异常作了一些初步的、定性的分析与解释。对八丈台的S_q异常则认为难以解释,值得深入研究。     

基于Swarm卫星GPS/TEC观测的顶部电离层三维电子密度分布的层析法研究

利用两颗伴飞的Swarm A/C卫星搭载的双频GPS接收机获取的TEC数据,在两个卫星轨道平面同时对顶部电离层电子密度进行层析成像,实现对顶部电离层电子密度的三维观测.为了能够重现扰动期间电离层电子密度的空间变化特征,在正则化求解过程中,我们引入了水平矩阵H和垂直矩阵V刻画电子密度的空间变化特征,引入整体约束矩阵C以调节不同空间对电子密度相对变化的权重.数值验证结果表明我们的算法对常见的观测误差具有较强的包容性,反演计算出的电子密度平均偏差优于10%.在不同地磁活动条件下,与第三方观测数据的对比,验证了本文反演算法的可靠性.实测数据反演结果表明我们的算法不仅能够较好地重现顶部电离层子午向百公里级别的不规则结构,还能有效分辨纬向相隔~150 km的两个卫星轨道平面的电子密度差异.     

改进的约束变密度界面反演策略及其应用

计算密度分界面的起伏变化在区域地质构造研究和石油矿产资源勘探中具有重要意义.已有密度界面反演方法更多侧重约束变密度界面反演算法,而对约束信息的准确性、研究区横向密度变化往往考虑不足,影响了最终反演结果的可信度.本文在变密度界面正反演算法基础之上,结合实际需求,提出已知深度信息约束下的变密度界面反演策略.该策略主要包括变密度约束反演算法、已知深度约束信息校验和分区变密度模型三个方面.其中反演算法提供了带已知深度约束信息的密度界面迭代反演方法;约束信息校验用于评估约束信息精度,通过调差降低约束信息的系统误差;在反演过程中引入水平密度分区以应对不同构造背景密度界面模型,提高反演结果的可信度.最后将本策略应用于南海莫霍面深度反演计算中,结果显示借助已知约束信息,利用分区密度模型能够获得更为可信的深度反演结果,验证了该策略的正确性.     

天山造山带及邻区孕震应力场及断裂带活动性特征的数值实验分析

天山造山带构造环境复杂,活动断裂带和强震分布广泛,且主要分布于阿尔泰山、天山、西昆仑—帕米尔弧形构造带上,尤以天山地区最为集中.迄今为止,天山造山带地区的主要断裂带的活动特征与孕震应力场特征之间的动力学机理尚未有清晰的认识.本文以GPS等实际观测数据为约束,建立有限元数值模型,计算了研究区域地壳形变、应力/应变积累速率、弹性应变能密度以及库仑应力变化率等关键因素.模拟计算结果显示地表速度场与研究区域实际GPS观测值基本一致,且主要断裂带上弹性应变能密度分布与实际地震活动性也基本吻合,验证了数值模型和结果的可靠性.结合最新的观测和数值模拟结果分析发现,研究区的断层和地震活动性主要受控于近南北向的主压应力,与主要观测特征相一致.同时,帕米尔高原北部边界带—塔什库尔干断裂（TKF）、天山造山带南边界的东侧—迈丹断裂（MDF）、兴地断裂（XDF）库仑应力增大明显,在未来强震发生的可能性较高,应予密切关注.     

基于改进全变差正则化的GPR多尺度全波形双参数同步反演

针对探地雷达（GPR）双参数全波形反演中电导率反演精度差、双参数存在串扰现象、反演计算量大、易陷入局部极值等问题.作者将具有多参数调节功能的L-BFGS算法引入到GPR时间域全波形反演中,它避免了对Hessian矩阵的直接存储与精确求解,减小了存储量和计算量.结合参数调节因子的选取,有效减小了同步反演时介电常数与电导率的串扰影响,在不降低介电常数反演精度的前提下,提高电导率参数的反演精度.通过在反演目标函数中加载改进全变差正则化方法,提高了反演的稳定性,使目标体边缘轮廓更加清晰.首先以简单模型为例,对比了单尺度反演与多尺度串行反演策略的优劣,说明多尺度串行反演有利于逐步搜索全局最优解;而开展参数调节因子的选取实验,说明合适的参数调节因子可以有效改善介质电导率的反演精度;测试了不同正则化的反演效果,表明改进全变差正则化能提高反演稳定性,显著降低模型重构误差.最后,分别对含噪合成数据和实测数据进行了反演测试,说明本文提出的多尺度、双参数反演具有较强的鲁棒性,能提供更丰富的信息约束,重构图像界面清晰、反演效果好.     

构造块体的相对运动和应变

在确定构造块体在球面上相对运动参数的基础上,进一步给出了描述构造块体内的应变变化及参数和确定这些参数的方法。讨论了将相对运动和应变变化综合起来以更充分地描述地壳运动和变形的模型。结合相对运动参数,提出了一种逐渐趋近的变形分析方法。     

Bjerhammar球面单层密度时间变化的实用解算理论

研究了Ｂｊｅｒｈａｍｍａｒ球面单层密度时变实用解算理论和稳定算法。从利用重复大地测量数据联合确定Ｂｊｅｒｈａｍｍａｒ球面单层密度时变基本理论出发,推导了不同Ｂｊｅｒｈａｍｍａｒ球面上ｑ＊ｔ和δＷ＊ｔ之间的一般解式;分析研究了ｑ＊ｔ的迭代求解过程和收敛特性,得出其收敛稳定性很大程度上取决于Ｒ值的结论,同时给出了界限Ｒ的关系式;提出了解求ｑ＊ｔ的分步迭代算法,该算法从理论上解决了以前算法存在的收敛稳定问题（Ｒ较小时）,并指出：地表观测点ｑｔ的确定应遵循匹配处理准则,即各观测点数据应蕴含相同物理过程。最后给出了ｑ＊ｔ的精度评估公式。