The hierarchical tessellation model and its use in spatial analysis

The hierarchid tessellation model belongs to a class of spatial data models based on the recursive decomposition of space. The quadtree is one such tessellation and is characterized by square cells and a 1:4 decomposition ratio. To relax these constraints in the tessellation, a generalized hierarchical tessellation data model, called Adaptive Recursive Tessellations (ART), has been proposed. ART increases flexibility in the tessellation by the use of rectangular cells and variable decomposition ratios. In ART, users can specify cell sizes which are intuitively meaningful to their applications, or which can reflect the scales of data. ART is implemented in a data structure called Adaptive Recursive Run-Encoding (ARRE), which is a variant of two-dimensional run-encoding whose running path can vary with the different tessellation structures incorporated in an ART model. Given the recognition of the benefits of implementing statistical spatial analysis in GIS, the use of hierarchical tessellation models such as ART in spatial analysis is discussed. Three examples are introduced to show how ART can: (1) be applied to solve the quadrat size problem in quadrat analysis of point patterns; (2) act as the data model in the variable resolution block kriging technique for geostatistical data to reduce variation in kriging error; and (3) facilitate the evaluation of spatial autocorrelation for area data at multiple map resolutions via the construction of a connectivity matrix for calculating spatial autocorrelation indices based on ARRE.     

Delineating agricultural field boundaries from TM imagery using dyadic wavelet transforms

Agricultural field boundary information is important in the study of agricultural crops by remote sensing. In this paper, the dyadic wavelet transform is adopted to extract field boundaries directly from TM imagery. A field boundary is a set of locally connected pixels characterized by their abrupt spectral intensity variation (singularity) on the image. Preliminary extraction is made on the near-infrared Band 4 (0.76–0.90 μm) image. It is found that most of the boundaries are well characterized by the wavelet modulus maxima curves, particularly for fully vegetated fields, while better results are obtained from the NDVI image for cropped field boundaries. To extract additional boundaries while suppressing unwanted noisy edges, the wavelet transform is then applied to the Wetness image calculated by the Tasselled Cap Transform. A technique is implemented to integrate the edges extracted from each source. The combined results show that the majority of the field boundaries are delineated. The use of a single date image fails to delineate some of the fallow field boundaries and boundaries between two cropped fields with similar spectral properties. Although techniques have been applied to reduce noise, still great difficulties are encountered when trying to delineate edges that are field boundaries from all the detected singularities. Manual editing is still required to trace the remaining boundaries to obtain a complete boundary map. It is recommended that a multi-temporal data set be used to obtain a complete delineation with a fully automated algorithm.     

Response of GPS occultation signals to atmospheric gravity waves and retrieval of gravity wave parameters

We show that the amplitude of the Global Positioning System (GPS) signals in the radio occultation (RO) experiments is an indicator of the activity of the gravity waves (GW) in the atmosphere. The amplitude of the GPS RO signals is more sensitive to the atmospheric wave structures than is the phase. Early investigations used only the phase of the GPS occultation signals for statistical investigation of the GW activity in the height interval 10–40 km on a global scale. In this study, we use the polarization equations and Hilbert transform to find the 1-D GW radio image in the atmosphere by analyzing the amplitude of the RO signal. The radio image, also called the GW portrait, consists of the phase and amplitude of the GW as functions of height. We demonstrate the potential of this method using the amplitude data from GPS/Meteorology (GPS/MET) and satellite mission Challenge Mini-satellite Payload (CHAMP) RO events. The GW activity is nonuniformly distributed with the main contribution associated with the tropopause and the secondary maximums related to the GW breaking regions. Using our method we find the vertical profiles of the horizontal wind perturbations and its vertical gradient associated with the GW influence. The estimated values of the horizontal wind perturbations are in fairly good agreement with radiosonde data. The horizontal wind perturbations v(h) are ±1 to ±5 m s with vertical gradients dv/dh ±0.5 to ±15 m s km at height 10–40 km. The height dependence of the GW vertical wavelength was inferred through the differentiation of the GW phase. Analysis of this dependence using the dispersion relationship for the GW gives the estimation of the projection of the horizontal background wind velocity on the direction of the GW propagation. For the event considered, the magnitude of this projection changes between 1.5 and 10 m s at heights of 10–40 km. We conclude that the amplitude of the GPS occultation signals contain important information about the wave processes in the atmosphere on a global scale.     

巴黎能力建设委员会助力全球气候治理

在全球气候治理中,能力建设涉及减缓、适应等多个方面.相对于发达国家,发展中国家在应对气候变化能力方面存在更多挑战,也亟待支持.巴黎能力建设委员会(PCCB)是《联合国气候变化框架公约》(以下简称公约)下属能力建设议题唯一的常设机构.在《巴黎协定》的授权下, PCCB于2016年成立,致力于解决发展中国家在实施和加强应对...     

Past,present and future vegetation-cloud feedbacks in the Amazon Basin

To begin exploring the underlying mechanisms that couple vegetation to cloud formation processes, we derive the lifting condensation level (LCL) to estimate cumulus cloud base height. Using a fully coupled land–ocean–atmosphere general circulation model (HadCM3LC), we investigate Amazonian forest feedbacks on cloud formation over three geological periods; modern-day (a.d. 1970–1990), the last glacial maximum (LGM; 21 kya), and under a future climate scenario (IS92a; a.d. 2070–2090). Results indicate that for both past and future climate scenarios, LCL is higher relative to modern-day. Statistical analyses indicate that the 800 m increase in LCL during the LGM is related primarily to the drier atmosphere promoted by lower tropical sea surface temperatures. In contrast, the predicted 1,000 m increase in LCL in the future scenario is the result of a large increase in surface temperature and reduced vegetation cover.     

Controlling factors of the δ 11B-pH proxy and its research direction

Significant boron isotope fractionation occurs in nature (?70 ‰ to +75 ‰) due to the high geochemical reactivity of boron and the large relative mass difference between ¹⁰B and ¹¹B. Since the 1990s, reconstruction of ancient seawater pH using the isotopic composition of boron in bio-carbonates (δ ¹¹B_carb), and then calculation of the past pCO₂ have become important issues for the international isotope geochemistry community, and are called the δ ¹¹B-pH proxy. Although many achievements have been made by this proxy, various aspects of boron systematics require rigorous evaluation. Based on the previous researches, mechanism of boron isotope fractionation, variation of boron isotope (δ ¹¹B) in nature (especially in bio-carbonates) and controlling factors of the δ ¹¹B-pH proxy, such as the dissociation constant of B(OH)₃ in seawater (pK_a), the δ ¹¹B of seawater (δ ¹¹B_SW), the boron isotopic fractionation factor between B(OH) ₄ ^? and B(OH)₃ (α _4–3), and the incorporated species of boron into bio-carbonates, are reviewed in detail and the research directions of this proxy are proposed. Generally, the controversy about pK_a, δ ¹¹B_sw, and α _4–3 is relatively less, but whether boron incorporated into bio-carbonates only in the form of B(OH) ₄ ^? remains doubtful. In the future, it is required that the physicochemical processes that control boron incorporation into carbonates be rigorously characterized and that the related chemical and isotopic fractionation be quantified. It is also necessary and important to establish a “best-fit empirically equation” between δ ¹¹B_carb and pH of seawater based on the precipitation experiments of inorganic or culture experiments of corals or foraminifera. In addition, extended application of the δ ¹¹B-pH proxy to the earlier part of the Phanerozoic relying on the Brachiopods is worthy of studying. Like other geochemical indicators, there are limiting factors of δ ¹¹B; however, it remains a very powerful tool in the reconstruction of past seawater pH at present.     

华北东部地幔改造作用和置换作用:单斜辉石激光探针研究

对山东蒙阴古生代金伯利岩中石榴石方辉橄榄岩、河南鹤壁新生代玄武岩中 铬尖晶石方辉橄榄岩、山东栖霞新生代玄武岩中央晶石二辉橄榄岩和山东山旺新生代 玄武岩中尖晶石（富ＣｐＸ）二辉橄榄岩单斜辉石的常量、微量元素的研究和对比表明：华 北地块古生代稳定存在的克拉通岩石圈地幔在中、新生代其东部被具大洋地幔性质的 显生宙地幔不均匀地发生了置换作用,残存于南北重力梯度带（岩石圈减薄过渡带）上 的鹤壁浅部克拉通地幔除保存有置换作用前改造作用的丰富信息外,也包含有其深部 被置换后的地幔改造作用信息;郯庐断裂带东部栖霞的新生地幔仍保留有较多的大洋 地幔成分特征,而郯庐断裂带内山旺的新生地幔受后续地幔改造作用的影响有向大陆 地幔演化的趋向．发生于华北地块东部中、新生代的岩石圈巨厚减薄作用与地幔的置换 作用有关,由渐进的富ＣＯ_２流体或碳酸岩熔体对岩石圈地幔的交代改造作用以及由此 所引起的熔体抽取作用是完成新生地幔对古老地幔置换作用的关键．郯庐断裂带等岩 石圈深大断裂为地幔改造作用和地幔置换作用提供了良好的通道作用．     

人工侧向冻结条件下土的冻结试验

在人工侧向冻结、正弦和恒温2种冻结温度、开放和封闭条件下,对含水率为32 %的重塑粉质黏土进行了4组冻结试验。结果表明,补水条件相同时侧向正弦冻结试样的冻结速度比封闭正弦冻结试验的缓慢,当其冻结状态进入稳态之后,试样的整体温度场随着侧向冻结温度的正弦变化而出现相应的波动,但径向各段的温度梯度的相位滞后于正弦温度相位半个周期;正弦冻结试样与恒温冻结试样在同一条件下的冻胀量相近,封闭体系下自上而下轴向冻结试样的冻胀量高于同条件下侧向冻结试样的冻胀量;人工冻结工程中应考虑大气环境温度变化（如昼夜温差）及冷能侧向输入土体内时对冻土温度场的变化的影响。     

西天山喇嘛苏铜矿成矿斑岩的岩石学、地球化学特征及成因探讨

位于西天山别珍套-科古琴晚古生代岛弧西段的喇嘛苏铜矿床是区内最大的铜矿床,与成矿作用有着密切关系的斑岩体为英云闪长斑岩、花岗闪长斑岩,是同源岩浆分异演化的产物,且花岗闪长斑岩可能属于岩浆演化晚期的产物。本区成矿斑岩的主量、微量元素和Sr-Nd同位素地球化学特征表明,其富集大离子亲石元素,而相对亏损高场强元素,出现了较为明显的Ta、Nb负异常,初始锶同位素ISr和εNd(t=390Ma)值分别为0.7072～0.7076和-0.32～0.17,显示壳幔混合源的特征,利用Sr和Nd同位素估算其源区物质约有50%来源于地壳。岩石地球化学特征指示了其为典型钙碱性火山弧花岗岩,暗示其形成于大陆弧环境。结合区域地质背景,推测本区成矿斑岩是在洋壳俯冲作用下发生部分熔融,交代原先的地幔楔,并混合了部分下地壳的物质,经历分离结晶作用的产物,其形成可能与晚古生代准噶尔洋板块向南的俯冲作用有关。结合东西天山的成矿斑岩的地球化学特征对比研究,岩浆源区的差别可能导致不同类型斑岩型矿床的形成,斑岩型铜矿床的形成较斑岩型钼矿床可能有更少的地壳物质贡献。