Constraint inversion Algorithm of Lidar Equationfor Deriving Aerosol Optical Property

1.IntroductionLidarisincreasinglyusedinmeasuringtheopticalpropertiesofaerosolsandcloudssincethe1960's.Ithasbeenshownthatthedeterminationoftheextinctioncoefficientasafunctionofheightisfeasible(Klett,1980;Fernald,1984;LuDarenetal.,1977;Spinhirneetal.,...     

大连初冬一次辐射平流雾天气过程分析

利用大连机场地面观测资料、Micaps系统下常规资料、探空资料和NCEP/NCAR全球再分析资料,从天气形势和背景、探空资料分析和物理量诊断方面,对2009年11月30日-12月2日发生在大连地区持续性大雾天气过程做了详细分析。结果表明,本次持续性雾过程属于辐射平流雾,是在稳定的大尺度天气背景下形成的。探空资料表明,大雾发生过程中,边界层内出现一层逆温和多层逆温;边界层内近地层的逆温和充沛的水汽条件对雾的形成和长时间的维持起着重要的作用。热力结构分析表明,温度日较差大表明地面辐射冷却对本次大雾过程具有明显的作用;低层持续的弱暖平流输入,有利于近地层逆温的建立和维持。动力场结构分析表明,在中低层,大雾发生前期和维持时期,存在弱的辐合上升运动;在大雾消散期,存在明显的辐散下沉运动。水汽条件分析表明,增湿和冷却使此次大雾过程中水汽达到饱和状态产生凝结,在大雾过程的前期,存在弱的水汽辐合;在大雾消散期,存在水汽辐散。     

天津机场地区冻雨天气分析

利用气象常规资料，对天津机场地区1979～2002年出现的3次冻雨天气进行对比分析，找出了冻雨发生的天气形势特点和大气垂直结构特征：地面受低压控制，高空处在槽前，有明显的增温现象，槽前西南暖湿气流强盛，存在低空急流；850hPa到地面存在逆温层；700～500hPa之间有冰晶层，冰晶层内的温度为-10～-14℃：850hPa附近为暖层，暖层底部的高度在1．0～2．0km左右，暖层的厚度大约为1．0～2．0km，暖层内的温度大约为0～2℃；从地面到1．0～2．0km的高度存在冷层，冷层内温度大约为0～-2℃。以2002年12月13～14日的冻雨天气为代表个例，利用客观分析产品进行诊断分析，总结出冻雨天气预报的着眼点，为冻雨天气预报提供参考。     

利用伴随法反演风应力系数

基于无结构网格有限体积法的正模型和伴随数据同化模型,建立了风应力系数反演表达式,将风应力系数的确定变成依赖于区域内观测值的数值迭代过程,实现了风应力系数确定的自动化,并成功进行了概化河口一系列风应力系数反演的伴随同化实验,反演得到的结果精度较高,实验结果表明：1）伴随法反演风应力系数是合理可行的;2）充分的观测数据能够提高反演的精度;3）风速在10 m/s或15 m/s的情况下,对于风应力系数反演的结果影响不大。     

正则化滤波在反演地表光谱反照率中的应用

利用有限的卫星观测数据,通过反演地物二向性反射分布函数(BRDF)模型提取复杂地表的真实信息,是一个不适定的地球科学反演问题。为了克服离散不适定性带来的困难,在分析引起模型解不稳定的原因和数据误差传递机制的基础上,使用构建滤波函数的方法实现BRDF模型的正则化约束反演。实验结果表明,正则化滤波算法与MODIS AMBRALS算法精度相当,完全适用于具有充足观测和稀疏观测情况下的地物参数反演。     

CAP方法反演震源机制的误差分析:以胶东半岛两次显著中等地震为例

利用区域波形数据使用CAP方法反演中强地震的震源机制正逐渐得到广泛应用.本文以胶东半岛近期发生的两次显著中等地震为例,讨论了使用CAP方法反演震源机制时的误差估计,展示了反演结果的不确定性分析过程.2013年11月23日和2014年1月7日在山东莱州和乳山分别发生了M4.6和M4.3级中等地震,两次事件均造成了较大影响.我们基于CAP方法,使用自助抽样(bootstrap)技术多次重复反演过程,得到大样本量的震源机制解数据;基于这些数据,使用粒子群算法和聚类分析技术给出了优化解,估计了震源机制解的误差范围,并利用震源机制解的P、T轴给出了震源球上的概率密度分布.     

多条人工地震测深剖面资料联合反演首都圈三维地壳结构

为了对多条二维剖面资料进行统一分析和解释,本文考虑了各剖面交点处的约束条件,通过建立统一的偏导矩阵,将各剖面资料联立求解.利用该方法对中国地震局地球物理勘探中心在首都圈内6条人工地震宽角反射/折射剖面资料进行了重新处理和解释,得到该地区的地壳三维速度结构和构造特征.在此基础上,结合其他地球物理研究成果,确定了壳内存在的4处低速层（体）的空间范围及断裂的走向.     

2018年9月8日墨江地震及周边地区构造应力场特征分析

2018年9月8日,云南省墨江县发生MS5.9地震并伴随一系列余震,探究该地震周围的应力场对于理解该地震的发生机制和后续地震的发展趋势具有着重要的参考意义.本研究收集了震源及其邻区中前人研究和Global CMT所给出的震源机制资料,对该地区进行了构造应力场反演,并同时利用反演得到的应力张量模拟墨江地区的震源机制解表现.结果表明:(1)在应力轴整体分布上,自西向东σ1轴(压轴)从NNE-SSW向逐渐转向NNW-SSE向,σ3轴(张轴)从WWN-EES向逐渐转向WWS-EEN向,张轴呈弧形分布,压轴呈放射状分布.(2)在应力轴倾伏角上,研究区域内的压应力轴和张应力轴倾伏角都比较小,即两轴均接近水平.(3)R值分布大体是在东南部相较于西北部大,结合当地地质背景分析得到,物质逃逸自西北向东南呈逐渐变缓的趋势.(4)利用反演得到的应力张量和应力状态计算墨江地震震源区的相对剪应力和相对正应力大小.由此推测,墨江地震恰好发生在相对剪切应力值和相对正应力正值最大的节面上.从而可以确定墨江地震的发震节面的基本参数:走向216.32°,倾角86.91°,滑动角0.27°,相对剪应力值0.9,相对正应力值0.3.本研究为此次墨江地震的发震背景和地震动力学研究提供了基础性资料.     

深渊着陆器坐底位置的精确测量和反演计算

深渊着陆器在下潜过程中,由于受到海流的影响,其坐底位置会偏离投放点,偏离量和偏离方向很难估计,而声学水下定位方法在深渊海区不太适用.通过在海面进行气枪十字放炮作业,利用炮点准确位置和时间信息,以及海底地震仪(Ocean Bottom Seismometer,OBS)记录的直达水波到时,对深渊着陆器的坐底位置进行精确测量可以提供有效定位数据.该方法对采集得到的数据进行预处理,使用蒙托卡洛方法结合最小二乘法,对万泉号着陆器实际坐底位置进行了反演校正,得到三个潜次的位置校正偏移量分别为211m、178m、861m,偏移方向各不相同,校正精度为±20m,比以往的定位精度大为提高.此外,该方法测量确定了三个海底标志物的精确位置,为未来万米无人或载人潜器提供了可靠的水下参考坐标.     

全球地磁感应测深数据三维反演

全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构.