Using GIS-based distributed soil loss modeling and morphometric analysis to prioritize watershed for soil conservation in Bago river basin of Lower Myanmar

Bago River is an important river in Myanmar. Although shorter than other rivers, it has its own river system, and people along the river rely heavily on it for their daily lives. The upper part of the watershed has changed rapidly from closed forest to open forest land in the 1990s. Since the recent degradation of the forest environment, annual flooding has become worse during the rainy season in Bago City. This paper aims at determining soil conservation prioritization of watershed based on soil loss due to erosion and morphometric analysis in the Bago Watershed by integrating remote sensing and geographic information system (GIS) techniques. In this study, soil erosion of the Bago watershed was determined using the Universal Soil Loss Equation. Such factormaps as rainfall, soil erodibility, slope length gradient, and crop management were compiled as input parameters for the modeling; and the soil loss from 26 sub-watersheds were estimated. Then, the soil erosion maps of the Bago watershed for 2005 were developed. The resulting Soil Loss Tolerance Map could be utilized in developing watershed management planning, forestry management planning, etc.     

Fast numerical solution of the linearized Molodensky problem

 When standard boundary element methods (BEM) are used in order to solve the linearized vector Molodensky problem we are confronted with two problems: (1) the absence of O(|x|⁻²) terms in the decay condition is not taken into account, since the single-layer ansatz, which is commonly used as representation of the disturbing potential, is of the order O(|x|⁻¹) as x→∞. This implies that the standard theory of Galerkin BEM is not applicable since the injectivity of the integral operator fails; (2) the N×N stiffness matrix is dense, with N typically of the order 10⁵. Without fast algorithms, which provide suitable approximations to the stiffness matrix by a sparse one with O(N(logN) ^s ), s≥0, non-zero elements, high-resolution global gravity field recovery is not feasible. Solutions to both problems are proposed. (1) A proper variational formulation taking the decay condition into account is based on some closed subspace of co-dimension 3 of the space of square integrable functions on the boundary surface. Instead of imposing the constraints directly on the boundary element trial space, they are incorporated into a variational formulation by penalization with a Lagrange multiplier. The conforming discretization yields an augmented linear system of equations of dimension N+3×N+3. The penalty term guarantees the well-posedness of the problem, and gives precise information about the incompatibility of the data. (2) Since the upper left submatrix of dimension N×N of the augmented system is the stiffness matrix of the standard BEM, the approach allows all techniques to be used to generate sparse approximations to the stiffness matrix, such as wavelets, fast multipole methods, panel clustering etc., without any modification. A combination of panel clustering and fast multipole method is used in order to solve the augmented linear system of equations in O(N) operations. The method is based on an approximation of the kernel function of the integral operator by a degenerate kernel in the far field, which is provided by a multipole expansion of the kernel function. Numerical experiments show that the fast algorithm is superior to the standard BEM algorithm in terms of CPU time by about three orders of magnitude for N=65 538 unknowns. Similar holds for the storage requirements. About 30 iterations are necessary in order to solve the linear system of equations using the generalized minimum residual method (GMRES). The number of iterations is almost independent of the number of unknowns, which indicates good conditioning of the system matrix. Received: 16 October 1999 / Accepted: 28 February 2001     

Thermodynamic modeling of binary CH4–CO2 fluid inclusions

An equation of state (EOS) explicit in Helmholtz free energy has been improved to calculate the PVTx and vapor–liquid phase equilibrium properties of CH₄–CO₂ fluid mixture. This EOS, where four mixing parameters are used, is based on highly accurate EOSs recommended by NIST for pure components (CH₄ and CO₂) and contains a simple generalized departure function presented by Lemmon and Jacobsen (1999). Comparison with experimental data available indicates that the EOS can calculate both vapor–liquid phase equilibrium and volumetric properties of this binary fluid system with accuracy close to that of experimental data up to high temperature and pressure within full range of composition. The EOS of CH₄–CO₂ fluid, together with the updated Gibbs free energy model of solid CO₂ (dry ice), is applied to calculate the CH₄ content (x_CH4) and molar volume (V_m) of the CH₄–CO₂ fluid inclusion based on the assumption that the volume of an inclusion keeps constant during heating and cooling. $V_{\text{m}}-x_{{\text{CH}}_4}$ diagrams are presented, which describe phase transitions involving vapor, liquid and CO₂ solid phases of CH₄–CO₂ fluid inclusions. Isochores of CH₄–CO₂ inclusions at given $x_{{\text{CH}}_4}$ and V_m can be easily calculated from the improved EOS.     

Equation of state constraints from neutron stars

Recent measurements of thermal radiation from neutron stars have suggested a rather broad range of radiation radii ( ). Sources in M13 and Omega Cen imply R _∞∼12–14 km, but X7 in 47 Tuc implies R _∞∼16–20 km and RX J1856-3754 R _∞>17 km. If these measurements are all correct, only a limited selection of EOS’s could be consistent with them, but a broad range of neutron star masses (up to 2 M_⊙) would also be necessary. The surviving equations of state are incompatible with significant softening above nuclear saturation densities, such as would occur with Boson condensates, a low-density quark-hadron transition, or hyperons. Other potential constraints, such as from QPO’s, radio pulsar mass and moment of inertia measurements, and neutron star cooling, are compared. US DOE Grant DE-FG02-87ER-40317.     

Equation of state of iron–silicon alloys to megabar pressure

The stability and pressure–volume equation of state of iron–silicon alloys, Fe-8.7 wt% Si and Fe-17.8 wt% Si, have been investigated using diamond-anvil cell techniques up to 196 and 124 GPa, respectively. Angular–dispersive X-ray diffractions of iron–silicon alloys were measured at room temperature using monochromatic synchrotron radiation and an imaging plate (IP). A bcc–Fe-8.7 wt% Si transformed to hcp structure at around 1636 GPa. The high-pressure phase of Fe-8.7 wt% Si with hexagonal close-packed (hcp) structure was found to be stable up to 196 GPa and no phase transition of bcc–Fe-17.8 wt% Si was observed up to 124 GPa. The pressure–volume data were fitted to a third-order Birch–Murnaghan equation of state (BM EOS) with zero–pressure parameters: V₀=22.2(8) ų, K₀=198(9) GPa, and K₀=4.7(3) for hcp–Fe-8.7 wt% Si and V₀=179.41(45) ų, K₀=207(15) GPa and K₀=5.1(6) for Fe-17.8 wt% Si. The density and bulk sound velocity of hcp–Fe-8.7 wt% Si indicate that the inner core could contain 3–5 wt% Si.     

Thermal equation of state of majorite with MORB composition

In situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press at SPring-8 on majoritic garnet synthesized from natural mid-ocean ridge basalt (MORB), whose chemical composition is close to the average of oceanic crust, at 19 GPa and 2200 K. Pressure-volume-temperature data were collected using a newly developed high-pressure cell assembly to 21 GPa and 1273 K. Data were fit to the high-temperature Birch-Murnaghan equation of state, with fixed values for the ambient cell volume (V₀ = 1574.14(4) Å³) and the pressure derivative of the isothermal bulk modulus (K^′T = 4). This yielded an isothermal bulk modulus of K_T0 = 173(1) GPa, a temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.022(5) GPa K⁻¹, and a volumetric coefficient of thermal expansivity α = a + bT with values of a = 2.0(3) × 10⁻⁵ K⁻¹ and b = 1.0(5) × 10⁻⁸ K⁻². The derived thermoelastic parameters are very similar to those of pyrope. The density of subducted oceanic crust compared to pyrolitic mantle at the conditions in Earth's transition zone (410-660 km depth) was calculated using these results and previously reported thermoelastic parameters for MORB and pyrolite mineral assembledges. These calculations show that oceanic crust is denser than pyrolitic mantle throughout the mantle transition zone along a normal geotherm, and the density difference is insensitive to temperature at the pressures in lower part of the transition zone.     

Stability and P–V–T equation of state of KAlSi3O8-hollandite determined by in situ X-ray observations and implications for dynamics of subducted continental crust material

In situ X-ray diffraction measurements of KAlSi₃O₈-hollandite (K-hollandite) were performed at pressures of 15–27 GPa and temperatures of 300–1,800 K using a Kawai-type apparatus. Unit-cell volumes obtained at various pressure and temperature conditions in a series of measurements were fitted to the high-temperature Birch-Murnaghan equation of state and a complete set of thermoelastic parameters was obtained with an assumed K′_300,0=4. The determined parameters are V _300,0=237.6(2) Å³, K _300,0=183(3) GPa, (?K _T,0/?T) _P =?0.033(2) GPa K^?1, a ₀=3.32(5)×10^?5 K^?1, and b ₀=1.09(1)×10^?8 K^?2, where a ₀ and b ₀ are coefficients describing the zero-pressure thermal expansion: α _T,0 = a ₀ + b ₀ T. We observed broadening and splitting of diffraction peaks of K-hollandite at pressures of 20–23 GPa and temperatures of 300–1,000 K. We attribute this to the phase transitions from hollandite to hollandite II that is an unquenchable high-pressure phase recently found. We determined the phase boundary to be P (GPa)=16.6 + 0.007 T (K). Using the equation of state parameters of K-hollandite determined in the present study, we calculated a density profile of a hypothetical continental crust (HCC), which consists only of K-hollandite, majorite garnet, and stishovite with 1:1:1 ratio in volume. Density of HCC is higher than the surrounding mantle by about 0.2 g cm^?3 in the mantle transition zone while this relation is reversed below 660-km depth and HCC becomes less dense than the surrounding mantle by about 0.15 g cm^?3 in the uppermost lower mantle. Thus the 660-km seismic discontinuity can be a barrier to prevent the transportation of subducted continental crust materials to the lower mantle and the subducted continental crust may reside at the bottom of the mantle transition zone.     

优化偏最小二乘光度法同时测定铑铱钯的研究

采用局部加权、随机矩阵法等步骤,使建立的校准模型能充分反映未知样品的数字结构特征,达到优化的目的,从而改善偏最小二乘法在处理实际未知样品数据时的解析能力,以提高计算结果的精度。用该优化偏最小二乘法对模拟地质样品中痕量贵金属Ｒｈ、Ｉｒ、Ｐｄ进行多组分光度测定,其分析结果的相对误差均小于１０％;分别测定９．５μｇ／Ｌ、３５．０μｇ／Ｌ、６６．０μｇ／Ｌ含量的Ｒｈ、Ｉｒ和Ｐｄ７次,其标准偏差Ｒｈ为０．２;Ｉｒ为１．１;Ｐｄ为２．０。