Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden's update and 1D derivatives

Dense profile-oriented resistivity data allows for 2D and 3D inversions. However, huge amounts of data make it practically impossible to do full 2D or 3D inversions on a routine basis. Therefore, a number of approximations have been suggested over the years to speed up computations. We suggest using a combination of Broyden's update on the Jacobian matrix with derivatives calculated using a 1D formulation on a parameterized 2D model of locally 1D layered models. The approximations bring down the effective number of 2D forward responses to a minimum, which again gives us the ability to invert very large sections. Broyden's update is not as useful with a parameterized problem as is the case with a smooth minimum structure problem that has been the usual application. 1D derivatives, however, seem to be very effective when initiating a full 2D solution with Broyden's update. We compare the different methods using two different kinds of data on two synthetic models and on two field examples. The most effective and reliable optimization combines 1D derivatives with a full 2D solution and Broyden's update. When using Broyden's update the Jacobian matrix needs to be reset every once in a while. We do this whenever the difference in data residual from the previous iteration is less than 5%. This combined inversion method reduces the computation time approximately a factor of 3 without losing model resolution.     

广东省东部地区的S波速度结构

利用广东省数字遥测台网记录到的2000年11月16日发生在所罗门群岛的两次Ms7．7和Ms7．8地震及2000年11月18日新爱尔兰地区的Ms7．3地震面波波形数据，测定了在广东地区汕头和韶关之间路径上的相速度频散曲线并反演得到了对应的地壳横波速度结构。结果表明该地区的地壳厚度为30km，上地壳横波速度变化较大，中下地壳相对均匀，没有发现地壳内部的低速层。     

Current states of stress in the northern Andes as indicated by focal mechanisms of earthquakes

Systematic inversion of double couple focal mechanisms of shallow earthquakes in the northern Andes reveals relatively homogeneous patterns of crustal stress in three main regions. The first region, presently under the influence of the Caribbean plate, includes the northern segment of the Eastern Cordillera of Colombia and the western flank of the Central Cordillera (north of 4°N). It is characterized by WNW–ESE compression of dominantly reverse type that deflects to NW–SE in the Merida Andes of Venezuela, where it becomes mainly strike–slip in type. A major bend of the Eastern thrust front of the Eastern Cordillera, near its junction with the Merida Andes, coincides with a local deflection of the stress regime (SW–NE compression), suggesting local accommodation of the thrust belt to a rigid indenter in this area. The second region includes the SW Pacific coast of Colombia and Ecuador, currently under the influence of the Nazca plate. In this area, approximately E–W compression is mainly reverse in type. It deflects to WSW–ENE in the northern Andes south of 4°N, where it is accommodated by right-lateral displacement of the Romeral fault complex and the Eastern front of the northern Andes. The third, and most complex, region is the area of the triple junction between the South American, Nazca and Caribbean plates. It reveals two major stress regimes, both mainly strike–slip in type. The first regime involves SW–NE compression related to the interaction between the Nazca and Caribbean plates and the Panama micro-plate, typically accommodated in an E–W left-lateral shear zone. The second regime involves NW–SE compression, mainly related to the interaction between the Caribbean plate and the North Andes block which induces left-lateral displacement on the Uramita and Romeral faults north of 4°N.Deep seismicity (about 150–170 km) concentrates in the Bucaramanga nest and Cauca Valley areas. The inversion reveals a rather homogeneous attitude of the minimum stress axis, which dips towards the E. This extension is consistent with the present plunge of the Nazca and Caribbean slabs, suggesting that a broken slab may be torn under gravitational stresses in the Bucaramanga nest. This model is compatible with current blocking of the subduction in the western northern Andes, inhibiting the eastward displacement of slabs, which are forced to break and sink in to the asthenosphere under their own weight.     

Structure de socle,sismostratigraphie et héritage structural au cours du rifting au niveau de la marge d'Ifni/Tan-Tan (Maroc sud-occidental)

Seismic reflection profiles from the Ifni/Tan-Tan Atlantic margin of southern Morocco, interpreted in the light of well data and field geology from the Western Anti-Atlas, allowed us to establish the seismostratigraphic framework of the syn-rift series and to reveal (i) a compressional structural style in the pre-Triassic basement similar to that established in the adjacent outcropping onshore basement but with an opposed western vergence, (ii) the importance of inherited anterior structures in the formation of Triassic-Liassic rift structures and (iii) an east–west propagation of these rift structures. To cite this article: N. AbouAli et al., C. R. Geoscience 337 (2005).     

刘公岛海洋牧场底层海水溶解氧浓度的变化特征

由于缺乏长期观测资料,前人对山东半岛邻近海域海水溶解氧的时间变化和空间分布特征的研究较少。本文基于威海刘公岛海洋牧场于2016年7月20日至2017年3月14日期间,利用生态环境实时在线观测系统获得的底层海水的温度、盐度、水深、溶解氧数据,分析了该牧场海水溶解氧浓度的时间变化特征及其影响因素,并探讨了低氧灾害发生的可能性。结果表明在观测期间,该牧场海水溶解氧浓度以季节变化为主,冬季最大、夏季最小,其中2月份平均值最高,约为10.86mg/L,8月份平均值最低,约为5.91mg/L。同时海水溶解氧浓度也存在显著的小时变化和日变化,且变化幅度于8月份最大、3月份最小。影响海水溶解氧浓度变化的主要因素是海水温度,溶解氧浓度随着温度的季节性变化而变化。夏季,水体分层会使溶解氧浓度发生大幅度的降低,大风过程对于溶解氧浓度也有一定的影响,通过打破夏季的季节性温跃层使水体发生垂向混合从而为海底提供氧气,但大风过程之后的几天会出现溶解氧浓度降低的现象。本次研究发现刘公岛海洋牧场在观测期间不存在低氧现象。     

GEOPHYSICAL EXPLORATION

正20141074Bao Xijie(Research Institute of Exploration and Development,Daqing Oilfield Company,PetroChina,Daqing 163712,China)Gather Optimal Processing and Application Effect of Prestack AVA Instantaneous Inversion     

The impact of northern Indian Ocean rivers on the Bay of Bengal using NEMO global ocean model

The effect of river runoff over the northern Indian Ocean(NIO) especially over the Bay of Bengal(Bo B) has been studied using global Nucleus for European Modelling of the Ocean(NEMO). Two sensitivity experiments, with and without river runoff are conducted and the influence of river runoff on the Indian Ocean hydrography,stratification and circulation features are studied. It is found that due to river runoff surface salinity over the northern Bo B decreases by more than 5 and the East India Coastal Current strengthens by 2 cm/s during post monsoon season. The fresh river water reaches up to 15°N in the Bo B and is the main cause for low salinity there.Sea surface temperature in the northwestern Bo B increases by more than 0.2℃ due to the river runoff in summer monsoon while surface cooling upto 0.2℃ is seen in north-west part of Bo B in winter season. The seasonal mixed layer depth in the region is found to be dependent on river runoff. The effect of vertical shear and Brunt Vaisala frequency on stratification is also examined. The ocean water becomes highly stratified up to 3 035 m due to the river runoff. It is found that the energy required for mixing is high in the northern and coastal Bo B.     

湖北省冻雨天气特征及其判定方法研究

利用常规地面气象站多时次气象要素和天气现象观测资料,以及武汉、宜昌和恩施的探空资料,对1998—2014年湖北地区冻雨天气的主要环流背景进行了分析,归纳了冻雨发生发展过程的气象要素和大气层结特征,给出了基于正负能量面积的冻雨发生定量化判断方法。结果表明:乌山阻高型和贝湖阻高型是导致湖北省出现冻雨天气的两种主要天气形势;2005年和2月是冻雨主要发生的年份和月份;共存在3种典型的气温层结演变特征,其中"暖雨-冻雨-固态"的气温层结演变最易导致持续时间长、影响范围大的冻雨天气;基于正负能量面积的冻雨判定方法为:3层层结时,正能量面积(A_(SP))小于80℃·hPa,负能量面积(A_(SN))小于400℃·hPa,且正负能量关系为(5.71A_(SP)-257.14)≤A_(SN)≤(6.25A_(SP)+200);2层层结时,ASP为350~650℃·hPa,A_(SN)为200~400℃·hPa,且正负能量关系为(A_(SP)-350)≤A_(SN)≤(ASP-100)。     

季节性冰封热融浅湖水温原位观测及其分层特征

为探究季节性冰封浅湖热力学特征,于2010年10月至2013年7月对高原腹地一典型热融湖塘冰层生消、水/冰温及气象条件开展原位观测,分析了水温分布时间变化、温跃层以及冰生消对水温结构的影响。结果发现:冰面升华显著,贯穿整个冰期;水温日变化、季节变化和垂直结构受气温、大气辐射、风速、冰生消和湖底沉积层热贡献影响显著;在"无冰期-结冰前-冰生长期-冰融化期-融化后-无冰期"年循环过程中水温垂直结构分别呈现出"分层-翻转-逆温分层-逆温与正温共存-翻转-分层"的循环过程。分层期水温结构仅由上部混合层和温跃层构成,且偶因强风搅动而全湖翻转混合。可见,相比大中型湖泊,季节性冰封浅湖热力学结构差异显著。     

Study on the Global Carbon Assimilation System Based on Multisource Remote Sensing Data

The Paris agreement signed in April, 2016 aims to balance global anthropogenic carbon emissions and terrestrial carbon sinks by the middle of the 21^st century. To fulfill this goal, it is necessary to calculate carbon fluxes of different regions reliably. The global carbon assimilation system is an effective technique for achieving this goal. The Ministry of Science and Technology of China supports the project entitled as study on the global carbon assimilation system based on multisource remote sensing data through the national key research and development programs for global change and adaptation during the thirteen-five period. This project will develop synergic inversion techniques for retrieving key parameters of biological and atmospheric cycles and for assimilating multisource remote sensing and ground based data. Then, the high resolution global carbon assimilation system coupled with an ecological model will be constructed. This system is able to assimilate jointly multisource observation data and to optimize key model parameters, photosynthesis and respiration carbon fluxes of global terrestrial ecosystems, and anthropogenic carbon emission fluxes of key regions. This system will be used to study quantitatively the spatial and temporal patterns of carbon fluxes of global terrestrial ecosystems and anthropogenic carbon emission fluxes of key regions and to identify the mechanisms driving the global terrestrial carbon sinks and sources. The outputs of this study will be helpful for the fulfillment of the key research and development programs for global change and adaptation and provide valuable data and technical support for the decision-making in China.