汶川MS8.0级地震余震分布及周边区域P波三维速度结构研究

利用川滇地区长期积累的地震走时观测资料和汶川地震余震观测资料对汶川地震震源区及周边区域地壳和上地幔P波三维速度结构进行了研究.结果表明,浅部P波速度分布与地表地质之间具有很好的对应关系.龙门山断裂带在20 km以上深度表现为高速异常带,彭灌杂岩体和宝兴杂岩体为局部高速异常区.龙门山断裂带中上地壳的局部高速异常体对汶川地震的余震分布具有明显的控制作用.在余震带南端,余震全部发生在与宝兴杂岩体对应的高速异常体的东北侧;在余震带的中段,与彭灌杂岩体对应的高速异常体在一定程度上控制了余震的分布;在余震带的东北端,宁强-勉县一带的高速异常体可能阻止了余震进一步向东北扩展.龙门山断裂带中上地壳的P波高速异常表明介质具有相对较高的强度,在青藏高原物质向东挤出过程中起到了较强的阻挡作用,有利于深部能量积累.在30 km深度之下,扬子地块具有明显的高速特征,其前缘随深度增加向青藏高原方向扩展,在下地壳和上地幔顶部已达到龙门山断裂带以西.     

龙门山中南段地壳上地幔三维密度结构

基于高精度布格重力异常资料,以川滇地区P波速度三维层析成像结果为约束建立初始模型,采用预优共轭梯度(Preconditional Conjugate Gradiem,PCG)反演方法得到了龙门山断裂带中南段的地壳上地幔(深度范围0～65km)三维密度结构(网格间距为10km(横向)×10km(纵向)×5km(深度))。密度成像结果表明:龙门山断裂带中南段两侧地壳密度结构存在明显差异,四川盆地有约10km厚的低密度沉积层,松潘-甘孜块体因沉积层较薄,且部分地区有基岩出露,上地壳表现为高密度结构;松潘-甘孜块体中、下地壳有大范围低密度层分布,介质强度明显低于高密度的四川盆地,青藏高原东移物质受到四川盆地阻挡后更易于在低密度的一侧发生挤压形变及隆升,从而形成龙门山逆冲推覆构造带;龙门山断裂带内部在地壳结构上具有明显的分段特征,表现为沿着龙门山断裂带地壳密度变化不连续,以汶川地震和芦山地震震中为界,形成多个高、低密度异常区;同时,结合地震精定位结果分析,汶川地震及其余震多分布于壳内中央断裂带西侧高密度体内,芦山地震及其余震则集中在地壳密度变化梯级带附近并偏向高密度体一侧。四川盆地下地壳密度较高,其前缘随深度增加向青藏高原方向扩展,在上地幔顶部接近龙门山断裂带以西。松潘-甘孜块体中、下地壳虽然有一定规模的低密度体分布,但其连通性差,在平面上多形成局部低密度异常区,是否存在下地壳流仍无法给出明确的证据。     

汶川Ms8．0地震余震序列重新定位及其地震构造研究

综合利用川西流动地震台阵观测数据和震后应急地震观测台站的震相数据,采用双差地震定位方法对汶川地震的余震序列进行了精确重新定位,并对汶川地震的地震构造进行了深入研究.其结果显示,汶川地震序列从彭灌杂岩南缘开始破裂,主震及其余震破裂带长约350 km,在大部分区域宽约20～30 km,其宽度和空间形态沿破裂带显示了强烈的分段和非均匀特征.坚硬的彭灌杂岩对余震的非均匀性分布和汶川地震复杂的破裂过程起到了重要的控制作用.以松潘-甘孜地块中地壳低速层顶部为底边界,余震主要分布在4～24 km深度范围内的龙门山东缘上地壳高速层内.余震深度分布削面清晰地显示了映秀-北川断裂和灌县江油断裂以及汶川-茂汶断裂在20～22 km深度合并为剪切带的特征.小鱼洞到理县方向存在一条长度超过60 km的垂直于龙门山走向的余震分布条带,综合震源机制解和地震破裂过程的研究结果,我们推测,这是坚硬的彭灌杂岩体底部在长期应力积累作用下发生破裂的反映,并成为汶川地震释放出巨大能量的主要原因.     

龙门山断裂带南段地壳一维P波速度结构

基于2009年1月1日至2013年5月6日四川地震台网、重庆地震台网记录的龙门山断裂带南段587个地震的5 012个P波到时数据,利用最小一维速度模型方法反演了龙门山断裂带南段地壳一维P波速度模型及台站校正值,并将其应用到龙门山地区地震重定位中。结果表明,台站校正值表征出龙门山断裂带南段地表速度结构的横向不均匀性,青藏高原的彭灌杂岩体及宝兴杂岩体在近地表表现为高速异常,而四川盆地的第四纪沉积表现为低速异常。重定位后地震震源在北西向的剖面上呈明显的条带状并向北西倾斜,该地震带与宝兴杂岩下方的滑脱带延伸趋势一致。此外,该地震带上方分布着一条反冲地震带,两地震带呈"y"型分布,这可能是宝兴杂岩上方的岩层为调节逆冲过程受阻而产生的反冲运动所致。     

川滇地区地壳上地幔三维速度结构研究

根据云南和四川地震台网174个台站记录的4625个区域地震初至Ｐ波和Ｓ波走时资料，并结合其它深部地球物理资料，确定了川滇地区地壳上地幔三维速度结构．在上地壳速度异常分布中，四川盆地为正异常，川西高原为负异常，龙门山断裂带为正、负异常的边界．龙门山断裂、鲜水河断裂以及红河断裂等，在下地壳和上地幔的速度异常中仍显示出构造分界特征，说明它们可能穿透了莫霍界面．腾冲火山区和攀西构造带在50km深度上呈现负速度异常，与上地幔温度和物质组成的差异相联系．川滇地区地壳结构的总体特征是:地壳和上地幔的低平均速度，地壳厚度变化剧烈，地壳和（或）上地幔存在高导层、高热流值．这些同印度板块与欧亚板块碰撞的构造背景有关．川滇菱形块体在地壳内总体上为正常或正异常速度，而其边界的深大走滑断裂存在负速度异常，它有助于地壳块体沿断裂的侧向挤出．在主要的地震带上，中下地壳的负速度异常与地震活动性相关．多数强烈地震发生在具有正速度异常或正常速度分布的上中地壳深度上，而其下方则通常是负速度异常带．      

小江断裂带周边地区三维P波速度结构及其构造意义

作为青藏高原的东南边界,小江断裂带在高原物质的侧向逃逸中发挥着重要的作用.本文利用流动地震台阵及固定台站的走时观测资料,对小江断裂带及周边区域的壳幔三维P波速度结构进行了研究.结果表明,在中上地壳,小江断裂带内部主要为低速异常,其东侧主要为高速异常.在中下地壳,小江断裂带中部为低速异常,北部和南部主要为高速异常,其中北部的高速异常可延伸到地表附近,南部的高速异常可一直延伸到上地幔.我们推测,小江断裂带中部的低速异常与深部热作用有关;北部的高速异常可能是晚古生代地幔柱活动导致大量基性和超基性幔源物质侵入地壳引起的,它的存在对青藏高原物质向南逃逸起到了一定的阻挡作用,可能是导致川滇活动块体北部次级块体快速抬升的重要因素;南部顶界面向北倾斜的高速异常体对川滇活动块体向南滑移起到了进一步的阻挡作用,导致其上覆的中上地壳低速异常区发生较强的变形和强烈的地震活动,同时在上地幔深度范围起到了稳定的作用,使其南部区域的介质受青藏高原物质向南挤出的影响明显减小.     

基于密集台阵研究龙门山断裂带南段地震空段的地震活动性

2016年12月—2018年4月间布设于汶川、芦山地震之间地震空段的密集监测台阵(LmsSGA)提供了密集的观测数据.通过拾取地震走时、初始定位,计算地方震级,得到了完备性震级为0级的地震目录.更加完备的地震目录为地震空段及周围地震活动的时空分布特征和孕震风险性评估提供了丰富的信息.重定位结果显示地震主要集中于龙门山断裂带深度为5～20km的孕震层内.地震活动频繁的汶川、芦山主震区,震源的空间分布模式与其早期余震相似,说明两次大地震的区域仍处于缓慢的应力调整阶段.青藏高原物质东向挤出受宝兴、彭灌杂岩阻挡,在两个杂岩体西北侧地震活动频繁.地震活动性分布显示汶川—茂县、映秀—北川断裂上存在一个清晰的长约30km,宽约20km的地震活动"空白"区域,与其下方因部分熔融而产生的低速体分布一致,我们推测熔融体的加温作用是导致空段内极低的地震活动性的主要原因.监测时段内仍观测到降雨变化率和地震数量呈反相关关系,再次证实了汶川—芦山地震间地震空段及邻区内季节性降雨对地震活动性存在一定调节作用.综合分析S波速度模型、历史强震活动及b值,我们推断地震空段东部的彭灌断裂中段及周围部分隐伏断层存在发生强震的风险.     

龙门山断裂带精细速度结构的双差层析成像研究

利用川西流动地震台阵、汶川地震震后应急台网记录到的P波到时资料,对2008年5月至2008年10月期间发生的汶川地震余震序列应用双差层析成像方法进行了地震震源和三维P波速度结构的联合反演.结果显示,联合反演获得的地震重定位结果与基于一维地壳参考模型的双差定位方法结果相近;研究区15 km以上速度结构与地表断裂分布密切相关,20 km以下深度呈现北东向和北西向交错结构.汶川地震破裂带南段龙门山断裂带之间上地壳呈现高速异常,速度结构的非均匀变化是控制余震分布和主震破裂传播的主要因素;联合反演结果给出了小鱼洞-理县方向存在隐伏断裂的速度结构证据,同时发现,破裂带北东段可能沿新发断裂扩展;结果确认了汶川地震起始段的高角度逆冲断裂特征,也确认了前山断裂和中央断裂在约20 km深度合并到脆韧转换带的特征.     

龙门山和安宁河断裂带周边地区的地壳三维P波速度结构

利用2008—2021年区域地震台网及部分流动地震台阵的近震走时观测资料，采用双差层析成像方法获得了龙门山断裂带和安宁河断裂带周边地区的高分辨率地壳三维P波速度结构。结果显示：四川盆地浅部呈明显的低速异常，在前陆盆地低速异常可一直延伸至约15 km深，推测与该地区接收大量来自青藏高原东缘的沉积物有关；在龙门山断裂带附近，中上地壳存在一些与断裂带平行的高速异常体，它们较好地揭示了该地区穹隆体和逆冲推覆体深浅部的分布特征，例如，位于龙门山断裂带中部的彭灌穹隆体下延深度可达15 km左右，而与之相邻的雪隆山穹隆体下延深度不超过5 km。在安宁河断裂带与大凉山断裂带之间存在震源深度超过20 km的微震活动带，其分布形态与P波高速异常的分布形态基本一致，推测高速异常体的存在可能导致地壳内部脆韧性转换带深度的增加，并且受安宁河与大凉山交会部位强烈变形的影响，该地区出现深度较大的地震活动。上述结果为认识该地区的深部地质构造和地震活动机制提供了新的依据。     

汶川MS8.0地震余震序列重新定位及其地震构造研究

综合利用川西流动地震台阵观测数据和震后应急地震观测台站的震相数据,采用双差地震定位方法对汶川地震的余震序列进行了精确重新定位,并对汶川地震的地震构造进行了深入研究.其结果显示,汶川地震序列从彭灌杂岩南缘开始破裂,主震及其余震破裂带长约350 km,在大部分区域宽约20~30 km,其宽度和空间形态沿破裂带显示了强烈的分段和非均匀特征.坚硬的彭灌杂岩对余震的非均匀性分布和汶川地震复杂的破裂过程起到了重要的控制作用.以松潘—甘孜地块中地壳低速层顶部为底边界,余震主要分布在4~24 km深度范围内的龙门山东缘上地壳高速层内.余震深度分布剖面清晰地显示了映秀—北川断裂和灌县—江油断裂以及汶川—茂汶断裂在20~22 km深度合并为剪切带的特征.小鱼洞到理县方向存在一条长度超过60 km的垂直于龙门山走向的余震分布条带,综合震源机制解和地震破裂过程的研究结果,我们推测,这是坚硬的彭灌杂岩体底部在长期应力积累作用下发生破裂的反映,并成为汶川地震释放出巨大能量的主要原因.     

Temporal variations and scaling of streamflow and baseflow and their nitrate-nitrogen concentrations and loads

The patterns of temporal variations of precipitation (P), streamflow (SF) and baseflow (BF) as well as their nitrate-nitrogen (nitrate) concentrations (C) and loads (L) from a long-term record (28 years) in the Raccoon River, Iowa, were analyzed using variogram and spectral analyses. The daily P is random but scaling may exist in the daily SF and BF with a possible break point in the scaling at about 18 days and 45 days, respectively. The nitrate concentrations and loads are shown to have a half-year cycle while daily P, SF, and BF have a one-year cycle. Furthermore, there may be a low-frequency cycle of 6–8 years in C. The power spectra of C and L in both SF and BF exhibit fractal 1/f scaling with two characteristic frequencies of half-year and one-year, and are fitted well with the spectrum of the gamma distribution. The nitrate input to SF and BF at the Raccoon watershed seems likely to be a white noise process superimposed on another process with a half-year and one-year cycle.     

燕郊等测点迁移优化与地磁观测研究

为优化地磁观测条件，开展了燕郊、夏垫等测点迁移工作；按照测眯迁移原则与实施技术方案，完整地收集并整理了地质构造与地球物理等方面的基本资料；进行了野外实地勘察，磁场梯度的测量，确定了新测点，在新老测点上进行了较长时间的地磁场对比观测；应用多种方法分析研究了地磁对比观测资料，结果表明，新老测点与有关测点的地磁变化具有良好的一致性，并得到了新老测点之间的地磁数据的按点差。     

东、南洞庭湖的径流、泥沙特征及冲淤规律

通过实地调查并对1957年以来水文、泥沙观测资料做系统分析和计算,探讨东、南洞庭湖出、入湖水量、沙量的年际和年内变化特征,以及长江下荆江段裁弯对湖区径流和泥沙的影响。提出了湖区泥沙汛淤枯冲的变化规律及水位升降与湖区泥沙冲淤的关系;论证了丰、平、枯年湖区淤积严重,面积日益缩小对径流的调节作用正在减弱。     

GIS与防灾减灾

地理信息系统(GIS)在我国已得到广泛的应用。在简要介绍G1S的发展历程、特点和发展趋势的基础上，重点从防灾减灾的角度介绍了GIS在气象灾害、地震灾害、地质灾害等领域的应用成果。     

Assessing Groundwater Depletion and Dynamics Using GRACE and InSAR: Potential and Limitations

In the last decade, remote sensing of the temporal variation of ground level and gravity has improved our understanding of groundwater dynamics and storage. Mass changes are measured by GRACE (Gravity Recovery and Climate Experiment) satellites, whereas ground deformation is measured by processing synthetic aperture radar satellites data using the InSAR (Interferometry of Synthetic Aperture Radar) techniques. Both methods are complementary and offer different sensitivities to aquifer system processes. GRACE is sensitive to mass changes over large spatial scales (more than 100,000 km²). As such, it fails in providing groundwater storage change estimates at local or regional scales relevant to most aquifer systems, and at which most groundwater management schemes are applied. However, InSAR measures ground displacement due to aquifer response to fluid‐pressure changes. InSAR applications to groundwater depletion assessments are limited to aquifer systems susceptible to measurable deformation. Furthermore, the inversion of InSAR‐derived displacement maps into volume of depleted groundwater storage (both reversible and largely irreversible) is confounded by vertical and horizontal variability of sediment compressibility. During the last decade, both techniques have shown increasing interest in the scientific community to complement available in situ observations where they are insufficient. In this review, we present the theoretical and conceptual bases of each method, and present idealized scenarios to highlight the potential benefits and challenges of combining these techniques to remotely assess groundwater storage changes and other aspects of the dynamics of aquifer systems.     

Correlation of foreshocks and aftershocks and asperities

A close correlation in spatial distribution of local seismic activity and energy release patterns before and after the 1979 Petatlan, Mexico earthquake suggests heterogeneity within the fault plane of this major low-angle thrust event associated with subduction along the Middle America Trench. A simple two-asperity model is proposed to account for the complexity. Foreshocks and aftershocks of the neighboring 1981 Playa Azul earthquake showed a similar pattern. As both events occurred at the junction of the Orozco Fracture Zone and the Middle America Trench, we speculate that the observed complex fault plane is caused by subduction of the rugged ocean floor of the Orozco Fracture Zone. Short-term precursory seismicity prior to the Petatlan earthquake can be explained by using the asperity model and migration of a slip front from the south-east to the north-west across the main shock source region.