Interpretation of a two-point flux stencil for skew parallelogram grids

Control-volume formulations for elliptic equations often use two-point flux stencils, even for skew grids. Any two-point flux stencil may be interpreted as a multipoint flux stencil. This yields a definition of the permeability (or conductivity) tensor. Formulas for calculating the permeability tensor, based on the user-specified quantities in the two-point flux stencil, are given. Numerical test examples demonstrate the validity of the derivation.      

长江中下游旱涝中期预报方法及其业务应用

介绍了双层多因子综合相似的长江中下游旱涝中期预报模型。该模型引入了具有明确动力学和热力学意义的地转西风动量经向输送和经向温度梯度诊断量，同时强调了西太平洋副热带高压对长江中下游夏季降水的关键作用，并考虑了对流层中、低层500hPa高度和850hPa温度的大尺度环流背景场。与过去业务应用的单层相似预报模型相比较，不仅物理意义和天气学含义明确，而且更符合对具有复杂动力和热力机制的降水预报的考虑。5年的预报试验和业务应用结果表明，该方法对长江中下游地区旱涝的中期趋势预报有较好的预报能力。     

祁连山疏勒河上游多年冻土区高寒草甸土壤CO2通量特征

研究青藏高原多年冻土区高寒草甸土壤CO₂通量有助于准确估算该区域的土壤CO₂排放, 对认识高原土壤碳循环及其对全球气候变化的响应具有重要意义. 利用静态箱-气相色谱法和LI-8100土壤CO₂通量自动测量系统对疏勒河上游多年冻土区高寒草甸土壤CO₂通量进行了定期观测, 结合气象和土壤环境因子进行了分析. 结果表明: 整个观测期高寒草甸土壤表现为CO₂的源, 土壤CO₂通量的日变化范围为2.52～532.81 mg·m^-2·h^-1. 土壤CO₂年排放总量为1 429.88 g·m^-2, 年均通量为163.23 mg·m^-2·h^-1; 其中, CO₂通量与空气温度和相对湿度、活动层表层2 cm、10 cm、20 cm、30 cm 土壤温度、含水量和盐分均显著相关. 2 cm土壤温度、空气温度和总辐射、空气温度、2 cm土壤盐分分别是影响活动层表层2 cm土壤完全融化期、冻结过程期、完全冻结期、融化过程期土壤CO₂通量的最重要因子. 在完全融化期、冻结过程期和整个观测期, 拟合最佳的温度因子变化分别能够解释土壤CO₂通量变化的72.0%、82.0%和38.0%, 对应的Q₁₀值分别为1.93、6.62和2.09. 冻融期(含融化过程期和冻结过程期)和完全冻结期的土壤CO₂排放量分别占年排放总量的15.35%和11.04%, 在年排放总量估算中不容忽视.     

On the role of the Asian monsoon in the angular momentum and kinetic energy balances of the tropics

The balance conditions of relative angular momentum and time-mean kinetic energy and their annual variations are studied for the Northern Hemisphere tropical belt. The belt is divided into two roughly equal size parts, the monsoon and the extramonsoon regions. The data used consist of all available daily rawinsonde reports from the world areological network for the two 5-year periods 1958–63 and 1968–73.In winter, the trade winds in the monsoon and extramonsoon regions are both sources of westerly relative angular momentum for the middle latitude circulation. However, it is found that the angular momentum gained in the extramonsoon region of the Tropics is mostly destroyed by a net southward flow of mass in that region, and becomes regenerated in the monsoon region by a net northward flow of mass there. This excess of angular momentum together with the angular momentum picked up locally in the monsoon region is almost all exported across its northern boundary. It is further found that in winter the Tropics are also an important source of mean kinetic energy for middle latitudes. Again almost all export of kinetic energy was found to take place across the northern boundary of the monsoon sector. Most of this energy must be generated through the pressure gradient term inside the monsoon region itself, the transformation from transient eddy kinetic energy being very small. The proper evaluation of the pressure gradient appears to be the main stumbling block in the present study, preventing us from estimating the generation and thereby, as a residual, the frictional dissipation in the two regions.In summer, the extramonsoon region remains a source of angular momentum, but the monsoon region with its surface westerlies acts as a sink, leading to a sharp reduction (and even a midsummer reversal) of the export into middle latitudes. Also the export of mean kinetic energy almost vanishes in summer, except for a small southward transfer across the equator. The calculations for two 5-year periods give very similar estimates and thereby show the reliability of the results.Parts of this paper were presented at the International Symposium on Monsoons, March 7–12, 1977 in New Delhi, India.     

Performance of a flux splitting scheme when solving the single-phase pressure equation discretized by MPFA

In this article we present a series of tests to study how well suited the TPFA coefficient matrix is as a preconditioner for the MPFA discrete system of equations in an iterative solver, using a flux splitting method. These tests have been conducted for single-phase flow for a wide range of anisotropy, heterogeneity, and grid skewness (mainly parallelogram grids). We use the K-orthogonal part of the MPFA transmissibilities for a parallelogram grid to govern the TPFA transmissibilities. The convergence of the flux splitting method is for each test case measured by the spectral radius of the iteration matrix.     

山地冰川流动模型探讨

近30 a来冰川动力学模型有了快速发展, 在南极、 格林兰冰盖预测中取得一系列重要成果, 对山地冰川的研究也初见端倪. 从冰川流动的力学过程出发, 利用本构方程、 理想冰川假设、 浅冰层近似(Shallow ice approximation)假设完整地推导了理想冰川流动的物理过程, 揭示了冰川流动的机理, 建立了气候变化和冰川自身重力引起的理想冰川物质和能量再分配的温度耦合三维流动模型. 结合山地冰川的冰床形态, 将理想冰川与实际冰川相结合, 使理想冰川流动模型更好地近似山地冰川的流动.     

Using climate to relate water discharge and area in modern and ancient catchments

Models relating sediment supply to catchment properties are important in order to use the geological record to deduce landscape evolution and interplay between tectonics and climate. Water discharge (Q _w) is an important factor in the widely used ‘BQART ’ model, which relates sediment load to a set of measurable catchment parameters. Although many of the factors in this equation may be independently estimated with some degree of certainty in ancient systems, water discharge (Q _w) certainly cannot. An analysis of a world database of modern catchments with 1255 entries shows that the commonly applied equation relating catchment area (A ) to water discharge (Q _w = 0·075A^0·8) does not predict water discharge from catchment area well in many cases (R ² = 0·5 and an error spanning about three orders of magnitude). This is because the method does not incorporate the effect of arid and wet climate on river water discharge. The inclusion of climate data into such estimations is an opportunity to refine these estimates, because generalized estimates of palaeoclimate can often be deduced on the basis of sedimentological data such as palaeosol types, mineralogy and palaeohydraulics. This paper investigates how the relationship between catchment area and river discharge varies with four runoff categories (arid, semi‐arid, humid and wet), which are recognizable in the geological record, and modifies the coefficient and exponent of the above‐mentioned equation according to these classes. This modified model yields improved results in relating discharge to catchment area (R ² = 0·95 and error spanning one order of magnitude) when core, outcrop or regional palaeoclimate reconstruction data are available in non‐arid systems. Arid systems have an inherently variable water discharge, and catchment area is less important as a control due to downstream losses. The model here is sufficient for many geological applications and makes it possible to include variations in catchment humidity in mass‐flux estimates in ancient settings.     

DNDC模型对我国旱地N2O释放的拟合对比分析

本文根据贵州省玉米－油菜轮作田、休耕地和大豆－冬小麦轮作田的N2O释放通量及其影响因子季节变化的是观测结果，对比分析了DNDC模型计算与三试验田田间观测在N2O释放通量及其主要影响因子季节变化上的拟合程度。结果表明，对N2O释放通量、土壤NO3^-和NH4^ 含量、土壤湿度和土温，模型计算结果与田间观测具有相似的季节变化模式和相近的数值，模型较好地反映了上述参数的动态变化过程。据此，推算出贵州省农业土壤N2O年释放量，并分析了气候变化和农业活动对释放量的潜在影响。     

Some aspects on the variations in depositional flux of excess Thorium-230 in the Central Indian basin during Late Quaternary

Examined in this paper is the tentative history of the depositional flux of²³⁰Th_xs (the unsupported fraction of²³⁰Th scavenged from the overlying water column), for the Late Quaternary period from a sediment core of the Central Indian Basin (CIB). The measured depositional flux of²³⁰Th_xs is found substantially higher than that of the possible theoretical flux from the overlying water column. Historical records, reconstructed from the²³⁰Th_xs chronology suggests that the depositional flux has varied considerably with time, reflecting an enhanced scavenging during the Holocene and the preceding interglacial periods whereas, comparatively lower flux than the predicted one occurred during the Last Glacial Maximum (LGM) period. The average ratio of the measured depositional flux to that of the predicted flux from the overlying water column, indicates that the core site acts as a sink for²³⁰Th_xs and based on the existence of bottom current activity; the²³⁰Th_xs could be the result of focusing of younger sediments. The depositional index (Di) has also been calculated to quantify the extent of lateral supply throughout the core with time. The estimated (Di) suggests that bottom focusing and feeble deposition and/or winnowing processes had occurred and that the former was most prevalent during the Holocene and the preceding interglacials, whereas the latter was observed at the LGM period.     

Inversion of multiple thermal indicators: Quantitative methods of determining paleoheat flux and geological parameters. II. Theoretical development for chemical, physical, and geological parameters

Using the logic for quantitative inversion of present-day downhole thermal indicators, the inversion procedure can be used to determine unknown, or poorly known, chemical and physical parameters as well as other geological quantities of interest which impact on burial history and thermal history of an evolving sedimentary basin. Some such quantities are: amount of erosion and timing of unconformities, paleo-overpressuring, stratigraphic age, timing of igneous intrusion and insertion temperature, overthrust timing and frictional heating, fault and slump timing, effects due to emplacement of a radiocative layer, and salt emplacement and dissolution timing. Combining a priori unknown values of these chemical/physical and geological parameters with unknown (a priori) paleoheat flux variations, a theoretical scheme, called thermal indicator tomography, is developed for the systematic determination of all parameters at the same time and on the same footing. Case histories will be discussed in subsequent papers in this series.