Implementation of Tidal Constituent Interpolation Method for the Israeli Coastline

Many ship-borne geodetic surveys at sea, such as Global Navigation Satellite System (GNSS)-based sea surface height (SSH) observation, acoustic profiling of the bottom, and others, deal with a dynamic topography which undergoes several changes during the survey campaign (e.g., changes in tide, salinity and currents). Those changes affect the measurements and may causes for some variations in the results. There are several methods for tidal variations correction, being the most dominant phenomena, such as tidal zoning, tidal constituent interpolation or ocean tidal models. In this study, we have implemented the tidal constituent interpolation method for the Israeli coastline in order to assess its quality and determine whether it is suitable for use in this particular region. This paper depicts the interpolation method, discusses some difficulties in the implementation for the Israeli coast and presents results from exemplary processing. In addition, we compare the results to those obtained using global and regional tidal models.     

介质密度反演偏导矩阵的精确计算方法

实现反演偏导矩阵的计算是基于导数最优化反演方法的关键,然而目前的地震反演几乎都是基于Zoeppritz方程近似实现的,使计算精度和适应范围受到限制.本文利用Zoeppritz方程建立了反射系数对地层介质密度比偏导方程,导出了Zoeppritz方程矩阵元对介质密度比的导数.通过求解偏导方程获得了反射系数对介质密度比偏导数的精确计算(考虑了速度中含介质密度的问题).利用数值算例分析了反射系数对介质密度比偏导数的变化特点.本文采用直接解法求解偏导矩阵方程组,获得了快的计算速度和高的计算精度,为实现地层介质密度反演(包括大角度反演)提供了偏导矩阵的计算方法.     

北京市房价与地价的动态关系——基于结构方程模型的实证分析

城市房价与地价之间的关系错综复杂,不仅受多种因素的交织影响,相互之间也存在动态关系。研究房价与地价关系的传统方法,如Granger因果检验和回归分析等,无法刻画房价与地价之间多维的网络状关系,相比之下,结构方程模型能同时处理多个内生潜变量,且不受观测指标共线性的影响,为刻画地价与房价的交互作用提供了新的工具。从住房与土地市场的供需传导机制出发,推导出房价与地价的结构模型,以北京市为例,运用2003-2013年居住用地价格和2014年在售楼盘价格,与北京市GIS电子地图相匹配,提取商服中心可达性、公共交通可达性、道路可达性、商服繁华度、设施便利性等解释变量,构建地价与房价结构方程模型,分析二者之间的结构关系。     

三维复电阻率电磁场正演研究

将Cole—Cole模型表示的地下异常体各剖分单元复电阻率,引入到体积分方程求解均匀大地三维电磁场计算中,实现电偶源地面激发、接收复电阻率三维电磁场正演模拟。     

Hydrodynamics and Thermodynamics of Newtonian Stars

We investigated spherically symmetric solution for nonrelativistic cosmological fluid equations and thermodynamic equation of state for Newtonian stars. It was shown that the assumption of a polytropic equation, , at the center of the star only suffices to integrate the equations explicitly. Our exact solution yields many fruitful results such as stellar stability, spherical oscillation and collapses of stars. Pressure, temperature, and density profiles inside stars were obtained. Central densities, pressures and temperatures of the Newtonian stars such as Sun, Jupiter and Saturn were also calculated. Collapse and expansion mechanism was explained by the heat transfer mechanism inside star. The upper bound value of white dwarf mass obtained by the Newtonian cosmological fluid equations turns out to be comparable to the static limit of Chandrasekhar one. Motion of the Universe was also discussed within the framework of Newtonian mechanics. Our calculation results without considering nuclear reactions inside stars may be applicable to the formation of protostars.     

转换波反射系数近似公式的比较及分析

多分量地震资料解释的首要问题就是多分量的标定和对比,其中最关键的是转换波反射系数的求取及子波的选取,然而Zoeppntz方程在数学表达上极其复杂,在物理意义上也不直观,因此很难得到直接的应用.为此,在详细介绍五种近似公式的基础上,利用六种模型分析了随入射角的变化各公式的适应性,最后得出结论,Aki和Richard的近似公式适应性普遍比较好.     

http://www.sciencedirect.com/science/article/pii/S1674987111001034

A comprehensive methodology that integrates Revised Universal Soil Loss Equation（RUSLE） model and Geographic Information System（GIS） techniques was adopted to determine the soil erosion vulnerability of a forested mountainous sub-watershed in Kerala,India.The spatial pattern of annual soil erosion rate was obtained by integrating geo-environmental variables in a raster based GIS method.GIS data layers including,rainfall erosivity（R）,soil erodability（K）,slope length and steepness（LS）,cover management （C） and conservation practice（P） factors were computed to determine their effects on average annual soil loss in the area.The resultant map of annual soil erosion shows a maximum soil loss of 17.73 t h^-1 y^-1 with a close relation to grass land areas,degraded forests and deciduous forests on the steep side-slopes（with high LS ）.The spatial erosion maps generated with RUSLE method and GIS can serve as effective inputs in deriving strategies for land planning and management in the environmentally sensitive mountainous areas.     

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop–soil system: Validation and application

Models for water transfer in the crop–soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop–soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop–soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration.     

Comparison of two different equations of state for application of carbon dioxide sequestration

We suggest that different equations of state (EOS) algorithms can and frequently will provide very different predictions of CO₂ migration following injection for sequestration. Rather than carry out an exhaustive examination of all EOS algorithms available, we elected to evaluate this general hypothesis by making detailed comparisons of simulation results of two very common EOS algorithms. We simulated and compared CO₂ migration patterns using two fundamentally different EOS algorithms – Modified Redlich-Kwong EOS (MRKEOS) and Span and Wagner EOS (SWEOS). In general, the predictions of thermophysical properties for both algorithms are close, except for a contrast in the predicted fugacity coefficient of CO₂, which subsequently propagates to a contrast in predicted solubility in water/brine. Typically, MRKEOS underestimates solubility of CO₂ compared to both SWEOS and experimental solubility data. In simulations of CO₂ migration, dissolution rates of separate-phase CO₂ predicted from the two EOS algorithms were significantly different, even for small contrasts in predicted fluid properties from EOS algorithms, resulting in markedly different migration patterns.     

Thermal equation of state of magnesite to 32 GPa and 2073 K

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.