北京大学陆面过程模式PKULM(Peking University Land Model)介绍及检验

陆面过程模式是气候模式和天气模式的核心组成部分之一.在土壤-植被-大气耦合模式(Soil-Plant-Atmosphere Model, SPAM)的基础上,发展了新一代北京大学陆面过程模式PKULM(Peking University Land Model).本文首先介绍了PKULM的辐射传输、湍流输送、光合作用、土壤水热输送等过程的参数化方案;采用隐式迭代计算框架,发展并应用了一个快速的线性方程组求解算法,提高了模式计算稳定性;提出并使用了二分搜索算法计算气孔阻抗,避免了CLM(Community Land Model)等使用的迭代方法在干旱区不稳定的情况,提高了模式的适用性;采用水势为基础的土壤水分扩散方程,使模式能够模拟土壤饱和区的水分输送过程,为进一步与水文过程模式耦合奠定了基础;还发展了一个地表积水与径流过程的机理模型,提高了模式对地表水分平衡过程的模拟能力;最后,使用"中国西北干旱区陆-气相互作用观测试验"平凉站的资料对模式进行了检验并与NOAH(National Center for Environmental Prediction, Oregon State University, Air Force, and Hydrology Lab model)陆面过程模式的模拟结果进行了比较,结果表明PKULM能够较好地模拟西北半干旱区农田下垫面地气交换过程.     

Diagnosis of river basins as CO2 sources or sinks subject to sediment movement

Sediment movement during erosion, transport and deposition greatly affects the ecosystem of river basins. However, there is presently no consensus as to whether particular river basins act as carbon dioxide (CO₂) sources or sinks related to these processes. This paper introduces a rule‐of‐thumb coordinate system based on sediment delivery ratio (SDR) and soil humin content (SHC) in order to evaluate the net effect of soil erosion, sediment transport and deposition on CO₂ flux in river basins. The SDR–SHC system delineates CO₂ source and sink areas, and further divides the sink into strong and weak areas according to the world‐average line. The Yellow River Basin, most severely suffering soil erosion in the world, only appears to be a weak erosion‐induced CO₂ sink in this system. The average annual CO₂ sequestration is ~0·235 Mt from 1960 to 2008, a relatively small value considering its 3·1% contribution to the World's sediment discharge. The temporal analysis shows that the Yellow River Basin was once a source in the 1960s, but changed its role to become a weak sink in the past 40 years due to both anthropogenic and climatic influences. The spatial analysis identifies the middle sub‐basin as the main source region, and the lower as the main sink. For comparison, sediment‐movement‐related CO₂ fluxes of eight other major basins in four continents are examined. It is found that the six basins considered in the Northern Hemisphere appear to be sinks, while the other two in the Southern Hemisphere act as sources. Copyright © 2012 John Wiley & Sons, Ltd.     

Difference between static and dynamic angle of repose of uniform sediment grains

In the investigation of sediment transport, it is necessary to differentiate various definitions of angle of repose (AoR) available in the literature. The static AoR, composed of upper and lower angle of slope, forms just before and after slope instability, while the dynamic AoR can be observed when sediment grains are moving continuously down an inclined plane. In the present study, a series of laboratory experiments was conducted to measure static and dynamic angle of repose for uniform natural sedi-ments with median diameter of 0.28–4.38 mm. The results show that the different slope angles have different characteristics. The upper and dynamic AoR increase slightly with increasing grain diameter, while the lower AoR is not sensitive to changes in sediment size and may assume a constant value. The average of the upper and lower AoR is equivalent to the dynamic AoR, and the difference between them increases with increasing grain diameter. The present study suggests that the different angles of repose should be treated with caution when applying in investigations of bedload transport, dune migration and local scour development.     

Hydropedology of a mildly‐arid loess covered area,southern Israel

Extensive loess covered areas characterize the mildly arid areas of western Israel, where average annual rainfall is 280 mm. Hydrological data point to a peculiar hydrological behavior of the ephemeral streams. The frequency of sporadic flash floods is very high. However, even in extreme rain events peak discharges are extremely low. Hydrographs are usually characterized by very steep rising and falling limbs, representative of saturated areas, extending over a limited part of the watershed. Following this observation we advanced the hypothesis that storm channel runoff originated in the channel itself, with negligible contribution from the adjoining hillslopes. The study was based on two complementary approaches. The hydrological approach was based on the detailed analysis of rainfall–runoff relationships in a small watershed (11 km²) and on the analysis of the hydrological characteristics of the drainage network. The second approach was based on the toposequence concept. Several boreholes were dug along a hillslope 400 m long. Chemical data obtained show no significant difference in the downslope direction. Similar results were also obtained for the particle size distribution and soil moisture content. Data obtained perfectly fit the concept of ‘Partial Area Contribution’ as it presents an extreme case of hydrological discontinuity at the hillslope–channel interface. The lack of pedological trends in the downslope direction is an additional indication of the limited connectivity between the hillslopes and the adjoining channel. The limited connectivity is attributed to the prevalence of low rain intensities in the study area. The present study is also relevant to our understanding of pedological processes in dryland areas. The high frequency of intermittent low intensity rainstorms limits runoff generation and flow distances, and casts doubt on the general application of the toposequence approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

The role of eddies in solute transport and recovery in rock fractures: Implication for groundwater remediation

A better understanding of solute transport and retention mechanism in rock fractures has been challenging due to difficulty in their direct observations in microscale rough‐walled fractures. Six representative troughs in a rough‐walled fracture were selected for microscale observations of eddy formation with increasing flow velocity and its effect on spatiotemporal changes of solute concentration. This experimental study was enabled by a microscale visualization technique of micro particle image velocimetry. With increasing flow velocity (Re ≤ 2.86), no eddies were generated, and solutes along the main streamlines transported rapidly, whereas those near the wall moved slowly. A larger amount of solutes remained trapped at all troughs at Re = 2.86 than Re < 1. For Re = 8.57, weak eddies started to be developed at the troughs on the lee side, which little contributed to overall solute flushing in the fracture. Accordingly, a large of amount of water was needed for solute flushing. The flow condition of 1 < Re < 10, before a full development of eddies, was least favourable in terms of time and amount of remediation fluid required to reach a target concentration. After large eddies were fully developed at troughs on the lee side for Re = 17.13, solutes were substantially reduced by eddies with less amount of water. Fully developed eddies were found to enhance solute transport and recovery, as opposed to a general consensus that eddies trap and delay solutes. Direct inflow into troughs on the stoss side also made a great contribution to solute flushing out of the troughs. This study indicates that fully developed eddies or strong inflows at troughs are highly possible to form for Re > 10 and this flow range could be favourable for efficient remediation.     

Using stable isotopes to understand seasonal and interannual dynamics in moisture sources and atmospheric circulation in precipitation

Stable isotopes in precipitation are useful tracers to strengthen understanding of climate change and hydrological processes. In this study, the moisture sources of 190 precipitation events in Beijing were analysed using the Hybrid Single‐particle Lagrangian Integrated Trajectory model, based on which we studied the relation between variations in precipitation δ¹⁸O and dynamics in moisture sources and atmospheric circulation in seasonal and interannual timescales. Categorization of 7 groups of moisture sources was performed, among which oceanic moisture sources presented lower δ¹⁸O in precipitation than continental moisture sources. The results show that seasonal variations of precipitation δ¹⁸O were caused by changes of moisture sources. In summer, moisture from proximal oceans dominated vapour transport to Beijing due to increasing monsoon strength and resulted in a relatively small variation of precipitation δ¹⁸O. At the interannual timescale, the variations of δ¹⁸O in summer precipitation were related to dynamics in oceanic moistures, showing depleted values when the contribution of oceanic moistures, especially the proportion of long‐distance oceanic moisture, was high. Further analysis indicated that changes of oceanic moisture sources were controlled by the strength of summer monsoons. These findings address the complexity of moisture sources in midlatitude monsoon areas and suggest that isotopic signals in precipitation have the potential to deduce changes in moisture sources and atmospheric circulation and can therefore serve for palaeoclimate reconstruction.     

Pathways and distances of fluid flow during low-grade metamorphism: evidence from pyrite deposits of the Cameros Basin, Spain

New geochemical and sulphur isotopic data are presented for a number of pyrite deposits from the Late Jurassic–Early Cretaceous Cameros Basin, Spain. The deposits were formed at, or close to, the peak of metamorphism and are always related to sandstone units in the mainly metapelite sequence. Iron remained immobile and conservative, pyrite iron being derived by sulphidation of chlorite in the host metapelites. Reduced sulphur, however, was supplied from two external sources: thermochemical reduction of sulphate and release of S during metamorphism of sedimentary sulphides. These sources provided isotopically heavy and light S, respectively, with variation in pyrite isotopic composition between different deposits resulting from differences in their relative importance at each site. During metamorphism, the sandstone units acted as aquifers, carrying the sulphidic pore waters to locations where permeability provided by syn-depositional fractures on a scale of 0.5–5  m allowed its interaction with the metapelites. Transport distances for sulphide during metamorphism were of the order of hundreds of metres.     

离心实验在污染物迁移研究中的应用

离心实验模拟以其可以获得与原型一致的应力水平并且能够大大缩短原型历时而被应用于污染物迁移实验研究。在离心相似理论及离心实验模拟的相似基础上,重点回顾了离心机在饱和带水分、非饱和带水分、保守性溶质、NAPLs、重金属、核素迁移以及污染场地修复方面的实验应用。最后讨论了离心模拟中土壤预制和加速度选择的问题,并简单介绍了离心监测方法。可以得出结论:离心机能够成功用于各类物质迁移的实验研究中,离心实验模拟能够为理论和数值等分析方法快速提供真实可靠的参数依据,但离心实验的理论基础和监测方法需要进一步完善,此外也应积极开展更接近实际情况下的离心实验模拟研究。