东亚副热带高空急流区中尺度扰动分布特征及可能机理

东亚副热带高空急流强度变化和天气气候密切相关,本文利用WRF模式输出的高时空分辨率模拟资料研究了东亚副热带高空急流区的中尺度扰动特征,并结合动力学理论,揭示了急流区中尺度扰动产生的可能机制.研究表明,急流轴南侧更容易出现水平尺度为几十公里的高频扰动,这些扰动的时空分布具有波动特征.对高空急流区中出现中尺度扰动区域的拉格朗日Rossby数、Richardson数以及绝对涡度的计算发现,高空急流轴南侧中尺度扰动出现的物理机制与非地转平衡流的不稳定发展有关,并且高空急流强度的大尺度整体变化与急流区中尺度扰动变化的累积效应有关.因此,开展高空急流强度变化规律研究不能忽视其内部中尺度动力过程的作用.     

La Coste地震计的原理及其应用

本文讨论了旋转型La coste地震计的原理和各种悬挂方式及其周期表达式。指出消除“拐腿”效应和振幅效应的最好办法是采用零长弹簧（l=0）并使绞链簧片的弹率c₁=-M′（0）否则,地震计将表现为正“拐腿”效应、正振幅效应或反效应,它们严重地阻碍着固有周期的增长。绞链簧片弹率c₁固定后,地震计只能工作在对应的周期。列举了满足La coste悬挂条件时各种因素变化对周期和零点的影响。     

声波测井时裸眼井中弹性波的传播

本文分析了简化模型下声波测井时裸眼井（内充流体一般是泥浆）中弹性波的传播,此波场可看成是轴向传播有不同波矢的分波的叠加。分析表明,对于井外只有单一均匀岩层且无限深的井,仅对如（18）式所示的特征频率成分才存在具有轴向传播最大相速度v_zm=c₂的所谓“首波”;而v_x > c₂及c₂ > v_z>b₂的分波却不存在,此处c₂和b₂分别是井外岩层的纵波和横波速度。这正解释了文献[3]中所报道的裸眼井的“过滤效应”。分析还表明,对于井外具有多层岩层的井,只在一定条件下才存在与上述特性相近的“首波”。本文结果对声速测井显然是重要的。     

热带海洋和大气中地形Rossby波和Rossby波的耦合不稳定

当大尺度背景场存在赤道急流时, 由相应不均匀的温跃层(海洋)和高度场(大气)激发出的地形Rossby波和由β效应激发出的Rossby波, 在一定条件下, 通过相互作用后可产生一类新的不稳定, 称为地形Rossby波和Rossby波的耦合不稳定. 讨论了这类波系在ENSO发展中可能起的作用.     

横向各向同性地层(VTI)井孔声弹效应对弯曲波的影响

本文发展了建立在地层参考状态为各向同性介质假定下的现行井孔声弹性理论,就井外为横向各向同性面与井轴垂直的、具有9个独立三阶弹性模量的横向各向同性介质(VTI井况),水平面内受双轴应力作用下给出了一个简洁的与井内压力、应力差、应力和以及多极源偏振方位角有关的井孔弯曲波声弹公式,并且导出了平面纵、横波速度的声弹公式.数值考察了弯曲波速度之改变量的灵敏系数随频率的变化、受井外水平双轴应力作用时两种偏振的偶极弯曲波频散曲线以及对应不同方位径向偏振的横波速度.研究结果表明弯曲波声弹公式与5个二阶弹性系数以及7个独立的三阶弹性模量有关;而且由内压引起的井孔弯曲波声弹性公式中的三阶弹性模量仅与6个独立的弹性模量有关.横向各向同性介质井孔弯曲波速度的交叉现象仍是判断地应力存在的标志;一个重要的认识是受双轴应力作用的弯曲波速度变化在低频区主要与c₁₄₄和c₁₅₅两个三阶弹性模量有关,而且此认识与径向偏振的平面横波一致.在缺乏足够实验条件情况下,对VTI情况,以c₁₄₄,c₁₅₅和c₁₂₃三个独立的量进行测量,然后可暂不考虑三阶弹性模量的各向异性,建立简化的应力反演公式.反之,如果已知地层的地应力信息,由简化的声弹公式可以反演三阶弹性模量c₁₄₄,c₁₅₅和c₁₂₃.     

基于准地转垂直运动方程对Q矢量的进一步探讨

Q矢量是诊断大气大尺度垂直运动的重要物理量.在f平面近似下,从原始运动方程组出发可以推导出Q矢量以及相关的半地转Q矢量、非地转Q矢量和湿Q矢量等表达式,但这种传统的推导方法没有考虑大气准地转运动的属性.本文利用数学物理变换,直接从准地转ω方程出发推导出Q矢量及其相关表达式,从而克服了传统方法在大气准地转运动属性表达方面...     

三维大地电磁自适应L1范数正则化反演

常规三维大地电磁反演的正则项为L₂范数,它以电阻率空间分布函数处处光滑为模型期望,弱化了算法对电性突变界面的分辨能力.本文实现了正则项为L₁范数的三维大地电磁反演算法,让模型空间梯度向量更有机会取得稀疏解,在充分正则的迭代下能够有效突出模型真实电性界面.为避免L₁范数零点不可导带来的求解困难,使用迭代重加权最小二乘法把原问题转换为一系列L₂正则子问题迭代求解.每个子问题的极小方法使用改进型拟牛顿法,其下降方向既能保证正则项海塞矩阵的精确性,又能允许反演过程随迭代灵活更新正则因子.使用比值法或分段衰减法自适应更新正则因子以避免迭代早期陷入奇异解,从而提升反演收敛的稳定性并降低初始模型依赖度.合成的无噪数据反演表明L₁正则算法的模型恢复效果优于L₂正则;不同噪声水平的合成数据反演表明本文的算法具有稳健性;实测数据反演对比表明在合理的正则因子调整策略下,L₁正则反演结果的模型分辨率优于L₂正则.另外,不同初始模型的反演...     

武汉和阿德莱德大气太阴N2潮及M2潮的对比分析

利用阿德莱德（35°S,138°E）和武汉（30.6°N,114.5°E）的流星雷达观测数据首次给出了80~100 km高度上周期为12.66太阳时的大气太阴半日潮汐（N₂潮）的季节、高度、年度变化及其与周期为12.42太阳时的大气太阴半日潮汐（M₂潮）的对比分析.分析结果表明：武汉和阿德莱德的N₂潮和M₂潮均有明显的季节、高度和年度变化.N₂潮与M₂潮的幅度比值大于其引力势之比0.191,在某些年份的不同季节和高度上,N₂潮的幅度甚至大于M₂潮的幅度.大多幅度之比接近或超过N₂潮和M₂潮引力势之比的2倍.中低热层的大气太阴N₂潮汐值得关注.     

跨赤道惯性急流的多平衡态

发展了Anderson和Moore跨赤道惯性急流理论, 应用非线性等值浅水模式, 根据与急流相联系的有旋大尺度气流(大气)或洋流(海洋)特征定出新的普适函数, 从而求出跨赤道惯性急流区中沿流线的位势涡度和能量守恒式, 由于最后的控制方程是非线性的, 因此可以有条件的存在跨赤道的多平衡态的惯性急流解. 通过急流区外不同的大尺度气流或洋流的特征分析, 讨论了跨赤道惯性急流区中可能存在平衡态的条件, 并给出多平衡态的解例.     

基于卫星测高交叉点的海洋表面地转流速度

在流体静力平衡状态下,海洋Coriolis力和压力梯度平衡就形成地转流,世界上大多数海流都近似为地转流.本文利用卫星测高交叉点方法计算海洋表面地转流速度,分析了利用测高交叉点计算地转流速度的不确定性,上升和下降弧段的海面倾斜在分辨率50 km上可以达到10^-7量级,才可能获得优于10 cm/s的地转流速度.在低纬度或者纬度接近卫星轨道倾角的地区,由交叉点方法计算的地转流速度精度低于中纬度地区.以中国台湾东部黑潮为试验区,利用最新的中国台湾周边海域大地水准面模型参考场计算高精度的大地水准面高,利用TOPEX/Poseidon和Jason-1的GDR数据(2002~2005年)计算海面高,然后计算交叉点的动力高,确定交叉点的地转流速度,结果与中国台湾NCOR(National Center for Ocean Research)的流速基本一致.     

海洋运动的尺度分析

本文利用尺度分析理论详尽地论述了旋转海洋中各类运动的特性及其相应的控制方程组,指出现实海洋中可能存在五种大尺度运动,并讨论了存在于这种运动中流场和质量场的适应过程,针对不同类型的运动分别分析了升降流的性质以及层化和摩擦对升降流的影响,估计了地球旋转对海水可压缩性的影响以及非静力平衡在各类海洋运动中的作用。     

动态PPP定位中周跳自动探测与处理策略

周跳的自动探测与有效处理是动态PPP定位中关键和棘手的问题.本文提出利用L₁、L₂、L₄、L₅历元间差分,并联合doppler积分综合探测与定位周跳.对准确判定发生周跳的卫星,不进行周跳的固定和修复,而是将其作为新升起的卫星重新估计模糊度.这一处理策略避免了周跳修复不准对整个定位精度的影响,当保证不少于4颗卫星连续跟踪时,PPP动态定位的精度是稳定的.实验结果表明：在相对比较恶劣的动态环境下,动态PPP定位精度可以达到0.1~0.2 m水平.     

THE ACTIVITY OF WESTERN LISHAN FAULT SINCE THE LATE PLEISTOCENE

Along the northern piedmont of Mt. Lishan, the characteristics and locations of the active normal Lishan fault in west of Huaqing Pool provide important evidences for determining the seismotectonic environment, seismic stability evaluation of engineering in the eastern Weihe Basin. After reviewing the results from high-density resistivity method, seismic profile data, geological drillhole section and trenching in west of the Huaqing Pool, it is found that the strike of western normal Lishan Fault changes from EW direction at the eastern part to the direction of N60°W, and the fault consists of two branches, dipping NE with a high dip angle of~75°. The artificial shallow seismic profile data reveals that the attitude of strata near Lishan Fault mainly dips to south, which is presumed to be related to the southward tilt movement of Mt. Lishan since the Cenozoic. The section of geological drillhole reveals that since the late middle Pleistocene, the displacement of the paleo-soil layer S₂ is about 10m. And the maximum displacement of western Lishan Fault recorded in the paleo-soil layer S₁ reaches 7.8m since the late Pleistocene. In addition, evidences from trench profile show that the western Lishan Fault was active at least 3 times since Malan loess deposition with ¹⁴ C dating age(32 170±530)Cal a BP. The multiple activities of the Lishan Fault result in a total displacement about 3.0m in the Malan loess layer L₁. The latest activity of the western Lishan Fault produced a displacement of about 0.9m in the early Holocene loess layer L₀((8 630±20)Cal a BP)and caused obvious tensile cracks in the Holocene dark leoss layer S₀((4 390±20)Cal a BP). Briefly, we have obtained a vertical movement rate of about 0.11~0.19mm/a since the Holocene((8 630±20)Cal a BP)in the western extension of the Lishan Fault, the recurrence interval of earthquakes on the fault is about(10.7±0.5)ka, and the co-seismic surface rupture in a single event is inferred to be about 0.9m.     

Palaeomagnetism of the loess/palaeosol sequence in Viatovo (NE Bulgaria) in the Danube basin

The results of a palaeomagnetic investigation of a 27 m thick loess/palaeosol sequence in Viatovo (NE Bulgaria) are presented in this paper. The sequence consists of topsoil S₀, seven loess horizons (L₁–L₇) and six interbedded palaeosols (S₁–S₆) overlying a red clay (terra rossa) complex. Magnetic viscosity experiments, IRM acquisition, AMS analysis and NRM stepwise alternating and thermal demagnetisation experiments of pilot samples were implemented for precise determination of the characteristic remanence and construction of a reliable magnetostratigraphical scheme. Analysis of IRM acquisition curves using the expectation – maximization algorithm of Heslop et al. [Heslop, D., Dekkers, M., Kruiver, P., van Oorschot, H., 2002. Analysis of isothermal remanent magnetization acquisition curves using the expectation – maximization algorithm. Geophys. J. Int., 148, 58–64] suggests that the best fitting is obtained by three coercivity components. Component 1 corresponds to SD maghemite/magnetite, while component 2 is probably related to the presence of oxidised detrital magnetites. The third component shows varying coercivities depending on the degree of pedogenic alteration of the samples and probably reflects the presence of detrital magnetite grains oxidised at different degree.

The relevance of the Viatovo section as a key representative sequence for the loess cover in the Danube basin is confirmed by the presence of geomagnetic polarity changes in the lower part of the sequence. The youngest one recorded in the seventh loess unit L₇ can be identified as corresponding to the Matuyama/Brunhes palaeomagnetic polarity transition. Two normal magnetozones were found in the red clay complex, probably corresponding to the Jaramillo and Olduvai subchronozones of the Matuyama chron.     

用宏观方法测定震源深度的量版

根据常用公式 h=Δ_i/10^{（I₀-I_i）/s}-1^1/2, 式中h为震源深度,Δ_i为烈度为I_i的等震綫半径,I₀为震中烈度,S为一系数;取（1）式的对数得 logh=logΔ_i-1/2log[10^{（I₀-I_i）/s}-1], 按（2）式可以作成一量版,以同时測定h和s。利用这个量版測定了19个中国地震的s和h,結合文献[4]的資料,指出中国东部的s系数比西部的偏低;且当深度加大时,s系数加大。采用文献[1,5]的資料測定了61个地震的s系数,結果表明s的数值随深度的增加而加大,占与低速层的关系并不明显。     

On utility of triangular C-grid type discretization for numerical modeling of large-scale ocean flows

Ocean circulation models based on triangular C-grid discretization are frequently employed to simulate coastal ocean dynamics on unstructured meshes. It is shown that on time and space scales dominated by slow geostrophic dynamics, this discretization tends to exhibit checkerboard noise in the field of horizontal velocity divergence and vertical velocity, respectively. The noise is linked to the geometry of triangular C-grid and is amplified in regimes that are close to geostrophic balance through the particular structure of the Coriolis operator. It can be partly suppressed in some cases but remains a problem in a general case and makes the triangular C-grid a suboptimal choice for large-scale ocean modeling.     

Can biomass burning produce a globally significant carbon-isotope excursion in the sedimentary record?

Negative carbon-isotope excursions have been comprehensively studied in the stratigraphic record but the discussion of causal mechanisms has largely overlooked the potential role of biomass burning. The carbon-isotopic ratios (δ¹³C) of vegetation, soil organic matter and peat are significantly lower than atmospheric carbon dioxide (CO₂), and thereby provide a source of low ¹³C CO₂ when combusted. In this study, the potential role of biomass burning to generate negative carbon isotope excursions associated with greenhouse climates is modeled. Results indicate that major peat combustion sustained for 1000 yr increases atmospheric CO₂ from 2.5× present atmospheric levels (PAL) to 4.6× PAL, and yields a pronounced negative δ¹³C excursion in the atmosphere ( 2.4‰), vegetation ( 2.4‰) and the surface ocean ( 1.2‰), but not for the deep ocean ( 0.9‰). Release of CO₂ initiates a short-term warming of the atmosphere (up to 14.4 °C, with a duration of 1628 yr), which is consistent with the magnitude and length of an observed Toarcian excursion event. These results indicate that peat combustion is a plausible mechanism for driving negative δ¹³C excursions in the rock record, even during times of elevated pCO₂.     

与东北冷涡相伴的高空急流诱发平流层重力波的数值模拟研究

使用中尺度数值模式WRF-ARW,针对2010年6月发生在中国东北地区一例伴随对流层高空西风急流（位于~9 km高度）演变过程出现的平流层重力波活动特征开展了数值模拟. 事件发生期间,对流层区域环流处在一个东北冷涡系统的控制之下. 模拟结果再现了该东北冷涡的发展和维持过程,以及与之相伴的高空急流的特征. 模拟结果揭示出在急流区域上空的平流层中存在显著重力波活动现象. 分析结果显示,重力波活动与急流存在紧密联系,在水平方向上,重力波呈显著的二维结构,出现在急流出口区上部并逆背景流向西传播. 功率谱分析结果表明盛行波动具有~700 km水平尺度、9~12 h时间尺度以及4~5 km垂直波长. 由于急流的存在,造成其与平流层中下部之间存在显著的水平风速垂直切变,与切变相伴的耗散使得上传的重力波动量通量数值随着高度升高而递减. 同时,在18~20 km高度间出现的西风-东风转换带极大地抑制了波动在垂直方向的传播,形成显著动量通量沉积效应. 估算结果表明,在11~20 km高度之间,这种效应的整体作用相当于对该层背景流施加强度为0.86 m·s^-1·day^-1的动力阻曳.     

基于二级近似离散复镜像法的低频电磁格林函数计算

本文讨论了利用二级近似离散复镜像法实现低频格林函数快速、精确的计算.通过数值分析,给出了在低频电磁场计算中该方法近似参数的选取原则.(1)二级近似的积分区间里谱格林函数采样个数N为十倍的近似多项式的个数M;(2)总的积分区间L2大小约为40/r,r为收发距离;(3)两个积分区间的分界值L1为总的积分区间大小L2与第二个...