Geophysical implications of Izu–Bonin mantle wedge hydration from chemical geodynamic modeling

Using two-dimensional dynamic models of the Northern Izu–Bonin (NIB) subduction zone, we show that a particular localized low-viscosity (η_LV =  3.3 × 10¹⁹ − 4.0 × 10²⁰ Pa s), low-density (Δρ ∼ −10 kg/m³ relative to ambient mantle) geometry within the wedge is required to match surface observations of topography, gravity, and geoid anomalies. The hydration structure resulting in this low-viscosity, low-density geometry develops due to fluid release into the wedge within a depth interval from 150 to 350 km and is consistent with results from coupled geochemical and geodynamic modeling of the NIB subduction system and from previous uncoupled models of the wedge beneath the Japan arcs. The source of the fluids can be either subducting lithospheric serpentinite or stable hydrous phases in the wedge such as serpentine or chlorite. On the basis of this modeling, predictions can be made as to the specific low-viscosity geometries associated with geophysical surface observables for other subduction zones based on regional subduction parameters such as subducting slab age.     

国家基础地理信息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.     

水平面代替水准面对两井贯通的影响

理论结合实际,用详实的数据说明了在两井贯通时,把水准面当作水平面处理的界限,避免了盲目应用而产生严重后果.     

新疆全区域似大地水准面的建立

简述了大地水准面精化的方法、原理及新疆似大地水准面确定的过程。     

High-precision geoid determination in small areas: A case study in Doñana National Park (Spain)

Doñana National Park is an area of approximately 500 km² located on the SW coast of Spain that shows one of the greatest geoid gradients on the entire Iberian Peninsula, due to its peculiar tectonic characteristics. So, it is necessary to elaborate an accurate geoid model that can be used with GPS for precise surveying, since the existing ones are insufficient, due to their poor resolution and their limited adaptation to a small area with such a strong gradient. The least squares prediction method was tested in order to obtain the undulation from GPS/orthometric points. The results obtained were unsatisfactory because of the strong geoid gradient. In order to improve accuracy the remove-restore technique was used. Global geopotential model EIGEN-CG01C and a Digital Elevation Model (DEM) with a 25 × 25 m resolution and an accuracy better than 3 m were used. Thus, the final geometrical geoid obtained reaches the precision required by other disciplines (3 cm in any point within the Park). Particularly, the geoid model has allowed for the acquisition of a precision DEM that is essential to formulate a hydrodynamic model for the Doñana marsh functions.     

Contributions of terrestrial and GRACE data to the study of the secular geoid changes in North America

This paper tests and discusses different statistical methods for modelling secular rates of change of the geoid in North America. In particular, we use the method of principal component/empirical orthogonal functions (PC/EOF) analysis to model the geoid rates from Gravity Recovery and Climate Experiment (GRACE) satellite data. As demonstrated, the PC/EOF analysis is useful for studying the contributions from different signals (mainly residual hydrology signals and leakage effects) to the GRACE-derived geoid rates. The PC/EOF analysis leads to smaller geoid rates compared to the conventional least-squares fitting of a trend and annual and semi-annual cycles to the time series of the spherical harmonic coefficients. This is because we filter out particular spatiotemporal modes of the regional geoid changes.We apply the method of least-squares collocation with parameters to combine terrestrial data (GPS vertical velocities from the Canadian Base Network and terrestrial gravity rates from the Canadian Gravity Standardization Net) with the GRACE-derived vertical motion to obtain again the geoid rates. The combined model has a peak geoid rate of 1.4 mm/year in the southeastern area of Hudson Bay contrary to the GRACE-derived geoid rates that show a large peak of 1.6–1.7 mm/year west of Hudson Bay. We demonstrate that the terrestrial data, which have a longer time span than the GRACE data, are important for constraining the GRACE-derived secular signal in the areas that are well sampled by the data.     

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

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

迈入新千年的大地测量学——第22届IUGG大会有关大地测量部分的技术总结

本文较系统地介绍了１９９９年７月在英国伯明翰举行的第２２届ＩＵＧＧ大会反映出来的大地测量的最新进展情况。全文共分４个部分：卫星定位技术的新进展,地球重力场的研究进展,地球动力学的新进展及大地测量理论与方法新进展。对其中若干领域,还提出我国今后发展的建议。