An Adaptive Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.     

Dehydration experiments on natural omphacites: qualitative and quantitative characterization by various spectroscopic methods

A series of natural omphacites from a wide range of P, T occurrences were investigated by electron microprobe (EMP), infrared (IR)-, Mössbauer (MS)- and optical spectroscopy in the UV/VIS spectral range (UV/VIS), secondary ion mass spectrometry (SIMS) and single crystal structure refinement by X-ray diffraction (XRD) to study the influence of hydrogen loss on valence state and site occupancies of iron. In accordance with literature data we found Fe²⁺ at M1 as well as at M2, and in a first approach assigned Fe³⁺ to M1, as indicated by MS and XRD results. Hydrogen content of three of our omphacite samples were measured by SIMS. In combination with IR spectroscopy we determined an absorption coefficient: ε _i,tot = 65,000 ± 3,000 lmol_H2O ^?1 cm^?2. Using this new ε _i,tot value, we obtained water concentrations ranging from 60 to 700 ppm H₂O (by weight). Hydrogen loss was simulated by stepwise heating the most water rich samples in air up to 800°C. After heat treatment the samples were analyzed again by IR, MS, UV/VIS, and XRD. Depending on the type of the OH defect, the grade of dehydration with increasing temperature is significantly different. In samples relatively poor in Fe³⁺ (<0.1 Fe³⁺ pfu), hydrogen associated with vacancies at M2 (OH bands around 3,450 cm^?1) starts to leave the structure at about 550°C and is completely gone at 780°C. Hydrogen associated with Al³⁺ at the tetrahedral site (OH bands around 3,525 cm^?1, Koch-Müller et al., Am Mineral, 89:921–931, 2004) remains completely unaffected by heat treatment up to 700°C. But all hydrogen vanished at about 775°C. However, this is different for a more Fe³⁺-rich sample (0.2 Fe³⁺ pfu). Its IR spectrum is characterized by a very intense OH band at 3,515 cm^?1 plus shoulder at 3,450 cm^?1. We assign this intense high-energy band to vibrations of an OH dipole associated with Fe³⁺ at M1 and a vacancy either at M1 or M2. OH release during heating is positively correlated with decrease in Fe²⁺ and combined with increase in Fe³⁺. That dehydration is correlated with oxidation of Fe²⁺ is indirectly confirmed by annealing of one sample in a gas mixing furnace at 700°C under reducing conditions keeping almost constant OH^? content and giving no indication of Fe²⁺-oxidation. Obtained data indicate that in samples with a relatively high concentration of Fe²⁺ at M2 and low-water concentrations, i.e., at a ratio of Fe^{2+ M2}/H > 10 dehydration occurs by iron oxidation of Fe²⁺ exclusively at the M2 site following the reaction: \( {\left[ {{\text{Fe}}^{{{\text{2 + [ M2]}}}}{\text{OH}}^{ - } } \right]} = {\left[ {{\text{Fe}}^{{{\text{3 + [ M2]}}}} {\text{O}}^{{{\text{2}} - }} } \right]} + {\text{1/2}}\;{\text{H}}_{{\text{2}}} \uparrow . \) In samples having relatively low concentration of Fe²⁺ at M2 but high-water concentrations, i.e., ratio of Fe^{2+ M2}/H < 5.0 dehydration occurs through oxidation of Fe²⁺ at M1.     

多源三维模型融合技术在三调城镇村内部细化调查中的应用研究

城镇村内部土地利用细化调查是第三次全国土调查中的重要组成部分,全面细化和完善的城镇村庄内部土地利用调查成果,有助于掌握翔实准确的全国国土利用现状和自然资源变化情况。本文基于多源三维模型融合技术,生成高精度高空地面、室内室外一体化三维模型,将其应用于城镇村内部土地利用细化调查工作中,修正城镇村内部土地利用现状图斑边界和地类属性等内容,最终形成准确的城镇村内部土地利用现状成果。利用融合后的三维模型上开展城镇村内部土地利用细化调查是一种有效的调查模式,可为三调高精度高效率调查工作提供技术支撑。     

似大地水准面拟合中已知点的选择

似大地水准面的拟合是测绘中的一个热点问题,其中已知点的选择对似大地水准面拟合结果有显著影响。文中探讨似大地水准面拟合中已知点的选择问题,分析比较传统手工方法、全组合法与逐点剔除法选择已知点,并用实际算例验证逐点剔除法的优点,所得结论对实际工程应用具有一定的参考价值。     

基于波形相似特征的微地震事件初至拾取及全局校正

微地震事件初至拾取是井下微地震监测数据处理的关键步骤之一.初至误差的存在会使微地震震源定位结果产生较大偏差,进而影响后续的压裂裂缝解释.通常初至拾取过程对所有的微地震事件选择相同的特征函数并采用一致的拾取参数进行统一处理,然而当事件的能量、震源机制、传播路径以及背景噪声等存在明显差异时,所得初至拾取结果差别显著.为了提高微地震事件初至拾取标准一致性,本文提出基于波形相似特征的初至拾取及全局校正方法.该方法首先利用互相关函数对每个事件内的各道记录进行时差校正,得到初始初至信息并形成叠加道,再对所有事件的叠加道进行全局互相关得到事件间初至相对校正量,最终初至结果可以通过各个事件的初始初至信息与其相对校正量相加得到.方法将所有微地震事件初至结果作为一个整体处理,从而能够克服常规方法初至拾取标准一致性差的缺陷.实际资料处理结果表明,相比于常规方法,该方法可以有效提高事件初至拾取和定位结果的一致性.     

Grid-size dependence of Cauchy boundary conditions used to simulate stream–aquifer interactions

This work examines the simulation of stream–aquifer interactions as grids are refined vertically and horizontally and suggests that traditional methods for calculating conductance can produce inappropriate values when the grid size is changed. Instead, different grid resolutions require different estimated values. Grid refinement strategies considered include global refinement of the entire model and local refinement of part of the stream. Three methods of calculating the conductance of the Cauchy boundary conditions are investigated. Single- and multi-layer models with narrow and wide streams produced stream leakages that differ by as much as 122% as the grid is refined. Similar results occur for globally and locally refined grids, but the latter required as little as one-quarter the computer execution time and memory and thus are useful for addressing some scale issues of stream–aquifer interactions. Results suggest that existing grid-size criteria for simulating stream–aquifer interactions are useful for one-layer models, but inadequate for three-dimensional models. The grid dependence of the conductance terms suggests that values for refined models using, for example, finite difference or finite-element methods, cannot be determined from previous coarse-grid models or field measurements. Our examples demonstrate the need for a method of obtaining conductances that can be translated to different grid resolutions and provide definitive test cases for investigating alternative conductance formulations.     

国家基础地理信息1:5万数据库内容完善的需求调查

通过对国内外基础地理信息的发展进行探讨,在深入研究我国现有1∶5万地形要素的内容和属性信息的基础上,设计调查问卷,其内容涉及国家基础地理信息在各专业部门/行业的用途、使用目的、数据精度和数据更新周期以及1∶5万数据库的要素内容、属性信息等,采用小组座谈、网络等形式进行了用户需求调查。     

A high order Godunov scheme with constrained transport and adaptive mesh refinement for astrophysical and geophysical MHD

We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying conservation law which, in our formulation, results in a “finite-surface” scheme for the induction equation. This naturally leads to the well-known “constrained transport” method, with additional continuity requirement on the magnetic field representation. The second ingredient in the MUSCL scheme is the predictor step that ensures second order accuracy both in space and time. We explore specific constraints that the mathematical properties of the induction equations place on this predictor step, showing that three possible variants can be considered. We show that the most aggressive formulations reach the same level of accuracy than the other ones, at a lower computational cost. More interestingly, these schemes are compatible with the Adaptive Mesh Refinement (AMR) framework. It has been implemented in the AMR code RAMSES. It offers a novel and efficient implementation of a second order scheme for the induction equation. The scheme is then adaptated to solve for the full MHD equations using the same methodology. Through a series of test problems, we illustrate the performances of this new code using two different MHD Riemann solvers (Lax–Friedrich and Roe) and the need of the Adaptive Mesh Refinement capabilities in some cases. Finally, we show its versatility by applying it to the ABC dynamo problem and to the collapse of a magnetized cloud core.     

Geostatistical approaches to refinement of digital elevation data

Data refinement refers to the processes by which a dataset’s resolution, in particular, the spatial one, is refined, and is thus synonymous to spatial downscaling. Spatial resolution indicates measurement scale and can be seen as an index for regular data support. As a type of change of scale, data refinement is useful for many scenarios where spatial scales of existing data, desired analyses, or specific applications need to be made commensurate and refined. As spatial data are related to certain data support, they can be conceived of as support-specific realizations of random fields, suggesting that multivariate geostatistics should be explored for refining datasets from their coarser-resolution versions to the finer-resolution ones. In this paper, geostatistical methods for downscaling are described, and were implemented using GTOPO30 data and sampled Shuttle Radar Topography Mission data at a site in northwest China, with the latter’s majority grid cells used as surrogate reference data. It was found that proper structural modeling is important for achieving increased accuracy in data refinement; here, structural modeling can be done through proper decomposition of elevation fields into trends and residuals and thereafter. It was confirmed that effects of semantic differences on data refinement can be reduced through properly estimating and incorporating biases in local means.     

自适应自然单元法研究——自适应细化

对h型自适应自然单元法的自适应细化方案进行了初步研究。在ZZ误差分析的基础上实现了节点的自动加密,使得随着自适应细化的进行,求解误差减小,而且误差分布趋于均匀。在数值分析中,主要有两种误差来源--插值误差和积分误差。随着节点的加密,Delaunay三角形的尺寸随之减小,三角形内的应力场趋于线性分布,那么插值误差和积分误差也都会随之减小。因而,h型自适应分析可以同时减小上述两种误差而达到不断提高求解精度的目的。由于自然单元法求解依赖于求解域内离散节点的Voronoi结构,建议的细化方案中新节点的引入只需局部调整Delaunay结构,算法的实现极为容易,程序实现简单、高效。研究表明,建议的自适应方案是可行的,自然单元法特别适合进行h型自适应分析。