Automatic mapping of lunar landforms using DEM-derived geomorphometric parameters

Developing approaches to automate the analysis of the massive amounts of data sent back from the Moon will generate significant benefits for the field of lunar geomorphology. In this paper, we outline an automated method for mapping lunar landforms that is based on digital terrain analysis. An iterative self-organizing (ISO) cluster unsupervised classification enables the automatic mapping of landforms via a series of input raster bands that utilize six geomorphometric parameters. These parameters divide landforms into a number of spatially extended, topographically homogeneous segments that exhibit similar terrain attributes and neighborhood properties. To illustrate the applicability of our approach, we apply it to three representative test sites on the Moon, automatically presenting our results as a thematic landform map. We also quantitatively evaluated this approach using a series of confusion matrices, achieving overall accuracies as high as 83.34% and Kappa coefficients (K) as high as 0.77. An immediate version of our algorithm can also be applied for automatically mapping large-scale lunar landforms and for the quantitative comparison of lunar surface morphologies.     

树轮稳定同位素记录:进展、问题及展望

树轮稳定同位素比率(δ13 C、δ18O 和δD)变化与树木生长季的气候和环境条件密切相关。控制实验表明, 树轮稳定同位素比率变化可较好地记录树木叶片气孔活动、水分来源、水分利用效率等相关信息, 已成为古气候和全球变化生态学研究的重要代用指标和手段之一。树轮稳定同位素记录能够较好地保留低频气候环境变化信息, 与年轮宽度相结合, 可有效地进行区域古气候重建和全球变化生态学方面的研究。文章综述了树轮稳定同位素的分馏机理及其发展, 系统阐述和评价了利用树轮稳定同位素指标进行古气候和古环境的相关研究成果, 并指出树轮稳定同位素研究在中国未来亟须着重发展的方向。     

Automated mapping of rock slope geometry,kinematics and stability with RSS-GIS

A GIS-implemented, deterministic approach for the automated spatial evaluation of geometrical and kinematical properties of rock slope terrains is presented. Based on spatially distributed directional information on planar geological fabrics and DEM-derived topographic attribute data, the internal geometry of rock slopes can be characterized on a grid cell basis. For such computations, different approaches for the analysis and regionalization of available structural directional information applicable in specific tectonic settings are demonstrated and implemented in a GIS environment. Simple kinematical testing procedures based on feasibility criteria can be conducted on a pixel basis to determine which failure mechanisms are likely to occur at particular terrain locations. In combination with hydraulic and strength data on geological discontinuities, scenario-based rock slope stability evaluations can be performed. For conceptual investigations on rock slope failure processes, a GIS-based specification tool for a 2-D distinct element code (UDEC) was designed to operate with the GIS-encoded spatially distributed rock slope data. The concepts of the proposed methodology for rock slope hazard assessments are demonstrated at three different test sites in Germany.     

起伏地形下我国太阳直接辐射空间制图

建立了起伏地形下太阳直接辐射分布式计算模型 ;成功地解决了起伏地形中地形相互遮蔽对太阳直接辐射影响的难题 ;采用数据集群技术 ,探讨了不同数据集情况下太阳直接辐射计算模式的时空有效性 ;以 1km× 1km分辨率的DEM数据作为地形的综合反映 ,完成了我国 1km× 1km分辨率各月气候平均太阳直接辐射空间分布制图     

GPS网平差观测量选取的理论分析

由多台GPS接收机在不同时段进行静态观测而组成的观测网形,当采用高精度的科研软件与常用商用软件解算基线时,其基线向量的方差阵不同,前者的方差阵结构是完全的、严密的。同时,在进行GPS网平差时,选取GPS基线的方式及数量的不同,构造的网形也就不同,采用不同的数学模型其得到的平差结果和精度指标也不相同。文章先从理论上进行分析,后用示例进行比较论证。     

应用NOAA/AVHRR数据测算局地水稻种植面积方法研究

本文采用在模糊监督分类中增加迭代过程的方法－－模糊监督分类一迭代法,在分解混合像元的基础上,利用ＡＶＨＲＲ数据求算水稻种植面积。根据稻田与旱地存在温度差异的特点,在分析ＡＶＨＲＲ数据统计特征的基础上,增加了第３和４两通道参加模糊监督分类,从而增加了分类橡元值矢量维数,增强了对水稻的鉴别能力。反采用的迭代法是收敛具有效的,经本方法输出的水稻种植面积百分含量图与实际水稻分布十分吻合,测算出的面积值与Ｔ     

Singular value decomposition and cluster analysis as regression diagnostics tools for geodetic applications

It is well known that high-leverage observations significantly affect the estimation of parameters. In geodetic literature, mainly redundancy numbers are used for the detection of single high-leverage observations or of single redundant observations. In this paper a further objective method for the detection of groups of important and less important (and thus redundant) observations is developed. In addition, the parameters which are predominantly affected by these groups of observations are identified. This method thus complements other diagnostics tools, such as, e.g., multiple row diagnostics methods as described in statistical literature (see, e.g., Belsley et al. in Regression diagnostics: identifying influential data and sources of collinearity. Wiley, New York, 1980). The method proposed in this paper is based on geometric aspects of adjustment theory and uses the singular value decomposition of the design matrix of an adjustment problem together with cluster analysis methods for regression diagnostics. It can be applied to any geodetic adjustment problem and can be used for the detection of (groups of) observations that significantly affect the estimated parameters or that are of negligible impact. One of the advantages of the proposed method is the improvement of the reliability of observation plans and thus the reduction of the impact of individual observations (and outliers) on the estimated parameters. This is of particular importance for the very long baseline interferometry technique which serves as an application example of the regression diagnostics tool.     

The determination of potential difference by the joint application of measured and synthetical gravity data: a case study in Hungary

In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10⁻³mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.     

Characterization of diurnal cycles in ZTD from a decade of global GPS observations

The diurnal cycle of the tropospheric zenith total delay (ZTD) is one of the most obvious signals for the various physical processes relating to climate change on a short time scale. However, the observation of such ZTD oscillations on a global scale with traditional techniques (e.g. radiosondes) is restricted due to limitations in spatial and temporal resolution. Nowadays, the International GNSS Service (IGS) provides an important data source for investigating the diurnal and semidiurnal cycles of ZTD and related climatic signals. In this paper, 10 years of ZTD data from 1997 to 2007 with a 2-hour temporal resolution are derived from global positioning system (GPS) observations taken at 151 globally distributed IGS reference stations. These time series are used to investigate diurnal and semidiurnal oscillations. Significant diurnal and semidiurnal oscillations of ZTD are found for all GPS stations used in this study. The diurnal cycles (24 hours period) have amplitudes between 0.2 and 10.9 mm with an uncertainty of about 0.5 mm and the semidiurnal cycles (12 h period) have amplitudes between 0.1 and 4.3 mm with an uncertainty of about 0.2 mm. The larger amplitudes of the diurnal and semidiurnal ZTD cycles are observed in the low-latitude equatorial areas. The peak times of the diurnal cycles spread over the whole day, while the peak value of the semidiurnal cycles occurs typically about local noon. These GPS-derived diurnal and semidiurnal ZTD signals are similar with the surface pressure tides derived from surface synoptic pressure observations, indicating that atmospheric tides are the main driver of the diurnal and semidiurnal ZTD variations.     

Preprocessing of gravity gradients at the GOCE high-level processing facility

One of the products derived from the gravity field and steady-state ocean circulation explorer (GOCE) observations are the gravity gradients. These gravity gradients are provided in the gradiometer reference frame (GRF) and are calibrated in-flight using satellite shaking and star sensor data. To use these gravity gradients for application in Earth scienes and gravity field analysis, additional preprocessing needs to be done, including corrections for temporal gravity field signals to isolate the static gravity field part, screening for outliers, calibration by comparison with existing external gravity field information and error assessment. The temporal gravity gradient corrections consist of tidal and nontidal corrections. These are all generally below the gravity gradient error level, which is predicted to show a 1/f behaviour for low frequencies. In the outlier detection, the 1/f error is compensated for by subtracting a local median from the data, while the data error is assessed using the median absolute deviation. The local median acts as a high-pass filter and it is robust as is the median absolute deviation. Three different methods have been implemented for the calibration of the gravity gradients. All three methods use a high-pass filter to compensate for the 1/f gravity gradient error. The baseline method uses state-of-the-art global gravity field models and the most accurate results are obtained if star sensor misalignments are estimated along with the calibration parameters. A second calibration method uses GOCE GPS data to estimate a low-degree gravity field model as well as gravity gradient scale factors. Both methods allow to estimate gravity gradient scale factors down to the 10⁻³ level. The third calibration method uses high accurate terrestrial gravity data in selected regions to validate the gravity gradient scale factors, focussing on the measurement band. Gravity gradient scale factors may be estimated down to the 10⁻² level with this method.