Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Numerical analysis of the symmetric methods

Aimed at the initial value problem of the particular second-order ordinary differential equations,y =f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.     

Low-low satellite-to-satellite tracking: a comparison between analytical linear orbit perturbation theory and numerical integration

Low-low satellite-to-satellite tracking (ll-SST) range-rate observations have been predicted by two methods: one based on a linear perturbation theory in combination with the Hill equations, and one based on solving the equations of motion of two low-flying satellites by numerical integration. The two methods produce almost equivalent Fourier spectra of the range-rate observations after properly taking into account a few resonant terms. For a typical GRACE-type configuration, where the two satellites trail each other at a distance of 300 km at an altitude of 460 km, and in the presence of the EGM96 gravity field model, complete to spherical harmonic degree and order 70, the agreement between the Fourier spectra is about 1 mm/s compared to a root-mean-square (RMS) value of more than 220 mm/s for the range-rate signal. The discrepancy of 1 mm/s can be reduced significantly when not taking into account perturbations caused by the J₂ term. Excluding the J₂ term, the agreement between the two methods improves to 0.4 mm/s compared to a RMS value of 6 mm/s for the range-rate signal. These values are 0.01 and 2.3 mm/s when ignoring the spectrum for frequencies below two cycles per orbital revolution, reducing the discrepancy even further to about 0.5% of the signal. The selected linear perturbation theory is thus capable of modeling gravity field induced range-rate observations with very high precision for a large part of the spectrum.     

Analysis and classification of hyperspectral data for mapping land degradation: An application in southern Spain

Desertification is a severe stage of land degradation, manifested by “desert-like” conditions in dryland areas. Climatic conditions together with geomorphologic processes help to mould desert-like soil surface features in arid zones. The identification of these soil features serves as a useful input for understanding the desertification process and land degradation as a whole. In the present study, imaging spectrometer data were used to detect and map desert-like surface features. Absorption feature parameters in the spectral region between 0.4 and 2.5 μm wavelengths were analysed and correlated with soil properties, such as soil colour, soil salinity, gypsum content, etc. Soil groupings were made based on their similarities and their spectral reflectance curves were studied. Distinct differences in the reflectance curves throughout the spectrum were exhibited between groups. Although the samples belonging to the same group shared common properties, the curves still showed differences within the same group.Characteristic reflectance curves of soil surface features were derived from spectral measurements both in the field and in the laboratory, and mean reflectance values derived from image pixels representing known features. Linear unmixing and spectral angle matching techniques were applied to assess their suitability in mapping surface features for land degradation studies. The study showed that linear unmixing provided more realistic results for mapping “desert-like” surface features than the spectral angle matching technique.     

Application of Linear Thermodynamics to the Atmospheric System. Part Ⅰ: Linear Phenomenological Relations and Thermodynamic Property of the Atmospheric System

从一般的热力学原理或其它自然原理对唯象关系所强加的限制，能够演绎出大气系统的一系列热力学性质。利用非平衡态线性热力学导出了湍流K闭合理论中湍流交换系数同唯象系数的关系，从理论上证明大气系统热量湍流输送同水泡之间存在交叉耦合，还导出了湍流强度同速度和位温梯度的关系，从而证明速度和位温空间分布的非均匀性是湍流之源。并证明湍流强度定理，不可压缩气体和各向同性湍流大气中，湍流强度正比于速度与位温梯度的标积。进而证明大气涡旋定理，位温的切变将导致涡旋运动或各种环流运动，速度涡度等于速度同位温相对梯度的矢积。展现了线性热力学在大气系统的应用前景。     

固溶体矿物的性质对计算河南西峡红柱石矿形成条件的影响

自然界的大多数造岩矿的是由两种或多种组分混合而成的固溶体。因此，在矿物，岩石和矿床等地质体的平衡热力学研究中不可避免地要涉及固溶体矿物热力学地为的问题。河南省西峡县杨乃沟红柱石矿泥质片岩中仅出现红柱石而没有蓝晶石和夕线石，红柱石矿物包裹体测温结果为420－600℃。通过矿物共生平衡p-t曲线所计算出的片岩形成温度和压力条件也应位于上述岩石学和岩相学等特征所限定的范围内。以此为基础，通过构筑矿物共生p-t平衡曲线的方法，对片岩中若干固溶体矿物理想性／非理想性问题进行了研究。其中斜长石必须以非理想固溶体处理，而石榴子石，黑云母和白云母视作理想固溶体的计算结果同片岩的岩石学和岩相学等特征相一致。     

高温高压下电解质溶液热力学性质的研究进展

高温高压下电解质溶液热力学性质的研究不仅对探讨地球内部流体的作用具有重要意义，而且在很多实际工业过程中也都有广泛应用。本综述了近10年来高温高压下电解质溶液热力学参数的实验测量和理论计算的进展，实验设备还有待进一步的改进和完善，以获得更高温度和压力下的数据，理论计算不仅有待于实验的发展，也与溶液的结构理论密切相关。     

An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

In situ methods for studying adsorbed phases at the solid/liquid interface: microcalorimetry and ellipsometryLes méthodes d'étude in situ des phases adsorbées : microcalorimétrie et ellipsométrie

The applications of calorimetric and ellipsometric methods to the understanding of adsorption phenomena occurring at the solid–liquid interface are presented with examples in the field of surfactant adsorption. The various experimental approaches are presented and compared, i.e., immersion calorimetry, titration calorimetry and liquid flow calorimetry. It is shown how they can be used to follow, in situ, the building steps of the surfactant adsorbed layer. On the other hand, ellipsometry allows the adsorption isotherm, the adsorption kinetics as well as the thickness of the adsorbed layer to be determined on well-defined flat surfaces. To cite this article: R. Denoyel, C. R. Geoscience 334 (2002) 689–702.     

Dynamic processes of a simple linear dune—a study in the Taklimakan Sand Sea, China

Samples of dune sands, surveys of the morphology and field measurements of wind velocity and direction of a simple linear dune in Taklimakan Sand Sea show that the airflow and sand flux vary with the change of wind direction on the dune surface. Decrease of the airflow stress on the lee flank does not result in much decrease of the sand flux because of the low threshold shear velocities and the airflow conditions. There are no significant relations between the sand flux on the lee flank and the angle of incidence of the airflow. The low threshold shear velocities and the maintenance of the sand flux at the lee flank are the main mechanisms keeping the linear shape of the dunes. Measurements of the sand flux shows that it reaches a maximum on the crest of the dune. The grain size of the transported sands has some differences compared to that of the dune surface. The sands transported are finer than that on the dune surface, but better sorted under the influence of the medium to low wind activity. The field experiment results exhibit that it is possible for the dunes to be shaped as linear dunes during the processes of accumulation and elongation.