Identification and genesis of the Late Jurassic volcanic rocks in Longquan area, southern Zhejiang Province
-
摘要:
浙江晚中生代大规模火山活动始于晚侏罗世还是早白垩世，长期以来存在争议。本文对浙南龙泉地区火山岩进行了锆石U-Pb定年，获得一批高精度的年代学数据，其SHRIMP锆石U-Pb和LA-ICP-MS锆石U-Pb年龄集中于163~145 Ma，确认浙江晚中生代火山活动始于晚侏罗世。地球化学特征和Sr-Nd同位素特征表明这套晚侏罗世火山岩具有高硅、高钾、低磷、贫铁镁的特征，属于过铝质高钾钙碱性岩石系列。由于晚中生代太平洋板块的俯冲作用，导致了下部地壳物质（基底变质岩）广泛熔融，形成了本区晚侏罗世火山岩浆活动。为了便于表达、对比，笔者建议新建“晚侏罗世黄茅尖群”地层单位。
Abstract:There has been controversy for a long time that the large-scale volcanic activity of Late Mesozioc in Zhejiang Province began with Late Jurassic or Early Cretaceous? Through zircon U-Pb geochronologic study of volcanic rocks in Longquan area of southern Zhejiang Province, a batch of high precision chronologic data were obtained. These zircon U-Pb ages are between 163 Ma to 145 Ma which were confirmed by SHRIMP and LA-ICP-MS. These data indicate that large-scale volcanic avtivity in Zhejiang Province began with Late Jurassic. According to geochemical characteristics and Sr-Nd isotope analysis, these Late Jurassic volcanic rocks belong to peraluminous and high potassium Ca-alkaline series characterized by high SiO2, K2O but low P2O5, FeOT and MgO. Due to the subduction of the Pacific plate in Late Mesozoic, the lower crust material (basement metamorphic rocks) was extensively melted and the late Jurassic volcanic-magmatic rocks in this area were formed. In addition, the authors propose defining a new rock-stratigraphic unit named" the Late Jurassic Huangmaojian Group"so as to better express and contrast.
-
1. 引言
机载高光谱成像是一种新兴技术，它是一种快速、无损的检测技术，具有光谱分辨率高、波段多和图谱合一的特点。自从1982年美国航天局喷气推进实验室(Jet Propulsion Laboratory，JPL)研制出第一台机载成像光谱仪(AIS-1)并获取遥感飞行数据以来，高光谱遥感的出现为遥感技术带来一场革命。而后包括中国在内的许多国家都研制成功了一系列高光谱遥感设备，开展对地观测和空间目标探测任务，搭载于卫星平台、车载平台(月球车、火星车等)和航空平台。进入21世纪以来，机载高光谱应用发展出现井喷态势，高光谱仪器性能的提高、功能的日益完善和高光谱数据快速、准确、自动化的分析处理，使得机载高光谱日益得到行业重视，广泛应用于地质调查、环境监测、精准农业、海洋观测、土地质量和城市规划等领域(张宗贵等，2000；甘甫平，2004；王润生等，2010)。
作为光谱细分遥感成像技术，光谱成像仪研制朝着高光谱、高空间分辨率和微型化轻量化发展。目前国外商业化的光谱成像仪有3种：光机扫描型、推扫型和凝视型(郑志忠等，2017；郑志忠等, 2020a, b)。以澳大利亚的HyMap为代表的光机扫描型机载成像仪，成为国内外高光谱行业应用的主打仪器。推扫型成像光谱仪有加拿大的CASI，芬兰的SPECIM、挪威的AISA和美国的Headwall成像光谱仪(郑志忠等, 2020b, c)。凝视型成像光谱仪有美国CRI的液晶可调谐成像光谱仪Varispec。国内小型化的机载高光谱仪器依赖进口，价格昂贵，且高性能的高光谱设备进口困难，限制了高光谱技术在国内的普及和推广应用。
近些年来，在可见光-近红外(400~1000 nm)和短波红外(1000~2500 nm)波段范围的高光谱成像仪研制方面，国内取得长足进步，优势单位包括上海技术物理研究所、长春光机所和中国地质调查局南京地质调查中心等。从20世纪80年代开始，上海技物所研制出OMIS高光谱成像系统，深圳中达瑞和科技有限公司研制了液晶可调谐滤波器(LCTF)，基于LCTC的成像光谱仪SHIS等。中国地质调查局南京地质调查中心在国家重大仪器开发专项和国家重点研发计划支持下，成功研制了以推扫式小型化成像光谱仪HMS400/1000为核心的轻型机载高光谱成像系统SSMap(魏桂华等，2019)，已在江苏、安徽、浙江等地进行了生态地质环境调查应用示范，取得了高质量的机载遥感数据，获得了一批重要调查成果。
2. 机载高光谱基本原理及其应用领域
2.1 基本原理
机载高光谱波段通常工作在400~2500 nm(0.4~2.5 μm)范围内，由可见光和近红外两个谱段组成，可见光谱段为400~780 nm(0.4~0.78 μm)，主要用于区分岩心颜色、植被、水体等；近红外光谱的波长在780~2500 nm(0.78~2.5 μm)，主要用于区分含羟基的物质，对此波段的近红外光谱产生吸收的官能团主要是含氢基团，包括C-H (甲基、亚甲基、甲氧基、羧基、方基等)、羟基O-H、巯基S-H、氨基N-H等，它们的合频和一级倍频位于780~2500 nm波段。由于含有羟基物质晶格中原子间的化学键的弯曲、伸缩和电子能级跃迁吸收某些区域的近红外光谱，所以根据某些官能团的在近红外区域的这种特征吸收光谱可以区分不同的物质(修连存等，2007)。
2.2 应用领域
随着机载高光谱遥感技术的发展，数据处理方法成熟，广泛应用于地质调查、环境监测、海洋生态评价、土地质量评估、精准农业、林业生态和城市规划领域(Berger et al., 2003)。在地质调查领域中，利用机载高光谱遥感具有空间分辨率高，能够精确反映地物目标这一特质，可以区分层状硅酸盐中单矿物(黏土矿物，绿泥石，蛇纹石等)，含羟基之硅酸盐矿物(绿帘石，闪石等)，硫酸盐矿物(明矾石，黄铁钾矾，石膏等)，碳酸盐矿物(方解石，白云石等)；而矿物种类揭示矿化作用过程中热液蚀变体系结晶时的温度、压力和化学环境，矿物中的离子交换反映在羟基的波长移动上，矿物的含量反映在羟基的吸光度上，这些参数揭示地质事件和成矿规律(修连存等，2009)。在环境监测领域，可迅速识别植物、建筑、水质、道路、湿地、土壤等分布与面积，通过建模，评价水土环境质量、农田作物长势和病虫害等。
高光谱遥感具有图谱合一的特点(图 1)，避免了以往航摄带来的弊端，它不仅能够用图像表达目标，还可以根据光谱识别目标，利用高光谱遥感技术，能够高效地对长江经济带生态环境、经济发展和城镇布局进行评价，为政府规划和治理提供技术服务。
3. 机载高光谱成像系统SSMap研发及参数指标
3.1 系统研发
SSMap采用中国地质调查局南京地质调查中心自主研发的VRIN可见光-近红外高光谱成像仪与SWIR短波红外高光谱成像仪(图 2)进行机载航空高光谱测量系统集成(图 3~图 5)，并研发了配套的软件对机载高光谱遥感数据进行预处理与信息提取。
![]()
图 2 VNIR可见光-近红外及SWIR短波红外成像光谱仪Figure 2. Visible-near-infrared and short-wave infrared imaging spectrometer![]()
图 3 有人机系统组成(a)和无人机系统组成(b)1—姿态传感器；2—短波红外成像光谱仪；3—可见光成像光谱仪；4—稳定平台；5—诺瓦泰惯性导航仪；6—支撑架；7—三轴稳定平台；8—短波红外成像光谱仪Figure 3. The composition of the manned vehicle system (a) and the composition of the UAV system (b)1- Attitude sensor; 2- Short- wave infrared imaging spectrometer; 3- Visible light imaging spectrometer; 4- Stable platform; 5- Novartis inertial navigator; 6-Support frame; 7-Three-axis stabilized platform; 8-Short-wave infrared imaging spectrometer![]()
图 4 直升机(a)运输机(b)固定翼无人机(c)和旋翼无人机(d)飞行平台Figure 4. Helicopter (a) transport aircraft (b) fixed-wing UAV (c) and rotor UAV (d) flying platform![]()
图 5 有人机平台安装(a)和无人机平台安装(b)Figure 5. The installation of the manned vehicle system (a) and the UAV system (b)系统由成像光谱仪、三轴稳定平台、惯性导航仪和计算机控制与采集模块等多个部分组成。
3.2 主要技术参数
表 1 小型成像光谱仪参数Table 1. Parameters of the small-size imaging spectrometer
表 2 惯导定位精度参数Table 2. The INS positioning accuracy parameters
3.3 数据处理软件研发
机载高光谱数据具有数据量大、处理过程复杂等特点，因此需要高性能的服务器电脑和专业的处理软件。硬件部分由多台高性能戴尔服务器电脑组成的运算系统，其服务器操作系统性能稳定，计算能力强，使数据处理工作得到了可靠保障。机载高光谱数据处理采用NovAtel的IE8.8软件进行Pos数据后差分解算，采用自主研发的HyperPic软件进行光谱数据处理，包括飞行规划、大气校正、几何校正、辐射定标和反射率反演等功能(Bierwirth et al., 2002)。
4. 机载高光谱数据采集流程与方法
4.1 系统安装及检查
飞行前将仪器设备安装在固定翼飞机、直升机或者无人机等平台上。检查系统外观是否有意有磕、碰和摔等情况。检查测量系统及其配套设备、地面测量基站等工作状态是否正常，确保系统各个部分组装牢固安全。电源、电缆连接正确无误，各仪器设备读数正常。进行光谱成像仪的数据采集测试，保证设备工作正常。进行惯性导航仪检查，确保GPS导航定位接收信号良好，通讯系统保持畅通。对系统进行干扰测试，保证数据收录正常。
4.2 野外数据采集
航空数据采集作业流程包括现场踏勘、飞行准备、飞行测量和地物测量。
飞行前，工作小组需对测区边缘和测区内进行现场踏勘，熟悉测区的地形地貌，对于地形复杂不容易达到的测区，通过地形图进行了解，确定地面测量基站和像控点的布设选址。通过调研，确定最佳机场，制定应急迫降方案。
飞行准备包括：收集气象资料，了解适航天气；收集测区内的地图并扫描成电子版；地图裁切，输入角坐标，进行地图定位，作为飞行导航底图；由专业机组人员对飞机及成像系统进行飞行前例行检查，检查飞机安全状况，成像系统的各接口及仪器性能稳定性进行调试监测和评估。
在飞行测量中，飞机引擎启动后，接通电源，启动仪器设备，检查设备状态是否正常，并按照设计要求，设置相关技术参数。起飞后，打开镜头窗口遮挡板。进入测线前2 km开始记录数据，并于飞出测线2 km后停止记录数据，检查记录状态，监测仪器状态，根据信号强弱调整相关参数。操作人员要时刻监测航迹、航高、航速是否符合要求，同时监测气象条件(云量、能见度)，填写飞行记录表。如不符合设计要求，及时进行补飞，补飞按原设计要求进行。飞行结束后，根据机上定位设备的需要，在静止状态下收集满足定位要求的数据后，关闭仪器，备份数据，为下一架次测量存储做好准备。
地面数据的获取主要包括GPS地面基站测量、地面光谱采样、水样采集等三个部分，在航飞数据采集的同时，地面同步进行测量工作。
4.3 GPS地面基站测量
使用地面卫星导航定位基站数据做后差分计算，在飞行前30 min打开地面基站，在卫星导航定位基站搜索卫星成功后(至少8颗卫星)，采用静态测量模式记录数据，数据记录频率高于1 Hz。飞行结束半小时后停止记录，关闭地面基站，导出基站数据并备份。
4.4 光谱采样及采样点定位
地面光谱采样采用宽谱段地物波谱仪(350~2500 nm)，采样点定位采用中海达iRTK2型号设备和江苏Cors定位系统，定位精度达厘米级。作业需要获取航空组规划好的飞行路线图，根据飞行路线图，设计地面定标路线。选择地面定标路线要求兼顾飞行路线条件与地面地物条件，每条航线范围内有至少一个地面采样点，每个采样点采样的地物种类数量视具体情况确定，在确保时间足够的情况下，尽量多采地物光谱。在地面道路条件允许的情况下，尽量使定标路线垂直或近似垂直地斜交航线，保证每条航线地面采样时间与飞行时间相差不大。在符合以上基本条件的地面路线中，选取沿途地物种类丰富，房屋建筑物遮挡较少，便于采样作业的地面路线。飞行测量过程中，地面光谱采样工作同步进行，并记录采集点位置、时间和对应的架次及航带，测量完成后，备份光谱数据和定位数据(荆平平等，2017；孙家波等，2018；蔡祥等，2021)。图 6为采集的几种典型地物地面光谱曲线。
![]()
图 6 几种典型地物地面光谱曲线1—草地；2—槐树；3—沥青路；4—滩涂；5—水泥路；6—小麦Figure 6. Spectrum curves of several typical ground features1- Grassland; 2- Locust tree; 3- Asphalt road; 4- Beach; 5- Cement road; 6-Wheat4.5 水体样采集及现场测试
采样前要根据检测项目的性质和采样方法要求，选择适宜材质的盛水容器和采样器，并充分清洗干净(李鑫星等，2021)。此外，还要准备好适当的交通工具。采样容器一般使用塑料容器或玻璃容器。塑料容器用作测定金属、放射性元素和其他无机物的盛水容器，玻璃容器则用作测定有机物和微生物等的盛水容器。有些项目要求提前加入保存剂，不同项目加入保存剂不同，根据保存剂的不同分别准备盛水容器。同时还准备一些采样的工具，比如水桶，采样器，现场使用的实验设备：pH计、温度计等。在采集河流、湖泊等地表水时，可以使用桶、瓶等容器直接采取。一般将其沉至水面下0.3~0.5 m处采集。在桥上等地方采样时，可将系着绳子的聚乙烯桶或带着坠子的采样瓶投于水中汲水，要注意不能混入漂浮于水面上的其他物质。采集深层水样时，可以使用带重锤的采样器。测定溶解氧的水样，要使用专门的双瓶采样器采样。湖泊、水库的采样点应该设置在湖水的主要出入口、中心区、沿湖泊(水库)水流方向、滞留区及湖边城市水源区(罗代洪等，2004；甘露等，2010；崔齐等，2018)。大部分水样在当天采样工作完成后，立即送往专业分析机构进行测试，测试的指标见表 3。
表 3 水常规化验指标(水质参数)一览Table 3. List of water routine test indexes
5. 高光谱数据处理与信息提取
5.1 高光谱数据预处理
高光谱机载平台获取数据时，由于平台高度及光谱仪视场角限制，导致数据获取幅宽有限，需多架次飞行获取的多航带数据才能完整覆盖研究区域(王瑞军等，2020)。因此机载高光谱数据预处理包括辐射校正、几何校正和图像拼接等工作，以获得在空间和辐射能量探测上的精准匹配的完整研究区域影像。最后通过同步获取的地物波谱数据验证影像的辐射精度，以进一步检验成像数据质量。
5.1.1 辐射定标
由于成像光谱仪输出的数据并不具有物理意义，因此通过辐射定标将成像光谱仪的DN值转为有物理意义的辐射亮度值，使遥感信息定量化。采用积分球系统(图 7)完成成像光谱仪的辐射定标。
通过对比积分球能级数据即对应功率、积分时间下输入的辐射亮度值与光谱仪输出DN值之间的关系，逐波段进行线性拟合，完成光谱仪的辐射定标。
在辐射定标前需要对积分球辐亮度数据进行线性插值，得到光谱仪的各波段所对应的辐亮度值。影像DN值与辐亮度值之间存在着线性关系，通过求解高光谱影像的各波段增益系数以及偏移量，完成影像数据的辐射定标(Moamen et al., 2021)。图 8和图 9为可见光影像和短波红外影像辐射定标后的光谱变化。
![]()
图 9 短波红外数据辐射定标结果Figure 9. Radiometric calibration results of short-wave infrared data5.1.2 反射率转换
在经过辐射定标后，求得辐亮值数据，进一步使用经验线性法进行反射率转换，参与反射率转换的地面同步光谱数达到324条，平均每一条航带2条以上的同步地物波谱。
基于经验模型的反射率反演算法，主要是应用图像像元自身的灰度值，不考虑大气、地物属性以及遥感平台特征，通过对象元灰度值进行求均值或回归运算，在统计意义上获得图像的反射率值(孙雨等，2015)。
经验线性法(Empirical Line, EL)已经在遥感定标和反射率反演中被广泛使用，主要依靠图像象元值和图像背景信息之间的拟合关系(Ben-Dor et al., 2004)完成地表反射率的反演计算。经验线性法基于如下简化的公式：
(1) 其中：DNb为给定象元在波段b的数字量化值；ρ(λ)为实测地物象元在波段b所在的波长的地表反射率值；Ab为由于传输和仪器本身所导致的倍数增益系数；Bb为由于大气程辐射和仪器导致的偏移。
该方法在应用中, 一般会实测地面两个或多个定标点的地面反射光谱数值(reflectance spectra)，然后计算遥感图像上对应象元点的平均辐射光谱值, 通过使用统计学中的回归运算，得到公式(1)中的增益值Ab和偏移值Bb，最后应用该公式，对整幅遥感图像进行反射率的计算。其中某一地物波谱反射率与地面波谱反射率对比如图 10。
![]()
图 10 槐树反射率和ASD数据对比1—影像光谱—槐树；2—ASD光谱—槐树Figure 10. Comparison of locust tree reflectance and ASD data1-Image spectrum-locust tree; 2-ASD spectrum-locust tree5.1.3 几何校正
首先进行IMU/GPS数据处理，Pos数据采用诺瓦泰公司的IE8.8软件进行处理(图 11)。该软件一般可以采用精密单点定位PPP和基站两种方式解算POS数据，本次采用基站后差分处理方式解算POS数据，精密计算每行数据于曝光时刻的机载GPS天线相位中心的WGS84框架坐标和姿态参数。
在完成POS数据的解算后进行插值处理，得到了高光谱影像各个扫描行所对应的外方位信息，可以通过建立旋转矩阵进行坐标转换，将导航坐标转换为投影坐标。具体坐标的转换过程如下：成图坐标系(m)→导航坐标系(g)→IMU坐标系(b)→传感器坐标系(c)→像空间坐标系(i)，由成图坐标系转换至像空间坐标系的旋转矩阵表达形式如下：
(2) 式中旋转矩阵可由已知条件观测获取，通过矩阵变换求得三个外方位元素姿态参数[(ω, ϕ, κ)]，三个外方位元素位置参数[(Xs, Ys, Zs)]需要考虑IMU坐标系中的坐标位置偏移[(Xt, Yt, Zt)]进行修正。
(3) 在求得所有外方位元素[(ω, ϕ, κ, Xs, Ys, Zs)]后，根据共线条件方程，采用直接法几何校正计算像点坐标[(x, y)]对应的地面点坐标[(X, Y, Z)]，并赋予其像点对应的灰度值。坐标变换公式为：
(4) 粗校正后的影像，精度还存在一定的误差，本实验利用采集的像片控制点对所有有效的像片均进行微分纠正。按照航带分区域纠正，最后通过镶嵌过程，得到最终高光谱影像。
5.2 地物分类与信息处理
工作中采用可见光数据287个波段进行地物的精细分类。光谱分辨率和空间分辨率分别为2.42 nm和1.0 m。地面光谱采样主要包含各类农作物和树木、草地、道路材质等68个地物类别，根据工区内地物存在的实际情况，仅选取了农作物、树木、草、水生植物、道路材质、建筑、裸露地表等七大类别，共计46种地物种类进行精细分类，具体如表 4。工区内的其他地物由于数量较少或者光谱干扰严重无法区分，如构树、枣叶、茶树、银杏、豇豆、葡萄、芒草、花岗岩地砖等，其中建筑类地物结合POI信息进行了人工解译划分，大致划分为11类用地信息，具体为：一般工厂、加油站、化工企业、医院诊所、商业购物、学校、居民住宅、广场公园、建筑零时用地、码头、行政机关等(李志忠等，2021)。
表 4 地物精细分类明细Table 4. Detailed classification of ground features
由于高光谱图像波段数目多，各波段间具有较强的相关性，主成分分析(PCA)方法对高光谱数据进行降维和去除噪声波段。支持向量机方法对高光谱遥感图像进行分类，可实现图像的分类识别，分类步骤及流程见图 12。
![]()
图 12 基于PCA和SVM的高光谱遥感影像分类流程图Figure 12. Flow chart of hyperspectral remote sensing image classification based on PCA and SVM主成分分析(PCA，Principal Components Analysis)是一种简化数据集的技术。它是一个线性变换。这个变换把数据变换到一个新的坐标系统中，使得任何数据投影的第一大方差在第一个坐标(称为第一主成分)上，第二大方差在第二个坐标(第二主成分)上，依次类推。主成分分析经常用减少数据集的维数，同时保持数据集对方差贡献最大的特征。
持向量机的基本思想是寻找一个分类超平面，使得训练样本中的两类样本能被分开，支持向量机(Support Vector Machine，SVM)作为一种最新的也是最有效的统计学习方法，近年来成为模式识别与机器学习领域一个新的研究热点。SVM具有小样本学习、抗噪声性能、学习效率高与推广性好的优点，能够用于解决空间信息处理分析领域的遥感影像处理。遥感图象分析与处理是SVM应用一个热门的研究方向。一些主要应用如土地利用分类，目前针对高光谱数据进行分类还是仅仅局限于传统的分类方法，不但运算速度慢，分类精度低，而且出现了严重的huges现象。而在高光谱遥感分类中SVM具有明显的优越性，因此SVM应用被归纳为高光谱遥感分类最重要的进展之一。通过飞行数据，研究了支持向量机在高光谱遥感图像分类中的应用，建立了基于支持向量机的高光谱遥感图像分类模型。
6. 机载高光谱遥感应用示范
6.1 沿江地物、水体和岸线稳定性探测
采用机载高光谱遥感技术，对长江下游干流镇江—扬州段沿江两岸岸线外5 km范围内进行生态地质环境遥感调查，总面积约670 km2，飞行测线长度2500 km。快速识别、提取工作区中的各项环境指标和影响因子，进行生态地质环境评价。飞行区域规划如图 13所示。
6.1.1 研究目标
在采集的高光谱数据基础上，识别和提取土地、植被、湿地、岸段、滩涂等自然资源要素，进行综合研究，实施生态地质环境评价，给出工作区内各项功能区的开发利用建议；结合水体采样化验分析和数学建模方式，对工作区内的水体进行半定量的评价，识别和查找污染源头，给出水产养殖用水、灌溉用水建议；利用光谱特征，圈定容易产生冲蚀崩塌的岸段；调查了解江滩湿地的植物种群分布及生长特点(苏龙飞等，2021)。
6.1.2 研究成果
(1) 地物精细分类
地物精细分类面积为670 km2，分类地物共计46类，具体如表 5。分类结果(局部)见图 14~图 16。其中分类精度最大90.2%，最小68.69%，平均分类精度达到80.57%，分类精度和指标符合要求。
表 5 地物精细分类统计Table 5. Statistics of fine classification of ground features
![]()
图 14 地物分类结果相对百分比图1—水稻；2—青菜；3—菱角；4—荷叶；5—大豆；6—上海青；7—花生；8—油菜；9—芝麻；10—萝卜；11—草莓；12—荞麦；13—桑树；14-杨树；15—女贞；16—冬青；17—栾树；18—松柏；19—香樟树；20-夹竹桃；21—杉树；22—桦树；23—竹子；24—罗汉松；25—玉兰；26—柽柳；27—海桐；28—松树；29—李子树；30—天竺葵；31—榉树；32—绿化草地；33—杂生草丛；34—芦苇；35—水体；36—藻；37—浮萍；38—水葫芦；39—一般建筑；40—工矿企业建筑；41—种植大棚；42—滩涂；43—裸土；44—沥青；45—水泥；46—未分类地物Figure 14. The relative percentage of ground feature classification results1- Rice; 2- Green vegetables; 3- Water chestnut; 4- Lotus leaf; 5- Soybean; 6- Shanghai green; 7- Peanut; 8- Rape; 9- Sesame; 10- Radish; 11- Strawberry; 12- Buckwheat; 13- Mulberry; 14- Poplar; 15- Privet; 16- Holly; 17- Luan tree; 18- Pine and cypress; 19- Camphor tree; 20-Oleander; 21-Fir tree; 22-Birch tree; 23-Bamboo; 24- Podocarpus; 25- Magnolia; 26- Tamarix; 27- Pitrus; 28- Pine; 29- Plum; 30- Geranium; 31- Zelkova; 32- Green grassland; 33- Weedy grass; 34-Reed; 35-Water body; 36-Algae; 37-Duckweed; 38- Water hyacinth; 39- General construction; 40- Industrial and mining enterprise construction; 41- Planting greenhouse; 42- Tidal flat; 43-Bare soil; 44-Asphalt; 45-Cement; 46-Unclassified ground features![]()
图 15 高光谱遥感功能区分布图Figure 15. Distribution of hyperspectral remote sensing functional areas![]()
图 16 高光谱遥感地物精细分类图Figure 16. Fine classification of hyperspectral remote sensing features(2) 地表水水质评价
本次高光谱遥感探测利用高光谱数据与地面实际水样测试分析数据，建立相关性模型，根据《地表水环境质量标准》(GB3838-2002)，对工作区内除长江以外的地表水体进行水质反演评价，评价指标包括总磷、总氮、氨氮、化学需氧量、五日生化需氧量、溶解氧、叶绿素7个指标。评价结果表明：Ⅲ类水占比最少，仅占0.03%；Ⅳ类水占比13.12%；Ⅴ类水体占比最多，为58.13%；劣Ⅴ类水体占比28.72%，主要集中在A、B、C、D、E五块区域(图 17~图 20)。其中，A区位于工作区西北侧，区东侧为工厂聚集区，北部为居民生活区，南部为大棚种植、养殖集中区，此处水质较差可能是因为工厂活动和施肥因素的影响；B区位于工作区西侧，区域内有塑料包装、建筑工程、机械配件等工厂，有较为密集的居民生活区，此处水质较差可能是受到工厂施工、居民生活用水排放的影响；C区位于工作区中部，区域内有钢铁厂、鞋厂等工厂，初步判断水体受到了工厂施工环境的影响；D区位于工作区北侧偏西，区内有河流，水质较差区域主要分布在居民生活区的周围，可能受到居民生活的影响；E区位于工作区东侧长江南岸，区内为大片水体种植、养殖区，养殖施肥导致此处水质较差。
![]()
图 17 高光谱遥感地表水综合等级划分图1—工作区范围；2—I类水；3—II类水；4—III类水；5—IV类水；6—V类水；7—劣V类水Figure 17. The comprehensive grade classification of surface water by hyperspectral remote sensing1-Working area scope; 2-Class I water; 3-Class II water; 4-Class III water; 5-Class IV water; 6-Class V water; 7-Inferior class V water![]()
图 18 高光谱遥感地表水总氮(N)元素浓度反演分级图1—工作区范围；2—I类水；3—II类水；4—III类水；5—IV类水；6—V类水；7—劣V类水Figure 18. Inversion and classification of total nitrogen (N) element concentration in surface water by hyperspectral remote sensing1-Working area scope; 2-Class I water; 3-Class II water; 4-Class III water; 5-Class IV water; 6-Class V water; 7-Inferior class V water![]()
图 19 高光谱遥感地表水总磷(P)元素浓度反演分级图1—工作区范围；2—I类水；3—II类水；4—III类水；5—IV类水；6—V类水；7—劣V类水Figure 19. Inversion and classification of total phosphorus (P) element concentration in surface water by hyperspectral remote sensing1-Working area scope; 2-Class I water; 3- Class II water; 4- Class III water; 5- Class IV water; 6- Class V water; 7-Inferior class V water![]()
图 20 高光谱遥感地表水五日生化需氧量(BOD5)反演分级图1—工作区范围；2—I类水；3—II类水；4—III类水；5—IV类水；6—V类水；7—劣V类水Figure 20. Five-day biochemical oxygen demand (BOD5) inversion classification of surface water by hyperspectral remote sensing1-Working area scope; 2-Class I water; 3- Class II water; 4- Class III water; 5- Class IV water; 6- Class V water; 7-Inferior class V water(3) 江滩湿地评价
对研究区内主要三处江滩湿地植被进行了分类，湿地区1位于镇江市焦山风景区以北，面积约16.84 km2。湿地区2位于共青团农场以南，面积约10.03 km2。湿地区3位于扬中市寿字圩以北，面积约6.13 km2。
分类结果表明(图 21~图 23)，湿地均存在不同程度的开发利用，开发利用的主要内容为水产养殖、苗圃、水稻。其中湿地1的水产养殖水域面积占比达22.32%，湿地2的水稻种植面积占比达18.89%，湿地3的水稻种植和苗圃积植面积加起来达到了42.17%，由此可见，这些人工开发肯定会对湿地生态环境产生不利影响。
![]()
图 21 湿地区1地物精细分类结果图1—水；2—水泥路面；3—水杉；4—柳树；5—蒿草；6—芦苇；7—水草；8—绿藻；9—泥路；10—毛草Figure 21. Results of fine classification of ground features in wet area 1#1-Water; 2-Cement road; 3-Metasequoia; 4-Willow; 5-Wormwood; 6-Reed; 7-Aquatic; 8-Green algae; 9-Dirt road; 10-Grass![]()
图 22 湿地区2地物精细分类结果图1—水；2—水泥路面；3—裸土；4—水稻；5—泥路；6—房屋；7—冲积物；8—芦苇；9—水草；10—柳树；11—香樟；12—绿藻；13—其他硬化地表；14—桦树；15—杨树；16—毛草；17—蒿草Figure 22. Results of fine classification of ground features in wet area 2#1- Water; 2- Cement road; 3- Bare soil; 4- Rice; 5- Dirt road; 6-House; 7- Alluvium; 8- reed; 9- Aquatic; 10- Willow; 11- Camphor; 12 - Green algae; 13- Other hardened ground; 14- Birch; 15- Poplar; 16- Grass; 17- Wormwood![]()
图 23 湿地区3地物精细分类结果图1—水；2—绿藻；3—房屋；4—水泥路面；5—泥路；6—砂石；7—未分类；8—冲积物；9—裸土；10—水稻；11—水草；12—芦苇；13—杨树；14—毛草Figure 23. Results of fine classification of ground features in wet area 3#1- Water; 2-Green algae; 3- House; 4- Cement road; 5- Dirt road; 6-Sand gravel; 7-Unclassified; 8 alluvium; 9- Bare soil; 10- Rice; 11- Aquatic; 12-Reed; 13-Poplar; 14-Grass(4) 长江岸段稳定性评价
利用高光谱数据的精细分类结果，在GIS(地理信息系统)软件支持下，对研究区内的江岸稳定性进行评价，评价因子包括岸线材质、岸前滩涂、岸前滩涂植被、岸后植被等。评价结果(图 24)表明，相对稳定岸段多分布于堤坝、沿江公路、港口、公园；相对稳定性较差岸段主要是自然岸段，主要分布于镇江市江心洲南北两岸、扬洲江都区自茅家圩往北直至王家桥一带、扬中市寿字圩以北环江路附近。
![]()
图 24 岸基稳定评价1—相对稳定性较差；2—稳定性一般；3—相对稳定Figure 24. Shore-based stability evaluation1-Relatively poor stability; 2-Generally stable; 3-Relatively stable6.2 土壤污染探测
研究区位于江苏省镇江市丹阳东北部地区，地貌以平原为主，地势较平坦，局部有较小范围的剥蚀残丘，最高峰为123 m，其余的地区为平原，海拔高度一般在7 m。飞行面积310 km2，飞行测线总长度为740 km(图 25)。
6.2.1 研究目标
通过机载高光谱技术采集耕地土壤高光谱影像，并进行同步野外地物波普测量，进行耕地土壤地面检测和高光谱传感器数据融合(杨思睿等，2018)。采集耕地土壤具有代表性的实验样品，实验室测定其全氮含量、有机质含量、水分含量、土壤肥力以及重金属污损情况等，与光谱特性分析结果对比，确定耕地土壤生态参数及其变化与光谱特性基本关系，建立耕地土壤影响因子快速监测的高光谱反演模型。
6.2.2 研究成果
通过反演结果可以看出，丹阳地区元素As的污染区域较多，主要集中研究区域的西北、西南和东部地区；元素Cr主要集中在中北部地区，而元素Cu则集中在中南部地区(图 26~ 图 28)。根据其分布可以发现，除了相关性较高的Hg和Ni外，其他污染元素分布并没有一定的相关性，分布区域各不相同。
![]()
图 26 丹阳地区元素As高光谱反演结果图Figure 26. Hyperspectral inversion results of element As in the Danyang area![]()
图 27 丹阳地区元素Cr高光谱反演结果图Figure 27. Hyperspectral inversion results of element Cr in the Danyang area![]()
图 28 丹阳地区元素Cu高光谱反演结果图Figure 28. Hyperspectral inversion results of element Cu in the Danyang area通过机载高光谱反演影像与实际采集样品分析结果插值对比分析，两者在数值趋势上基本一致(图 29)，说明通过高光谱影像进行污染元素含量反演具有较强的说服力。丹阳地区土壤重金属污染元素主要为As、Cr、Cu。
![]()
图 29 丹阳地区土壤元素As实测插值分析图Figure 29. Interpolation analysis diagram of measured soil element As in the Danyang area实地调研后发现，图中污染元素含量较高的区域基本都位于工业区，其中，与汽车行业的相关工业园对污染元素Cu的影响较大，而一些家具、橡胶加工行业对污染元素Cr和As的影响较大。
6.3 生态环境调查
6.3.1 研究目标
通过机载高光谱成像技术对安徽省池州市贵池区进行生态环境调查。通过采集的高光谱影像、同步水体样采集测试与实际土壤样品分析，反演得出藻类分布、土壤重金属污染分布及矿山开采区的粉尘污染分布信息。
6.3.2 研究成果
利用航空高光谱数据对安徽贵池秋浦河沿岸地区生态环境调查发现，水体藻类有着特征波段，通过反演能够有效提取到藻类信息(图 30)，精度可达74.33%以上。
由机载高光谱反演的影像与实际采集土壤样品分析结果插值对比分析，两者结果基本一致(图 31，图 32)，进一步说明了机载高光谱探测土壤污染的有效性。
![]()
图 31 机载高光谱Cr元素反演影像(上)与实际采集土壤样品分析结果(下)对比Figure 31. Comparison of the airborne hyperspectral Cr element inversion image (above) and the actual collected soil samples (below)![]()
图 32 机载高光谱Cu元素反演影像(上)与实际采集土壤样品分析结果(下)对比Figure 32. Comparison of the airborne hyperspectral Cu element inversion image (above) and the actual collected soil samples (below)在贵池某矿山开采区通过机载高光谱反演影像可以清晰识别出矿山粉尘污染分布(图 33)。
![]()
图 33 高光谱遥感识别贵池某矿山粉尘(红色区域)污染分布Figure 33. Hyperspectral remote sensing to identify the distribution of dust pollution(red area)in some mine in Guichi7. 结论
高光谱成像技术将成像技术与光谱技术相结合，探测目标的二维几何空间及一维光谱信息，获取高光谱分辨率的连续、窄波段的图像数据。随着机载高光谱遥感技术的发展，数据处理方法日益成熟，这项技术被广泛应用于地质调查、环境监测、海洋生态评价、土地质量评估、精准农业、林业生态和城市规划等领域。中国地质调查局南京地质调查中心在国家重大仪器开发专项和国家重点研发计划支持下，成功研制出以推扫式小型化成像光谱仪HMS400/1000为核心的轻型机载高光谱成像系统SSMap，建立了机载高光谱遥感综合调查技术方法和作业流程，形成了一套可靠的地物分类方法和水土污染光谱定量反演模型。该系统具有快捷、高效和高空间分辨率的特点，已在江苏、安徽、浙江等地进行了生态地质环境调查应用示范，总飞行面积逾4500 km2，获得了一批重要成果。研究表明，在生态环境调查、水体监测、土地利用和土地质量评价等方面具有广阔的应用前景。具有自主知识产权的国产化机载高光谱成像仪降低了采购成本，推动了机载高光谱技术普及，服务了国家重大区域发展战略，具有显著的经济和社会效益。
-
![]()
图 1 研究区地质简图
1—第一旋回第一亚旋回火山岩；2—第一旋回第二亚旋回火山岩；3—第一旋回第三亚旋回火山岩；4—第二旋回火山岩；5—第三旋回火山岩；6—早侏罗世枫坪组；7—潜火山岩；8—前泥盆纪变质岩、侵入岩；9—溢流相流纹岩/石泡流纹岩；10—火山颈相集块角砾熔岩；11—空落相流纹质含角砾晶玻屑凝灰岩；12—流纹质玻屑凝灰岩；13—碎屑流相流纹质晶屑熔结凝灰岩；14—喷发沉积相凝灰质砂岩；15—角度不整合界线；16—断层界线；17—同位素采样位置及编号；18—剖面位置
Figure 1. Geological map of the study area
1-Volcanic rocks of Ⅰ1 cycle; 2- Volcanic rocks of Ⅰ2 cycle; 3- Volcanic rocks of Ⅰ3 cycle; 4- Volcanic rocks of Ⅱ cycle; 5- Volcanic rocks of Ⅲ cycle; 6-Fengping Formation of early Jurassic; 7-Subvolcanic rock; 8-Pre-Devonian metamorphic rock and intrusive rock; 9-Flooding rhyolite/ lithophysa rhyolite; 10-Volcanic neck facies breccia and conglomerate lava; 11-Rhyolitic breccias crystal glass tuff formed by Empty stacked phase; 12-Rhyolitic glass tuff formed by Empty stacked phase; 13-Rhyolitic crystal welded tuff formed from pyroclastic flow accumulation; 14- Tuffaceous sandstone formed from eruption-sedimentary; 15-Angular unconformity; 16-Fault; 17-Isotope sampling position; 18-Profile position
![]()
图 2 研究区火山岩柱状图简图
1—流纹质晶玻屑熔结凝灰岩；2—流纹质含角砾晶屑熔结凝灰岩；3—流纹质角砾熔结凝灰岩；4—流纹质晶玻屑凝灰岩；5—流纹质凝灰熔岩；6—（石泡）流纹岩；7—流纹质凝灰岩；8—凝灰质砂岩
Figure 2. Columnar section of the volcanic rocks in the study area
1-Rhyolitic crystal-vitric welded tuff; 2-Rhyolitic breccias crystal glass welded tuff; 3-Rhyolitic breccia welded tuff; 4-Rhyotaxitic crystal-vitric tuff; 5-rhyolitic tufflava; 6-Lithophysa rhyolite; 7- Rhyolitic tuff; 8-Tuffaceous sandstone
![]()
图 3 流纹质熔结凝灰岩（a，b）、流纹岩（c）及石泡流纹岩（d）显微照片（正交偏光）
qz—石英；kfs—钾长石
Figure 3. Microphotograph of rhyolitic welded tuff (a, b), rhyolite (c) and lithophysa rhyolite (d) (crossed nicols)
qz-Quartz; kfs-K-feldspar
![]()
图 5 研究区晚侏罗世火山岩锆石U-Pb年龄谐和图
Figure 5. Zircon U-Pb concordia diagrams of the volcanic rocks in the study area
![]()
图 6 研究区火山岩判别图解
a—火山岩TAS图解（Le Maitre，1989；Pc—苦橄玄武岩；B—玄武岩；O1—玄武安山岩；O2—安山岩；O3—英安岩；R—流纹岩；S1—粗面玄武岩；S2—玄武质粗面安山岩；S3-粗面安山岩；T—粗面岩、粗面英安岩；F—副长石岩；U1—碱玄岩、碧玄岩；U2—响岩质碱玄岩；U3-碱玄质响岩；Ph—响岩；Ir—Irvine分界线，上方为碱性，下方为亚碱性）；b—A/CNK-A/NK图解（Rickwood，1989；Le Maitre et al., 1989）；c-FAM图解（Irvine, 1971；TH—拉斑玄武岩系列；T—钙碱性火成岩系列）；d—SiO2-K2O图解（据Rickwood，1989；Le Maitre et al., 1989）
Figure 6. Diagrams of the volcanic rocks in the study area
a-TAS plot for the volcanic rocks(after Le Maitre, 1989; Pc-Picrite basalt; B-Basalt; O1-Basaltic andesite; O2-Andesite; O3-Dacite; R-Rhyolite; S1-Trachybasalt; S2- Basaltic trachyandesite; S3-Trachyandesite; T-Trachyte, trachydacite; F-Feldspathoidite; U1-Tephrite, basanite; U2-Phonolitic tephrite; U3-Tephriphonolite; Ph-Pphonolite; Ir-Irvine boundary, upper is basic, lowerw is alkaline); b-A/CNK-A/NK plot (after Rickwood, 1989; Le Maitre et al., 1989); c-FAM plot (after Irvine, 1971; TH-Tholeiite series; T- Calc-alkaline series); d-SiO2-K2O plot (after Rickwood, 1989; Le Maitre et al., 1989)
![]()
图 7 火山岩微量元素原始地幔标准化蛛网图（a）和球粒陨石标准化稀土配分图（b）（原始地幔和球粒陨石标准值均引自Sun and McDonough, 1989）
Figure 7. Primitive mantle-normalized trace element spidergram (a) and chondrite-normalized REE patterns (b) of volcanic rocks (after Sun and McDonough, 1989)
![]()
图 8 火山岩构造环境判别图解
a—Yb-Th/Ta图解（Michael et al., 2000）；b—Ta/Yb-Th/Ta图解（Michael et al., 2000）；c—火山岩戈蒂尼—里特曼指数图解（A区—非造山带地区火山岩；B区—造山带地区火山岩（岛弧及活动大陆边缘）；C区—A区、B区派生的碱性、富碱岩；J—日本火山岩）；d—Rb/10-Hf-Ta*3三角图解（Wood，1980）；e—Rb-（Yb+Ta）图解；f—Rb-（Y+Yb）图解（VAG—火山岛弧；WPG—板内；syn—COLG-同碰撞；ORG—洋中脊）
Figure 8. Discrimination diagrams illustrating tectonic setting of the volcanic rocks
a-Yb-Th/Ta diagram (after Michael et al., 2000); b-Ta/Yb-Th/Ta diagram (after Michael et al., 2000); c-logσ-logτ diagram (A-Non-orogenic volcanic zone; B-Orogenic volcanic zone; C-Alkaline and alkaline-rich derived from A and B; J-Japanese volcanic rocks); e- Rb-(Yb+Ta) diagram; f- Rb-(Y+Yb)diagram (VAG-Volcanic arc; WPG-Within plate; syn-COLG-Syn-collisional; ORG-Mid-ocean ridge)
![]()
图 9 火山岩年龄-Isr值图
M—幔源型；MC—壳幔混合源；C-1—下部壳源型；C-2—上部壳源型
Figure 9. Age-Isr of volcanic rocks
M-Mantle source; MC- Crust and mantle source; C-1-Lower crust; C-2-Upper crust
表 2 研究区火山碎屑岩和潜火山岩锆石U-Pb分析结果
Table 2 Zircon U-Pb dating results of the pyroclastic rock and subvolcanic rock in the study area
下载: 导出CSV
表 4 研究区火山岩微量元素测试结果（10-6）
Table 4 Trace elements (10-6) of volcanic rocks in the study area
下载: 导出CSV
表 5 研究区火山岩稀土元素测试结果（10-6）
Table 5 REE abundances (10-6) of volcanic rocks in the study area
下载: 导出CSV
-
Atherton M P, Tarney J. 1979. Origin of granite batholiths:Geochemical evidence[J]. Mineralogical Society. Geochemistry Group:1-148. https://www.researchgate.net/publication/256159340_Origin_of_granite_batholiths_Geochemical_evidence
Andersen T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb[J]. Chemical Geology, 192(1/2):59-79. https://www.researchgate.net/publication/222924679_Correction_of_common_lead_in_U-Pb_analyses_that_do_not_report_204Pb
Bureau of Geology and Mineral Resources of Zhejiang Province. 1989. Regional Geology of Zhejiang Province[M]. Beijing:Geological Publishing House (in Chinese with English abstract).
Bureau of Geology and Mineral Resources of Zhejiang Province. 1996. Stratigraphy (Lithostratic) of Zhejiang Province[M]. Wuhan:China University of Geosciences Press (in Chinese with English abstract).
Bao Chaomin, Wang Kongzhong. 1999. Mesozoic volcanicsedimentary strata and tectonicmovements in Zhejiang[J].Regional Geology of China, 18(2):201-204(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZQYD902.013.htm
Chen Jiangfeng, Guo Xinsheng, Tang Jiafu, Zhou Taixi, 1999. Nd isotopic model ages:Implications of the growth of the continental crust of southeastern China[J]. Journal of Nanjing University(Natural sciences), 35(6):649-658(in Chinese with English abstract).
Chen Peirong, Kong Xinggong, Wang Yinxi, Ni Qisheng, Zhang Bangtong, Ling Hongfei. 1999. Rb-Sr Isoptopic dating and significance of Early Yanshanian bimodal volcanic-intrusive complex from south Jiangxi Province[J].Geological Journal of China Universities, 5(4):378-383(in Chinese with English abstract).
Chen Peiji. 2000. Comments on the Classification and Correlation of Non-marine Jurassic and Cretaceous of China[J]. Journal of Stratigraphy, 24(2):114-119(in Chinese with English abstract). http://cn.bing.com/academic/profile?id=bc7137f0445a9b8f89443540fca23f45&encoded=0&v=paper_preview&mkt=zh-cn
Ding Baoliang, Lan Shanxian, Wang Yingping. 1987. Discussion on the non-marine Juro-cretaceous boundary in Zhejiang, Fujian and Jiangxi Provinces[J]. Bull. Nanjing Inst. Geol. M.R., Chinese Acad. Geol. Sci., 8(3):43-55(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-HSDZ198703004.htm
Dong Shuwen, Zhang Yueqiao, Chen Xuanhua, Long Changxing, Wang Tao, Yang Zhenyu, Hu Jianmin.2008. The formation and deformational characteristics of East Asia multi-direction convergent tectonic system in Late Jurassic[J]. Acta Geoscientica Sinica, 29(3):306-317(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqxb200803005
Dai Hongwen. 2011. On the relation between Late Jurassic volcanic rock and lnstrusive rock in Panyang County of Fuan City, Fujian Province[J]. Geology of Fujian, 30(4):335-340(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-FJDZ201104011.htm
Duan Zheng, Xing Guangfu, Yu Minggang, Zhao Xilin, Jin Guodong, Chen Zhihong. 2013. Time sequence and geological process of Late Mesozoic volcanic activities in the area of Zhejiang-Fujian boundary[J]. Geological Review, 59(3):454-469(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzlp201303006
Gu Zhiwei. 2002. Review on the genus Ferganoconcha[J]. Acta Palaeontologica Sinica, 41(1):152-158(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gswxb200201016
Gu Zhiwei.2005. On the Mincheian movement[J]. Journal of Stratigraphy, 29(1):1-6(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dcxzz200501001
Irvine T N, Barager W R A. 1971. A guide to the chemical classification of the common volcanic rocks[J]. Canadian Journal of Earth Sciences, 8:523-548. doi: 10.1139/e71-055
Jiang Weisan. 1995. Discussion on the upper ane basal limite of Cretaceous in Zhejiang Province[J]. Geology of Zhejiang, 11(1):15-23(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZJDZ199501001.htm
Le Maitre R W. 1989. A Classification of Igneous Rocks and Glossary of Terms[M]. Oxford:Blackwell, 193.
Ludwig K R. 2003. User's manual for Isoplot 3.0:A geochronological toolkit for Microsoft Excel[J]. Berkeley Geochronology Center Special Publication, (4):1-70.
Liu Chun, Yan Jun, Song Chuanzhong, Ling Quanzhong, Peng Ge, Shi Lei, Liu Xiaoqiang. 2012. Geochronology and geochemistryof the volcanic rocks from Fanchang basin in the Middle-Lower Yangtze River:Petrogenesis and geological significance[J]. Acta Petrologica Sinica, 28(10):3228-3240(in Chinese with English bstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-YSXB201210015.htm
Ma Wuping. 1994. Division of the late mesozoic strata in Zhejiang[J]. Journal of Stratigraphy, 18(2):91-101(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400025465
Michael P G, Eva S S. 2000. From continents to island arcs:A geochemical index of tectonic setting for arc-related and withinplate felsic to intermediate volcanic rocks[J]. The Canadian Mineralogist, 38:1065-1073. doi: 10.2113/gscanmin.38.5.1065
Nie Yonghong, Li Shuguang, Ge Ningjie, Liu Deliang, Hu Xiongjian. 1995. Metamorphic features and ages of Garnet amphibolites in Beiyuan, Southwestern Zhejiang Province, China[J]. Acta Mineralogica Sinica, 15(4):442-450(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-KWXB199504013.htm
Pearce J A, Harris N B W, Tindle A G. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks[J]. J. Petrol, 25:956-983. doi: 10.1093/petrology/25.4.956
Rickwood P C. 1989. Boundary lines within petrologic diagrams which use oxides of major and minor elements[J]. Lithos, 22:247-263. doi: 10.1016/0024-4937(89)90028-5
Sun S S, McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes[C]//Saunders A D, Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42 (1): 313-345.
Song Biao, Zhang Yuhai, Wan Yusheng, Jian Ping. 2002.Mountmakingand procedure of the SHRIMP dating[J]. Geological Review, 48:26-30(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZLP2002S1006.htm
Shu Liangshu, Zhou Xinmin.2002. Late Mesozoic tectonism of Southeast China[J]. Geological Review, 48(3):249-260(in Chinese with English abstract). http://cn.bing.com/academic/profile?id=0d840a6cd9030d1c0c2b2f58503a6a69&encoded=0&v=paper_preview&mkt=zh-cn
Shu L S, Zhou X M, Deng P, Wang B, Jiang S Y, Yu J H, Zhao X X. 2009. Mesozoic tectonic evolution of the southeast China block:New insights from basin analysis[J]. Journal of Asian Earth Sciences, 34:376-391. doi: 10.1016/j.jseaes.2008.06.004
Tao Kuiyuan, Xing Guangfu, Yang Zhuliang, Mao Jianren, Zhao Yu, Xu Naizhen. 2000. Determination of and discussion on the ages of mesozoic volcanic rocks in Zhejiang-Comments on the argumentation of Lapierre et al.[J]. Geological Review, 46(1):14-21(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzlp200001003
Williams I S. 1998. U-Th-Pb geochronology by ion microprobe[C]//Mickibben M A, Shanks ⅢW C, Ridley W I(eds.).Applications of Microanalytical Techniques to Understanding Mineralizing Processes. Reviews in Economic Geology, 7: 1-35.
Wang Qinghua. 2001. Discussion of ages of the Jiande group and the Moshishan group in zhejiang[J]. Volcanologu & Mineral Resources, 22(3):163-169(in Chinese with English abstract). http://cn.bing.com/academic/profile?id=380665551d0a406f7f8b9b84f84771ee&encoded=0&v=paper_preview&mkt=zh-cn
Huang Jialong, Lu Qingdi, Zhang Zhengyi, Nie Tongchun, Zhou Guowu, Lin Min, Zheng Ping. 2008. Establishment of a neostratotype section for the Nanyuan Formation in the Yuanzhuang area, Xianyuan, Fujian, China, and revision of its age[J]. Geological Bulletin of China, 27(6):785-792(in Chinese with English abstract).
Wang Duixing, Gao Wanli, Li Chunlin, Wang Zongxiu, Zhao Zhidan. 2015. LA-ICP-MS zircon U-Pb geochronology, petrochemistry of the Late Jurassic granite porphyry in central Zhejiang Province and their geological significanc[J]. Geology in China, 42(6):1684-1699(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-DIZI201506003.htm
Wang Jiaen, Wang Zhen, Liu Yuandong, Liang He.2016. Geological age and structural environment of volcanic rocks from the Laocun Fromation(Jiande group)[J]. Journal of Stratigraphy, 40(3):262-271(in Chinese with English abstract). http://www.en.cnki.com.cn/Article_en/CJFDTotal-DCXZ201603004.htm
Xu Butai, Yu Yunwen. 1997. The isotopic ages of mesozoic volcanic rock in Zhejiang and their geological time[J]. Acta Geoscientia Sinica, 18:117-119(in Chinese with English abstract). http://cn.bing.com/academic/profile?id=2940b532b88f5bf2d35cfb5f5a8298b2&encoded=0&v=paper_preview&mkt=zh-cn
Xing Guangfu, Yang Zhuliang, Mao Jianren, Shu Liangshu, Shen Jialin, Sun Qianghui, Chen Rong, Tao Kuiyuan. 2002. Characteristics of early Jurassic ingnerous rocks on the continental margin of southeastern China and their tectonic significanse[J]. Geological Bulletin of China, 21(7):384-391(in Chinese with English abstract).
Xing Guangfu, Chen Rong, Yang Zhuliang, Shen Jialin, Zhou Yuzhang, Yu Minggang, Li Longming. 2004. Discussion on Mesozoic volcanic rocks in Fujian, Zhejiang and Jiangxi Provinces-Comments on Zheng et al.'s paper[J]. Geology of Fujian, 23(4):233-239(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fujdz200404010
Xing Guangfu, Lu Qingdi, Chen Rong, Zhang Zhengyi, Nie Tongchun, Li Longming, Huang Jialong, Lin Min. 2008. Study on the ending time of Late Mesozoic tectonic regime transition in South China-Comparing to the Yanshan area in North China[J]. Acta Geologica Sinica, 82(4):451-463(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dizhixb200804003
Yuan Zhongxin, Wu Liangshi, Zhang Zongqing, Ye Xiaojiang. 1991.The Sm-Nd, Rb-Sr isotopic age-dating of Mayuan Group in northern Fujian[J]. Acta Petrologica et Mineralogica, 10(2):127-132(in Chinese with English abstract).
Yu Yunwen, Xu Butai. 1999. Stratigraphical sequence and geochronology of the Upper Mesozoic volcano-sedimentary rock series in Zhejiang[J]. Journal of Stratigraphy, 23(2):136-145(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DCXZ902.007.htm
Yu Yunwen, Xu Butai, Chen Jiangfeng, Dong Chuanwan.2001. Nd isotopic systematics of the Late Mesozoic volcanic rocks from southeasterm Zhejiang Province, China:Implications for stratigraphic study[J]. GeologicaI JournaI of China Universities, 7(1):62-69(in Chinese with English abstract). http://cn.bing.com/academic/profile?id=5d9c13cd2103129ac32e04d8f329b022&encoded=0&v=paper_preview&mkt=zh-cn
Zhang H F, Gao S, Zhong Z Q, Zhang B R, Zhang L, Hu S H. 2002.Geochemical and Sr-Nd-Pb isotopic compositions of Cretaceous granitioids:Constraints on tectonic framework and crustal structure of the Dabieshan ultrahigh-pressure metamorphic belt, China[J]. Chemical Geology, 186(3/4):281-299.
Zhang Qi, Wang Yan, Li Chengdong, Wang Yuanlong, Jin Weijun, Jia Xiuqin.2006. Granite classification on the basis of Sr and Yb contents and ite implications[J]. Acta Prtrologica Sinica, 22(9):2249-2269(in Chinese with English abstract). http://cn.bing.com/academic/profile?id=889afca2e396020a7abd0b34f19953f7&encoded=0&v=paper_preview&mkt=zh-cn
Zhou Xinmin, Sun Tao, Shen Weizhou, Shu Liangshu, Niu Yaoling.2006. Petrogenesis of Mesozoic granitoids and volcanicrocks in South China:A response to tectonic evolution[J]. Episodes, 29(1):26-33. doi: 10.18814/epiiugs/2006/v29i1/004
Zhang Jiajing, Shi Guanghai, Tong Guaisheng, Zhang Zhiyu, Liu Hai, Wu Rongshi, Chen Lei. 2009. Geochemistry and geochronology of copper and polymetal-bearing volcanic rocks of the Erhuling Formation in Xujiadun, Zhejiang Province[J]. Acta Geologica Sinica, 83(6):791-799(in Chinese with English abstract). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dizhixb200906005
Zhang Qi, Jun Weijun, Li Chengdong, Wang Yuanlong. 2010.Revisiting the new classification of granitic rocks based on wholerock Sr and Yb contents:Index[J]. Acta Petrologica Sinica, 26(4):985-1015(in Chinese with English abstract).
Zheng Keli. 2012. New researching development of Late Mesozoicvolcano formation in Fujian Province[J]. Geology of Fujian, 31(4):325-335(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-FJDZ201204006.htm
包超民, 王孔忠. 1999.浙江中生代火山-沉积地层和构造运动[J].中国区域地质, 18(2):201-204. doi: 10.3969/j.issn.1671-2552.1999.02.014 陈江峰, 郭新生, 汤加富, 周泰禧. 1999.中国东南地壳增长与Nd同位素模式年龄[J].南京大学学报(自然科学), 35(6):649-658. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199901025350 陈培荣, 孔兴功, 王银喜, 倪琦生, 章邦桐, 凌洪飞. 1999.赣南燕山早期双峰式火山-侵入杂岩的Rb-Sr同位素定年及意义[J].高校地质学报, 5(4):378-383. http://www.cnki.com.cn/Article/CJFDTotal-GXDX199904002.htm 陈丕基. 2000.中国陆相侏罗、白垩系划分对比述评[J].地层学杂志, 24(2):114-119. doi: 10.3969/j.issn.0253-4959.2000.02.006 丁保良, 蓝善先, 汪迎平. 1987.试论浙闽赣地区非海相侏罗、白垩纪界线[J].中国地质科学院南京地质矿产研究所所刊, 8(3):43-55. http://www.cnki.com.cn/Article/CJFDTotal-HSDZ198703004.htm 董树文, 张岳桥, 陈宣华, 龙长兴, 王涛, 杨振宇, 胡建民. 2008.晚侏罗世东亚多向汇聚构造体系的形成与变形特征[J].地球学报, 29(3):306-317. doi: 10.3321/j.issn:1006-3021.2008.03.005 戴洪文. 2011.福建福安潘洋地区晚侏罗世火山岩与侵入岩的关系探讨[J].福建地质, 30(4):335-340. doi: 10.3969/j.issn.1001-3970.2011.04.010 段政, 邢光福, 余明刚, 赵希林, 靳国栋, 陈志洪. 2013.浙闽边界区晚中生代火山作用时序与过程分析[J].地质论评, 59(3):454-469. doi: 10.3969/j.issn.0371-5736.2013.03.006 顾知微. 2002.关于费尔干蚌Ferganoconcha[J].古生物学报, 41(1):152-158. doi: 10.3969/j.issn.0001-6616.2002.01.016 顾知微. 2005.论闽浙运动[J].地层学杂志, 29(1):1-6. doi: 10.3969/j.issn.0253-4959.2005.01.001 蒋维三. 1995.试论浙江白垩系的顶、底界[J].浙江地质, 11(1):15-23. http://www.cnki.com.cn/Article/CJFDTotal-ZJDZ199501001.htm 刘春, 闫峻, 宋传中, 李全忠, 彭戈, 史磊, 刘晓强. 2012.长江中下游繁昌盆地火山岩年代学和地球化学:岩石成因和地质意义[J].岩石学报, 28(10):3228-3240. http://www.cnki.com.cn/Article/CJFDTotal-YSXB201210015.htm 马武平. 1994.论浙江中生代晚期地层划分[J].地层学杂志, 18(2):91-101. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400025465 聂永红, 李曙光, 葛宁洁, 刘德良, 胡雄健. 1995.浙西南北园八都群石榴斜长角闪岩的变质特征及年代[J].矿物学报, 15(4):442-450. doi: 10.3321/j.issn:1000-4734.1995.04.014 宋彪, 张玉海, 万渝生, 简平.2002.锆石SHRIMP样品靶制作、年龄测定及有关现象讨论[J].地质评论, 48(增刊):26-30. http://www.cnki.com.cn/Article/CJFDTotal-DZLP2002S1006.htm 舒良树, 周新民. 2002.中国东南部晚中生代构造作用[J].地质论评, 48(3):249-260. doi: 10.3321/j.issn:0371-5736.2002.03.004 陶奎元, 邢光福, 杨祝良, 毛建仁, 赵宇, 许乃政. 2000.浙江中生代火山岩时代厘定和问题讨论兼评Lapierre等关于浙江中生代火山活动时代的论述[J].地质论评, 46(1):14-21. doi: 10.3321/j.issn:0371-5736.2000.01.003 汪庆华. 2001.试论浙江建德群和磨石山群时代[J].火山地质与矿产, 22(3):63-169. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hsdzykc200103002 皇家龙, 卢清地, 张正义, 聂童春, 周国武, 林敏, 郑平. 2008.福建仙游园庄地区南园组新层型剖面的建立及时代的重新厘定[J].地质通报, 27(6):785-792. doi: 10.3969/j.issn.1671-2552.2008.06.006 王对兴, 高万里, 李春麟, 王宗秀, 赵志丹. 2015.浙中地区晚侏罗世花岗斑岩LA-ICP-MS锆石U-Pb年龄、地球化学特征及其地质意义[J].地质学报, 42(6):1684-1699. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgdizhi201506003 王加恩, 王振, 刘远栋, 梁河. 2016.浙西建德群老村组火山岩地质时代及构造环境研究[J].地层学杂志, 40(3):262-271. http://xueshu.baidu.com/usercenter/paper/show?paperid=7952553ede0c1fa742fde44a957c68ec&site=xueshu_se&hitarticle=1 徐步台, 俞云文. 1997.浙江中生代火山岩同位素年龄及地层时代归属[J].地球学报, 18(增刊):117-119. http://www.cqvip.com/qk/98325X/1997A00/2883749.html 邢光福, 杨祝良, 毛建仁, 舒良树, 沈加林, 孙强辉, 陈荣, 陶奎元. 2002.东南大陆边缘早侏罗世火成岩特征及其构造意义[J].地质通报, 21(7):384-391. doi: 10.3969/j.issn.1671-2552.2002.07.004 邢光福, 陈荣, 杨祝良, 沈加林, 周宇章, 余明刚, 李龙明. 2004.闽浙赣中生代火山岩研究问题浅议——兼与郑承意等商榷[J].福建地质, 23(4):233-239. doi: 10.3969/j.issn.1001-3970.2004.04.010 邢光福, 卢清地, 陈荣, 张正义, 聂童春, 李龙明, 黄家龙, 林敏. 2008.华南晚中生代构造体制转折结束时限研究——兼与华北燕山地区对比[J].地质学报, 82(4):451-463. doi: 10.3321/j.issn:0001-5717.2008.04.003 袁忠信, 吴良士, 张宗清, 叶笑江. 1991.闽北麻源群Sm-Nd、Rb-Sr同位素年龄研究[J].岩石矿物学杂志, 10(2):127-132. http://www.cnki.com.cn/Article/CJFDTotal-YSKW199102003.htm 俞云文, 徐步台. 1999.浙江中生代晚期火山-沉积岩系层序和时代[J].地层学杂, 23(2):136-145. doi: 10.3969/j.issn.0253-4959.1999.02.008 俞云文, 徐步台, 陈江峰, 董传万. 2001.浙东南中生代晚期火山岩Nd同位素组成及其地层学意义[J].高校地质学报, 7(1):62-69. doi: 10.3969/j.issn.1006-7493.2001.01.007 浙江省地质矿产局. 1989.浙江省区域地质志[M].北京:地质出版社. 浙江省地质矿产局. 1996.浙江省岩石地层[M].武汉:中国地质大学出版社. 张旗, 王焰, 李承东, 王元龙, 金惟俊, 贾秀勤. 2006.花岗岩的Sr-Yb分类及其地质意义[J].岩石学报, 22(9):2249-2269. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ysxb98200609001 张家菁, 施光海, 童贵生, 张智宇, 刘海, 吴荣土, 陈磊. 2009.浙江徐家墩鹅湖岭组含铜多金属矿火山岩的地球化学与年代学[J].地质学报, 83(6):791-799. doi: 10.3321/j.issn:0001-5717.2009.06.005 张旗, 金惟俊, 李承东, 王元龙. 2010.再论花岗岩按照Sr-Yb的分类标志[J].岩石学报, 26(4):985-1015. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ysxb98201004001 郑克丽. 2012.福建晚中生代火山地层研究新进展[J].福建地质, 31(4):325-335. doi: 10.3969/j.issn.1001-3970.2012.04.005
计量
- 文章访问数: 2603
- HTML全文浏览量: 477
- PDF下载量: 4119
