沙漠下垫面地面温度及能量收支的BATS模式参数化改进

利用陆面过程模式BATS,引入地表发射率及两种大气发射率参数化方案,同时引入不同的地表粗糙度参数化方案,对比各种参数化方案对沙漠下垫面地面温度及能量收支的模拟状况。结果表明：采用Van Bavel等发展的地表发射率及Chung等发展的大气发射率方案可以明显改进地面温度及向上长波辐射的模拟,Chung等方案在夜间与正午的模拟效果更好,减小了1 ℃左右的地面温度模拟偏差,减小了10 W·m^-2左右的向上长波辐射模拟偏差。晴天地面温度及向上长波辐射的模拟结果优于阴天。利用Zhang等发展的裸土粗糙度参数化方案也会提高模式对地表感热通量模拟的准确性。     

一种利用历史DLG数据辅助优化数字正射影像镶嵌线生成的方法

为了提升大面积数字正射影像图生产的效率,介绍了一种利用历史DLG数据辅助优化数字正射影像镶嵌线生成的方法,通过对DLG数据处理,提取出区域具有地物的面层,生成避让区域,在此基础上再考虑镶嵌线尽可能穿过原始DLG中较少的面及色差较小的区域,使用贪婪蛇模型算法实现了优化镶嵌线的生成。通过生产试验,证明了该方法可部分减少镶嵌线编辑的工作量。     

基于移动-多面函数的高程异常拟合方法

将多面函数拟合法获得局部区域的最优结果作为移动二次曲面拟合时选择搜索半径的依据。然后,将获得的搜索半径作为最优拟合距离,运用移动二次曲面法实现了高程异常拟合,即移动-多面函数拟合方法。以某区域获取的高程异常作为实例,比较了3种方法的拟合精度。结果显示,移动二次曲面法、多面函数和移动多面函数的中误差分别为1.499、0.870和0.304 m,可见,移动-多面函数法比另外两种方法的拟合精度都高。     

利用接收函数反演青藏高原西部地壳S波速度结构

相对于宽阔的腹地,青藏高原西部南北向宽度仅约600km,却记录了印度和欧亚板块汇聚的深部过程及其响应.本文用22台宽频带流动地震台站在西缘构建了一条南北向探测剖面(~80°E,TW-80试验).利用接收函数反演剖面下方S波速度结构,综合西部已有的宽频带探测结果,分析认为:印度板块向北俯冲可能已到达班公湖—怒江缝合带附近,俯冲过程中下地壳发生榴辉岩化;喀拉昆仑断裂带、班公湖—怒江缝合带、阿尔金断裂带均为切穿地壳的深断裂,莫霍面发生错断;喀拉昆仑断裂带和龙木错断裂带之间的中上地壳没有发现连续的S波低速体,说明可能缺乏解耦层,支持青藏高原西部地壳为整体缩短增厚模式.     

利用重力/GPS联合观测数据计算地壳垂向构造应力的新方法

本文提出一种基于重力/GPS联合观测数据计算垂向构造应力的新方法.计算步骤如下:(1)通过重力/GPS联合观测数据计算布格重力异常;(2)依据布格重力异常数据推算莫霍面深度;(3)依据GPS观测数据,通过均衡理论计算均衡面深度;(4)依据莫霍面与均衡面之间剩余物质(壳幔物质密度差)所承受的附加浮力,计算地壳承载的垂向构造应力.本文利用上述构造应力新算法,计算了巴颜喀拉块体东边界及周边地区垂向构造应力分布,发现龙泉山断裂带以东地区垂向构造应力基本为零,龙泉山断裂带与龙门山断裂带之间地区垂向构造应力为正值,巴颜喀拉地块东部垂向构造应力为负值.鲜水河断裂带东南段周边蓄积了-40~-50 MPa的垂向构造应力,且梯度变化剧烈;松潘高原蓄积的垂向构造应力大约为-10~-20 MPa,相对较小.     

基于微动方法研究五大连池火山区尾山火山锥浅层剪切波速度结构

为了研究五大连池火山区尾山火山锥浅层三维波速结构特征,在尾山火山锥附近区域布设了无线地震检波器密集台阵,记录连续地震背景噪声波形数据.基于微动方法(拓展空间自相关方法)提取了台站间2~5 Hz频率范围的Rayleigh面波相速度频散曲线.利用面波层析成像方法反演获得2~5 Hz Rayleigh面波二维相速度图像,基于每一个网格节点的频散曲线,进一步反演获得了尾山火山锥附近区域地表至700m深度的三维剪切波速度结构.成像结果显示:在0~150m较浅深度,靠近尾山火山锥区域显示为相对高速异常,远离火山锥区域则显示为相对低速异常.而至150~700m较深深度,波速异常特征与浅部相反,靠近尾山火山锥的区域显示为相对低速异常,远离火山锥的区域显示为相对高速异常.在远离尾山火山锥区域,浅层的相对低速异常可能与松散沉积层有关,深部的高速异常则反映了结晶变质岩的影响.在靠近尾山火山锥区域,浅部的相对高速异常应该反映了出露地表的玄武岩,而深部的相对低速异常则可能反映了火山通道周围广泛发育的破碎裂隙结构及其火山喷发后孔隙流体填充的影响.     

基于交叉梯度约束的地震初至纵波与瑞雷面波联合反演

本文提出了一种初至纵波(P波)与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演(Windowed-Amplitude Waveform Inversion,w-AWI)可获得近地表横波(S波)速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构.     

青藏高原东南缘基于背景噪声的Rayleigh面波方位各向异性研究

青藏高原东南缘地区是现今地壳形变和地震活动最强烈的地区之一,也是研究青藏高原现今变形机制和构造演化规律的重要区域.本研究使用云南区域地震台网的55个宽频带地震台站连续地震背景噪声数据,采用双台站互相关方法获得Rayleigh(瑞利)面波经验格林函数,提取相速度频散曲线,反演得到云南地区周期5~34s范围内方位各向异性分布图像.反演结果揭示:短周期(5~12s)Rayleigh面波快波优势方向与区域断裂走向有很好的一致性,快波方向随着断裂走向的变化而变化.周期16~26s快波优势方向与反映上地壳特性的5~12s图像总体图像相似,但细节略有不同.其中,滇中块体内易门断裂和滇中块体内东侧的普渡河断裂附近,各向异性快波方向从NS向NW方向旋转;易门断裂以西呈NW向.这反映了青藏高原物质东流和川滇块体受到青藏块体的南东向挤压作用.周期30~34s范围的各向异性,滇缅泰块体和印支块体,快波优势方向为NS和NNW向;而在滇中块体内部,各向异性快波方向呈顺时针旋转变化,可能与青藏高原物质向东逃逸有关.本文还开展了与体波各向异性的对比分析,通过与近震S波分裂、Pms转换波分裂和远震SKS、PKS和SKKS(以后简称为XKS)分裂的对比研究,发现随着周期的增大,得到的快波优势方向与XKS剪切波快波偏振方向趋向一致,与地壳快剪切波偏振方向呈一定夹角.本研究认为,青藏高原东南缘地区壳幔各向异性具有不同的特征和形成机制.     

利用嫦娥二号微波辐射计数据的全月亮温制图

中国嫦娥绕月卫星在世界上首次搭载了微波辐射计,用于获取不同深度月壤的微波辐射亮温,而这正是月壤厚度研究和月壤表层参数信息反演的基础。基于嫦娥二号的微波辐射计2C级数据,建立4个频率通道(3.0、7.8、19.35、37.0GHz)在不同纬度下的亮温日变化模型,并通过该模型对亮温数据进行时角校正,得到了正午时刻的全月微波亮温图。结合CCD影像和DEM数据比较的结果表明,处理后的亮温数据很好地反映了月球的地形特征及月表的反射率信息。基于高频的亮温数据对月球亮温异常区域进行了分析,发现其都是哥白尼纪的新鲜撞击坑。     

风云三号C星微波湿温探测仪的定标和验证

风云三号C星(FY-3C)已经于2013年9月23日发射升空,其上装载的微波湿温探测仪(MWHTS)已于9月30日开机正常工作.MWHTS具有对大气温度和湿度垂直分布进行同步探测的能力.MWHTS为跨轨扫描式微波辐射计,在89~191GHz毫米波段内设置了十五个探测通道,其中包括118.75GHz氧气吸收线附近的8个大气温度探测通道,183.31GHz水汽吸收线附近的5个大气湿度探测通道,以及89GHz和150GHz两个窗区通道.设置在118.75GHz的一组毫米波探测通道是国际上业务卫星首次使用的大气探测通道,这组通道和183.31GHz通道对大气进行联合探测,将获得更加精细的大气温湿度垂直分布数据,为数值预报和气候研究提供丰富信息.为保证MWHTS观测资料的定量应用,对仪器性能和定标精度进行了在轨测试.利用MWHTS在轨正常工作后的三个月数据,对仪器在轨定标的基础数据:冷空和黑体计数值,黑体和仪器温度进行监测分析和质量检验,经过质量检验的在轨定标基础数据,结合发射前真空试验得到的非线性订正项在轨定标生成MWHTS观测亮温数据.评估MWHTS在轨辐射定标结果的精度和偏差特性使用了三种方法:1通过场地定标试验获取大气温湿廓线和地面温度等大气参数信息,结合微波逐线正演辐射传输模式MonoRTM(Monochromatic Radiative Transfer Model)模拟MWHTS的上行微波辐射亮温,与MWHTS实际观测结果进行对比分析;2两个通道特性一致的同类星载被动微波载荷同时观测同一目标,观测亮温的差异主要取决于两个载荷的定标系统偏差.选取美国SNPP上搭载的微波探测仪器ATMS作为MWHTS的参考载荷,基于SNO(simultaneous nadir overpass)技术,对两个仪器的观测亮温进行交叉比对,观测亮温时空匹配及均匀性检验的条件为:观测时间差异小于20min,观测像元中心距离小于3km,观测角度在星下点附近差异小于5°,观测像元周围3×3像元内的亮温标准差小于1K;3基于美国国家环境预测中心的全球数据同化系统GDAS(Global Data Assimilation System)数据,利用快速辐射传输模式CRTM(Community Radiative Transfer Model)对MWHTS各通道亮温进行正演模拟,模拟结果(O)和仪器实际观测的亮温(B)之间的差异记为"O-B",对偏差值"O-B"进行统计特征分析.仪器中心频率的变化、正演模式模拟精度和模式输入廓线自身的误差都会对"O-B"产生影响.但是对于首次使用的探测频点而言(如118.75GHz通道),由于国际上没有同类载荷可以进行交叉比对,借助于正演辐射传输模式计算得到"O-B"偏差的分析结果可以在一定程度上反映仪器整体定标情况.外场地定标试验结果显示除通道14外,其他14个通道的亮温差都在1.3K以内;与同类载荷ATMS的在轨观测进行直接交叉比对表明通道14与ATMS的亮温偏差最大,但中心频点一致的5个水汽探测通道的标准差都小于1K;将MWHTS观测结果和正演辐射传输模式模拟结果即"O-B"进行偏差分析显示,靠近118.75GHz吸收线中心的通道2—6"O-B"标准差小于0.5K,其他通道"O-B"标准差和ATMS相应通道的结果相当;MWHTS观测和模拟偏差随角度变化的研究表明通道1,7~13和15观测结果对角度有一定依赖性.