平缓地区地形湿度指数的计算方法

地形湿度指数( topographic wetness index) 可定量模拟流域内土壤水分的干湿状况, 在流域 的土壤及分布式水文模型等研究中具有重要的意义。但现有的地形湿度指数计算方法在应用于 地形平缓地区时会得到明显不合理的结果, 即在河谷地区内, 地形湿度指数仅在狭窄的汇水线上 数值较高, 而在汇水线以外的位置则阶跃式地变为异常低的地形湿度指数值。本文针对此问题对 地形湿度指数的计算方法提出改进: 以多流向算法MFD- fg 计算汇水面积, 相应地以最大下坡计 算地形湿度指数, 再基于一个正态分布函数对河谷平原地区内的地形湿度指数进行插值处理。应 用结果表明, 所得地形湿度指数的空间分布不但能合理地反映平缓地区坡面上的水分分布状况, 并且在河谷地区内地形湿度指数值也都比较高, 其空间分布呈平滑过渡, 因而整个研究区域的水 分分布状况得到了比较合理的反映。     

The evolution of sedimentary basins on a viscoelastic lithosphere: theory and examples

Summary. A systematic approach is suggested for modelling the development of sedimentary basins. The theory, which partitions basin formation into initiating and isostatic adjustment processes, is applicable to all modes of basin formation if these processes are linear, or can he represented with sufficient accuracy in an incrementally linear form.
The dynamics of regional isostatic adjustment are characterized by the Heaviside space-time Green functions for the response of elastic and viscoelastic (Maxwell) thin plate models of the lithosphere. It is shown, by convolving the Heaviside—Green functions with cylindrical surface loads, that the rate of isostatic adjustment on a viscoelastic lithosphere is a function of the wavelength of the surface load, long wavelengths being compensated most rapidly.
Six archetypal initiating processes for sedimentary basin development are presented. These processes are those responsible for the subsidence of the Earth's surface which creates a depression in which water and sediments collect. Isostatic amplification of subsidence by sediment and water loads is cast in the form of an integral equation with isostatic Heaviside—Green functions as kernel.
Specific examples, the basins that result from a graben initiating process, are compared with the largest scale structure of the North Sea Basin, a basin that is known to be underlain by a graben system. A model, in which a 50-km wide graben subsides exponentially with a time constant of 5 × 10⁷yr during the interval 180–100 Myr bp , is shown to be consistent with the largest scale structure of the North Sea Basin if the underlying lithosphere is viscoelastic with a flexural rigidity of ∼5 × 10²⁵ Nm and relaxation time constant ∼ 10⁶ yr.     

基于DEM的地形湿度指数及其应用研究进展

地形湿度指数以数字高程模型（DEM）为基础,综合考虑了地形和土壤特性对土壤水分分布的影响,在流域土壤水分空间分布的研究中具有重要的理论与应用价值。根据对汇流面积时空变化的分析,地形湿度指数可分为静态、半动态和动态地形湿度指数等类型。在利用该指数评价土壤水分空间分布状况时,需要考虑计算方法与DEM网格单元大小的影响及其普遍适应性等问题。同时,针对黄土高原干旱半干旱区实际的水分循环过程和产流机制与现有地形湿度指数假设条件的差异,提出了地形湿度指数在黄土高原地区的研究设想与应用途径,以期为该区的水土保持和植被恢复研究提供理论和实践支持。     

基于模糊聚类的新疆典型高寒草原土壤pH值空间制图

准确、高效地掌握草原土壤属性的空间分布能够为草地资源境管理提供基础信息和参考依据。相比于传统土壤调查方法,基于模糊逻辑的土壤—环境推理能够提高野外采样效率和预测制图精度,被广泛应用于数字土壤制图。但由于土壤自身的空间变异性及其与环境条件间的非线性,现有推理模型的稳定性较低,尚未在高寒草原区进行应用。选择新疆巴音布鲁克典型亚高山草原地区约4 km^2区域为研究区,以高程、坡度、坡向、沿剖面曲率、沿等高线曲率、地形湿度指数6个地形因子为土壤环境因子,采用模糊C均值聚类(Fuzzy C-means Clustering,FCM)方法对环境因子聚类,得到9个环境因子组合,并在隶属度值高的环境因子组合中心共设置18个典型点。运用土壤—环境推理方法模拟研究区表层土壤pH值空间分布,其变化范围在7.170~8.186之间。选取35个独立样本进行精度检验(均匀采样点16个,横截面采样点9个,垂直带采样点10个),模拟结果与实测值基本吻合,且基于模糊聚类和土壤—环境推理方法的模拟精度高于普通克里格法和反距离权重法。通过基于模糊逻辑和土壤—环境推理的数字土壤制图方法在小尺度区域的运用验证,结果表明基于典型点的采样方案能够快速、有效地对区域土壤属性进行空间模拟,该方法对于类似小尺度的研究区同样有效。     

DEM resolution dependencies of terrain attributes across a landscape

This paper documents resolution dependencies in terrain analysis and describes how they vary across landform location. Six terrain attributes were evaluated as a function of DEM resolution—slope, plan curvature, profile curvature, north–south slope orientation, east–west slope orientation, and topographic wetness index. The research highlights the effect of varying spatial resolution through a spatial sampling/resampling scheme while maintaining sets of indexed sample points at various resolutions. Tested sample points therefore coincide exactly between two directly compared resolutions in terms of their location and elevation value. An unsupervised landform classification procedure based on statistical clustering algorithms was employed to define landform classes in a reproducible manner. Correlation and regression analyses identified sensitive and consistent responses for each attribute as resolution was changed, although the tested terrain attributes responded in characteristically different ways. These responses displayed distinguishable patterns among various landform classes, a conclusion that was further verified by a series of two‐sample, two‐tailed t‐tests.     

Sediment transport to the Laptev Sea–hydrology and geochemistry of the Lena River

This study focuses on the fluvial sediment input to the Laptev Sea and concentrates on the hydrology of the Lena basin and the geochemistry of the suspended particulate material. The paper presents data on annual water discharge, sediment transport and seasonal variations of sediment transport. The data are based on daily measurements of hydrometeorological stations and additional analyses of the SPM concentrations carried out during expeditions from 1975 to 1981. Samples of the SPM collected during an expedition in 1994 were analysed for major, trace, and rare earth elements by ICP-OES and ICP-MS. Approximately 700 km³ freshwater and 27 times 10⁶ tons of sediment per year are supplied to the Laptev Sea by Siberian rivers, mainly by the Lena River. Due to the climatic situation of the drainage area, almost the entire material is transported between June and September. However, only a minor part of the sediments transported by the Lena River enters the Laptev Sea shelf through the main channels of the delta, while the rest is dispersed within the network of the Lena Delta. Because the Lena River drains a large basin of 2.5 times 10⁶ km², the chemical composition of the SPM shows a very uniform composition. In contrast, smaller rivers with more restricted catchment areas exhibit significant differences.     

Depositional and tectonic evolution of a supradetachment basin: 40Ar/39Ar geochronology of the Nova Formation, Panamint Range, California

The Nova Basin contains an upper Miocene to Pliocene supradetachment sedimentary succession that records the unroofing of the Panamint metamorphic core complex, west of Death Valley, California. Basin stratigraphy reflects the evolution of sedimentation processes from landslide emplacement during basin initiation to the development of alluvial fans composed of reworked, uplifted sections of the basin fill. ⁴⁰Ar/³⁹Ar geochronology of volcanic units in middle and lower parts of the sequence provide age control on the tectonic and depositional evolution of the basin and, more generally, insights regarding the rate of change of depositional environments in supradetachment basins. Our work, along with earlier research, indicate basin deposition from 11.38 Ma to 3.35 Ma. The data imply sedimentation rates, uncorrected for compaction, of ~100 m Myr⁻¹ in the lower, high-energy part to ~1000 m Myr⁻¹ in the middle part characterized by debris-flow fan deposition. The observed variation in sediment flux rate during basin evolution suggests that supradetachment basins have complex depositional histories involving rapid transitions in both the style and rate of sedimentation.     

Unstable Asia: active deformation of Siberia revealed by drainage shifts

Regional incision and lateral shifts of rivers in the West Siberian Basin and surrounding areas show the action of long wavelength surface tilting, directed away from the Urals and Central Asian mountains and towards the Siberian Craton. In the north of the basin, surface uplift of individual folds is recorded by local lateral drainage migration. Lateral slopes of river valleys vary in gradient from 0.001 to 0.0001, generally decreasing with increasing river discharge. As a result of this surface deformation significant drainage shifts are taking place in three of the longest and highest discharge river systems on Earth: the Yenisei, Ob' and Irtysh. The deformation is most plausibly caused by subtle faulting at depth, below the thick basin fill of Mesozoic and Lower Cenozoic sediments. Active deformation of western Siberia appears to represent a previously unrecognised, far-field effect of the India–Eurasia collision, up to ∼1500 km north of the limit of major seismicity and mountain building. It adds ∼2.5 × 10⁶ km² to the region deformed by the collision, which is an area greater than the Himalayas and Tibet combined. It is also an analogue for the formation of low-angle unconformities in terrestrial sedimentary basins on the periphery of other orogenic belts.     

Palaeoclimatic record from a loess–soil profile in northeastern Bulgaria—II. Correlation with global climatic events during the Pleistocene

Magnetic susceptibility ( χ ) variations and the behaviour of the ratio of susceptibility to saturation magnetization ( χ/J _s ) along the loess/palaeosol section at Koriten (NE Bulgaria) are used to deduce climatic changes during the Pleistocene in southeastern Europe. A good correlation of susceptibility variations with the astronomically tuned oxygen isotope record from ODP site 677 enables us to propose a more precise dating of the upper part of the Bulgarian loess complex. Close correspondence between susceptibility and δ ^{1 8} O records demonstrates the global significance of the palaeoclimatic signal recorded, although differences in relative amplitudes of χ and χ/J _s and δ ^{1 8} O create difficulties in making quantitative estimates of the climatic humidity in the past. The role of local factors affecting the palaeoclimatic mineral magnetic record deduced from the profile studied in Bulgaria is discussed.     

Surface vertical displacements, potential perturbations andgravity changes of a viscoelastic earth model induced by internal point dislocations

Using the viscoelastic correspondence principle, we utilize the surface coseismic spheroidal deformation fields (i.e. vertical displacements, potential perturbations and gravity changes) of SNREI earth models caused by four typical types of point dislocation, derived by Sun & Okubo (1993 ), to deduce the fundamental formulas for spheroidal fields relevant to viscoelastic earth models. In computations, we employ a strike-slip dislocation on a vertical plane buried at the bottom of the lithosphere to estimate the maximal viscous relaxation responses to this kind of source that possibly exist on the surface of the earth. We take the seismic moment as 10²²  N  m, which is characteristic of an average large earthquake. The numerical results demonstrate that, if we take the viscosity as 10¹⁹  Pa  s in the asthenosphere, and 10²¹  Pa  s in the other mantle layers, the rates of surface vertical displacements and gravity changes within about 2.5° for the 10 postseismic years are respectively 1.5–8.1  cm  yr⁻¹ and 4.0–14.9  μgal  yr⁻¹ : the viscous relaxation for this mantle viscosity profile proceeds much faster than for a constant mantle viscosity of 10²¹  Pa  s.     

基于GIS的侵蚀冲沟与地貌因子的关系

冲沟是多种外力作用的结果,展示为多种地貌现象。选择胶东半岛北岩子口小流域作为研究对象,以正射影像（DOM）和数字高程模型（DEM）为数据源,选择土地利用、地形湿度指数、坡度坡长因子、水流动力指数、坡向、平面曲率和地表粗糙度等7个地貌因子,计算小流域土壤侵蚀敏感性,研究侵蚀冲沟与地貌因子的关系。结果表明:冲沟多发生在坡度坡长较大、水流动力强、地表湿度较低的林草地区域,多分布在阳坡的凹面;冲沟大多分布在土壤侵蚀中等敏感性及以上的区域,约占总冲沟的96.4%;在土壤侵蚀中等敏感性及以上区域,实验权重值对验证区冲沟分布的响应精度为87.9%,与实际值93.3%相差较小,精度较高。研究结果可为区域内冲沟侵蚀的评价提供一定科学依据。     

地形指数的物理意义分析

地形指数模型 ( TOPMODEL)用地形指数在流域中的空间格局来确定流域饱和缺水量的空间分布和产流区的空间位置与范围 ,物理意义明确。本文介绍了地形指数的物理基础及其与土壤水分的关系。分析了地形指数空间变化与汇流面积 a及局地坡度空间变化的关系 ,a范围内有效汇流面积的变化 ,提出了一种确定 a上限值的方法。     

Approaches to Modelling the Surface Albedo of a High Arctic Glacier

Broadband surface albedo measurements, made during the summer melt season at three weather stations on John Evans Glacier (79°40 ' N, 74°00 ' W), varied strongly with the solar zenith angle ( θ _z ). Tests were carried out to assess the impact of diurnal variations in surface albedo on seasonal net shortwave radiation ( K ^* ) totals. Removing the diurnal signal from albedo measurements by daily averaging of hourly measurements, or by applying midday measurements to all hours of the day, changed K ^* by up to 16%. Ignoring measurements made at θ _z & 75°, to account for measurement (cosine) error at high θ _z , decreased K ^* by between 5 and 18%. Given the sensitivity of K ^* to diurnal patterns in surface albedo, experiments were carried out with two albedo models. One model accounted for albedo variations with θ _z and one did not. The model driven by θ _z , when implemented within a surface energy balance model for John Evans Glacier, produced better melt estimates. This suggests that diurnal variations in surface albedo should be accounted for in energy balance models of glacier melt.     

Role of pore pressure generation in sediment transport within half-grabens

Sediment transport and overpressure generation are coupled primary through the impact of effective stress on subsidence and compaction. Here, we use mathematical modeling to explore the interactions between groundwater flow and diffusion-controlled sediment transport within alluvial basins. Because of lateral variation in permeability, proximal basin facies will have pore pressure close to hydrostatic levels while distal fine-grained facies can reach near lithostatic levels. Lateral variation in pore pressure leads to differential compaction, which deforms basins in several ways. Differential compaction reduces basin size, bends isochron surfaces across the sand–clay interface, restricts basinward progradation of sand facies, and reduces the amplitude of oscillation in the lateral position of the sand–clay interface especially in the deepest part of the section even when temporal sediment supply are held constant. Overpressure generation was found to be sensitive to change in sediment supply in permeable basins (at least 10⁻¹⁷ m² in our model). We found that during basin evolution, temporal variations in overpressure and sediment supply fluctuations are not necessarily in phase with each other, especially in tight (low permeability) basins (<10⁻¹⁷ m² in our model).     

Physiographically Controlled Allometry of Specific Sediment Yield in the Canadian Cordillera: A Lake Sediment-based Approach

It is generally supposed that specific sediment yield declines as the drainage basin area increases, as part of the mobilized sediment becomes trapped in the downstream cascade of storage zones. In British Columbia, using fluvial suspended sediment load data, Church and Slaymaker (Nature 1989, Vol 337, pp. 452–454) have observed a pattern of increasing specific sediment yield at all spatial scales up to 3×10⁴km². This trend has been attributed to the dominance of secondary remobilization of Quaternary sediments over primary denudation of the land surface. Using a larger data set of lake sediment-based estimates of long-term sediment yield, sub-regional patterns of specific yield have been investigated for the Canadian Cordillera. Between spatial scales of 0.9 and 190 km² sediment yield trends are differentiated by physiography, as indicated by the variable allometry observed in the specific sediment yield–drainage basin area relations. Highest sediment yields were observed in the Coast Mountains where specific sediment yields conform to the regional pattern described by Church and Slaymaker. However, in flat-lying plateau and major valley areas specific sediment yield decreases with increasing drainage area, thus conforming to the conventional model of sediment delivery. In several other sub-regions of intermediate relief there were no significant relations between specific yield and drainage area. These results suggest that no single model of sediment yield is adequate to describe sediment transfer processes in the Canadian Cordillera at the sub-regional scale.     

Deep-water renewal in a Scottish fjord: temperature, salinity and oxygen isotopes

The sea lochs (fjords) of north-west Scotland are located in a region of Europe particularly well situated to monitor changes in westerly air streams. The moisture transported in these air streams has a profound effect on regional precipitation, freshwater run-off and, in turn, sea loch circulation. The gentle slope of the regional salinity:δ¹⁸O mixing-line, defined as 0.18 ‰ per salinity unit, suggests that the temperature: δ¹⁸O relationship may be readily resolved in these coastal waters. Deep-water renewal events, both observed and predicted from empirical models, in the bottom-waters of Loch Etive provide an opportunity to assess the temperature, salinity and δ¹⁸O relationship. Predicted changes in δ¹⁸O_calcite as a function of changing salinity (ΔS) and changing temperature (ΔT) during deep-water renewal events suggest that >80% fall above analytical detection limits. The theoretical likelihood of recording such renewal events in the "palaeoclimate" record appears to be promising, but temperature and salinity change during renewal events may have either sign. Scottish fjords, because of the relatively small impact which salinity has on δ¹⁸O_water, may provide useful study sites in palaeoclimate research, particularly where palaeotemperature is the primary record of interest.     

Soil-landscape modelling and spatial prediction of soil attributes

Abstract

Explicit and quantitative models for the spatial prediction of soil and landscape attributes are required for environmental modelling and management. In this study, advances in the spatial representation of hydrological and geomorphological processes using terrain analysis techniques are integrated with the development of a field sampling and soil-landscape model building strategy. Statistical models are developed using relationships between terrain attributes (plan curvature, compound topographic index, upslope mean plan curvature) and soil attributes (A horizon depth, Solum depth, E horizon presence/absence) in an area with uniform geology and geomorphic history. These techniques seem to provide appropriate methodologies for spatial prediction and understanding soil landscape processes.     

Regional variations of heat flow differences with depth in Alberta, Canada

Summary. Differences between estimated average heat flow values for the Mesozoic and Cenozoic formations ( Q ₁) and estimated average heat flow values for the Palaeozoic formations below the erosional unconformity ( Q ₂) are calculated for the Alberta part of the western Canadian sedimentary basin. Significant heat flow differences exist for these two intervals and the map of Δ Q = Q ₁– Q ₂ shows that Q ₂ is generally greater than Q ₁ in the western and south-western part of Alberta, while in the northern part of the province Q ₂ is generally less than Q ₁. The regional variations of Δ Q are large, with standard deviation of 26 mW m⁻² and average value –13.5 mW m⁻². A regional trend of Δ Q correlates with topographic relief and the hydraulic head variations in the basin. It is shown that there is a heat flow increase with depth in water recharge areas and a decrease in heat flow with depth in the low topographic elevation water discharge areas when comparing the average heat flow in Mesozoic + Cenozoic and Palaeozoic formations.     

Models of crustal flow in the India–Asia collision zone

Surface velocities in parts of the India–Asia collision zone are compared to velocities calculated from equations describing fluid flow driven by topographically produced pressure gradients. A good agreement is found if the viscosity of the crust is ∼10²⁰ Pa s in southern Tibet and ∼10²² Pa s in the area between the Eastern Syntaxis and the Szechwan Basin. The lower boundary condition of the flow changes between these two areas, with a stress-free lower boundary in the area between the Szechwan basin and the Eastern Syntaxis, and a horizontally rigid but vertically deformable boundary where strong Indian lithospheric material underlies southern Tibet. Deformation maps for olivine, diopside and anorthite show our findings to be consistent with laboratory measurements of the rheology of minerals. Gravitationally driven flow is also suggested to be taking place in the Indo–Burman Ranges, with a viscosity of ∼10¹⁹–10²⁰ Pa s. Flow in both southern Tibet and the Indo–Burman Ranges provides an explanation for the formation of the geometry of the Eastern Himalayan Syntaxis. The majority of the normal faulting earthquakes in the Tibetan Plateau occur in the area of southern Tibet which we model as gravitationally spreading over the Indian shield.