Two‐dimensional controlled‐source electromagnetic interferometry by multidimensional deconvolution: spatial sampling aspects

We use numerically modelled data sets to investigate the sensitivity of electromagnetic interferometry by multidimensional deconvolution to spatial receiver sampling. Interferometry by multidimensional deconvolution retrieves the reflection response below the receivers after decomposition of the fields into upward and downward decaying fields and deconvolving the upward decaying field by the downward decaying field. Thereby the medium above the receiver level is replaced with a homogeneous half‐space, the sources are redatumed to the receiver level and the direct field is removed. Consequently, in a marine setting the retrieved reflection response is independent of any effect of the water layer and the air above. A drawback of interferometry by multidimensional deconvolution is a possibly unstable matrix inversion, which is necessary to retrieve the reflection response. Additionally, in order to correctly separate the upward and the downward decaying fields, the electromagnetic fields need to be sampled properly. We show that the largest possible receiver spacing depends on two parameters: the vertical distance between the source and the receivers and the length of the source. The receiver spacing should not exceed the larger of these two parameters. Besides these two parameters, the presence of inhomogeneities close to the receivers may also require a dense receiver sampling. We show that by using the synthetic aperture concept, an elongated source can be created from conventionally acquired data in order to overcome these strict sampling criteria. Finally, we show that interferometry may work under real‐world conditions with random noise and receiver orientation and positioning errors.     

Glacial paleoceanography off the mouth of the Mekong River, southern South China Sea, during the last 500 ka

We have analyzed core MD01-2392, ∼ 360 km east of the Mekong River mouth in the South China Sea (SCS). Over the past 500 ka, planktonic foraminiferal oxygen-isotopic values are consistently lighter than northern SCS and open-ocean records by up to 0.5‰, indicating the influence of run-off from the Mekong River during both glacial and interglacial periods. Carbonate content is higher during interglacials; sedimentation rates were higher during glacials. Increased sedimentation rates since 30 ka imply increased run-off during the last glacial maximum and Holocene Period. Contrary to general experience, in which it is classed as a warm species for temperature estimates, the thermocline-dwelling species Pulleniatina obliquiloculata increased its numbers during glacial periods. This implies an estuarine circulation and even brackish-water caps during glacial periods, reinforcing the sense of strong run-off. In an overall decline of warm water, the thermocline shoaled stepwise, with rapid rises across the glacial terminations. We infer that the southern SCS was opened to an influx of Indian Ocean waters through southern passages at those times of rising sea levels.     

Three-dimensional Q structure in Jiashi earthquake region of Xinjiang

3-D S-waveQ structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, and a rough configuration of deep crustal faults in the earthquake region is presented. First, amplitude spectra of S-waves are extracted from 450 carefully-chosen earthquake records, called observed amplitude spectra. Then, after instrumental and site effect correction, theoretical amplitude spectra are made to fit observed amplitude spectra with nonlinear damped least-squares method to get the observed travel time overQ, provided that earthquake sources conform to Brune’s disk dislocation model. Finally, by 3-D ray tracing method, theoretical travel time overQ is made to fit observed travel time overQ with nonlinear damped least-squares method. In the course of fitting, the velocity model, which is obtained by 3-D travel time tomography, remains unchanged, while onlyQ model is modified. When fitting came to the given accuracy, the ultimateQ model is obtained. The result shows that an NE-trending lowQ zone exists at the depths of 10–18 km, and an NW-trending lowQ zone exists at the depths of 12–18 km. These roughly coincide with the NE-trending and the NW-trending low velocity zones revealed by other scientists. The difference is that the lowQ zones have a wider range than the low velocity zones. Foundation item: Joint Seismological Science Foundation of China (957-07-414) and State Key Basic Research Development and Programming Project (95-13-02-02). Contribution No. RCEG200105, Research Center of Exploration Geophysics, China Seismological Bureau.     

小入射角下海浪后向散射系数不对称性和各向异性分析验证

基于北京遥测技术研究所机载波谱仪实测数据以及ECMWF提供的风场数据,分析了不同风速条件下小入射角海浪平均后向散射系数不对称性与各向异性随入射角的变化特性,并对导致这一现象的原因进行了初步分析。结果表明:小入射角下,海浪平均后向散射系数随入射角增大呈现递减趋势并在18°时达到某一稳定值,之后基本维持不变;海浪后向散射系数不对称性与各向异性均随入射角增大呈现递增趋势,同样在18°增大至某一稳定值,之后基本维持不变;低风速条件下,海浪不对称性和各向异性与风速有关,且风速越大不对称性与各向异性越明显。     

内蒙古呼和浩特-包头地区地震目录的完整性分析

以内蒙古中部的呼和浩特-包头地区(40°～42°N,108°～114°E)作为研究区域,收集整理2001年1月1日至2010年4月30日ML≥1.0地震作为研究对象,按地震样本数目滑动扫描并计算最小完整性震级Mc,绘制Mc的时序变化曲线.分析认为Mc在时间上的变化主要反映台网监测能力和地震活动规律,但台站布局、地震信噪比和人为因素也会引起Mc的短期或不连续性的变化.同时在Rydelek等提出的2个假设条件下,对呼包地区的地震目录完整性进行了RS检验,其结果与上述扫描计算结果一致,也与根据台网监测能力所得到的震级范围吻合,证明本文的计算结果能够比较真实的反映2001年以来呼包地区地震目录在各时段内的完整性.     

金活动态提取剂提取-电感耦合等离子体质谱法测定深穿透地球化学样品中的金

金的地球化学勘查基于金的准确测定,地球化学样品中金含量通常处于ng/g水平,需先进行分离富集,再采用电感耦合等离子体质谱法(ICP-MS)或石墨炉原子吸收光谱法(GFAAS)进行测定。当前,隐伏矿床勘查是地球化学探测技术的发展前沿,金活动态提取技术是寻找隐伏金矿的有效手段之一。相比于全量分析,金的活动态含量更低,需要解决选择性提取、高效预富集与准确测定等一系列难题。本文采用柠檬酸铵与土壤中黏土矿物及次生矿物作用促使吸附和可交换组分的金进入提取液,以硫脲和硫代硫酸钠络合金使活动态金向提取液中扩散,达到选择性提取的目的,建立了提取液中金的预富集及ICP-MS测定方法。实验确定的分析条件为:采用5g/L柠檬酸铵-2g/L硫脲-5g/L硫代硫酸钠为提取剂,提取时间24h,在酸性硫脲介质下用活性炭富集金,金吸附率可达89.6%～109.2%,灰化解吸温度为650～700℃。本方法检出限为0.05ng/g,相对标准偏差(RSD)为9.4%～10.2%,加标回收率为91.2%～93.4%。与已报道的硫酸铁-硫脲-硫代硫酸钠溶液提取再GFAAS测定的方法相比,本方法具有检出限低、测试线性范围宽、测试速度快的优势;应用于森林覆盖区黑龙江东安金矿区地球化学探测试验,金活动态异常与隐伏金矿位置一致。     

辽西北票(金-羊)盆地北票组层序地层及沉积体系特征

为揭示北票(金-羊)盆地内北票组层序地层特征及沉积体系分布规律,应用层序地层学及沉积学理论和方法,综合利用野外露头、剖面等资料,建立了北票(金-羊)盆地北票组的层序地层格架,确定研究区的沉积相类型及分布规律。研究结果表明,北票(金-羊)盆地北票组可划分为2个三级层序,其内部可进一步划分为7个准层序组。在层序内识别出冲积扇、辫状河三角洲、湖泊3个沉积相。明确了沉积相的分布规律,向上形成由冲积扇到辫状河三角洲,再到湖相,最后是辫状河三角洲的沉积体系;平面上呈近北西—南东向展布的沉积体系,整体构成2个水进、水退的沉积体系。古地貌影响了底部冲积扇相沉积物的沉积特征。气候不仅影响了沉积物类型,且和湖平面变化共同控制了沉积物的演化规律。     

Sea Level Prediction in the Yellow Sea From Satellite Altimetry With a Combined Least Squares-Neural Network Approach

Accessible high-quality observation datasets and proper modeling process are critically required to accurately predict sea level rise in coastal areas. This study focuses on developing and validating a combined least squares-neural network approach applicable to the short-term prediction of sea level variations in the Yellow Sea, where the periodic terms and linear trend of sea level change are fitted and extrapolated using the least squares model, while the prediction of the residual terms is performed by several different types of artificial neural networks. The input and output data used are the sea level anomalies (SLA) time series in the Yellow Sea from 1993 to 2016 derived from ERS-1/2, Topex/Poseidon, Jason-1/2, and Envisat satellite altimetry missions. Tests of different neural network architectures and learning algorithms are performed to assess their applicability for predicting the residuals of SLA time series. Different neural networks satisfactorily provide reliable results and the root mean square errors of the predictions from the proposed combined approach are less than 2?cm and correlation coefficients between the observed and predicted SLA are up to 0.87. Results prove the reliability of the combined least squares-neural network approach on the short-term prediction of sea level variability close to the coast.     

长江口特大洪水浑浊带悬沙粒度分布特征

2020年长江流域发生了历史第二大洪水,大通站洪峰流量达到84 500 m3/s。本文基于2020年7月长江口特大洪水期间最大浑浊带多站位的水沙观测数据,重点分析了悬沙粒度组分的时空分布特征,并与常态水文条件下的粒度数据进行对比。结果表明：（1）最大浑浊带悬沙垂向平均中值粒径为10.4μm,变化范围为6～27μm,以黏性细颗粒泥沙为主;其中核心区南槽、北槽及北港的中值粒径分别为8.4μm、7.6μm和8.5μm,过渡区分别为7.2μm、16.4μm和14.5μm。（2）悬沙中值粒径垂向分布受不同组分影响,核心区底层中值粒径为8.8～9.6μm;底层黏土含量在28%～31%之间,粉砂含量在61%～64%之间,中值粒径主要受黏土及粉砂组分影响;过渡区北港和北槽垂向平均砂组分高达19%,南槽砂组分平均仅占5%,中值粒径主要受砂组分影响。（3）对比2013年洪季浑浊带数据,2020年粒径整体增大5.4μm,核心区黏土含量相较2013年减少12.7%,砂增加6.3%;过渡区北槽与北港平均粒径增大10μm。     

2017年九寨沟MS7.0地震对周围地区的静态应力影响

基于九寨沟M_S7.0地震的破裂模型及均匀弹性半空间模型,本文计算了该地震在周围主要活动断层上产生的库仑应力变化、在周围地区产生的应力场和位移场和同震库仑应力变化对余震的触发.结果表明：（1）九寨沟地震造成虎牙断裂中段库仑破裂应力有较大增加,已经超过0.01 MPa的阈值,虎牙断裂北段、塔藏断裂中段和岷江断裂北段北部的库仑破裂应力有较大降低,因此尤其要注意虎牙断裂中段的危险性.（2）水平面应力场在该地震震中东西两侧增加（拉张）,张应力起主要作用.在震中南北两侧降低（压缩）,压应力起主要作用.从水平主压和主张应力方向来看,均呈现出条形磁铁的磁场形态.从剖面上的应力场来看,在上盘的面膨胀区域内,大部分点的主张应力方向与地表是垂直的,在其他区域内,主张应力和主压应力均以震中为中心,向外呈辐射状.（3）从地表水平位移场来看,震中东西两侧物质朝震中位置汇聚,南北两侧物质向外流出,在震中处的最大水平位移量达43 mm.从地表垂直位移场来看,震中南北两侧出现明显的隆升,隆升最大值达56.8 mm.震中东西两侧出现明显的沉降,沉降最大值达74.5 mm.从剖面的位移场来看,九寨沟地震为左旋走滑地震,且有一定的正断成分.由分析可以推测该断层破裂在大致22~26 km的深度上就截止了.并推测下盘物质运动的动力来自震源北东东方向（四川块体）深度在6~30 km的下盘下层物质,上盘物质运动的动力来自震源北西西方向（巴颜喀拉块体）深度在0~6 km的上盘上层物质.（4）通过计算不同深度上主震对余震的触发作用可知,主震后的最大余震受到了主震的触发作用,多数其他余震也受到主震的触发作用.主震的发生促使了库仑应力增加地区余震的发生,抑制了一部分库仑应力减少地区余震的发生.