STUDIES ON TERRACES AND RIVER SYSTEM CHANGES OF THE UPPER YELLOW RIVER

Ｉ.ＴｈｅＴｅｒｒａｃｅｓｏｆｔｈｅＵｐｐｅｒＹｅｌｌｏｗＲｉｖｅｒＴｈｅｔｅｒｒａｃｅｓｏｆｔｈｅｕｐｐｅｒＹｅｌｌｏｗＲｉｖｅｒａｒｅｍｏｓｔｌｙｄｉｓｐｅｒｓｅｄａｌｏｎｇｔｈｅｗｉｄｅｒｃｏｕｒｓｅｓｉｎｓｅｖｅｒａｌＭｅｓｏｚｏｉｃＣｅｎｏｚｏｉｃｂａｓｉｎｓ(Ｆｉｇｕｒｅ1).1.ＴｅｒｒａｃｅｓｂｅｔｗｅｅｎＮｇｏｒｉｎｇＬａｋｅａｎｄＧｙａｒｉｎｇＬａｋｅＴｈｅｓｕｒｆａｃｅｌａｙｅｒｓｏｆｔｈｅｓｅｄｉｍｅｎｔｉｎｔｈｅｗｉｄｅａｎｄｓｈａｌｌｏｗｂａｓｉｎｓｗｈｅｒｅｔｈｅＹｅｌｌ…     

Mass-front velocity of dry granular flows influenced by the angle of the slope to the runout plane and particle size gradation

The mass-front velocities of granular flows results from the joint action of particle size gradations and the underlying surfaces.However,because of the complexity of friction during flow movement,details such as the slope-toe impedance effects and momentum-transfer mechanisms have not been completely explained by theoretical analyses,numerical simulations,or field investigations.To study the mass-front velocity of dry granular flows influenced by the angle of the slope to the runout plane and particle size gradations we conducted model experiments that recorded the motion of rapid and long-runout rockslides or avalanches.Flume tests were conducted using slope angles of 25°,35°,45°,and 55° and three particle size gradations.The resulting mass-front motions consisted of three stages:acceleration,velocity maintenance,and deceleration.The existing methods of velocity prediction could not explain the slowing effect of the slope toe or the momentum-transfer steady velocity stage.When the slope angle increased from 25° to 55°,the mass-front velocities dropped significantly to between 44.4% and59.6% of the peak velocities and energy lossesincreased from 69.1% to 83.7% of the initial,respectively.The velocity maintenance stages occurred after the slope-toe and mass-front velocity fluctuations.During this stage,travel distances increased as the angles increased,but the average velocity was greatest at 45°.At a slope angle of 45°,as the median particle size increased,energy loss around the slope toe decreased,the efficiency of momentum transfer increased,and the distance of the velocity maintenance stage increased.We presented an improved average velocity formula for granular flow and a geometrical model of the energy along the flow line.     

Practical Incorporation of Multivariate Parameter Uncertainty in Geostatistical Resource Modeling

There is a need to bridge theory and practice for incorporating parameter uncertainty in geostatistical simulation modeling workflows. Simulation workflows are a standard practice in natural resource and recovery modeling, but the incorporation of multivariate parameter uncertainty into those workflows is challenging. However, the objectives can be met without considerable extra effort and programming. The sampling distributions of statistics comprise the core theoretical notion with the addition of the spatial degrees of freedom to account for the redundancy in the spatially correlated data. Prior parameter uncertainty is estimated from multivariate spatial resampling. Simulation-based transfer of prior parameter uncertainty results in posterior distributions which are updated by data conditioning and the model domain extents and configuration. The results are theoretically tractable and practical to achieve, providing realistic assessments of uncertainty by accounting for large-scale parameter uncertainty, which is often the most important component impacting a project. A simulation-based multivariate workflow demonstrates joint modeling of intrinsic shale properties and uncertainty in estimated ultimate recovery in a shale gas project. The multivariate workflow accounts for joint prior parameter uncertainty given the current well locations and results in posterior estimates on global distributions of all modeled properties. This is achieved by transferring the joint prior parameter uncertainty through conditional simulations.     

An Analytical Model For Threshold Crossing Rates Of Concentration Fluctuations In Dispersing Plumes

A simple analytical model is developed for the meanupcrossing rate of plume concentration fluctuations assuming that thisprocess can be well approximated by a lognormal process. The resultingexpression requires only the specification of the in-plume fluctuationintensity and in-plume Taylor micro-time scale and, hence, does notexplicitly involve the joint probability density function of theconcentration and its derivative. The analytical model provides agood fit to some field measurements of the mean upcrossing rate ina dispersing plume.     

Study on uncertainty of inversion and resolution limit in the presence of noise

This study analyzed the uncertainty of inversion and the resolution limit in the presence of noise by means of statistical experiments. The exhaustive method is adopted to obtain the global optimal solution in each experiment. We found that even with small level of noise, solutions fluctuate in a large range for the thin bed. The distribution of solutions in the presence of noise is closely related to the spread of the cost function in the absence of noise. As a result, the area of a certain neighborhood around the true solution on the spread of the cost function in the absence of noise is used to evaluate the uncertainty of inversion and the resolution limit in the presence of noise. In the case that the SNR (signal-to-noise ratio) is 5 in this study, solutions focus around the true solution with a very small uncertainty only when the bed thickness is greater than the reciprocal of the double predominant frequency of the convoluting wavelet.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

地层证据--对大别造山带汇聚历史的制约

大别造山带的汇聚演化历史一直存在不同的认识,来自扬子陆块和华北陆块显生宙以来的地层记录以及北淮阳地区的地层记录可能为两个陆块的汇聚模式提供制约。通过对扬子陆块与华北陆块显生宙以来岩相古地理进行比较,其中3个阶段相似, 2个阶段差别很大。以浊积岩发育作为俯冲阶段的标志,磨拉石建造出现作为闭合的标志,显生宙以来,大别造山带曾经历了两次俯冲和两次汇聚,俯冲时间分别是晚奥陶世—早志留世和早三叠世,汇聚的时间分别为中 晚泥盆世和中三叠世,早期的汇聚具有软碰撞特征,以上认知得到了北淮阳地区地层、古生物和沉积记录的支持,同时也与大别山造山带榴辉岩的同位素年龄数据基本吻合。高压 超高压变质岩在晚三叠世处于快速折返阶段,早侏罗世以前到达地表,其平均折返隆升速率为39～66mm/a。     

Reconstruction of cropland spatial patterns and its spatiotemporal changes over the 20th century on the Songnen Plain,Northeast China

We initially estimated the cropland area at county level using local historical documents for the Songnen Plain (SNP) in the 1910s and 1930s. We then allocated this cropland area to grid cells with a size of 1 km × 1 km, using a range of cultivation possibilities from high to low; this was based on topography and minimum distances to rivers, settlements, and traffic lines. Cropland areas for the 1950s were obtained from the Land Use Map of Northeast China, and map vectorization was performed with ArcGIS technology. Cropland areas for the 1970s, 1980s, 1990s, 2000s, and 2010s were retrieved from Landsat images. We found that the cropland areas were 4.92 × 10⁴ km² and 7.60 × 10⁴ km², accounting for 22.8% and 35.2% of the total area of the SNP in the 1910s and 1930s, respectively, which increased to 13.14 × 10⁴ km², accounting for 60.9% in the 2010s. The cropland increased at a rate of 1.18 × 10⁴ km² per decade from the 1910s to 1970s while it was merely 0.285 × 10⁴ km² per decade from the 1970s to 2010s. From the 1910s to 1930s, new cultivation mainly occurred in the central SNP while, from the 1930s to 1970s, it was mainly over the western and northern parts. This spatially explicit reconstruction could be offered as primary data for studying the effects of changes in human-induced land cover based on climate change over the last century.     

CGE模型在地理学综合研究中的应用

地理学研究以综合性和区域性为特色。提高地表过程研究中人文与自然要素综合分析的深度与广度,寻找学科交叉的综合途径是地理学中具有挑战性的课题。可计算一般均衡模型(Computable General Equilibrium,CGE) 是环境与经济系统之间互动影响关系模拟的有力工具,在地理学综合研究中得到广泛应用。本文概括了CGE模型在资源环境的可持续性、气候变化与资源税、灾害与综合风险管理、污染物的区域环境过程等多个方面的应用,从模拟主体、模型约束、模拟方法以及过程与结果表达4 个角度,总结了CGE模型与地理学综合研究之间的契合点,并指出与地理学其他模型的耦合以及情景模拟的地理学依据是CGE模型下一步的发展方向。