SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

地震折射波的虚拟射线理论及其应用

本文提出了一种地震折射液的虚拟射线理论,根据该理论,只要已知地表层的速度,即可直接由折射波信息提取地震参数,从而可实现折射界面的反演。 文中通过实际介质模型的计算机实验结果,验证了该理论的正确性。     

引入差比关系法分析西北太平洋TOPEX/POSEIDON卫星高度计测高数据

TOPEX/POSEIDON(T/P)卫星高度计数据信息中存在周期成分混淆问题.对其中的一类混淆引入差比关系方法对混淆的分潮进行分离.卫星轨道交叉点资料包括升轨和降轨资料,资料量比沿轨点资料多1倍,经分析发现:在已有为期6a多的观测资料时间序列中,在沿轨处混淆的分潮如K₁和SSA在交叉点处不再混淆,可以直接分离.因此首先对交叉点资料进行调和分析.然后由交叉点的分析结果得到分潮间的差比关系,处理到相近的沿轨点处,从而得到沿轨点的调和常数.用引入差比关系方法,对西北太平洋海区6a多的T/P卫星高度计资料进行了潮汐分析,并与沿岸及岛屿验潮站资料进行了比较,所得结果较满意.     

高斯投影公式的机助推导

随着GPS应用的普及,高斯投影在测量中出现的机会越来越多。许多文献上都给出了推导过程,这个推导过程的运算量是非常大的,而且计算精度难于提高。经过研究,一种易于编程的递推公式被发现,虽然这个公式在计算量上没有减少,但是可以较好地解决精度问题。     

A framework of region-based spatial relations for non-overlapping features and its application in object based image analysis

Object based image analysis (OBIA) is an approach increasingly used in classifying high spatial resolution remote sensing images. Object based image classifiers first segment an image into objects (or image segments), and then classify these objects based on their attributes and spatial relations. Numerous algorithms exist for the first step of the OBIA process, i.e. image segmentation. However, less research has been conducted on the object classification part of OBIA, in particular the spatial relations between objects that are commonly used to construct rules for classifying image objects and refining classification results. In this paper, we establish a context where objects are areal (not points or lines) and non-overlapping (we call this “single-valued” space), and propose a framework of binary spatial relations between segmented objects to aid in object classification. In this framework, scale-dependent “line-like objects” and “point-like objects” are identified from areal objects based on their shapes. Generally, disjoint and meet are the only two possible topological relations between two non-overlapping areal objects. However, a number of quasi- topological relations can be defined when the shapes of the objects involved are considered. Some of these relations are fuzzy and thus quantitatively defined. In addition, we define the concepts of line-like objects (e.g. roads) and point-like objects (e.g. wells), and develop the relations between two line-like objects or two point-like objects. For completeness, cardinal direction relations and distance relations are also introduced in the proposed context. Finally, we implement the framework to extract roads and moving vehicles from an aerial photo. The promising results suggest that our methods can be a valuable tool in defining rules for object based image analysis.     

Environmental tracers as indicators of karst conduits in groundwater in South Dakota,USA

Environmental tracers sampled from the carbonate Madison aquifer on the eastern flank of the Black Hills, South Dakota, USA indicated the approximate locations of four major karst conduits. Contamination issues are a major concern because these conduits are characterized by direct connections to sinking streams, high groundwater velocities, and proximity to public water supplies. Objectives of the study were to estimate approximate conduit locations and assess possible anthropogenic influences associated with conduits. Anomalies of young groundwater based on chlorofluorocarbons (CFCs), tritium, and electrical conductivity (EC) indicated fast moving, focused flow and thus the likely presence of conduits. δ¹⁸O was useful for determining sources of recharge for each conduit, and nitrate was a useful tracer for assessing flow paths for anthropogenic influences. Two of the four conduits terminate at or near a large spring complex. CFC apparent ages ranged from 15 years near conduits to >50 years in other areas. Nitrate-N concentrations >0.4 mg/L in groundwater were associated with each of the four conduits compared with concentrations ranging from <0.1 to 0.4 mg/L in other areas. These higher nitrate-N concentrations probably do not result from sinking streams but rather from other areas of infiltration.     

Earthquake hazard of the northern parts of the Bohemian Massif and Western Carpathians

The new procedure of earthquake hazard evaluation developed by Kijko and Sellevoll is tested and applied for the border region of Czechoslovakia and Poland. The new method differs from the conventional approach. It incorporates the uncertainty of earthquake magnitudes, and accepts mixed data containing only large historical events and recent, complete catalogues. Seismic hazard has been calculated for nine regions determined in the border area. In the investigated area, data of historical catalogues are uncertain or, in many cases, the epicentral intensities are unknown. Thus, a number of assumptions have to be adopted in data preparation of catalogues since the year 1200. The calculated values of parameters b in the Gutenberg-Richter frequency-intensity relation as well as the return periods, seem to be reasonable and are generally confirmed by the results obtained from catalogues for the last 80–130 years.     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

Monte Carlo simulation of mixing in Ca3Fe2Ge3O12–Ca4Ge4O12 garnets and implications for the thermodynamic stability of pyrope–majorite solid solution

Static lattice energy calculations, based on empirical pair potentials, were performed for a large set of structures differing in the arrangement of octahedral cations within the garnet 2 × 2 × 2 supercell. The compositions of these structures varied between Ca₃Fe₂Ge₃O₁₂ and Ca₄Ge₄O₁₂. The energies were cluster expanded using pair and quaternary terms. The derived ordering constants were used to constrain Monte Carlo simulations of temperature-dependent mixing properties in the ranges of 1,073–3,673 K and 0–10 GPa. The free energies of mixing were calculated using the method of thermodynamic integration. The calculations predict a wide miscibility gap between Fe-rich (cubic) and Fe-pure (tetragonal) garnets consistent with recent experimental observations of Iezzi et al. (Phys Chem Miner 32:197–207, 2005). It is shown that the miscibility gap arises due to a very strong cation ordering at the Fe-pure composition, driven by the charge difference between Ca²⁺ and Ge⁴⁺ cations. The structural and thermodynamic analogies between Ca–Ge and Mg–Si systems suggest that a similar miscibility gap should exist between pyrope and Mg–Si-majorite.     

Use of electrical resistivity methods for detecting subsurface fresh and saline water and delineating their interfacial configuration: a case study of the eastern Dead Sea coastal aquifers, Jordan

The alluvial aquifer is the primary source of groundwater along the eastern Dead Sea shoreline, Jordan. Over the last 20 years, salinity has risen in some existing wells and several new wells have encountered brackish water in areas thought to contain fresh water. A good linear correlation exists between the water resistivity and the chloride concentration of groundwater and shows that the salinity is the most important factor controlling resistivity. Two-dimensional electrical tomography (ET) integrated with geoelectrical soundings were employed to delineate different water-bearing formations and the configuration of the interface between them. The present hydrological system and the related brines and interfaces are controlled by the Dead Sea base level, presently at 410 m b.s.l. Resistivity measurements show a dominant trend of decreasing resistivity (thus increasing salinity) with depth and westward towards the Dead Sea. Accordingly, three zones with different resistivity values were detected, corresponding to three different water-bearing formations: (1) strata saturated with fresh to slightly brackish groundwater; (2) a transition zone of brine mixed with fresh to brackish groundwater; (3) a water-bearing formation containing Dead Sea brine. In addition, a low resistivity unit containing brine was detected above the 1955 Dead Sea base level, which was interpreted as having remained unflushed by infiltrating rain.