晚喜山期以来四川盆地构造-热演化模拟

四川盆地位于扬子板块西缘,是我国重要的含油气盆地之一.25 Ma以来的晚喜山期是四川盆地构造-热演化的重要时期,此时,盆地大部分区域受到挤压处于隆升剥蚀的构造动力学环境.本文采用有限元数值模拟方法,基于二维瞬态热传导(含平流项)的基本方程,并引入修正的Airy均衡理论模型,通过覆盖全盆地的八条剖面模拟研究了晚喜山期以来四川盆地的构造-热演化特征,且利用现今大地热流对模拟结果进行了有效约束.模拟结果显示在晚喜山期(~25 Ma)川中地区地表热流较高,为60~64 mW/m²;川西南地区次之,为60~62 mW/m²;川东北地区最低,为50~54 mW/m².该期基底热流,也是川中隆起区热流高,川东北强烈剥蚀区热流低.热流的空间分布特征揭示了四川盆地深部动力学机制.四川盆地晚喜山期以来,抬升剥蚀作用降低了其地表热流和基底热流,其降低幅度与对应的剥蚀速率相关,即剥蚀速率越大,这种降低作用越明显.     

利用气象地温资料反演大地热流

利用气象台记录的深度为0.8m,1.6m和3.2m的月平均地温资料,进行Fourier频谱分析,得到定常波和子谐波（变化波）两部分.同时取地表以下为水平分层的均匀双层介质,建立数学模型,根据线性热传导和最优化理论对定常波以及子谐波的振幅谱和相位谱进行联合反演,计算气象台地表附近的定常地温梯度以及地表土壤的热扩散系数等物性参数,并据此计算该地的大地热流值.初步计算结果表明,上述方法得到的某些地区的大地热流值与传统的通过钻孔测量法得到的大地热流值符合得相当好.     

重力球面外部校正技术

山区重力勘探中，重力外部校正是重力勘探料处理和解释的前提，其内容包括地形校正、中间层校正等。本文在分析常规外部校正方法的基础上，提出了一套适合于山区重力勘探的球面地形校正及有限球壳中间层校正方法，并严格推出了球面地形校正公式。对中国南方ZJJ地区重力资料处理结果表明，该方法提高了地形校正的精度，使山区重力资料品质达到或接近平原区的水平。     

南海珠江口盆地钻井BHT温度校正及现今地温场特征

温度是盆地现今地温场研究最重要的基础参数之一,针对钻井井底温度(BHT)一般低于地层真实温度的情况,本文对珠江口盆地16口钻井的BHT数据进行了校正,并根据校正后的温度数据采用最小二乘法新增计算了16个地温梯度及大地热流数据.结合前人发表数据,本文绘制了珠江口盆地最新地温梯度等值线图及大地热流等值线图.除此之外,系统收集整理了盆地钻井大地热流、生热率以及热导率等参数,利用一维热传导方程计算了盆地1~5km深度处的温度,并绘制了盆地深部温度等值线图.在此基础上,分析了盆地现今地温场特征及其成因,以期为珠江口盆地地球动力学背景研究及油气资源评价与勘探提供地热学参数.     

青藏高原东南缘热流估算及与地震活动相关性分析

青藏高原东南缘地区内部构造运动强烈,是地热资源发育与地震事件频发的活动地区.大地热流记录了发生在地球深部各种作用过程的热学信息,可以作为地质构造活动和地震活动研究的有效约束,但是大范围的热流数据测量很难实现,因此,本文根据居里面深度结合放射性元素分布等计算了青藏高原东南缘的大地热流分布.首先,通过地表放射性元素的分布计算出地表产热量的分布,然后,利用相关地热参数之间的关系迭代计算出该地区地壳上下层的热导率分布,最终估算出地表热流及地下不同深度处热流值的分布.本文结果表明:(1)青藏高原东南缘的大地热流位于44～108mW·m~(-2)之间,平均75mW·m~(-2),符合研究地区西南高、东北低的背景趋势,地壳内部热流值随深度的增加而降低.大部分地区地表热流异常与实际地热带分布相吻合,如川西、藏东南与滇西地区等地为地热高值区,川东和楚雄等地为热流低值区.(2)结合其他地球物理探测结果,总结了地壳内部热流与地震事件的联系:在地热梯度带地区,当两侧地层在一定深度范围内存在明显物性差异时,地震事件高发.     

中国大陆地区大地热流数据汇编(第三版)

大地热流是地球内部热作用过程最直接的地表显示，其中蕴涵着丰富的地质、地球物理和地球动力学信息.我国大陆地区已定期进行了6次热流汇编，其中第一、二次汇编的热流数据已正式公布(第一、二版)，此后4次汇编仅发表了统计结果.而未再公布汇编数据.本文基于第六次热流数据汇编，将自第二版以来新增的450个热流数据汇编成中国大陆地区大地热流数据第三版.新版热流数据的统计结果表明，中国大陆地区实测热流值变化于23-319mW/m2，平均63±24.2mW/m2;剔除与地表热异常相关的数据后，热流值变化范围为30-140mW/m2，平均61±15.5mW/m2.     

复杂山区光学遥感反射率计算模型

提出一种复杂山区光学遥感反射率计算模型,以有效消除地形影响,计算地表反射率．通过推导地形影响下的方向反射,结合辐射传输方程,考虑地表反射目标接收的太阳直射辐射为方向．方向二向反射,反射目标接收的大气漫散射和周围地形背景反射辐射为半球．方向二向反射的综合过程,发展了基于方向反射的地形影响消除和地表反射率计算模型．利用GOMS模型模拟数据及真实光学遥感数据对模型应用能力进行了评价,结果表明此模型具有有效的地形消除能力以及大气校正能力,模型的建立为后续山区定量遥感的应用提供技术支持．     

云南深部热流研究

阐述了云南大地热流分布特征，并从地壳结构入手探讨了地壳、上地幔的生热率垂向分布，进而研究了云南的深部热流结构．研究结果表明，垂向上，上地壳所产生的热流量最大，中地壳次之，下地壳最小；横向上，各构造单元深部热流结构差别较大，全区地幔热流变化范围是１９～５７ｍＷ／ｍ２．     

利用气象地温资料计算区域地球表面平均热流

利用最新收集到的全国范围内177个气象测点,时间跨度最长达51年的浅层地表长时间温度序列,通过频谱分析及数值求解的方法去除温度序列中受大气影响的高频变化部分,得到反映来自地球内部热信息的稳定温度梯度,进而求得地表平均热流,并分析其在多大程度上能够反映大地热流场,希望在丰富大地热流数据库,及对我国地震、旱涝灾害的预测方面给予启示.     

磷灰石裂变径迹分析在恢复盆地沉降抬升史中的应用──以鄂尔多斯盆地为例

对鄂尔多斯盆地磷灰石裂变径迹资料深入分析表明．最迟２３Ｍａ以来盆地发生了一期由于快速抬升剥蚀引起的冷却事件．盆地东部以９５ｍ／Ｍａ的速率抬升，造成约２０００ｍ的剥蚀量；而盆地西部则以５６ｍ／Ｍａ的速率抬升，导致了约１０００ｍ的剥蚀量．盆地东、西部的差异抬升剥蚀导致了盆地现今微微西倾的构造面貌．这一抬升剥蚀事件是印度板块与欧亚板块碰撞引起亚洲构造运动形式以挤压为主，转换为中新世以来以地壳增厚为主的结果．Ｋ－Ａｒ年龄和镜质体反射率资料分析表明，盆地在１７０－１６０Ｍａ（中侏罗末）曾发生一期热事件，使古地温梯度达５７℃／ｋｍ，古热流值达９６－１０９ｍｗ／ｍ^２．     

Formation,uplift and dissection of planation surfaces at passive continental margins – a new approach

The usefulness of large‐scale, low‐relief, high‐level landscapes as markers of uplift events has become a subject of disagreement among geomorphologists. We argue that the formation of low‐relief surfaces over areas of large extent and cutting across bedrock of different age and resistance must have been guided by distinct base levels. In the absence of other options the most likely base level is sea level. We have analysed West Greenland landscapes in a recent study by combining the cooling history from apatite fission‐track analysis (AFTA) data with the denudation history from landscape analysis and the stratigraphic record. An important difference between our approach and that of classical geomorphology is that we now have the ability to document when thick sections of rocks have been deposited and then removed. The present‐day high‐level plateau in West Greenland is the remnant of a planation surface that was formed by denudation that lasted c. 20 million years during which up to 1 km of cover was removed after maximum burial at the Eocene–Oligocene transition. Here we present additional AFTA data to show that the planation surface is the end‐product of Cenozoic denudation even in basement areas and argue that Phanerozoic sediments – most likely of Cretaceous–Palaeogene age – must have been present prior to denudation. The planation surface was offset by reactivated faults and uplifted to present‐day altitudes of up to 2 km. The uplift occurred in two late Neogene phases that caused incision of valleys below the planation surface and their subsequent uplift. We therefore find that the elevated and deeply dissected plateau is evidence of episodic post‐rift uplift that took place millions of years after cessation of sea‐floor spreading west of Greenland. We suggest that other margins with similar morphology may also be characterized by episodic post‐rift uplift unrelated to the processes of rifting and continental separation, rather than being permanently uplifted since the time of rifting, as is commonly assumed. Copyright © 2009 John Wiley & Sons, Ltd.     

Geomorphological development of the eastern margin of the Australian Craton

The denudation chronologies of five uplands of contrasting geological structure located at the eastern margin of the Australian Craton are examined. They are the Isa Highlands, MacDonnell Ranges and Flinders Ranges (fold mountains, orogenic setting); the Arcoona Plateau (dissected plateau, platform setting); and the Gawler Ranges (massif of old silicic volcanics, cratonic setting). In each, surfaces of Mesozoic age, many of them exhumed and of pre-Cretaceous age, are preserved. Each also appears to have been uplifted recurrently. Each was either overwhelmed or bordered by the Early Cretaceous (Neocomian-Aptian) sea. Tectonism associated with the break-up of Gondwana probably allowed this important marine transgression. Thereafter, thalassostatic and erosional/depositional isostatic responses have maintained the uplands as uplands and the intervening basins as negative topographic and structural units; the pattern of topography has been constant for the last 60–100 Ma.     

Numerical calculations of the perturbation of normal heat flow by topographic and local conductivity structure

A numerical method is used to investigate the effect of topographic and local thermal conductivity anomalies on near-surface heat flow for two-dimensional models. Heat flow associated with a sloping topographic structure is calculated. Also, the effects of a fault structure associated with the sloping topography are considered. Vertical and horizontal heat flow components are calculated alone; the surface of the earth as well as throughout the whole region of interest. The results indicate that surface heat flow is substantially affected by topographic relief and the horizontal heat flow component associated with topographic features can be large. Also, regional heat flow is greatly perturbed by local thermal conductivity anomalies and the effect of topographic features may be considerably modified by the subsurface structure.     

Denudation intensity and control in the Chinese Tian Shan: new constraints from mass balance on catchment‐alluvial fan systems

Tectonics and climate are usually seen as the main controlling factors of denudation rates, which seem to rise with the tectonic activity and to decrease when the climate becomes drier. However, the low denudation rates observed in semi‐arid to arid contexts are generally measured on orogenic plateaus where the respective influence of the flat relief and the dry climate cannot really be unravelled. The Chinese Tian Shan was chosen as a case study. In the northern piedmont of this mountain range, a series of well‐preserved Quaternary alluvial fans offer the opportunity to perform a mass balance study at the scale of several catchment areas and several hundreds of thousands of years. Based on a geometrical reconstruction of these fans, the volumes of sediments exported out of 10 drainage basins during the Middle–Late Pleistocene (from ～300 to ～12 kyr) and the Holocene (from ～12 kyr to present) have been estimated. From these volumes, an average denudation rate of ～135 m/Myr was determined in the Tian Shan Range for the last 300 kyr. In agreement with other mass balances performed in the same area, the typical denudation intensity of the northern Tian Shan is thus of a few hundred meters per million years at most, regardless of the space and time scales considered. From a comparison with denudation rates in other mountain ranges throughout the world, we suggest that a dry climate can dramatically limit the denudation intensity even in active orogenic systems with a high topographic gradient like the Tian Shan. As a result, the time required to reach equilibrium between denudation and rock uplift in these systems could be extremely long (i.e. of more than several million years). Copyright © 2016 John Wiley & Sons, Ltd.     

Mathematical simulation of geotemperature and heat flow patterns

Geotemperature and heat flow patterns in a large-scale Meso-Cenozoic basin such as the North China Basin are strongly affected by the relief of the basement, and controlled by the contrast of thermal conductivity between basement rock and sedimentary cover. Usually, heat flow observed at the surface of a basement uplift is greater than that of a basement depression. Calculation revealed, that the ratio of the former and the latter is determined by the uplifted height (H) of the bed-rock roof of the basement and the thickness (h) of the sedimentary cover. The relief of the basement also disturbs the geotemperature and, hence, the heat flow patterns at shallow depth. Consequently, the more or less “uniform” one dimensional heat flow from the deep interior of the Earth becomes two dimensional at shallow depth with great lateral and vertical variations. The extent of the disturbed zone is also controlled by the contrast of the thermal conductivity between basement rock and sedimentary cover as well as the uplifted heigh (H) of the bed-rock roof of the basement. Numerical computation demonstrated that the disturbed depth (Ze) is usually about 3–6 times of the uplifted height (H) of a basement uplift.     

Effects of terrain smoothing on topographic shielding correction factors for cosmogenic nuclide‐derived estimates of basin‐averaged denudation rates

Estimation of spatially averaged denudation rates from cosmogenic nuclide concentrations in sediments depends on the surface production rates, the scaling methods of cosmic ray intensities, and the correction algorithms for skyline, snow and vegetation shielding used to calculate terrestrial cosmogenic nuclide production. While the calculation of surface nuclide production and application of latitude, altitude and palaeointensity scaling algorithms are subjects of active research, the importance of additional correction for shielding by topographic obstructions, snow and vegetation is the subject of ongoing debate. The derivation of an additional correction factor for skyline shielding for large areas is still problematic. One important issue that has yet to be addressed is the effect of the accuracy and resolution of terrain representation by a digital elevation model (DEM) on topographic shielding correction factors. Topographic metrics scale with the resolution of the elevation data, and terrain smoothing has a potentially large effect on the correction of terrestrial cosmogenic nuclide production rates for skyline shielding. For rough, high‐relief landscapes, the effect of terrain smoothing can easily exceed analytical errors, and should be taken into account. Here we demonstrate the effect of terrain smoothing on topographic shielding correction factors for various topographic settings, and introduce an empirical model for the estimation of topographic shielding factors based on landscape metrics. Copyright © 2008 John Wiley and Sons, Ltd.     

Tortoises and hares: dissolution,erosion and isostasy in landscape evolution

Denudation mechanisms differ fundamentally between limestone and silicate rock types, which are subject to very different rate thresholds and enhancers/inhibitors. Silicates are removed largely by erosion, the mechanical entrainment and transport of particles. This is a relatively high energy, and highly episodic, process which occurs only when a minimum threshold ?ow velocity is exceeded; it is inhibited by vegetation cover and favoured by strongly seasonal runoff. Limestone is removed largely by chemical dissolution at a rate directly proportional to runoff. Dissolution is a relatively low energy process that can occur at any ?ow velocity or in static water; in general it is enhanced by vegetation cover and non‐seasonality of runoff. These contrasting factors in the denudation of silicates versus limestone can produce strikingly uneven rates of surface lowering across a landscape, sometimes akin to the well known ‘tortoise and hare race’, where the slow and steady denudation of limestones may in the long term exceed the sometimes rapid, but often localized and episodic, erosion of silicates. Prolonged exposure of limestone to a humid temperate climate in a tectonically stable environment produces low‐relief corrosion plains in which limestone uplands are anomalous and, in most instances, due to recent unroo?ng from beneath a siliciclastic cover. In a highly seasonal or semi‐arid climate almost the exact inverse may develop, with ‘?ashy’ runoff and sparse vegetation favouring erosion rather than dissolution. Even under a constant humid climate progressive unroo?ng of a thick limestone unit within folded siliciclastics may lead to a topographic inversion over time, with the limestone outcrop always forming a topographic low ?anked by siliciclastic uplands. Valleys will be initiated on anticlinal crests, where the limestone is ?rst unroofed, but progressive lowering of the limestone causes these valleys to migrate to their ?nal position in the synclinal troughs. In humid climates isostatic compensation in response to slow, but continuous, denudation of extensive limestone outcrops may be a signi?cant factor in the development of relief on adjacent, more slowly eroding, silicate outcrops. Copyright © 2004 John Wiley & Sons, Ltd.     

APPLICATION OF TOPOGRAPHIC SLOPE AND ELEVATION VARIATION COEFFICIENT IN IDENTIFYING THE MOTUO ACTIVE FAULT ZONE

The eastern Himalaya syntaxis is located at the southeastern end of the Qinghai-Tibet Plateau and is the area where the Eurasian plate collides and converges with the Indian plate. The Namjabawa is the highest peak in the eastern section of the Himalayas, and the Yarlung Zangbo River gorge is around the Namjabawa Peak. The NE-striking Aniqiao Fault with right-lateral strike-slip is the eastern boundary fault of the Namjabawa syntaxis. Motuo Fault is in the east of and parallel to the Aniqiao Fault, distributing along the valley of the Yarlung Zangbo River. The section of Yarlung Zangbo River valley at the eastern side of the Namjabawa area is located in the southern foothills of the Himalayas and belongs to the subtropical humid climate zone with dense tropical rainforest vegetation. Dense vegetation, large terrain elevation difference, strong endogenetic and exogenic forces, and abundant valley deposition bring enormous difficulty to the research on active faults in this area. Since 1990s, surface morphology can be quantitatively expressed by digital elevation models as the rapid development of remote sensing technology. Geomorphic types and their characteristics can be quantified by geomorphological parameters which are extracted from DEM data, describing geomorphologic evolution and tectonic activity. But to date, researches based on quantitative geomorphic parameters are mainly focus on the differential uplift of regional blocks. In the study and mapping of active faults, surface traces of active faults are acquired by visual interpretation of remote sensing images. It has not been reported to identify the location of active faults via the change of quantitative geomorphic parameters. The distribution map of topographic elevation variation coefficient is suitable to reflect the regional erosion cutting and topographic relief, and the places with higher topographic elevation variation coefficient are more strongly eroded. In this paper, we attempt to identify the active faults and explore their distribution in the Yarlung Zangbo Gorge in the east of the Namjabawa Peak based on the application of two quantitative geomorphic parameters, namely, the topographic slope and the elevation variation coefficient. Using the DEM data of 30m resolution, two quantitative geomorphic parameters of topographic slope and elevation variation coefficient in Namjabawa and its surrounding areas were obtained on the ArcGIS software platform. On the topographic slope distribution map, the slope of the eastern and western banks of the Yarlung Zangbo River near Motuo is steep with a slope angle of more than 30°. Under the background of steep terrain, there are gentle slope belts of 5°~25° distributing intermittently and NE-striking. On the distribution map of topographic elevation variation coefficient, the elevation variation coefficient of the Yarlung Zangbo River near Motuo is greater than 0.9. On the background of the high topographic fluctuation area, it develops gently topographic undulating belts with elevation variation coefficient of 0.2~0.9. The belts are intermittently distributed and northeastern trending. Through the field geological and geomorphological investigation and trench excavation, it is found that the abnormal strips of the above-mentioned geomorphological parameters are the locations where the active faults pass. The above results show that the quantitative analysis of the topographic slope and the coefficient of variation of elevation can help us find active faults in areas with large terrain slope, serious vegetation coverage and high denudation intensity.     

Quantifying river incision into low‐relief surfaces using local and catchment‐wide 10Be denudation rates

Relief generation in non‐glaciated regions is largely controlled by river incision into bedrock but datable fluvial terraces that allow quantifying incision rates are not always present. Here we suggest a new method to determine river incision rates in regions where low‐relief surfaces are dissected by streams. The approach consists of three steps and requires the ¹⁰Be concentrations of a stream sediment sample and a regolith sample from the low‐relief surface. In the first step, the spatial distribution of ¹⁰Be surface concentrations in the given catchment is modelled by assuming that denudation rates are controlled by the local hillslope angles. The slope–denudation rate relation for this catchment is then quantified by adjusting the relation between slope angle and denudation rate until the average ¹⁰Be concentration in the model is equal to the one measured in the stream sediment sample. In the second step, curved swath profiles are used to measure hillslope angles adjacent to the main river channel. Third, the mean slope angle derived from these swath profiles and the slope–denudation relation are used to quantify the river incision rate (assuming that the incision rate equals the denudation rate on adjacent hillslopes). We apply our approach to two study areas in southern Tibet and central Europe (Black Forest). In both regions, local ¹⁰Be denudation rates on flat parts of the incised low‐relief surface are lower than catchment‐wide denudation rates. As the latter integrate across the entire landscape, river incision rates must exceed these spatially averaged denudation rates. Our approach yields river incision rates between ~15 and ~30 m/Ma for the Tibetan study area and incision rates of ~70 to ~100 m/Ma in the Black Forest. Taking the lowering of the low‐relief surfaces into account suggests that relief in the two study areas increases at rates of 10–20 and 40–70 m/Ma, respectively. Copyright © 2018 John Wiley & Sons, Ltd.     

Grain size-dependent 10Be concentrations in alluvial stream sediment of the Huasco Valley,a semi-arid Andes region

Terrestrial cosmogenic nuclide concentrations in sediment are used to quantify mean denudation rates in catchments. This article explores the differences between the ¹⁰Be concentration in fine (sand) and in coarse (1–3 or 5–10 cm pebbles) river sediment. Sand and pebbles were sampled at four locations in the Huasco Valley, in the arid Chilean Andes. Sand has ¹⁰Be concentrations between 4.8 and 8.3·10⁵ at g⁻¹, while pebbles have smaller concentrations between 2.2 and 3.3·10⁵ at g⁻¹. It appears that the different concentrations, systematically measured between sand and pebbles, are the result of different denudation rates, linked with the geomorphologic processes that originated them. We propose that the ¹⁰Be concentrations in sand are determined by the mean denudation rate of all of the geomorphologic processes taking place in the catchment, including debris flow processes as well as slower processes such as hill slope diffusion. In contrast, the concentrations in pebbles are probably related to debris flows occurring in steep slopes. The mean denudation rates calculated in the catchment are between 30 and 50 m/Myr, while the denudation rates associated with debris flow are between 59 and 81 m/Myr. These denudation rates are consistent with those calculated using different methods, such as measuring eroded volumes.