稀疏反演求解基追踪降噪问题的地震谱分解方法

谱分解描述为一线性反演问题,由于该问题欠定性,需使用稀疏反演算法,然后将该问题转化为基追踪降噪问题(BPDN),引入谱投影梯度(SPGL1)稀疏反演算法提高反演谱分解(ISD)分辨率,并进一步研究其潜在优势。结果表明:基于基追踪降噪问题的反演谱分解方法 (ISD-BPDN)在烃类检测中有更高的时频分辨率,分层更准确,精细地指示了烃类的存在。     

南海文昌地区内波振幅反演研究

The field experiment is conducted from April 16,2005 to July 20,2005 at Wenchang area east of Hainan Island(19°35'N,112°E) of China.Internal wave packets are observed frequently with thermistor chains during the experiment.Meanwhile,internal waves are also detected from a synthetic aperture radar(SAR) image on June 19,2005 and several other moderate-resolution imaging spectroradiometer(MODIS) images near a mooring position.The distance between the positive and negative peaks induced by the internal wave can be obtained from satellite images.Combined with remote sensing images and in situ data,a new method to inverse the amplitude of the internal wave is proposed based on a corrected nonlinear Schr?dinger(NLS) equation.Two relationships are given between the peak-to-peak distance and the characteristic wavelength of the internal wave for different nonlinear and dispersion coefficients.Based on the satellite images,the amplitude inversion of the internal waves are carried out with the NLS equation as well as the Kd V equation.The calculated amplitudes of the NLS equation are close to the observation amplitude which promise the NLS equation a reliable method.     

Parameterization of sheared entrainment in a well-developed CBL. Part I: Evaluation of the scheme through large-eddy simulations

The entrainment flux ratio A _e and the inversion layer (IL) thickness are two key parameters in a mixed layer model. A _e is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat flux. The IL is the layer between the mixed layer and the free atmosphere. In this study, a parameterization of A _e is derived from the TKE budget in the firstorder model for a well-developed CBL under the condition of linearly sheared geostrophic velocity with a zero value at the surface. It is also appropriate for a CBL under the condition of geostrophic velocity remaining constant with height. LESs are conducted under the above two conditions to determine the coefficients in the parameterization scheme. Results suggest that about 43% of the shear-produced TKE in the IL is available for entrainment, while the shear-produced TKE in the mixed layer and surface layer have little effect on entrainment. Based on this scheme, a new scale of convective turbulence velocity is proposed and applied to parameterize the IL thickness. The LES outputs for the CBLs under the condition of linearly sheared geostrophic velocity with a non-zero surface value are used to verify the performance of the parameterization scheme. It is found that the parameterized A _e and IL thickness agree well with the LES outputs.     

地表温度反演的算法综述

地表温度是研究地球表面的重要参数,被称为地表的皮肤温度。反演地表温度对自然灾害监测、城市热岛等有着重要意义,利用热红外遥感技术是反演地表温度的主要方式。本文阐述了目前主要用于反演地表温度的各算法特性、优缺点和适用情况,对各种算法进行对比分析,并指出反演地表温度存在的问题及解决问题的方向,并对反演地表温度发展趋势进行了展望。     

雾霾天气个例气象条件对比分析

应用常规观测资料、NCEP 1°×1°再分析资料和L波段探空资料从环流形势、扩散条件和边界层特征3个方面对2013年两次雾、霾天气个例进行对比分析,结果表明:500hPa西北气流冷平流、地面弱风场、垂直速度呈弱上升-下沉的垂直分层特点和逆温是两次雾、霾天气出现和维持的共同特征。地面西北风、850hPa弱冷平流、近地层浅薄的接地逆温(100~200m)和湿层与霾天气对应,地面偏东风、850hPa暖平流、925hPa以下深厚的悬浮逆温(400m)和湿层与雾天气对应,霾过程较雾过程逆温强度强,上升运动高度高。消散时雾较霾下沉运动中心高度低,强度弱;霾消散时接地逆温特征变化不大,雾消散时悬浮逆温有底部抬升和大气稳定层结向中性层结转变的变化特征;但均有下沉气流接地、垂直风切变较强和高层低露点干空气下传到地面的特点。     

红外临边探测发展现状

临边探测器垂直分辨率高,辐射信号的变化主要受大气参数影响,且辐射校准后几乎不受背景辐射影响,对痕量气体可探测性强,探测高度约为10—100 km,对臭氧层的大气遥感监测具有特殊意义。红外临边传感器经过40余年的发展过程,在通道数目、信噪比、可探测目标种类、冷冻技术、在轨寿命方面都有明显改进。本文以平流层O_3及其相关的痕量气体为中心,按照传感器搭载的卫星平台及卫星发射时间的先后顺序,比较了国际上已有的LRIR、LIMS、SAMS、ISAMS、CLAES、MIPAS、HIRDLS、TES共8个星载红外临边探测器的特性、探测目标、技术改进、性能提升等;介绍了红外临边探测大气成分的反演原理;总结了红外临边探测器在痕量气体,平流层云与气溶胶,重力波与极涡,非局部热力学平衡效应4个方面的应用成果。最后对红外临边探测的前景趋势提出了思考。     

琉球沟弧盆系莫氏面深度的三维广义反演

本文利用三维广义反演技术对重力资料进行反演，计算出了琉球沟弧盆系地区的莫氏面深度。由于考虑了分辨率和方差之间的最佳折衷问题，反演结果不论是在精度上还是在分辨能力方面都有了明显改进。分析结果也表明，重力反演得到的莫氏面深度同李昭兴等人用地震折射波法得到的莫氏面深度有较大的差别。文中地产生这种差别的可能原因提出了自己的见解。     

Morphology and crustal structure of a small transform fault along the Mid-Atlantic Ridge: The Atlantis Fracture Zone

The Atlantis Fracture Zone (30° N) is one of the smallest transform faults along the Mid-Atlantic Ridge with a spatial offset of 70 km and an age offset of ~ 6 Ma. The morphology of the Atlantis Fracture Zone is typical of that of slow-slipping transforms. The transform valley is 15–20 km wide and 2–4 km deep. The locus of strike-slip deformation is confined to a narrow band a few kilometers wide. Terrain created at the outside corners of the transform is characterized by ridges which curve toward the ridge-transform intersections and depressions which resemble nodal basins. Hooked ridges are not observed on the transform side of the ridge-transform intersections. Results of the three-dimensional inversion of the surface magnetic field over our survey area suggest that accretionary processes are sufficiently organized within 3–4 km of the transform fault to produce lineated magnetic anomalies. The magnetization solution further documents a 15-km, westward relocation of the axis of accretion immediately south of the transform about 0.25 Ma ago. The Atlantis Transform is associated with a band of high mantle Bouguer anomalies, suggesting the presence of high densities in the crust and/or mantle along the transform, or anomalously thin crust beneath the transform. Assuming that all the mantle Bouguer anomalies are due to crustal thickness variations, we calculate that the crust may be 2–3 km thinner than a reference 6-km thickness beneath the transform valley, and 2–3 km thicker beneath the mid-points of the spreading segments which bound the transform. Our results indicate that crustal thinning is not uniform along the strike of the fracture zone. Based on studies of the state of compensation of the transform, we conclude that the depth anomaly associated with the fracture zone valley is not compensated everywhere by thin crust. Instead, the regional relationship between bathymetry and gravity is best explained by compensation with an elastic plate with an effective thickness of ~ 4 km or greater. However, the remaining isostatic anomalies indicate that there are large variations away from this simple model which are likely due to variations in crustal thickness and density near the transform.     

热带西太平洋海温异常对大气影响的数值模拟

用IAP-GCM二层楼式,模拟1月份西太平洋暖池区海温增暖对1-8月份大气环流的影响,初步探讨暖池区海温异常对我国夏季降水的影响过程。指出北半球500hPa环流场对暖地区1月份海温增暖的响应存在着PNA和EAP波列,在冬季较为明显。其能量传播特征表示为南北半球的中纬度地区基本上是向东传播,在热带低纬度地区存在着准2个月的低频振荡。     

Improvements in navigation using SeaBeam crossing errors

Adjustments to satellite constrained navigation are required to match SeaBeam bathymetric data at track crossings due to errors in dead reckoning and inaccuracies in satellite fixes. By shifting one of the SeaBeam swaths involved in a track crossing relative to the other and calculating the sum of the squares of the differences in bathymetry within the area of overlapping coverage, we map a two-dimensional error surface whose minimum corresponds to the best estimate of the correction to navigation required at the crossing point. Estimates of the covariance of this correction are derived from the error surface. We employ the curve fitting technique of Tarantola and Valette (1982) to invert for a smooth correction function to a starting model of the position of the ship as a function of time. This technique incorporates formal errors assigned to dead reckoning, satellite fixes, and the shifts required to match bathymetric swaths at crossing points in a simultaneous inversion for the correction function for all tracks within the study area.In a test of the method in a study area on the southern Mid-Atlantic Ridge, a data set involving two cruises, 30 days of SeaBeam data, and 753 track crossings, we found that crossing SeaBeam swaths can potentially resolve the relative position of the ship on the two tracks to within 30 to 70 m. The inversion procedure yielded a much better constrained navigation function and much improved match of bathymetry. The final model of the navigation fit crossing shifts about as well as satellite data (with respect to their assigned data errors) with the RMS value of the crossing shifts decreasing from 1200 m in the original satellite-constrained navigation to 200 m in the final solution. However, the potential resolution of position using SeaBeam swaths was not fully achieved in the solution because there are systematic bathymetric artifacts in SeaBeam data, multiple local minima in the error surfaces in highly lineated topography, inadequate dead reckoning data, occasional bad satellite fixes, and limitations on the short period corrections allowed in the model.