各向异性标准化方差计算重磁源边界

在重磁源边界定位方法中,传统的梯度方法易受干扰的影响使计算的边界混乱,而且在弱异常处由于叠加异常的影响很难识别场源边界.本文首先利用坐标旋转构造了各向异性高斯函数,提出了各向异性标准化方差计算重磁源边界的方法.理论分析与模型实验详细阐明了该方法的数学含义,并通过干扰分析验证了方法的稳定性与有效性,结合中扬子地区航磁异常处理实例,表明通过计算各向异性标准化方差,可以有效地确定重磁源边界,尤其对微弱异常能够实现较好的边界定位,计算的边界分辨率高且信息丰富,有利于资料的综合解释.     

基于Slepian方法和地面重力观测确定时变重力场模型: 以2011—2013年华北地区数据为例

时变重力场是研究地球系统内部物质运动和时空演化过程的有效途径.目前广泛使用的GRACE时变重力场模型受限于其空间分辨率（约400 km）,难以探测较小空间尺度的重力变化.本文首次尝试利用Slepian局部谱分析方法和多期地面重力观测确定更高空间分辨率的时变重力场模型.Slepian方法通过构建研究区域内的正交基函数,将信号能量集中在研究区域内部,是构建球面局部重力场模型的理想方法.本文根据Slepian方法的特点给出了区域重力场建模及参数优化的步骤,以我国华北地区为例,基于2011—2013多期地面观测确定了区域时变重力场模型,并与同区域由Slepian方法和GRACE卫星数据确定的重力变化进行了对比分析.结果表明：（1）贝叶斯信息量准则可作为确定Slepian展开最佳截断数的有效手段;（2）基于研究区域内现有重复测点数据,能够恢复120阶时变重力场,空间分辨率（半波长）约150 km;（3）2011—2013年间研究区域内GRACE估计结果与120阶地面结果在时空分布的显著趋势上存在较好的对应,证明了本文利用Slepian方法和地面观测所得时变重力场模型的可靠性.本文研究结果可为区域重力场建模提供新的参考,也可为华北地区水资源变化监测、构造活动分析以及地震风险性评估等研究提供高分辨率的时变重力场模型支撑.

     

基于重磁震资料的南海新生代盆地分布综合研究

作为西太平洋最大的边缘海,南海分布有30多个新生代沉积盆地,其蕴含着丰富的油气资源.但由于资料的限制,南海存在不同区域盆地研究程度不同,不同区域盆地面积差别较大,部分盆地只是坳陷而没有达到盆地的级别以及盆地外围可能存在凹陷等问题.南海新生代盆地分布问题制约了其油气分布规律、储量等基础地质问题的研究.本文以地震剖面数据为约束,以重力资料为主、辅以磁力资料,研究了南海新生代盆地分布及构造区划.通过提取新生代盆地及其构造单元引起的重力异常,结合地震剖面等资料反演了新生界底界面深度及新生界厚度.在充分调研已有盆地和构造单元划分方案的基础上,根据南海的地质及地球物理特征,确定了盆地及构造单元划分标准.以新生界厚度为基础并结合重、磁、震、地质等资料,进行地质-重磁震联合解释,将南海原有的36个盆地重新划分为24个盆地,盆地总面积扩大了约15万km2.研究表明,南海新生代盆地沉积层厚度在1.5~16 km之间,有6个北东东/北东向沉积坳陷带、2个近南北向沉积坳陷带以及1个三角沉积坳陷区;盆地展布方向主要为北东和北东东向,其次为北西和近南北向,呈现"南三北三"的分布特征.     

东北地区典型湖沼沉积物溶解态有机质组成特征及来源解析：以库里泡为例

溶解态有机质(DOM)是湖泊沉积物的重要组分,解析DOM的组成和来源对于深入理解湖泊有机质生物地球化学循环及控制水体富营养化具有重要意义。 本研究于2021年5月和8月采集了大庆市库里泡表层及柱状沉积物、泡内水生物(浮游藻类、挺水植物和沉水植物)、岸边土壤、陆生C3和C4植物、禽畜粪便以及城镇污水等样品,分析了样品DOM的稳定同位素(δ¹³C 和 δ¹⁵N)组成及三维荧光光谱特征,并利用IsoSource软件计算了不同来源样品对沉积物DOM的贡献率。结果显示:(1)库里泡内表层沉积物DOM的δ¹³C和δ¹⁵N组成存在季节性差异,5月δ¹³C和δ¹⁵N均值分别为-25.54‰和9.02‰,8月分别为-26.81‰和8.40‰。(2)库里泡内柱状沉积物DOM的δ¹³C和δ¹⁵N组成在垂直方向存在差异,表层(0~3 cm)δ¹³C和δ¹⁵N均值分别为-26.58‰和9.04‰,深层(3~30 cm)均值分别为-25.40‰和10.61‰。(3)表层沉积物DOM的三维荧光组分存在季节性差异。5月以类腐殖质荧光组分为主,占比为87.89%,HIX和BIX分别为6.27和0.67;8月蛋白类荧光组分占比为49.58%,HIX和BIX指数均值分别为1.72和0.87。(4)5月表层沉积物DOM外源输入占比为61%,以土壤(21.40%)和城镇污水(18.08%)为主;而8月内源贡献占比稍高(55.10%),且以挺水植物为主(48.68%)。(5)柱状沉积物(0~30 cm)不同深度DOM来源组成相近,主要为城镇污水、沉水植物/藻和挺水植物,贡献率均值分别为42.13%、25.07%和18.53%。整体上,库里泡沉积物DOM来源主要受到流域内人类活动及气候特征的影响,本文研究结果有利于加深对我国东北地区湖泊沉积物有机质迁移转化及累积规律的认识。     

宜兴太湖防护林对农业面源污染的减源增汇作用浅析

为探讨太湖沿岸防护林宜兴段对农业面源污染物的控制效果,通过在宜兴市太湖防护林采集土壤样品,测定其中营养元素含量,结合周边农业生产施肥概况,初步分析了该区域防护林对农业面源污染的减源—增汇效应.结果表明,2009年10月份防护林地上部分固定的C、N、P和K分别为2439.48、18.34、2.48和13.42 g/m2;2011年5月份固定量分别为3114.12、23.42、3.16和17.13 g/m2;而2013年8月份固定量分别为4722.46、35.53、4.81和25.99 g/m2.可见C、N、P和K的固定量是逐年增加的,相比2009年,分别增加了93.6%、93.7%、94.0%和93.7%.结果说明宜兴段太湖沿岸防护林对农业面源污染的控制效果较显著,有助于太湖水污染的治理.     

基于先验信息约束的重磁三维交互反演建模技术——以铜陵矿集区为例

三维地质建模是实现深部矿产勘查突破的重要途径,其通过控矿地质体的三维建模,直观刻画控矿要素之间的空间、成因和演化关系,帮助理解成矿系统,开展深部找矿预测.当前三维建模主要采用地质资料构建,在缺少钻孔等已知资料的情况下,难于构建出可靠的三维地质模型.为了克服常规三维地质建模方法可信度低、精度差的缺点,本文将重磁交互反演技术引入到三维地质建模中,提出了基于先验信息约束,通过二度半剖面交互反演、三维物性反演联合修正的三维地质模型方法.采用该方法建立了铜陵矿集区的三维地质模型,并开展了深部找矿预测,取得以下主要认识:1)基于先验信息约束的重磁交互反演建模技术能大幅提高三维地质模型的可信度,是实现地下地质体"透明化"的重要途径;2)铜陵地区重要控矿地层(C-P-T2)主要分布于向斜区和火山岩覆盖区,深部岩浆岩条件优越,深部找矿工作应从传统的背斜隆起区转向向斜区和火山岩覆盖区以及部分凹陷区,这将极大拓展铜陵地区找矿空间;3)铜陵地区传统的五大矿田的侵入岩均具有复合岩体特征,特别是与辉石闪长岩共生的矿产以铁铜矿为主,而与金矿密切相关的侵入岩则多为中酸性岩体,这一推断如果成立,将对铜陵地区找矿突破产生重大影响;4)铜陵地区的推覆构造十分发育,木镇凹陷区的早古生代基底和南陵盆地北西缘的戴公山背斜就是一套区域性推覆构造的产物,三维地质模型还揭示铜陵地区侵入岩具有深、中、浅三重结构,这为铜陵隆起的推覆成因提供了证据,也指示了在铜陵之外再寻找类似铜陵矿集区的良好前景.     

羌塘盆地中央隆起带的重磁场证据及其构造意义

对于羌塘盆地是否存在横贯东西的中央隆起带,目前学术界仍有分歧.本文提供的最新高精度航空重、磁资料证实存在呈东西向贯通羌塘盆地的中央隆起带,并对该带的构造特征进行了精细刻画.隆起带受南北两侧深大断裂控制,其空间跨度（宽度）由西向东逐渐收敛,并被一组近南北向的隐伏断裂系切割、左滑错动.重、磁场资料还显示中央隆起带在双湖东、西两侧存在明显差异：西段基底大规模隆起,基岩深度一般在3~5 km以内,明显浅于南北羌塘坳陷7~15 km的基底埋深;东段基底隆起幅度明显降低,主要表现为潜伏的低隆起,其中双湖—雅曲段基底埋深5~7 km,雅曲—岗尼段基底埋深7~9 km;即中央隆起带基底自西向东"台阶状"降低,隆起的幅度和分布范围受到近南北向断裂控制.构造分层表明,与南羌塘地块相比,北羌塘地块的基底隆起幅度小、稳定性更好.南北羌塘基底地球物理属性的显著差异说明羌塘盆地并不存在统一的前寒武系变质基底,中央隆起带的形成应该与古特提斯洋关闭时形成的混杂岩带有关.     

Ecological–economic modeling as a tool for watershed management: A case study of Lake Qionghai watershed, China

This paper presents an ecological–economic model for a lake and its watershed systems. We describe the linkage between the watershed system and the lake aquatic ecosystem and the modeling process. The lake–watershed system was divided into six subsystems: social system, economic system, terrestrial ecosystem, lake water system, pollutant system, and lake aquatic ecosystem. The model equations were constructed based on five main assumptions. The Lake Qionghai watershed in southwestern China, which is undergoing rapid eutrophication, was used as a case study. The targeted goals for total phosphorus (TP) and chlorophyll a (Chl a) concentrations in the lake in 2015 are 0.025 and 10.0 mg m⁻³, respectively. We present two scenarios from 2004 to 2015 based on the ecological–economic model. In both scenarios, the TP and Chl a concentrations in the lake are predicted to increase under the effects of watershed pressures and the targeted goals cannot be met. The application of techniques to reduce pollutants loading and the corresponding pollutants reductions are reflected again in the constructed model. The model predicts that TP and Chl a concentrations will decrease to 0.024 and 7.71 mg m⁻³, respectively, which meet the targeted thresholds. The model results provide directions for local government management of watersheds and lake aquatic ecosystem restoration.     

青藏高原冈底斯带及邻区重磁三维反演及岩浆岩特征研究

岩浆岩在青藏高原的大陆动力学研究中有着重要的作用,它既是构造演化的记录,又是重要构造-岩浆-成矿带的指示.本文主要基于冈底斯带及邻区的地面重力和航磁数据,首先进行地质-地球物理先验信息约束下的重磁2.5维交互式反演,再将2.5维反演结果作为参考模型加入到三维反演计算中,得到地下三维密度和磁化率结构.结合岩浆岩密度、磁化率统计资料和岩浆岩地球化学成果,推断研究区基性岩、Ⅰ型花岗岩和S型花岗岩的三维分布图,得到如下结论：S型花岗岩主要分布在冈底斯东带和冈底斯弧背断隆带以北;北冈底斯的西部无明显的岩浆活动,而在其南侧和北侧,发现大量的隐伏基性岩和零散分布的Ⅰ型花岗岩;中生代Ⅰ型花岗岩在南冈底斯和冈底斯弧背断隆带广泛分布,且到新生代才出现大量的S型花岗岩.上述结果为中生代班公湖-怒江洋壳和新特提斯洋壳的双向剪刀式俯冲模式的观点提供了重要佐证,并认为班公湖-怒江洋壳在北冈底斯西部约84°E-88°E的范围内先后存在向北和向南俯冲的可能,北向羌塘地体下俯冲,南向冈底斯地体下俯冲.

     

Preliminary sensitivity study of eruption source parameters for operational volcanic ash cloud transport and dispersion models — A case study of the August 1992 eruption of the Crater Peak vent, Mount Spurr, Alaska

Ash clouds are one of the major hazards that result from volcanic eruptions. Once an eruption is reported, volcanic ash transport and dispersion (VATD) models are used to forecast the location of the ash cloud. These models require source parameters to describe the ash column for initialization. These parameters include: eruption cloud height and vertical distribution, particle size distribution, and start and end time of the eruption. Further, if downwind concentrations are needed, the eruption mass rate and/or volume of ash need to be known. Upon notification of an eruption, few constraints are typically available on many of these source parameters. Recently, scientists have defined classes of eruption types, each with a set of pre-defined eruption source parameters (ESP). We analyze the August 18, 1992 eruption of the Crater Peak vent at Mount Spurr, Alaska, which is the example case for the Medium Silicic eruption type. We have evaluated the sensitivity of two of the ESP – the grain size distribution (GSD) and the vertical distribution of ash – on the modeled ash cloud. HYSPLIT and Puff VATD models are used to simulate the ash clouds from the different sets of source parameters. We use satellite data, processed through the reverse absorption method, as reference for computing statistics that describe the modeled-to-observed comparison. With the grain size distribution, the three options chosen, (1) an estimated distribution based on past eruption studies, (2) a distribution with finer particles and (3) the National Oceanic and Atmospheric Administration HYSPLIT GSD, have little effect on the modeled ash cloud. For the initial vertical distribution, both linear (uniform concentration throughout the vertical column) and umbrella shapes were chosen. For HYSPLIT, the defined umbrella distribution (no ash below the umbrella), apparently underestimates the lower altitude portions of the ash cloud and as a result has a worse agreement with the satellite detected ash cloud compared to that with the linear vertical distribution for this particular eruption. The Puff model, with a Poisson function to represent the umbrella cloud, gave similar results as for a linear distribution, both having reasonable agreement with the satellite detected cloud. Further sensitivity studies of this eruption, as well as studies using the other source parameters, are needed.     

3D joint inversion of magnetotelluric and airborne tipper data: a case study from the Morrison porphyry Cu–Au–Mo deposit,British Columbia,Canada

Z‐axis tipper electromagnetic and broadband magnetotelluric data were used to determine three‐dimensional electrical resistivity models of the Morrison porphyry Cu–Au–Mo deposit in British Columbia. Z‐axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Ground‐based magnetotelluric surveys can achieve a greater exploration depth than Z‐axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z‐axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z‐axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z‐axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non‐uniform and relatively sparse spatial sampling. The joint Z‐axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z‐axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and ≥0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z‐axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z‐axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z‐axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.     

A new approach of dynamic monitoring of 5‐day snow cover extent and snow depth based on MODIS and AMSR‐E data from Northern Xinjiang region

Taking the Northern Xinjiang region as an example, we develop a snow depth model by using the Advanced Microwave Scanning Radiometer‐Earth Observing System (AMSR‐E) horizontal and vertical polarization brightness temperature difference data of 18 and 36 GHz bands and in situ snow depth measurements from 20 climatic stations during the snow seasons November–March) of 2002–2005. This article proposes a method to produce new 5‐day snow cover and snow depth images, using Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) daily snow cover products and AMSR‐E snow water equivalent and daily brightness temperature products. The results indicate that (1) the brightness temperature difference (Tb_18h–Tb_36h) provides the most accurate and precise prediction of snow depth; (2) the snow, land and overall classification accuracies of the new images are separately 89.2%, 77.7% and 87.2% and are much better than those of AMSR‐E or MODIS products (in all weather conditions) alone; (3) the snow classification accuracy increases as snow depth increases; and (4) snow accuracies for different land cover types vary as 88%, 92.3%, 79.7% and 80.1% for cropland, grassland, shrub, and urban and built‐up, respectively. We conclude that the new 5‐day snow cover–snow depth images can provide both accurate cloud‐free snow cover extent and the snow depth dynamics, which would lay a scientific basis for water management and prevention of snow‐related disasters in this dry and cold pastoral area. After validations of the algorithms over other regions with different snow and climate conditions, this method would also be used for monitoring snow cover and snow depth elsewhere in the world. Copyright © 2011 John Wiley & Sons, Ltd.     

长江中下游典型浅水湖泊表层水体微塑料时空分布特征——以湖北保安湖为例

以湖北保安湖为研究对象,分别于2021年4、7、10月和2022年1月采集表层水样,测定水体中微塑料丰度、粒径、形状、颜色和类型,探究浅水型湖泊表层水体中微塑料的时空分布特征及其影响因素。结果表明,保安湖表层水体微塑料的年平均丰度为(16.20±2.23)items/L。微塑料丰度呈现明显的季节性差异,其中夏季(采样时间7月)平均丰度最低,为(1.40±0.09)items/L。在所有微塑料颗粒中,粒径0.064～1 mm占比最大,为82.57%;黑色和无色微塑料占主导地位,占比分别为36.16%和21.31%;纤维状微塑料分布最广泛,占比达40.01%;聚乙烯(PE)和低密度聚乙烯(LDPE)是最主要的微塑料类型,两者之和占比达46.05%,其次为聚苯乙烯(PS,占比17.2%)、聚对苯二甲酸乙二醇酯(PET,占比8.33%)和聚乙酸乙酯(PVAC,占比8.31%)。统计分析表明,微塑料丰度与湖水水质无显著相关关系。该研究揭示了长江中下游地区典型浅水湖泊微塑料分布现状,为评估类似湖泊微塑料潜在污染风险提供了科学依据。     

南方红壤区长期水土流失综合治理显著减少河流输沙——以鄱阳湖流域赣江上游平江为例

长时间序列水沙数据分析有助于科学评价流域尺度水土保持减流减沙效应,为科学开展区域水土保持成效评估提供实现路径.本文选取南方红壤区水土保持工程最为集中、持续时间最长的典型流域——平江流域,基于1975-2014年的逐日降水量、流量、含沙量数据,综合应用预置白M-K检验法、Theil-Sen趋势度估计法、Pettitt检验...     

岩石磁组构对剩磁稳定性的制约探讨:以印支地块中生代碎屑岩和拉萨林周盆地设兴组红层为例

岩石磁组构因能提供磁性矿物晶体形状、排列方式等赋存信息而被广泛应用于判别岩石剩磁是否受到了后期构造应力的显著影响; 但常规岩石磁化率各向异性(AMS)是否能够准确限定岩石剩磁的稳定性, 目前尚无深入探讨.本文以印支地块Nakhon Thai盆地中生代Nam Phong、Phu Kradung和Phra Wihan组三套碎屑岩样品及拉萨地块林周盆地设兴组红层样品为例, 通过岩石磁化率组构和剩磁组构的对比分析发现, 尽管Nam Phong组绝大多数样品和Phu Kradung组全部样品的AMS组构显示其具有铅笔状至强劈理过渡型构造变形组构特征, 但高场等温剩磁各向异性(hf-AIR)显示其高矫顽力赤铁矿所携带的特征剩磁组构仍具有典型沉积组构特征, 表明其以赤铁矿为主的载磁矿物未遭受后期构造应力的显著影响, 仍然能够准确记录岩石形成时期的古地球磁场方向.另一方面发现有且仅有剩磁组构才是判别碎屑沉积岩特征剩磁是否遭受了后期构造应力影响的充分必要条件.也就是说, 如果岩石剩磁组构(如以赤铁矿为主要载磁矿物的岩石hf-AIR组构)指示其原始沉积组构已被构造组构显著叠加或取代, 则必然说明该岩石剩磁方向已受到构造应力的作用而发生了显著偏转, 因而不能直接用于构造演化和古地理重建等块体运动学研究.